Models
iron_accuracy_GAMM <- gamm4(brain ~ s(age, k=4) + sex + Overall_Accuracy + averageManualRating + roi, data=r2star, REML=TRUE, random=~(1|bblid))
Results
summary(iron_accuracy_GAMM$gam)
##
## Family: gaussian
## Link function: identity
##
## Formula:
## brain ~ s(age, k = 4) + sex + Overall_Accuracy + averageManualRating +
## roi
##
## Parametric coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 17.10228 0.32551 52.540 <2e-16 ***
## sex.L 0.03977 0.05106 0.779 0.436
## Overall_Accuracy 0.08633 0.04914 1.757 0.079 .
## averageManualRating -0.13376 0.16787 -0.797 0.426
## roicaudate -1.59725 0.07245 -22.047 <2e-16 ***
## roipallidum 5.30185 0.07242 73.207 <2e-16 ***
## roiputamen -0.10108 0.07236 -1.397 0.163
## ---
## 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.196 2.196 89.6 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## R-sq.(adj) = 0.669
## lmer.REML = 18378 Scale est. = 2.9424 n = 4493