如何在 R 中创建广义加性混合模型输出的汇总表

Question

我正在使用广义加性混合模型（GAMM，带有 mgcv 包），并且我正在尝试找到一种方法来在 R 中创建模型的日志就绪汇总表，并将它们打印为我的手稿的图形。我找到了使用 GLM、GLMM、LM 和 LMM 以及 R Markdown 格式来实现此目的的不同方法。但我找不到 GAM 或 GAMM 的任何内容。这可能吗？

谢谢！

Answer 1

s$s.table

只是一个普通矩阵，所以：

library(mgcv)
library(stargazer)

set.seed(0) 
dat <- gamSim(1, n = 200, scale = 2)

b <- gamm(y~s(x0)+s(x1)+s(x2)+s(x3),data=dat)

s <- summary(b$gam)
class(s$s.table)
## [1] "matrix" "array" 

stargazer(s$s.table)

给予

% Table created by stargazer v.5.2.3 by Marek Hlavac, Social Policy Institute. E-mail: marek.hlavac at gmail.com
% Date and time: Wed, Dec 04, 2024 - 2:22:31 PM
\begin{table}[!htbp] \centering 
  \caption{} 
  \label{} 
\begin{tabular}{@{\extracolsep{5pt}} ccccc} 
\\[-1.8ex]\hline 
\hline \\[-1.8ex] 
 & edf & Ref.df & F & p-value \\ 
\hline \\[-1.8ex] 
s(x0) & $2.489$ & $2.489$ & $5.266$ & $0.006$ \\ 
s(x1) & $2.209$ & $2.209$ & $56.511$ & $0$ \\ 
s(x2) & $7.305$ & $7.305$ & $36.985$ & $0$ \\ 
s(x3) & $1.000$ & $1.000$ & $0.212$ & $0.646$ \\ 
\hline \\[-1.8ex] 
\end{tabular} 
\end{table}