summary_gamm: Summarize a GAMM
In m-clark/gammit: Using Generalized Additive Mixed Models

Description Usage Arguments Details Value See Also Examples

View source: R/summary_gamm.R

Summarize a gam model in a clean mixed effects style.

1	summary_gamm(model, digits = 3)

`model`	The mgcv model
`digits`	number of digits to display

This displays the variance components and fixed effects from a gam model. Assumes an mgcv model of the form gam(... + s(g, bs='re')), but should work with just about any gam with a smooth term.

Invisibly returns a list with the variance components via extract_vc and fixed effects parameter, labeled vc and fe respectively.

extract_vc

library(mgcv)
library(lme4)

lmer_model <- lmer(Reaction ~ Days + (Days || Subject), data = sleepstudy)
ga_model <- gam(Reaction ~ Days + s(Subject, bs = "re") + s(Days, Subject, bs = "re"),
  data = sleepstudy,
  method = "REML"
)

summary(lmer_model)
summary_gamm(ga_model)