Models
iron_efficiency_GAMM <- gamm4(brain ~ s(age, k=4) + sex + Overall_Efficiency + averageManualRating + roi, data=r2star, REML=TRUE, random=~(1|bblid))
Results
summary(iron_efficiency_GAMM$gam)
##
## Family: gaussian
## Link function: identity
##
## Formula:
## brain ~ s(age, k = 4) + sex + Overall_Efficiency + averageManualRating +
## roi
##
## Parametric coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 17.05773 0.32404 52.642 <2e-16 ***
## sex.L 0.03517 0.05094 0.690 0.490
## Overall_Efficiency 0.07799 0.05122 1.523 0.128
## averageManualRating -0.11071 0.16693 -0.663 0.507
## roicaudate -1.59733 0.07245 -22.047 <2e-16 ***
## roipallidum 5.30186 0.07243 73.204 <2e-16 ***
## roiputamen -0.10111 0.07237 -1.397 0.162
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Approximate significance of smooth terms:
## edf Ref.df F p-value
## s(age) 2.14 2.14 88.49 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## R-sq.(adj) = 0.669
## lmer.REML = 18379 Scale est. = 2.9427 n = 4493