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