ACURA Theta Analysis V1.6
Russell Chan
8th March 2019
##Set working directory
setwd("C:/Users/ChanRW/OneDrive - Universiteit Twente/2020_Alpha_SeqLearn/Current Analysis")##Packages
library(lme4)
library(effects)
library(lattice)
library(car)
library(ggplot2)
library(knitr)
library(reshape2)
library(dplyr)
library(forcats)
library(DHARMa)
library(Hmisc)
library(qpcR)
library(phia)
library(lsmeans)
library(emmeans)
library(multcomp)
library(plotly)
library(stargazer)
library(lmerTest)##Import dataset
d.MU<-read.table("Pow_Data_DEC18_V3.4.CSV", sep = ",", header = T, stringsAsFactors = F)#Factors
#Create Factors
d.MU$ID <- factor(d.MU$ID)
d.MU$Group <- factor(d.MU$Group)
d.MU$RoI <- factor(d.MU$RoI, levels=c('F', 'C', 'P'))
d.MU$Session <- factor(d.MU$Session, levels=c('1', '2', '3'))
d.MU$Timepoint <-as.numeric(d.MU$Timepoint, levels=c('1', '2', '3', '4', '5', '6', '7', '8', '9', '10'))
#d.MU$TimepointC <- factor(d.MU$TimepointC)
d.MU$muTheta <- as.numeric(d.MU$muTheta)
d.MU$muAlpha <- as.numeric(d.MU$muAlpha)#Model0: Theta Random Intercepts: ID. Anti-conservative.
m.T0 <- lmer(muTheta ~ Group * Session * Timepoint * RoI + (1 | ID), data=d.MU, REML = FALSE)
Anova(m.T0)## Analysis of Deviance Table (Type II Wald chisquare tests)
##
## Response: muTheta
## Chisq Df Pr(>Chisq)
## Group 12.8853 2 0.001592 **
## Session 186.0409 2 < 2.2e-16 ***
## Timepoint 348.8268 1 < 2.2e-16 ***
## RoI 7361.9089 2 < 2.2e-16 ***
## Group:Session 536.6922 4 < 2.2e-16 ***
## Group:Timepoint 37.1507 2 8.567e-09 ***
## Session:Timepoint 53.3599 2 2.588e-12 ***
## Group:RoI 163.7957 4 < 2.2e-16 ***
## Session:RoI 38.0504 4 1.094e-07 ***
## Timepoint:RoI 324.7699 2 < 2.2e-16 ***
## Group:Session:Timepoint 61.3865 4 1.483e-12 ***
## Group:Session:RoI 102.7651 8 < 2.2e-16 ***
## Group:Timepoint:RoI 62.1296 4 1.035e-12 ***
## Session:Timepoint:RoI 4.3648 4 0.358878
## Group:Session:Timepoint:RoI 49.5454 8 4.997e-08 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1m.T0test <- as(m.T0,"lmerModLmerTest")
summary(m.T0test, ddf="Satterthwaite")## Linear mixed model fit by maximum likelihood . t-tests use Satterthwaite's
## method [lmerModLmerTest]
## Formula: muTheta ~ Group * Session * Timepoint * RoI + (1 | ID)
## Data: d.MU
##
## AIC BIC logLik deviance df.resid
## 213031.6 213525.2 -106459.8 212919.6 49694
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -4.1636 -0.6243 -0.1100 0.4882 10.6939
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 0.7697 0.8773
## Residual 4.2146 2.0530
## Number of obs: 49750, groups: ID, 29
##
## Fixed effects:
## Estimate Std. Error df t value
## (Intercept) 5.011e+00 3.277e-01 3.590e+01 15.293
## GroupMED1 9.082e-01 4.504e-01 3.590e+01 2.017
## GroupMED21 1.027e+00 4.231e-01 3.591e+01 2.426
## Session2 4.908e-01 1.495e-01 4.972e+04 3.283
## Session3 3.045e-01 1.495e-01 4.972e+04 2.037
## Timepoint 3.092e-02 1.704e-02 4.972e+04 1.815
## RoIC -1.150e+00 1.579e-01 4.972e+04 -7.286
## RoIP -2.418e+00 1.585e-01 4.972e+04 -15.253
## GroupMED1:Session2 1.468e-01 2.055e-01 4.972e+04 0.714
## GroupMED21:Session2 5.191e-02 1.949e-01 4.972e+04 0.266
## GroupMED1:Session3 -5.052e-01 2.055e-01 4.972e+04 -2.459
## GroupMED21:Session3 -1.167e-01 1.932e-01 4.972e+04 -0.604
## GroupMED1:Timepoint -2.320e-02 2.341e-02 4.972e+04 -0.991
## GroupMED21:Timepoint -3.694e-02 2.201e-02 4.972e+04 -1.678
## Session2:Timepoint -5.919e-02 2.409e-02 4.972e+04 -2.457
## Session3:Timepoint -4.385e-02 2.409e-02 4.972e+04 -1.820
## GroupMED1:RoIC 7.697e-01 2.168e-01 4.972e+04 3.550
## GroupMED21:RoIC -1.225e-01 2.037e-01 4.972e+04 -0.601
## GroupMED1:RoIP -8.014e-02 2.176e-01 4.972e+04 -0.368
## GroupMED21:RoIP -4.490e-01 2.043e-01 4.972e+04 -2.198
## Session2:RoIC -7.359e-01 2.233e-01 4.972e+04 -3.296
## Session3:RoIC -2.915e-01 2.231e-01 4.972e+04 -1.307
## Session2:RoIP -8.278e-01 2.237e-01 4.972e+04 -3.700
## Session3:RoIP 1.914e-03 2.235e-01 4.972e+04 0.009
## Timepoint:RoIC -4.908e-03 2.545e-02 4.972e+04 -0.193
## Timepoint:RoIP 1.366e-02 2.555e-02 4.972e+04 0.535
## GroupMED1:Session2:Timepoint 5.922e-03 3.311e-02 4.972e+04 0.179
## GroupMED21:Session2:Timepoint 5.056e-02 3.140e-02 4.972e+04 1.610
## GroupMED1:Session3:Timepoint 7.851e-02 3.311e-02 4.972e+04 2.371
## GroupMED21:Session3:Timepoint 1.344e-01 3.113e-02 4.972e+04 4.319
## GroupMED1:Session2:RoIC 1.763e-01 3.067e-01 4.972e+04 0.575
## GroupMED21:Session2:RoIC 6.739e-01 2.907e-01 4.972e+04 2.318
## GroupMED1:Session3:RoIC -5.631e-01 3.066e-01 4.972e+04 -1.837
## GroupMED21:Session3:RoIC 4.420e-01 2.880e-01 4.972e+04 1.535
## GroupMED1:Session2:RoIP 5.222e-01 3.072e-01 4.972e+04 1.700
## GroupMED21:Session2:RoIP 1.264e+00 2.911e-01 4.972e+04 4.342
## GroupMED1:Session3:RoIP -3.724e-01 3.070e-01 4.972e+04 -1.213
## GroupMED21:Session3:RoIP 4.996e-01 2.884e-01 4.972e+04 1.732
## GroupMED1:Timepoint:RoIC -2.434e-02 3.494e-02 4.972e+04 -0.697
## GroupMED21:Timepoint:RoIC 7.236e-02 3.284e-02 4.972e+04 2.204
## GroupMED1:Timepoint:RoIP 1.736e-01 3.507e-02 4.972e+04 4.949
## GroupMED21:Timepoint:RoIP 1.470e-01 3.293e-02 4.972e+04 4.464
## Session2:Timepoint:RoIC 8.172e-02 3.599e-02 4.972e+04 2.271
## Session3:Timepoint:RoIC 2.867e-02 3.595e-02 4.972e+04 0.797
## Session2:Timepoint:RoIP 1.395e-01 3.606e-02 4.972e+04 3.870
## Session3:Timepoint:RoIP 3.218e-02 3.602e-02 4.972e+04 0.893
## GroupMED1:Session2:Timepoint:RoIC -3.056e-02 4.943e-02 4.972e+04 -0.618
## GroupMED21:Session2:Timepoint:RoIC -1.058e-01 4.685e-02 4.972e+04 -2.259
## GroupMED1:Session3:Timepoint:RoIC 4.921e-02 4.941e-02 4.972e+04 0.996
## GroupMED21:Session3:Timepoint:RoIC -8.364e-02 4.641e-02 4.972e+04 -1.802
## GroupMED1:Session2:Timepoint:RoIP -6.759e-02 4.951e-02 4.972e+04 -1.365
## GroupMED21:Session2:Timepoint:RoIP -2.274e-01 4.691e-02 4.972e+04 -4.847
## GroupMED1:Session3:Timepoint:RoIP -6.435e-02 4.948e-02 4.972e+04 -1.301
## GroupMED21:Session3:Timepoint:RoIP -4.281e-02 4.648e-02 4.972e+04 -0.921
## Pr(>|t|)
## (Intercept) < 2e-16 ***
## GroupMED1 0.051271 .
## GroupMED21 0.020406 *
## Session2 0.001028 **
## Session3 0.041678 *
## Timepoint 0.069599 .
## RoIC 3.25e-13 ***
## RoIP < 2e-16 ***
## GroupMED1:Session2 0.475063
## GroupMED21:Session2 0.789968
## GroupMED1:Session3 0.013949 *
## GroupMED21:Session3 0.545580
## GroupMED1:Timepoint 0.321831
## GroupMED21:Timepoint 0.093274 .
## Session2:Timepoint 0.014024 *
## Session3:Timepoint 0.068776 .
## GroupMED1:RoIC 0.000385 ***
## GroupMED21:RoIC 0.547738
## GroupMED1:RoIP 0.712656
## GroupMED21:RoIP 0.027984 *
## Session2:RoIC 0.000983 ***
## Session3:RoIC 0.191355
## Session2:RoIP 0.000216 ***
## Session3:RoIP 0.993168
## Timepoint:RoIC 0.847045
## Timepoint:RoIP 0.592938
## GroupMED1:Session2:Timepoint 0.858068
## GroupMED21:Session2:Timepoint 0.107353
## GroupMED1:Session3:Timepoint 0.017751 *
## GroupMED21:Session3:Timepoint 1.57e-05 ***
## GroupMED1:Session2:RoIC 0.565386
## GroupMED21:Session2:RoIC 0.020441 *
## GroupMED1:Session3:RoIC 0.066245 .
## GroupMED21:Session3:RoIC 0.124815
## GroupMED1:Session2:RoIP 0.089164 .
## GroupMED21:Session2:RoIP 1.42e-05 ***
## GroupMED1:Session3:RoIP 0.225118
## GroupMED21:Session3:RoIP 0.083236 .
## GroupMED1:Timepoint:RoIC 0.485995
## GroupMED21:Timepoint:RoIC 0.027548 *
## GroupMED1:Timepoint:RoIP 7.49e-07 ***
## GroupMED21:Timepoint:RoIP 8.06e-06 ***
## Session2:Timepoint:RoIC 0.023164 *
## Session3:Timepoint:RoIC 0.425215
## Session2:Timepoint:RoIP 0.000109 ***
## Session3:Timepoint:RoIP 0.371658
## GroupMED1:Session2:Timepoint:RoIC 0.536471
## GroupMED21:Session2:Timepoint:RoIC 0.023879 *
## GroupMED1:Session3:Timepoint:RoIC 0.319302
## GroupMED21:Session3:Timepoint:RoIC 0.071533 .
## GroupMED1:Session2:Timepoint:RoIP 0.172185
## GroupMED21:Session2:Timepoint:RoIP 1.26e-06 ***
## GroupMED1:Session3:Timepoint:RoIP 0.193415
## GroupMED21:Session3:Timepoint:RoIP 0.357000
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1##
## Correlation matrix not shown by default, as p = 54 > 12.
## Use print(x, correlation=TRUE) or
## vcov(x) if you need it#Model1: Theta Random Intercepts: ID and session(timepoint). Slightly more conservative, but this might be the most accurate for our design since we are interested in Sessions and Timepoints muTheta changes as an effect of the Group.
m.T1 <- lmer(muTheta ~ Group * Session * Timepoint * RoI + (1 | ID) + (1 | Session:Timepoint), data=d.MU, REML = FALSE)
Anova(m.T1)## Analysis of Deviance Table (Type II Wald chisquare tests)
##
## Response: muTheta
## Chisq Df Pr(>Chisq)
## Group 12.8714 2 0.0016033 **
## Session 16.6476 2 0.0002427 ***
## Timepoint 30.8614 1 2.771e-08 ***
## RoI 7406.0471 2 < 2.2e-16 ***
## Group:Session 539.9095 4 < 2.2e-16 ***
## Group:Timepoint 37.3734 2 7.664e-09 ***
## Session:Timepoint 4.7705 2 0.0920669 .
## Group:RoI 164.7777 4 < 2.2e-16 ***
## Session:RoI 38.2785 4 9.817e-08 ***
## Timepoint:RoI 326.7170 2 < 2.2e-16 ***
## Group:Session:Timepoint 61.7546 4 1.241e-12 ***
## Group:Session:RoI 103.3812 8 < 2.2e-16 ***
## Group:Timepoint:RoI 62.5021 4 8.637e-13 ***
## Session:Timepoint:RoI 4.3910 4 0.3556695
## Group:Session:Timepoint:RoI 49.8425 8 4.382e-08 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1m.T1test <- as(m.T1,"lmerModLmerTest")
summary(m.T1test, ddf="Satterthwaite")## Linear mixed model fit by maximum likelihood . t-tests use Satterthwaite's
## method [lmerModLmerTest]
## Formula: muTheta ~ Group * Session * Timepoint * RoI + (1 | ID) + (1 |
## Session:Timepoint)
## Data: d.MU
##
## AIC BIC logLik deviance df.resid
## 212806.7 213309.1 -106346.3 212692.7 49693
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -4.1736 -0.6240 -0.1108 0.4866 10.8191
##
## Random effects:
## Groups Name Variance Std.Dev.
## Session:Timepoint (Intercept) 0.02593 0.1610
## ID (Intercept) 0.77053 0.8778
## Residual 4.18952 2.0468
## Number of obs: 49750, groups: Session:Timepoint, 30; ID, 29
##
## Fixed effects:
## Estimate Std. Error df t value
## (Intercept) 5.011e+00 3.457e-01 4.354e+01 14.495
## GroupMED1 9.082e-01 4.505e-01 3.585e+01 2.016
## GroupMED21 1.027e+00 4.232e-01 3.586e+01 2.426
## Session2 4.908e-01 2.154e-01 8.876e+01 2.278
## Session3 3.045e-01 2.154e-01 8.876e+01 1.413
## Timepoint 3.092e-02 2.455e-02 8.876e+01 1.259
## RoIC -1.150e+00 1.574e-01 4.969e+04 -7.307
## RoIP -2.418e+00 1.580e-01 4.969e+04 -15.298
## GroupMED1:Session2 1.468e-01 2.049e-01 4.969e+04 0.716
## GroupMED21:Session2 5.191e-02 1.943e-01 4.969e+04 0.267
## GroupMED1:Session3 -5.052e-01 2.049e-01 4.969e+04 -2.466
## GroupMED21:Session3 -1.167e-01 1.926e-01 4.969e+04 -0.606
## GroupMED1:Timepoint -2.320e-02 2.335e-02 4.969e+04 -0.994
## GroupMED21:Timepoint -3.694e-02 2.195e-02 4.969e+04 -1.683
## Session2:Timepoint -5.919e-02 3.472e-02 8.876e+01 -1.705
## Session3:Timepoint -4.385e-02 3.472e-02 8.876e+01 -1.263
## GroupMED1:RoIC 7.697e-01 2.162e-01 4.969e+04 3.561
## GroupMED21:RoIC -1.225e-01 2.031e-01 4.969e+04 -0.603
## GroupMED1:RoIP -8.014e-02 2.170e-01 4.969e+04 -0.369
## GroupMED21:RoIP -4.490e-01 2.037e-01 4.969e+04 -2.204
## Session2:RoIC -7.359e-01 2.226e-01 4.969e+04 -3.306
## Session3:RoIC -2.915e-01 2.224e-01 4.969e+04 -1.311
## Session2:RoIP -8.278e-01 2.231e-01 4.969e+04 -3.711
## Session3:RoIP 1.914e-03 2.228e-01 4.969e+04 0.009
## Timepoint:RoIC -4.908e-03 2.537e-02 4.969e+04 -0.193
## Timepoint:RoIP 1.366e-02 2.547e-02 4.969e+04 0.536
## GroupMED1:Session2:Timepoint 5.922e-03 3.302e-02 4.969e+04 0.179
## GroupMED21:Session2:Timepoint 5.056e-02 3.130e-02 4.969e+04 1.615
## GroupMED1:Session3:Timepoint 7.851e-02 3.302e-02 4.969e+04 2.378
## GroupMED21:Session3:Timepoint 1.344e-01 3.104e-02 4.969e+04 4.332
## GroupMED1:Session2:RoIC 1.763e-01 3.058e-01 4.969e+04 0.577
## GroupMED21:Session2:RoIC 6.739e-01 2.898e-01 4.969e+04 2.325
## GroupMED1:Session3:RoIC -5.631e-01 3.057e-01 4.969e+04 -1.842
## GroupMED21:Session3:RoIC 4.420e-01 2.871e-01 4.969e+04 1.539
## GroupMED1:Session2:RoIP 5.222e-01 3.063e-01 4.969e+04 1.705
## GroupMED21:Session2:RoIP 1.264e+00 2.902e-01 4.969e+04 4.355
## GroupMED1:Session3:RoIP -3.724e-01 3.061e-01 4.969e+04 -1.217
## GroupMED21:Session3:RoIP 4.996e-01 2.875e-01 4.969e+04 1.737
## GroupMED1:Timepoint:RoIC -2.434e-02 3.484e-02 4.969e+04 -0.699
## GroupMED21:Timepoint:RoIC 7.236e-02 3.274e-02 4.969e+04 2.210
## GroupMED1:Timepoint:RoIP 1.736e-01 3.497e-02 4.969e+04 4.964
## GroupMED21:Timepoint:RoIP 1.470e-01 3.283e-02 4.969e+04 4.477
## Session2:Timepoint:RoIC 8.172e-02 3.588e-02 4.969e+04 2.278
## Session3:Timepoint:RoIC 2.867e-02 3.585e-02 4.969e+04 0.800
## Session2:Timepoint:RoIP 1.395e-01 3.595e-02 4.969e+04 3.881
## Session3:Timepoint:RoIP 3.218e-02 3.592e-02 4.969e+04 0.896
## GroupMED1:Session2:Timepoint:RoIC -3.056e-02 4.929e-02 4.969e+04 -0.620
## GroupMED21:Session2:Timepoint:RoIC -1.058e-01 4.671e-02 4.969e+04 -2.266
## GroupMED1:Session3:Timepoint:RoIC 4.921e-02 4.926e-02 4.969e+04 0.999
## GroupMED21:Session3:Timepoint:RoIC -8.364e-02 4.627e-02 4.969e+04 -1.808
## GroupMED1:Session2:Timepoint:RoIP -6.759e-02 4.936e-02 4.969e+04 -1.369
## GroupMED21:Session2:Timepoint:RoIP -2.274e-01 4.677e-02 4.969e+04 -4.862
## GroupMED1:Session3:Timepoint:RoIP -6.435e-02 4.933e-02 4.969e+04 -1.304
## GroupMED21:Session3:Timepoint:RoIP -4.281e-02 4.634e-02 4.969e+04 -0.924
## Pr(>|t|)
## (Intercept) < 2e-16 ***
## GroupMED1 0.051323 .
## GroupMED21 0.020436 *
## Session2 0.025121 *
## Session3 0.161060
## Timepoint 0.211289
## RoIC 2.76e-13 ***
## RoIP < 2e-16 ***
## GroupMED1:Session2 0.473743
## GroupMED21:Session2 0.789348
## GroupMED1:Session3 0.013666 *
## GroupMED21:Session3 0.544378
## GroupMED1:Timepoint 0.320385
## GroupMED21:Timepoint 0.092298 .
## Session2:Timepoint 0.091743 .
## Session3:Timepoint 0.209952
## GroupMED1:RoIC 0.000370 ***
## GroupMED21:RoIC 0.546541
## GroupMED1:RoIP 0.711833
## GroupMED21:RoIP 0.027518 *
## Session2:RoIC 0.000949 ***
## Session3:RoIC 0.190029
## Session2:RoIP 0.000207 ***
## Session3:RoIP 0.993148
## Timepoint:RoIC 0.846593
## Timepoint:RoIP 0.591832
## GroupMED1:Session2:Timepoint 0.857648
## GroupMED21:Session2:Timepoint 0.106305
## GroupMED1:Session3:Timepoint 0.017413 *
## GroupMED21:Session3:Timepoint 1.48e-05 ***
## GroupMED1:Session2:RoIC 0.564224
## GroupMED21:Session2:RoIC 0.020067 *
## GroupMED1:Session3:RoIC 0.065437 .
## GroupMED21:Session3:RoIC 0.123690
## GroupMED1:Session2:RoIP 0.088211 .
## GroupMED21:Session2:RoIP 1.33e-05 ***
## GroupMED1:Session3:RoIP 0.223733
## GroupMED21:Session3:RoIP 0.082317 .
## GroupMED1:Timepoint:RoIC 0.484690
## GroupMED21:Timepoint:RoIC 0.027087 *
## GroupMED1:Timepoint:RoIP 6.94e-07 ***
## GroupMED21:Timepoint:RoIP 7.57e-06 ***
## Session2:Timepoint:RoIC 0.022755 *
## Session3:Timepoint:RoIC 0.423830
## Session2:Timepoint:RoIP 0.000104 ***
## Session3:Timepoint:RoIP 0.370228
## GroupMED1:Session2:Timepoint:RoIC 0.535252
## GroupMED21:Session2:Timepoint:RoIC 0.023461 *
## GroupMED1:Session3:Timepoint:RoIC 0.317856
## GroupMED21:Session3:Timepoint:RoIC 0.070689 .
## GroupMED1:Session2:Timepoint:RoIP 0.170905
## GroupMED21:Session2:Timepoint:RoIP 1.17e-06 ***
## GroupMED1:Session3:Timepoint:RoIP 0.192085
## GroupMED21:Session3:Timepoint:RoIP 0.355563
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1##
## Correlation matrix not shown by default, as p = 54 > 12.
## Use print(x, correlation=TRUE) or
## vcov(x) if you need it#Model2: Theta Random Slope: Session intercept:ID. Timepoint tracks the average change within the session. #Converged.
m.T2 <- lmer(muTheta ~ Group * Session * Timepoint * RoI + (1 + Session | ID), data=d.MU, REML= FALSE)
Anova(m.T2)## Analysis of Deviance Table (Type II Wald chisquare tests)
##
## Response: muTheta
## Chisq Df Pr(>Chisq)
## Group 13.2378 2 0.001335 **
## Session 2.6648 2 0.263843
## Timepoint 377.1543 1 < 2.2e-16 ***
## RoI 7950.2238 2 < 2.2e-16 ***
## Group:Session 8.0451 4 0.089940 .
## Group:Timepoint 40.1676 2 1.895e-09 ***
## Session:Timepoint 57.6932 2 2.965e-13 ***
## Group:RoI 177.2142 4 < 2.2e-16 ***
## Session:RoI 40.7356 4 3.049e-08 ***
## Timepoint:RoI 351.1437 2 < 2.2e-16 ***
## Group:Session:Timepoint 66.3716 4 1.323e-13 ***
## Group:Session:RoI 112.0423 8 < 2.2e-16 ***
## Group:Timepoint:RoI 67.1750 4 8.955e-14 ***
## Session:Timepoint:RoI 4.7193 4 0.317332
## Group:Session:Timepoint:RoI 53.5689 8 8.368e-09 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1m.T2test <- as(m.T2,"lmerModLmerTest")
summary(m.T2test, ddf="Satterthwaite")## Linear mixed model fit by maximum likelihood . t-tests use Satterthwaite's
## method [lmerModLmerTest]
## Formula: muTheta ~ Group * Session * Timepoint * RoI + (1 + Session | ID)
## Data: d.MU
##
## AIC BIC logLik deviance df.resid
## 209397.3 209935.0 -104637.7 209275.3 49689
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -4.7465 -0.6090 -0.1027 0.4825 11.3491
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## ID (Intercept) 0.6824 0.8261
## Session2 0.9483 0.9738 -0.28
## Session3 0.9783 0.9891 -0.16 0.52
## Residual 3.8981 1.9744
## Number of obs: 49750, groups: ID, 29
##
## Fixed effects:
## Estimate Std. Error df t value
## (Intercept) 5.011e+00 3.093e-01 3.574e+01 16.205
## GroupMED1 9.082e-01 4.250e-01 3.574e+01 2.137
## GroupMED21 1.028e+00 3.993e-01 3.575e+01 2.575
## Session2 4.908e-01 3.731e-01 3.822e+01 1.315
## Session3 3.045e-01 3.781e-01 3.856e+01 0.805
## Timepoint 3.091e-02 1.638e-02 4.966e+04 1.887
## RoIC -1.149e+00 1.518e-01 4.966e+04 -7.568
## RoIP -2.415e+00 1.524e-01 4.966e+04 -15.841
## GroupMED1:Session2 1.468e-01 5.128e-01 3.822e+01 0.286
## GroupMED21:Session2 3.968e-02 4.875e-01 3.870e+01 0.081
## GroupMED1:Session3 -5.052e-01 5.197e-01 3.856e+01 -0.972
## GroupMED21:Session3 -1.193e-01 4.882e-01 3.858e+01 -0.244
## GroupMED1:Timepoint -2.320e-02 2.252e-02 4.966e+04 -1.030
## GroupMED21:Timepoint -3.694e-02 2.117e-02 4.966e+04 -1.745
## Session2:Timepoint -5.919e-02 2.317e-02 4.966e+04 -2.555
## Session3:Timepoint -4.385e-02 2.317e-02 4.966e+04 -1.892
## GroupMED1:RoIC 7.686e-01 2.085e-01 4.966e+04 3.686
## GroupMED21:RoIC -1.251e-01 1.959e-01 4.966e+04 -0.639
## GroupMED1:RoIP -7.485e-02 2.093e-01 4.966e+04 -0.358
## GroupMED21:RoIP -4.514e-01 1.965e-01 4.966e+04 -2.297
## Session2:RoIC -7.378e-01 2.147e-01 4.966e+04 -3.436
## Session3:RoIC -2.926e-01 2.145e-01 4.966e+04 -1.364
## Session2:RoIP -8.310e-01 2.152e-01 4.966e+04 -3.862
## Session3:RoIP -8.205e-04 2.150e-01 4.966e+04 -0.004
## Timepoint:RoIC -4.908e-03 2.447e-02 4.966e+04 -0.201
## Timepoint:RoIP 1.366e-02 2.457e-02 4.966e+04 0.556
## GroupMED1:Session2:Timepoint 5.922e-03 3.185e-02 4.966e+04 0.186
## GroupMED21:Session2:Timepoint 5.056e-02 3.020e-02 4.966e+04 1.674
## GroupMED1:Session3:Timepoint 7.851e-02 3.185e-02 4.966e+04 2.465
## GroupMED21:Session3:Timepoint 1.344e-01 2.994e-02 4.966e+04 4.491
## GroupMED1:Session2:RoIC 1.731e-01 2.950e-01 4.966e+04 0.587
## GroupMED21:Session2:RoIC 6.773e-01 2.796e-01 4.966e+04 2.423
## GroupMED1:Session3:RoIC -5.625e-01 2.948e-01 4.966e+04 -1.908
## GroupMED21:Session3:RoIC 4.457e-01 2.770e-01 4.966e+04 1.609
## GroupMED1:Session2:RoIP 5.174e-01 2.954e-01 4.966e+04 1.751
## GroupMED21:Session2:RoIP 1.267e+00 2.800e-01 4.966e+04 4.525
## GroupMED1:Session3:RoIP -3.777e-01 2.953e-01 4.966e+04 -1.279
## GroupMED21:Session3:RoIP 5.030e-01 2.774e-01 4.966e+04 1.813
## GroupMED1:Timepoint:RoIC -2.434e-02 3.360e-02 4.966e+04 -0.724
## GroupMED21:Timepoint:RoIC 7.236e-02 3.158e-02 4.966e+04 2.291
## GroupMED1:Timepoint:RoIP 1.736e-01 3.373e-02 4.966e+04 5.146
## GroupMED21:Timepoint:RoIP 1.470e-01 3.167e-02 4.966e+04 4.642
## Session2:Timepoint:RoIC 8.172e-02 3.461e-02 4.966e+04 2.361
## Session3:Timepoint:RoIC 2.867e-02 3.458e-02 4.966e+04 0.829
## Session2:Timepoint:RoIP 1.395e-01 3.468e-02 4.966e+04 4.024
## Session3:Timepoint:RoIP 3.218e-02 3.464e-02 4.966e+04 0.929
## GroupMED1:Session2:Timepoint:RoIC -3.056e-02 4.754e-02 4.966e+04 -0.643
## GroupMED21:Session2:Timepoint:RoIC -1.058e-01 4.506e-02 4.966e+04 -2.349
## GroupMED1:Session3:Timepoint:RoIC 4.921e-02 4.752e-02 4.966e+04 1.036
## GroupMED21:Session3:Timepoint:RoIC -8.364e-02 4.463e-02 4.966e+04 -1.874
## GroupMED1:Session2:Timepoint:RoIP -6.759e-02 4.761e-02 4.966e+04 -1.420
## GroupMED21:Session2:Timepoint:RoIP -2.274e-01 4.512e-02 4.966e+04 -5.040
## GroupMED1:Session3:Timepoint:RoIP -6.435e-02 4.759e-02 4.966e+04 -1.352
## GroupMED21:Session3:Timepoint:RoIP -4.281e-02 4.470e-02 4.966e+04 -0.958
## Pr(>|t|)
## (Intercept) < 2e-16 ***
## GroupMED1 0.039535 *
## GroupMED21 0.014320 *
## Session2 0.196187
## Session3 0.425568
## Timepoint 0.059194 .
## RoIC 3.84e-14 ***
## RoIP < 2e-16 ***
## GroupMED1:Session2 0.776271
## GroupMED21:Session2 0.935557
## GroupMED1:Session3 0.337038
## GroupMED21:Session3 0.808253
## GroupMED1:Timepoint 0.302943
## GroupMED21:Timepoint 0.080952 .
## Session2:Timepoint 0.010636 *
## Session3:Timepoint 0.058446 .
## GroupMED1:RoIC 0.000228 ***
## GroupMED21:RoIC 0.523081
## GroupMED1:RoIP 0.720617
## GroupMED21:RoIP 0.021613 *
## Session2:RoIC 0.000592 ***
## Session3:RoIC 0.172653
## Session2:RoIP 0.000113 ***
## Session3:RoIP 0.996955
## Timepoint:RoIC 0.841035
## Timepoint:RoIP 0.578303
## GroupMED1:Session2:Timepoint 0.852481
## GroupMED21:Session2:Timepoint 0.094070 .
## GroupMED1:Session3:Timepoint 0.013696 *
## GroupMED21:Session3:Timepoint 7.12e-06 ***
## GroupMED1:Session2:RoIC 0.557259
## GroupMED21:Session2:RoIC 0.015405 *
## GroupMED1:Session3:RoIC 0.056396 .
## GroupMED21:Session3:RoIC 0.107584
## GroupMED1:Session2:RoIP 0.079890 .
## GroupMED21:Session2:RoIP 6.05e-06 ***
## GroupMED1:Session3:RoIP 0.200810
## GroupMED21:Session3:RoIP 0.069775 .
## GroupMED1:Timepoint:RoIC 0.468802
## GroupMED21:Timepoint:RoIC 0.021941 *
## GroupMED1:Timepoint:RoIP 2.67e-07 ***
## GroupMED21:Timepoint:RoIP 3.46e-06 ***
## Session2:Timepoint:RoIC 0.018220 *
## Session3:Timepoint:RoIC 0.407018
## Session2:Timepoint:RoIP 5.74e-05 ***
## Session3:Timepoint:RoIP 0.352921
## GroupMED1:Session2:Timepoint:RoIC 0.520374
## GroupMED21:Session2:Timepoint:RoIC 0.018823 *
## GroupMED1:Session3:Timepoint:RoIC 0.300417
## GroupMED21:Session3:Timepoint:RoIC 0.060955 .
## GroupMED1:Session2:Timepoint:RoIP 0.155734
## GroupMED21:Session2:Timepoint:RoIP 4.67e-07 ***
## GroupMED1:Session3:Timepoint:RoIP 0.176273
## GroupMED21:Session3:Timepoint:RoIP 0.338181
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1##
## Correlation matrix not shown by default, as p = 54 > 12.
## Use print(x, correlation=TRUE) or
## vcov(x) if you need it#Model2.5: Theta Random Slope:Timepoint intercept:ID. Session tracks the average change within the Timepoints. #Converged.
m.T2.5 <- lmer(muTheta ~ Group * Session * Timepoint * RoI + (1 + Timepoint | ID), data=d.MU, REML= FALSE)
Anova(m.T2.5)## Analysis of Deviance Table (Type II Wald chisquare tests)
##
## Response: muTheta
## Chisq Df Pr(>Chisq)
## Group 12.8792 2 0.001597 **
## Session 190.0786 2 < 2.2e-16 ***
## Timepoint 9.5311 1 0.002020 **
## RoI 7521.6245 2 < 2.2e-16 ***
## Group:Session 548.3342 4 < 2.2e-16 ***
## Group:Timepoint 1.0036 2 0.605440
## Session:Timepoint 53.7373 2 2.143e-12 ***
## Group:RoI 167.3493 4 < 2.2e-16 ***
## Session:RoI 38.8758 4 7.390e-08 ***
## Timepoint:RoI 332.4108 2 < 2.2e-16 ***
## Group:Session:Timepoint 62.8444 4 7.317e-13 ***
## Group:Session:RoI 104.9945 8 < 2.2e-16 ***
## Group:Timepoint:RoI 63.9815 4 4.217e-13 ***
## Session:Timepoint:RoI 4.4431 4 0.349353
## Group:Session:Timepoint:RoI 50.5542 8 3.197e-08 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1m.T2.5test <- as(m.T2.5,"lmerModLmerTest")
summary(m.T2.5test, ddf="Satterthwaite")## Linear mixed model fit by maximum likelihood . t-tests use Satterthwaite's
## method [lmerModLmerTest]
## Formula: muTheta ~ Group * Session * Timepoint * RoI + (1 + Timepoint |
## ID)
## Data: d.MU
##
## AIC BIC logLik deviance df.resid
## 212073.7 212584.9 -105978.8 211957.7 49692
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -4.1785 -0.6213 -0.1079 0.4842 10.5624
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## ID (Intercept) 1.09283 1.0454
## Timepoint 0.01071 0.1035 -0.54
## Residual 4.12514 2.0310
## Number of obs: 49750, groups: ID, 29
##
## Fixed effects:
## Estimate Std. Error df t value
## (Intercept) 5.011e+00 3.841e-01 3.315e+01 13.047
## GroupMED1 9.082e-01 5.279e-01 3.315e+01 1.720
## GroupMED21 1.028e+00 4.959e-01 3.316e+01 2.073
## Session2 4.908e-01 1.479e-01 4.969e+04 3.318
## Session3 3.045e-01 1.479e-01 4.969e+04 2.059
## Timepoint 3.091e-02 4.029e-02 4.041e+01 0.767
## RoIC -1.151e+00 1.562e-01 4.969e+04 -7.370
## RoIP -2.415e+00 1.568e-01 4.969e+04 -15.398
## GroupMED1:Session2 1.468e-01 2.033e-01 4.969e+04 0.722
## GroupMED21:Session2 5.584e-02 1.931e-01 4.970e+04 0.289
## GroupMED1:Session3 -5.052e-01 2.033e-01 4.969e+04 -2.485
## GroupMED21:Session3 -1.167e-01 1.911e-01 4.969e+04 -0.611
## GroupMED1:Timepoint -2.320e-02 5.537e-02 4.041e+01 -0.419
## GroupMED21:Timepoint -3.721e-02 5.202e-02 4.043e+01 -0.715
## Session2:Timepoint -5.919e-02 2.384e-02 4.969e+04 -2.483
## Session3:Timepoint -4.385e-02 2.384e-02 4.969e+04 -1.840
## GroupMED1:RoIC 7.707e-01 2.145e-01 4.969e+04 3.593
## GroupMED21:RoIC -1.230e-01 2.016e-01 4.969e+04 -0.610
## GroupMED1:RoIP -9.035e-02 2.153e-01 4.969e+04 -0.420
## GroupMED21:RoIP -4.516e-01 2.022e-01 4.969e+04 -2.234
## Session2:RoIC -7.337e-01 2.209e-01 4.969e+04 -3.321
## Session3:RoIC -2.905e-01 2.207e-01 4.969e+04 -1.316
## Session2:RoIP -8.316e-01 2.214e-01 4.969e+04 -3.757
## Session3:RoIP -9.644e-04 2.211e-01 4.969e+04 -0.004
## Timepoint:RoIC -4.732e-03 2.517e-02 4.969e+04 -0.188
## Timepoint:RoIP 1.313e-02 2.528e-02 4.969e+04 0.520
## GroupMED1:Session2:Timepoint 5.922e-03 3.276e-02 4.969e+04 0.181
## GroupMED21:Session2:Timepoint 4.985e-02 3.111e-02 4.970e+04 1.602
## GroupMED1:Session3:Timepoint 7.851e-02 3.276e-02 4.969e+04 2.396
## GroupMED21:Session3:Timepoint 1.344e-01 3.080e-02 4.969e+04 4.365
## GroupMED1:Session2:RoIC 1.694e-01 3.035e-01 4.969e+04 0.558
## GroupMED21:Session2:RoIC 6.731e-01 2.876e-01 4.969e+04 2.341
## GroupMED1:Session3:RoIC -5.602e-01 3.033e-01 4.969e+04 -1.847
## GroupMED21:Session3:RoIC 4.410e-01 2.849e-01 4.969e+04 1.548
## GroupMED1:Session2:RoIP 5.333e-01 3.039e-01 4.969e+04 1.755
## GroupMED21:Session2:RoIP 1.267e+00 2.880e-01 4.969e+04 4.401
## GroupMED1:Session3:RoIP -3.622e-01 3.037e-01 4.969e+04 -1.192
## GroupMED21:Session3:RoIP 5.006e-01 2.853e-01 4.969e+04 1.755
## GroupMED1:Timepoint:RoIC -2.452e-02 3.457e-02 4.969e+04 -0.709
## GroupMED21:Timepoint:RoIC 7.246e-02 3.249e-02 4.969e+04 2.230
## GroupMED1:Timepoint:RoIP 1.754e-01 3.470e-02 4.969e+04 5.055
## GroupMED21:Timepoint:RoIP 1.475e-01 3.258e-02 4.969e+04 4.526
## Session2:Timepoint:RoIC 8.131e-02 3.560e-02 4.969e+04 2.284
## Session3:Timepoint:RoIC 2.849e-02 3.557e-02 4.969e+04 0.801
## Session2:Timepoint:RoIP 1.402e-01 3.567e-02 4.969e+04 3.931
## Session3:Timepoint:RoIP 3.270e-02 3.564e-02 4.969e+04 0.918
## GroupMED1:Session2:Timepoint:RoIC -2.929e-02 4.891e-02 4.969e+04 -0.599
## GroupMED21:Session2:Timepoint:RoIC -1.057e-01 4.635e-02 4.969e+04 -2.281
## GroupMED1:Session3:Timepoint:RoIC 4.868e-02 4.888e-02 4.969e+04 0.996
## GroupMED21:Session3:Timepoint:RoIC -8.346e-02 4.592e-02 4.969e+04 -1.818
## GroupMED1:Session2:Timepoint:RoIP -6.962e-02 4.898e-02 4.969e+04 -1.421
## GroupMED21:Session2:Timepoint:RoIP -2.280e-01 4.642e-02 4.969e+04 -4.913
## GroupMED1:Session3:Timepoint:RoIP -6.621e-02 4.895e-02 4.969e+04 -1.352
## GroupMED21:Session3:Timepoint:RoIP -4.301e-02 4.598e-02 4.969e+04 -0.935
## Pr(>|t|)
## (Intercept) 1.31e-14 ***
## GroupMED1 0.094687 .
## GroupMED21 0.046024 *
## Session2 0.000906 ***
## Session3 0.039524 *
## Timepoint 0.447316
## RoIC 1.73e-13 ***
## RoIP < 2e-16 ***
## GroupMED1:Session2 0.470312
## GroupMED21:Session2 0.772420
## GroupMED1:Session3 0.012952 *
## GroupMED21:Session3 0.541254
## GroupMED1:Timepoint 0.677452
## GroupMED21:Timepoint 0.478461
## Session2:Timepoint 0.013023 *
## Session3:Timepoint 0.065839 .
## GroupMED1:RoIC 0.000327 ***
## GroupMED21:RoIC 0.541748
## GroupMED1:RoIP 0.674728
## GroupMED21:RoIP 0.025501 *
## Session2:RoIC 0.000897 ***
## Session3:RoIC 0.188078
## Session2:RoIP 0.000172 ***
## Session3:RoIP 0.996520
## Timepoint:RoIC 0.850891
## Timepoint:RoIP 0.603342
## GroupMED1:Session2:Timepoint 0.856553
## GroupMED21:Session2:Timepoint 0.109155
## GroupMED1:Session3:Timepoint 0.016560 *
## GroupMED21:Session3:Timepoint 1.27e-05 ***
## GroupMED1:Session2:RoIC 0.576796
## GroupMED21:Session2:RoIC 0.019258 *
## GroupMED1:Session3:RoIC 0.064731 .
## GroupMED21:Session3:RoIC 0.121615
## GroupMED1:Session2:RoIP 0.079277 .
## GroupMED21:Session2:RoIP 1.08e-05 ***
## GroupMED1:Session3:RoIP 0.233071
## GroupMED21:Session3:RoIP 0.079327 .
## GroupMED1:Timepoint:RoIC 0.478141
## GroupMED21:Timepoint:RoIC 0.025729 *
## GroupMED1:Timepoint:RoIP 4.31e-07 ***
## GroupMED21:Timepoint:RoIP 6.02e-06 ***
## Session2:Timepoint:RoIC 0.022386 *
## Session3:Timepoint:RoIC 0.423097
## Session2:Timepoint:RoIP 8.48e-05 ***
## Session3:Timepoint:RoIP 0.358801
## GroupMED1:Session2:Timepoint:RoIC 0.549245
## GroupMED21:Session2:Timepoint:RoIC 0.022581 *
## GroupMED1:Session3:Timepoint:RoIC 0.319268
## GroupMED21:Session3:Timepoint:RoIC 0.069112 .
## GroupMED1:Session2:Timepoint:RoIP 0.155204
## GroupMED21:Session2:Timepoint:RoIP 8.99e-07 ***
## GroupMED1:Session3:Timepoint:RoIP 0.176231
## GroupMED21:Session3:Timepoint:RoIP 0.349669
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1##
## Correlation matrix not shown by default, as p = 54 > 12.
## Use print(x, correlation=TRUE) or
## vcov(x) if you need it#Model3: Theta Random Slope: Timepoint+Session, intercept:ID. Session tracks the overall differences between sessions.
m.T3<-lmer(muTheta ~ Group * Session * Timepoint * RoI + (1 + Timepoint+Session | ID), data=d.MU, REML= FALSE)## Warning in checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
## Model failed to converge with max|grad| = 0.0105483 (tol = 0.002, component 1)Anova(m.T3)## Analysis of Deviance Table (Type II Wald chisquare tests)
##
## Response: muTheta
## Chisq Df Pr(>Chisq)
## Group 13.3204 2 0.001281 **
## Session 2.6544 2 0.265216
## Timepoint 9.5215 1 0.002031 **
## RoI 8137.6577 2 < 2.2e-16 ***
## Group:Session 8.0479 4 0.089838 .
## Group:Timepoint 1.0026 2 0.605739
## Session:Timepoint 58.1991 2 2.303e-13 ***
## Group:RoI 181.3913 4 < 2.2e-16 ***
## Session:RoI 41.6963 4 1.928e-08 ***
## Timepoint:RoI 360.0707 2 < 2.2e-16 ***
## Group:Session:Timepoint 68.0670 4 5.807e-14 ***
## Group:Session:RoI 114.6848 8 < 2.2e-16 ***
## Group:Timepoint:RoI 69.3073 4 3.178e-14 ***
## Session:Timepoint:RoI 4.8129 4 0.307037
## Group:Session:Timepoint:RoI 54.7607 8 4.914e-09 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1m.T3test <- as(m.T3,"lmerModLmerTest")
summary(m.T3test, ddf="Satterthwaite")## Linear mixed model fit by maximum likelihood . t-tests use Satterthwaite's
## method [lmerModLmerTest]
## Formula: muTheta ~ Group * Session * Timepoint * RoI + (1 + Timepoint +
## Session | ID)
## Data: d.MU
##
## AIC BIC logLik deviance df.resid
## 208354.3 208927.3 -104112.2 208224.3 49685
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -4.8136 -0.6101 -0.1007 0.4852 11.2245
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## ID (Intercept) 0.98744 0.9937
## Timepoint 0.01075 0.1037 -0.56
## Session2 0.94871 0.9740 -0.16 -0.13
## Session3 0.97960 0.9897 -0.18 0.08 0.52
## Residual 3.80830 1.9515
## Number of obs: 49750, groups: ID, 29
##
## Fixed effects:
## Estimate Std. Error df t value
## (Intercept) 5.011e+00 3.654e-01 3.289e+01 13.714
## GroupMED1 9.082e-01 5.022e-01 3.289e+01 1.808
## GroupMED21 1.030e+00 4.718e-01 3.290e+01 2.182
## Session2 4.908e-01 3.725e-01 3.797e+01 1.317
## Session3 3.045e-01 3.777e-01 3.824e+01 0.806
## Timepoint 3.092e-02 4.007e-02 3.940e+01 0.772
## RoIC -1.150e+00 1.501e-01 4.964e+04 -7.664
## RoIP -2.412e+00 1.507e-01 4.964e+04 -16.007
## GroupMED1:Session2 1.468e-01 5.120e-01 3.797e+01 0.287
## GroupMED21:Session2 4.546e-02 4.867e-01 3.852e+01 0.093
## GroupMED1:Session3 -5.052e-01 5.191e-01 3.824e+01 -0.973
## GroupMED21:Session3 -1.193e-01 4.876e-01 3.826e+01 -0.245
## GroupMED1:Timepoint -2.320e-02 5.507e-02 3.940e+01 -0.421
## GroupMED21:Timepoint -3.721e-02 5.173e-02 3.942e+01 -0.719
## Session2:Timepoint -5.919e-02 2.290e-02 4.964e+04 -2.584
## Session3:Timepoint -4.385e-02 2.290e-02 4.964e+04 -1.915
## GroupMED1:RoIC 7.696e-01 2.061e-01 4.964e+04 3.734
## GroupMED21:RoIC -1.256e-01 1.937e-01 4.964e+04 -0.649
## GroupMED1:RoIP -8.512e-02 2.069e-01 4.964e+04 -0.411
## GroupMED21:RoIP -4.540e-01 1.942e-01 4.964e+04 -2.337
## Session2:RoIC -7.355e-01 2.123e-01 4.964e+04 -3.465
## Session3:RoIC -2.916e-01 2.121e-01 4.964e+04 -1.375
## Session2:RoIP -8.349e-01 2.127e-01 4.964e+04 -3.925
## Session3:RoIP -3.715e-03 2.125e-01 4.964e+04 -0.017
## Timepoint:RoIC -4.732e-03 2.419e-02 4.964e+04 -0.196
## Timepoint:RoIP 1.313e-02 2.429e-02 4.964e+04 0.541
## GroupMED1:Session2:Timepoint 5.922e-03 3.148e-02 4.964e+04 0.188
## GroupMED21:Session2:Timepoint 4.985e-02 2.990e-02 4.965e+04 1.668
## GroupMED1:Session3:Timepoint 7.851e-02 3.148e-02 4.964e+04 2.494
## GroupMED21:Session3:Timepoint 1.344e-01 2.959e-02 4.964e+04 4.543
## GroupMED1:Session2:RoIC 1.661e-01 2.916e-01 4.964e+04 0.570
## GroupMED21:Session2:RoIC 6.766e-01 2.763e-01 4.964e+04 2.448
## GroupMED1:Session3:RoIC -5.597e-01 2.914e-01 4.964e+04 -1.921
## GroupMED21:Session3:RoIC 4.447e-01 2.737e-01 4.964e+04 1.624
## GroupMED1:Session2:RoIP 5.286e-01 2.920e-01 4.964e+04 1.810
## GroupMED21:Session2:RoIP 1.270e+00 2.767e-01 4.964e+04 4.591
## GroupMED1:Session3:RoIP -3.674e-01 2.919e-01 4.964e+04 -1.259
## GroupMED21:Session3:RoIP 5.040e-01 2.741e-01 4.964e+04 1.839
## GroupMED1:Timepoint:RoIC -2.452e-02 3.321e-02 4.964e+04 -0.738
## GroupMED21:Timepoint:RoIC 7.246e-02 3.121e-02 4.964e+04 2.321
## GroupMED1:Timepoint:RoIP 1.754e-01 3.334e-02 4.964e+04 5.262
## GroupMED21:Timepoint:RoIP 1.475e-01 3.130e-02 4.964e+04 4.711
## Session2:Timepoint:RoIC 8.131e-02 3.421e-02 4.964e+04 2.377
## Session3:Timepoint:RoIC 2.849e-02 3.417e-02 4.964e+04 0.834
## Session2:Timepoint:RoIP 1.402e-01 3.428e-02 4.964e+04 4.091
## Session3:Timepoint:RoIP 3.271e-02 3.424e-02 4.964e+04 0.955
## GroupMED1:Session2:Timepoint:RoIC -2.928e-02 4.699e-02 4.964e+04 -0.623
## GroupMED21:Session2:Timepoint:RoIC -1.057e-01 4.453e-02 4.964e+04 -2.374
## GroupMED1:Session3:Timepoint:RoIC 4.869e-02 4.697e-02 4.964e+04 1.037
## GroupMED21:Session3:Timepoint:RoIC -8.346e-02 4.412e-02 4.964e+04 -1.892
## GroupMED1:Session2:Timepoint:RoIP -6.963e-02 4.706e-02 4.964e+04 -1.480
## GroupMED21:Session2:Timepoint:RoIP -2.280e-01 4.460e-02 4.964e+04 -5.113
## GroupMED1:Session3:Timepoint:RoIP -6.622e-02 4.704e-02 4.964e+04 -1.408
## GroupMED21:Session3:Timepoint:RoIP -4.301e-02 4.418e-02 4.964e+04 -0.973
## Pr(>|t|)
## (Intercept) 3.69e-15 ***
## GroupMED1 0.079695 .
## GroupMED21 0.036333 *
## Session2 0.195567
## Session3 0.425099
## Timepoint 0.444981
## RoIC 1.84e-14 ***
## RoIP < 2e-16 ***
## GroupMED1:Session2 0.775947
## GroupMED21:Session2 0.926078
## GroupMED1:Session3 0.336558
## GroupMED21:Session3 0.808048
## GroupMED1:Timepoint 0.675858
## GroupMED21:Timepoint 0.476182
## Session2:Timepoint 0.009755 **
## Session3:Timepoint 0.055555 .
## GroupMED1:RoIC 0.000188 ***
## GroupMED21:RoIC 0.516522
## GroupMED1:RoIP 0.680717
## GroupMED21:RoIP 0.019433 *
## Session2:RoIC 0.000530 ***
## Session3:RoIC 0.169095
## Session2:RoIP 8.67e-05 ***
## Session3:RoIP 0.986051
## Timepoint:RoIC 0.844911
## Timepoint:RoIP 0.588732
## GroupMED1:Session2:Timepoint 0.850773
## GroupMED21:Session2:Timepoint 0.095411 .
## GroupMED1:Session3:Timepoint 0.012630 *
## GroupMED21:Session3:Timepoint 5.55e-06 ***
## GroupMED1:Session2:RoIC 0.568891
## GroupMED21:Session2:RoIC 0.014351 *
## GroupMED1:Session3:RoIC 0.054790 .
## GroupMED21:Session3:RoIC 0.104276
## GroupMED1:Session2:RoIP 0.070244 .
## GroupMED21:Session2:RoIP 4.43e-06 ***
## GroupMED1:Session3:RoIP 0.208029
## GroupMED21:Session3:RoIP 0.065984 .
## GroupMED1:Timepoint:RoIC 0.460372
## GroupMED21:Timepoint:RoIC 0.020276 *
## GroupMED1:Timepoint:RoIP 1.43e-07 ***
## GroupMED21:Timepoint:RoIP 2.47e-06 ***
## Session2:Timepoint:RoIC 0.017461 *
## Session3:Timepoint:RoIC 0.404451
## Session2:Timepoint:RoIP 4.30e-05 ***
## Session3:Timepoint:RoIP 0.339507
## GroupMED1:Session2:Timepoint:RoIC 0.533175
## GroupMED21:Session2:Timepoint:RoIC 0.017624 *
## GroupMED1:Session3:Timepoint:RoIC 0.299922
## GroupMED21:Session3:Timepoint:RoIC 0.058522 .
## GroupMED1:Session2:Timepoint:RoIP 0.138974
## GroupMED21:Session2:Timepoint:RoIP 3.18e-07 ***
## GroupMED1:Session3:Timepoint:RoIP 0.159172
## GroupMED21:Session3:Timepoint:RoIP 0.330351
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1##
## Correlation matrix not shown by default, as p = 54 > 12.
## Use print(x, correlation=TRUE) or
## vcov(x) if you need it## optimizer (nloptwrap) convergence code: 0 (OK)
## Model failed to converge with max|grad| = 0.0105483 (tol = 0.002, component 1)#Model4: Theta Random Slope: Timepoint*Session, intercept:ID. Interaction between the two tracks the difference in the changes between sessions.
m.T4<-lmer(muTheta ~ Group * Session * Timepoint * RoI + (1 + Timepoint*Session | ID), data=d.MU, REML= FALSE)## Warning in checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
## Model failed to converge with max|grad| = 0.0170338 (tol = 0.002, component 1)Anova(m.T4)## Analysis of Deviance Table (Type II Wald chisquare tests)
##
## Response: muTheta
## Chisq Df Pr(>Chisq)
## Group 11.7504 2 0.0028082 **
## Session 2.1096 2 0.3482621
## Timepoint 11.2940 1 0.0007776 ***
## RoI 8323.4255 2 < 2.2e-16 ***
## Group:Session 7.5529 4 0.1093995
## Group:Timepoint 1.6346 2 0.4416221
## Session:Timepoint 3.1151 2 0.2106498
## Group:RoI 185.5377 4 < 2.2e-16 ***
## Session:RoI 42.6478 4 1.224e-08 ***
## Timepoint:RoI 369.7188 2 < 2.2e-16 ***
## Group:Session:Timepoint 3.6086 4 0.4615548
## Group:Session:RoI 117.3101 8 < 2.2e-16 ***
## Group:Timepoint:RoI 71.1658 4 1.288e-14 ***
## Session:Timepoint:RoI 4.9274 4 0.2948239
## Group:Session:Timepoint:RoI 56.1413 8 2.648e-09 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1m.T4test <- as(m.T4,"lmerModLmerTest")
summary(m.T4test, ddf="Satterthwaite")## Linear mixed model fit by maximum likelihood . t-tests use Satterthwaite's
## method [lmerModLmerTest]
## Formula: muTheta ~ Group * Session * Timepoint * RoI + (1 + Timepoint *
## Session | ID)
## Data: d.MU
##
## AIC BIC logLik deviance df.resid
## 207423.0 208093.0 -103635.5 207271.0 49674
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -4.9268 -0.6013 -0.1029 0.4823 10.9848
##
## Random effects:
## Groups Name Variance Std.Dev. Corr
## ID (Intercept) 1.04116 1.0204
## Timepoint 0.02104 0.1451 -0.61
## Session2 1.58411 1.2586 -0.13 0.22
## Session3 1.55020 1.2451 -0.38 0.46 0.42
## Timepoint:Session2 0.03273 0.1809 0.16 -0.56 -0.65 -0.34
## Timepoint:Session3 0.03104 0.1762 0.41 -0.67 -0.22 -0.63 0.54
## Residual 3.72328 1.9296
## Number of obs: 49750, groups: ID, 29
##
## Fixed effects:
## Estimate Std. Error df t value
## (Intercept) 5.011e+00 3.742e-01 3.174e+01 13.393
## GroupMED1 9.082e-01 5.143e-01 3.174e+01 1.766
## GroupMED21 1.028e+00 4.831e-01 3.175e+01 2.127
## Session2 4.908e-01 4.666e-01 3.164e+01 1.052
## Session3 3.045e-01 4.621e-01 3.261e+01 0.659
## Timepoint 3.092e-02 5.373e-02 3.248e+01 0.575
## RoIC -1.151e+00 1.484e-01 4.958e+04 -7.757
## RoIP -2.412e+00 1.490e-01 4.958e+04 -16.186
## GroupMED1:Session2 1.468e-01 6.413e-01 3.164e+01 0.229
## GroupMED21:Session2 7.102e-02 6.111e-01 3.193e+01 0.116
## GroupMED1:Session3 -5.052e-01 6.351e-01 3.261e+01 -0.795
## GroupMED21:Session3 -1.149e-01 5.966e-01 3.262e+01 -0.193
## GroupMED1:Timepoint -2.320e-02 7.384e-02 3.248e+01 -0.314
## GroupMED21:Timepoint -3.689e-02 6.937e-02 3.249e+01 -0.532
## Session2:Timepoint -5.919e-02 6.785e-02 3.289e+01 -0.872
## Session3:Timepoint -4.385e-02 6.628e-02 3.365e+01 -0.662
## GroupMED1:RoIC 7.705e-01 2.038e-01 4.958e+04 3.781
## GroupMED21:RoIC -1.229e-01 1.915e-01 4.958e+04 -0.642
## GroupMED1:RoIP -9.860e-02 2.046e-01 4.958e+04 -0.482
## GroupMED21:RoIP -4.545e-01 1.921e-01 4.958e+04 -2.366
## Session2:RoIC -7.308e-01 2.099e-01 4.958e+04 -3.482
## Session3:RoIC -2.907e-01 2.097e-01 4.958e+04 -1.386
## Session2:RoIP -8.336e-01 2.103e-01 4.958e+04 -3.964
## Session3:RoIP -3.787e-03 2.101e-01 4.958e+04 -0.018
## Timepoint:RoIC -4.557e-03 2.392e-02 4.958e+04 -0.191
## Timepoint:RoIP 1.312e-02 2.402e-02 4.958e+04 0.546
## GroupMED1:Session2:Timepoint 5.922e-03 9.325e-02 3.289e+01 0.064
## GroupMED21:Session2:Timepoint 4.777e-02 8.852e-02 3.343e+01 0.540
## GroupMED1:Session3:Timepoint 7.851e-02 9.109e-02 3.365e+01 0.862
## GroupMED21:Session3:Timepoint 1.336e-01 8.558e-02 3.367e+01 1.562
## GroupMED1:Session2:RoIC 1.533e-01 2.883e-01 4.958e+04 0.532
## GroupMED21:Session2:RoIC 6.701e-01 2.732e-01 4.958e+04 2.453
## GroupMED1:Session3:RoIC -5.664e-01 2.882e-01 4.958e+04 -1.966
## GroupMED21:Session3:RoIC 4.394e-01 2.707e-01 4.958e+04 1.623
## GroupMED1:Session2:RoIP 5.408e-01 2.887e-01 4.958e+04 1.873
## GroupMED21:Session2:RoIP 1.269e+00 2.736e-01 4.958e+04 4.639
## GroupMED1:Session3:RoIP -3.540e-01 2.886e-01 4.958e+04 -1.227
## GroupMED21:Session3:RoIP 5.019e-01 2.711e-01 4.958e+04 1.851
## GroupMED1:Timepoint:RoIC -2.470e-02 3.284e-02 4.958e+04 -0.752
## GroupMED21:Timepoint:RoIC 7.196e-02 3.086e-02 4.958e+04 2.332
## GroupMED1:Timepoint:RoIP 1.779e-01 3.297e-02 4.958e+04 5.396
## GroupMED21:Timepoint:RoIP 1.476e-01 3.096e-02 4.958e+04 4.767
## Session2:Timepoint:RoIC 8.045e-02 3.383e-02 4.958e+04 2.379
## Session3:Timepoint:RoIC 2.832e-02 3.379e-02 4.958e+04 0.838
## Session2:Timepoint:RoIP 1.400e-01 3.390e-02 4.958e+04 4.131
## Session3:Timepoint:RoIP 3.272e-02 3.386e-02 4.958e+04 0.966
## GroupMED1:Session2:Timepoint:RoIC -2.695e-02 4.646e-02 4.958e+04 -0.580
## GroupMED21:Session2:Timepoint:RoIC -1.045e-01 4.403e-02 4.958e+04 -2.374
## GroupMED1:Session3:Timepoint:RoIC 4.991e-02 4.644e-02 4.958e+04 1.075
## GroupMED21:Session3:Timepoint:RoIC -8.249e-02 4.362e-02 4.958e+04 -1.891
## GroupMED1:Session2:Timepoint:RoIP -7.185e-02 4.653e-02 4.958e+04 -1.544
## GroupMED21:Session2:Timepoint:RoIP -2.279e-01 4.410e-02 4.958e+04 -5.168
## GroupMED1:Session3:Timepoint:RoIP -6.868e-02 4.651e-02 4.958e+04 -1.477
## GroupMED21:Session3:Timepoint:RoIP -4.261e-02 4.369e-02 4.958e+04 -0.975
## Pr(>|t|)
## (Intercept) 1.32e-14 ***
## GroupMED1 0.087015 .
## GroupMED21 0.041245 *
## Session2 0.300858
## Session3 0.514542
## Timepoint 0.568989
## RoIC 8.87e-15 ***
## RoIP < 2e-16 ***
## GroupMED1:Session2 0.820472
## GroupMED21:Session2 0.908203
## GroupMED1:Session3 0.432094
## GroupMED21:Session3 0.848480
## GroupMED1:Timepoint 0.755412
## GroupMED21:Timepoint 0.598429
## Session2:Timepoint 0.389315
## Session3:Timepoint 0.512749
## GroupMED1:RoIC 0.000156 ***
## GroupMED21:RoIC 0.520882
## GroupMED1:RoIP 0.629786
## GroupMED21:RoIP 0.017988 *
## Session2:RoIC 0.000498 ***
## Session3:RoIC 0.165681
## Session2:RoIP 7.39e-05 ***
## Session3:RoIP 0.985618
## Timepoint:RoIC 0.848898
## Timepoint:RoIP 0.584999
## GroupMED1:Session2:Timepoint 0.949749
## GroupMED21:Session2:Timepoint 0.592970
## GroupMED1:Session3:Timepoint 0.394874
## GroupMED21:Session3:Timepoint 0.127732
## GroupMED1:Session2:RoIC 0.594967
## GroupMED21:Session2:RoIC 0.014187 *
## GroupMED1:Session3:RoIC 0.049337 *
## GroupMED21:Session3:RoIC 0.104558
## GroupMED1:Session2:RoIP 0.061084 .
## GroupMED21:Session2:RoIP 3.50e-06 ***
## GroupMED1:Session3:RoIP 0.220000
## GroupMED21:Session3:RoIP 0.064130 .
## GroupMED1:Timepoint:RoIC 0.452094
## GroupMED21:Timepoint:RoIC 0.019726 *
## GroupMED1:Timepoint:RoIP 6.85e-08 ***
## GroupMED21:Timepoint:RoIP 1.88e-06 ***
## Session2:Timepoint:RoIC 0.017386 *
## Session3:Timepoint:RoIC 0.402041
## Session2:Timepoint:RoIP 3.63e-05 ***
## Session3:Timepoint:RoIP 0.333877
## GroupMED1:Session2:Timepoint:RoIC 0.561861
## GroupMED21:Session2:Timepoint:RoIC 0.017615 *
## GroupMED1:Session3:Timepoint:RoIC 0.282469
## GroupMED21:Session3:Timepoint:RoIC 0.058643 .
## GroupMED1:Session2:Timepoint:RoIP 0.122577
## GroupMED21:Session2:Timepoint:RoIP 2.38e-07 ***
## GroupMED1:Session3:Timepoint:RoIP 0.139751
## GroupMED21:Session3:Timepoint:RoIP 0.329392
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1##
## Correlation matrix not shown by default, as p = 54 > 12.
## Use print(x, correlation=TRUE) or
## vcov(x) if you need it## optimizer (nloptwrap) convergence code: 0 (OK)
## Model failed to converge with max|grad| = 0.0170338 (tol = 0.002, component 1)#ANOVA on model differences.
#Will they run effectively?
anova(m.T0, m.T1, m.T2, m.T2.5, m.T3, m.T4)## Data: d.MU
## Models:
## m.T0: muTheta ~ Group * Session * Timepoint * RoI + (1 | ID)
## m.T1: muTheta ~ Group * Session * Timepoint * RoI + (1 | ID) + (1 |
## m.T1: Session:Timepoint)
## m.T2.5: muTheta ~ Group * Session * Timepoint * RoI + (1 + Timepoint |
## m.T2.5: ID)
## m.T2: muTheta ~ Group * Session * Timepoint * RoI + (1 + Session |
## m.T2: ID)
## m.T3: muTheta ~ Group * Session * Timepoint * RoI + (1 + Timepoint +
## m.T3: Session | ID)
## m.T4: muTheta ~ Group * Session * Timepoint * RoI + (1 + Timepoint *
## m.T4: Session | ID)
## npar AIC BIC logLik deviance Chisq Df Pr(>Chisq)
## m.T0 56 213032 213525 -106460 212920
## m.T1 57 212807 213309 -106346 212693 226.93 1 < 2.2e-16 ***
## m.T2.5 58 212074 212585 -105979 211958 734.99 1 < 2.2e-16 ***
## m.T2 61 209397 209935 -104638 209275 2682.37 3 < 2.2e-16 ***
## m.T3 65 208354 208927 -104112 208224 1050.98 4 < 2.2e-16 ***
## m.T4 76 207423 208093 -103636 207271 953.31 11 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1#Modelfit checking.
AIC(m.T0, m.T1, m.T2, m.T2.5, m.T3, m.T4)## df AIC
## m.T0 56 213031.6
## m.T1 57 212806.7
## m.T2 61 209397.3
## m.T2.5 58 212073.7
## m.T3 65 208354.3
## m.T4 76 207423.0#All effects Model 0 to 4 (6 models). In all these, Timepoint is numeric.
#Model 0 Effect
ae.m.T0 <- allEffects(m.T0)
ae.m.df.T0 <- as.data.frame(ae.m.T0[[1]])
#Ordering RoI
ae.m.df.T0$RoI <- factor(ae.m.df.T0$RoI, levels=c('F', 'C', 'P'))
#Ordering Timepoint
ae.m.df.T0$Timepoint <- as.character(ae.m.df.T0$Timepoint)
ae.m.df.T0$Timepoint <- as.numeric(ae.m.df.T0$Timepoint)
#Model 1 Effects
ae.m.T1 <- allEffects(m.T1)
ae.m.df.T1 <- as.data.frame(ae.m.T1[[1]])
#Ordering RoI
ae.m.df.T1$RoI <- factor(ae.m.df.T1$RoI, levels=c('F', 'C', 'P'))
#Ordering Timepoint
ae.m.df.T1$Timepoint <- as.character(ae.m.df.T1$Timepoint)
ae.m.df.T1$Timepoint <- as.numeric(ae.m.df.T1$Timepoint)
#Model 2 Effects
ae.m.T2 <- allEffects(m.T2)
ae.m.df.T2 <- as.data.frame(ae.m.T2[[1]])
#Ordering RoI
ae.m.df.T2$RoI <- factor(ae.m.df.T2$RoI, levels=c('F', 'C', 'P'))
#Ordering Timepoint
ae.m.df.T2$Timepoint <- as.character(ae.m.df.T2$Timepoint)
ae.m.df.T2$Timepoint <- as.numeric(ae.m.df.T2$Timepoint)
#Model 2.5 Effects
ae.m.T2.5 <- allEffects(m.T2.5)
ae.m.df.T2.5 <- as.data.frame(ae.m.T2.5[[1]])
#Ordering RoI
ae.m.df.T2.5$RoI <- factor(ae.m.df.T2.5$RoI, levels=c('F', 'C', 'P'))
#Ordering Timepoint
ae.m.df.T2.5$Timepoint <- as.character(ae.m.df.T1$Timepoint)
ae.m.df.T2.5$Timepoint <- as.numeric(ae.m.df.T1$Timepoint)
#Model 3 Effects
ae.m.T3 <- allEffects(m.T3)
ae.m.df.T3 <- as.data.frame(ae.m.T3[[1]])
#Ordering RoI
ae.m.df.T3$RoI <- factor(ae.m.df.T3$RoI, levels=c('F', 'C', 'P'))
#Ordering Timepoint
ae.m.df.T3$Timepoint <- as.character(ae.m.df.T3$Timepoint)
ae.m.df.T3$Timepoint <- as.numeric(ae.m.df.T3$Timepoint)
#Model 4 Effects
ae.m.T4 <- allEffects(m.T4)
ae.m.df.T4 <- as.data.frame(ae.m.T4[[1]])
#Ordering RoI
ae.m.df.T4$RoI <- factor(ae.m.df.T4$RoI, levels=c('F', 'C', 'P'))
#Ordering Timepoint
ae.m.df.T4$Timepoint <- as.character(ae.m.df.T4$Timepoint)
ae.m.df.T4$Timepoint <- as.numeric(ae.m.df.T4$Timepoint)#Plotting models 0 to 4. They all look the same….
#Model 0 plot: Interacttions
ae.T0 <- ggplot(ae.m.df.T0, aes(x=Timepoint, y=fit, ymin=lower, ymax=upper))+
geom_point(aes(color = Group))+
geom_path(aes(x=Timepoint, y=fit, color=Group))+
geom_ribbon(aes(ymin=lower, ymax=upper, fill=Group), alpha=0.2)+
ylab("Mean Theta Power (n.u.)")+
scale_x_continuous(name="Timpoint", breaks=c(1,2,3,4,5,6,7,8,9,10) , limits=c(1,10))+
ggtitle("muTheta Group x Session x Timepoint x RoI")+
theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),panel.background = element_blank(), axis.line = element_line(colour = "black"))+
facet_grid(RoI ~ Session)
plot(ae.T0)
#Model 1 effects
ae.T1 <- ggplot(ae.m.df.T1, aes(x=Timepoint, y=fit, ymin=lower, ymax=upper))+
geom_point(aes(color = Group))+
geom_path(aes(x=Timepoint, y=fit, color=Group))+
geom_ribbon(aes(ymin=lower, ymax=upper, fill=Group), alpha=0.2)+
ylab("Mean Theta Power (n.u.)")+
scale_x_continuous(name="Timpoint", breaks=c(1,2,3,4,5,6,7,8,9,10) , limits=c(1,10))+
ggtitle("muTheta Group x Session x Timepoint x RoI")+
theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),panel.background = element_blank(), axis.line = element_line(colour = "black"))+
facet_grid(RoI ~ Session)
plot(ae.T1)
#Model 2 effects
ae.T2 <- ggplot(ae.m.df.T2, aes(x=Timepoint, y=fit, ymin=lower, ymax=upper))+
geom_point(aes(color = Group))+
geom_path(aes(x=Timepoint, y=fit, color=Group))+
geom_ribbon(aes(ymin=lower, ymax=upper, fill=Group), alpha=0.2)+
ylab("Mean Theta Power (n.u.)")+
scale_x_continuous(name="Timpoint", breaks=c(1,2,3,4,5,6,7,8,9,10) , limits=c(1,10))+
ggtitle("muTheta Group x Session x Timepoint x RoI")+
theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),panel.background = element_blank(), axis.line = element_line(colour = "black"))+
facet_grid(RoI ~ Session)
plot(ae.T2)
#Model 2.5 effects
ae.T2.5 <- ggplot(ae.m.df.T2.5, aes(x=Timepoint, y=fit, ymin=lower, ymax=upper))+
geom_point(aes(color = Group))+
geom_path(aes(x=Timepoint, y=fit, color=Group))+
geom_ribbon(aes(ymin=lower, ymax=upper, fill=Group), alpha=0.2)+
ylab("Mean Theta Power (n.u.)")+
scale_x_continuous(name="Timpoint", breaks=c(1,2,3,4,5,6,7,8,9,10) , limits=c(1,10))+
ggtitle("muTheta Group x Session x Timepoint x RoI")+
theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),panel.background = element_blank(), axis.line = element_line(colour = "black"))+
facet_grid(RoI ~ Session)
plot(ae.T2.5)
#Model 3 effects
ae.T3 <- ggplot(ae.m.df.T3, aes(x=Timepoint, y=fit, ymin=lower, ymax=upper))+
geom_point(aes(color = Group))+
geom_path(aes(x=Timepoint, y=fit, color=Group))+
geom_ribbon(aes(ymin=lower, ymax=upper, fill=Group), alpha=0.2)+
ylab("Mean Theta Power (n.u.)")+
scale_x_continuous(name="Timpoint", breaks=c(1,2,3,4,5,6,7,8,9,10) , limits=c(1,10))+
ggtitle("muTheta Group x Session x Timepoint x RoI")+
theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),panel.background = element_blank(), axis.line = element_line(colour = "black"))+
facet_grid(RoI ~ Session)
plot(ae.T3)
#Model 4 effects
ae.T4 <- ggplot(ae.m.df.T4, aes(x=Timepoint, y=fit, ymin=lower, ymax=upper))+
geom_point(aes(color = Group))+
geom_path(aes(x=Timepoint, y=fit, color=Group))+
geom_ribbon(aes(ymin=lower, ymax=upper, fill=Group), alpha=0.2)+
ylab("Mean Theta Power (n.u.)")+
scale_x_continuous(name="Timpoint", breaks=c(1,2,3,4,5,6,7,8,9,10) , limits=c(1,10))+
ggtitle("muTheta Group x Session x Timepoint x RoI")+
theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),panel.background = element_blank(), axis.line = element_line(colour = "black"))+
facet_grid(RoI ~ Session)
plot(ae.T4)
#Plotting, All the plots look the same…however we’ll make Timepoints a factor and analyse this as well. #Running Timepoint as a factor and then plotting
d.MU$Timepoint <-factor(d.MU$Timepoint, levels=c('1', '2', '3', '4', '5', '6', '7', '8', '9', '10'))
m.T5 <- lmer(muTheta ~ Group * Session * Timepoint * RoI + (1 | ID) + (1 | Session:Timepoint), data=d.MU, REML = FALSE)## boundary (singular) fit: see ?isSingularAnova(m.T5)## Analysis of Deviance Table (Type II Wald chisquare tests)
##
## Response: muTheta
## Chisq Df Pr(>Chisq)
## Group 12.885 2 0.0015921 **
## Session 189.257 2 < 2.2e-16 ***
## Timepoint 609.893 9 < 2.2e-16 ***
## RoI 7489.149 2 < 2.2e-16 ***
## Group:Session 545.967 4 < 2.2e-16 ***
## Group:Timepoint 182.402 18 < 2.2e-16 ***
## Session:Timepoint 129.978 18 < 2.2e-16 ***
## Group:RoI 166.627 4 < 2.2e-16 ***
## Session:RoI 38.708 4 8.004e-08 ***
## Timepoint:RoI 391.113 18 < 2.2e-16 ***
## Group:Session:Timepoint 220.101 36 < 2.2e-16 ***
## Group:Session:RoI 104.541 8 < 2.2e-16 ***
## Group:Timepoint:RoI 109.775 36 2.182e-09 ***
## Session:Timepoint:RoI 50.272 36 0.0574544 .
## Group:Session:Timepoint:RoI 124.439 72 0.0001239 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1m.T5test <- as(m.T5,"lmerModLmerTest")
summary(m.T5test, ddf="Satterthwaite")## Linear mixed model fit by maximum likelihood . t-tests use Satterthwaite's
## method [lmerModLmerTest]
## Formula: muTheta ~ Group * Session * Timepoint * RoI + (1 | ID) + (1 |
## Session:Timepoint)
## Data: d.MU
##
## AIC BIC logLik deviance df.resid
## 212613.6 215020.0 -106033.8 212067.6 49477
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -4.1063 -0.6222 -0.1092 0.4901 10.8782
##
## Random effects:
## Groups Name Variance Std.Dev.
## Session:Timepoint (Intercept) 0.0000 0.0000
## ID (Intercept) 0.7697 0.8773
## Residual 4.1430 2.0354
## Number of obs: 49750, groups: Session:Timepoint, 30; ID, 29
##
## Fixed effects:
## Estimate Std. Error df t value
## (Intercept) 4.909e+00 3.461e-01 4.465e+01 14.185
## GroupMED1 7.716e-01 4.756e-01 4.465e+01 1.622
## GroupMED21 9.840e-01 4.468e-01 4.467e+01 2.202
## Session2 6.497e-01 2.170e-01 4.972e+04 2.994
## Session3 6.476e-01 2.170e-01 4.972e+04 2.984
## Timepoint2 1.718e-01 2.170e-01 4.972e+04 0.792
## Timepoint3 3.014e-01 2.170e-01 4.972e+04 1.389
## Timepoint4 5.354e-01 2.170e-01 4.972e+04 2.467
## Timepoint5 -2.465e-02 2.170e-01 4.972e+04 -0.114
## Timepoint6 2.320e-01 2.170e-01 4.972e+04 1.069
## Timepoint7 5.183e-01 2.170e-01 4.972e+04 2.389
## Timepoint8 1.396e-01 2.170e-01 4.972e+04 0.643
## Timepoint9 3.829e-01 2.170e-01 4.972e+04 1.765
## Timepoint10 4.696e-01 2.170e-01 4.972e+04 2.164
## RoIC -1.300e+00 2.292e-01 4.972e+04 -5.674
## RoIP -2.564e+00 2.301e-01 4.972e+04 -11.146
## GroupMED1:Session2 8.678e-02 2.982e-01 4.972e+04 0.291
## GroupMED21:Session2 -3.419e-01 2.828e-01 4.972e+04 -1.209
## GroupMED1:Session3 -8.791e-01 2.982e-01 4.972e+04 -2.948
## GroupMED21:Session3 -4.183e-01 2.803e-01 4.972e+04 -1.492
## GroupMED1:Timepoint2 1.632e-01 2.982e-01 4.972e+04 0.547
## GroupMED21:Timepoint2 -1.066e-01 2.803e-01 4.972e+04 -0.380
## GroupMED1:Timepoint3 2.404e-01 2.982e-01 4.972e+04 0.806
## GroupMED21:Timepoint3 -2.478e-01 2.803e-01 4.972e+04 -0.884
## GroupMED1:Timepoint4 -3.034e-01 2.982e-01 4.972e+04 -1.017
## GroupMED21:Timepoint4 -5.123e-02 2.803e-01 4.972e+04 -0.183
## GroupMED1:Timepoint5 2.521e-01 2.982e-01 4.972e+04 0.845
## GroupMED21:Timepoint5 1.739e-01 2.803e-01 4.972e+04 0.620
## GroupMED1:Timepoint6 1.564e-01 2.982e-01 4.972e+04 0.524
## GroupMED21:Timepoint6 -4.374e-02 2.803e-01 4.972e+04 -0.156
## GroupMED1:Timepoint7 -1.910e-01 2.982e-01 4.972e+04 -0.640
## GroupMED21:Timepoint7 -3.584e-01 2.803e-01 4.972e+04 -1.278
## GroupMED1:Timepoint8 -7.099e-02 2.982e-01 4.972e+04 -0.238
## GroupMED21:Timepoint8 -2.569e-01 2.803e-01 4.972e+04 -0.917
## GroupMED1:Timepoint9 -1.967e-02 2.982e-01 4.972e+04 -0.066
## GroupMED21:Timepoint9 -3.933e-01 2.803e-01 4.972e+04 -1.403
## GroupMED1:Timepoint10 -1.369e-01 2.982e-01 4.972e+04 -0.459
## GroupMED21:Timepoint10 -3.227e-01 2.803e-01 4.972e+04 -1.151
## Session2:Timepoint2 -4.386e-01 3.069e-01 4.972e+04 -1.429
## Session3:Timepoint2 -3.774e-01 3.069e-01 4.972e+04 -1.230
## Session2:Timepoint3 -1.255e-01 3.069e-01 4.972e+04 -0.409
## Session3:Timepoint3 -5.836e-01 3.069e-01 4.972e+04 -1.902
## Session2:Timepoint4 -7.446e-01 3.069e-01 4.972e+04 -2.426
## Session3:Timepoint4 -9.319e-01 3.069e-01 4.972e+04 -3.037
## Session2:Timepoint5 -3.247e-01 3.069e-01 4.972e+04 -1.058
## Session3:Timepoint5 -4.205e-01 3.069e-01 4.972e+04 -1.370
## Session2:Timepoint6 -5.397e-01 3.069e-01 4.972e+04 -1.759
## Session3:Timepoint6 -9.301e-01 3.069e-01 4.972e+04 -3.031
## Session2:Timepoint7 -7.441e-01 3.069e-01 4.972e+04 -2.425
## Session3:Timepoint7 -1.163e+00 3.069e-01 4.972e+04 -3.790
## Session2:Timepoint8 -5.964e-01 3.069e-01 4.972e+04 -1.944
## Session3:Timepoint8 -2.294e-01 3.069e-01 4.972e+04 -0.748
## Session2:Timepoint9 -8.525e-01 3.069e-01 4.972e+04 -2.778
## Session3:Timepoint9 -2.076e-01 3.069e-01 4.972e+04 -0.677
## Session2:Timepoint10 -4.780e-01 3.069e-01 4.972e+04 -1.558
## Session3:Timepoint10 -9.990e-01 3.069e-01 4.972e+04 -3.256
## GroupMED1:RoIC 6.123e-01 3.147e-01 4.972e+04 1.946
## GroupMED21:RoIC -9.098e-03 2.957e-01 4.972e+04 -0.031
## GroupMED1:RoIP -1.561e-01 3.158e-01 4.972e+04 -0.494
## GroupMED21:RoIP 2.184e-02 2.966e-01 4.972e+04 0.074
## Session2:RoIC -6.488e-01 3.241e-01 4.972e+04 -2.002
## Session3:RoIC -4.194e-01 3.238e-01 4.972e+04 -1.295
## Session2:RoIP -5.547e-01 3.247e-01 4.972e+04 -1.708
## Session3:RoIP -3.247e-01 3.244e-01 4.972e+04 -1.001
## Timepoint2:RoIC 6.041e-02 3.241e-01 4.972e+04 0.186
## Timepoint3:RoIC 1.168e-01 3.241e-01 4.972e+04 0.360
## Timepoint4:RoIC 1.183e-01 3.241e-01 4.972e+04 0.365
## Timepoint5:RoIC 5.779e-01 3.241e-01 4.972e+04 1.783
## Timepoint6:RoIC 3.829e-01 3.241e-01 4.972e+04 1.182
## Timepoint7:RoIC -1.671e-01 3.241e-01 4.972e+04 -0.516
## Timepoint8:RoIC 6.708e-03 3.241e-01 4.972e+04 0.021
## Timepoint9:RoIC -1.249e-02 3.241e-01 4.972e+04 -0.039
## Timepoint10:RoIC 1.447e-01 3.241e-01 4.972e+04 0.446
## Timepoint2:RoIP 6.990e-02 3.253e-01 4.972e+04 0.215
## Timepoint3:RoIP 1.971e-01 3.253e-01 4.972e+04 0.606
## Timepoint4:RoIP 2.726e-02 3.253e-01 4.972e+04 0.084
## Timepoint5:RoIP 8.740e-01 3.253e-01 4.972e+04 2.686
## Timepoint6:RoIP 5.943e-01 3.253e-01 4.972e+04 1.827
## Timepoint7:RoIP -1.519e-01 3.253e-01 4.972e+04 -0.467
## Timepoint8:RoIP 2.030e-01 3.253e-01 4.972e+04 0.624
## Timepoint9:RoIP 4.127e-02 3.253e-01 4.972e+04 0.127
## Timepoint10:RoIP 3.602e-01 3.253e-01 4.972e+04 1.107
## GroupMED1:Session2:Timepoint2 -3.366e-02 4.217e-01 4.972e+04 -0.080
## GroupMED21:Session2:Timepoint2 3.317e-01 3.999e-01 4.972e+04 0.829
## GroupMED1:Session3:Timepoint2 4.613e-01 4.217e-01 4.972e+04 1.094
## GroupMED21:Session3:Timepoint2 6.217e-01 3.964e-01 4.972e+04 1.568
## GroupMED1:Session2:Timepoint3 -4.936e-01 4.217e-01 4.972e+04 -1.171
## GroupMED21:Session2:Timepoint3 3.328e-01 3.999e-01 4.972e+04 0.832
## GroupMED1:Session3:Timepoint3 5.219e-01 4.217e-01 4.972e+04 1.237
## GroupMED21:Session3:Timepoint3 5.860e-01 3.964e-01 4.972e+04 1.478
## GroupMED1:Session2:Timepoint4 6.310e-01 4.217e-01 4.972e+04 1.496
## GroupMED21:Session2:Timepoint4 1.240e+00 3.999e-01 4.972e+04 3.100
## GroupMED1:Session3:Timepoint4 1.305e+00 4.217e-01 4.972e+04 3.095
## GroupMED21:Session3:Timepoint4 1.096e+00 3.964e-01 4.972e+04 2.765
## GroupMED1:Session2:Timepoint5 -1.317e-01 4.217e-01 4.972e+04 -0.312
## GroupMED21:Session2:Timepoint5 8.482e-01 3.999e-01 4.972e+04 2.121
## GroupMED1:Session3:Timepoint5 5.858e-01 4.217e-01 4.972e+04 1.389
## GroupMED21:Session3:Timepoint5 1.029e+00 3.964e-01 4.972e+04 2.596
## GroupMED1:Session2:Timepoint6 3.876e-01 4.217e-01 4.972e+04 0.919
## GroupMED21:Session2:Timepoint6 7.915e-01 3.999e-01 4.972e+04 1.979
## GroupMED1:Session3:Timepoint6 1.294e+00 4.217e-01 4.972e+04 3.068
## GroupMED21:Session3:Timepoint6 1.488e+00 3.964e-01 4.972e+04 3.753
## GroupMED1:Session2:Timepoint7 7.697e-01 4.217e-01 4.972e+04 1.825
## GroupMED21:Session2:Timepoint7 1.196e+00 3.999e-01 4.972e+04 2.992
## GroupMED1:Session3:Timepoint7 1.508e+00 4.217e-01 4.972e+04 3.577
## GroupMED21:Session3:Timepoint7 1.894e+00 3.964e-01 4.972e+04 4.777
## GroupMED1:Session2:Timepoint8 1.493e-01 4.217e-01 4.972e+04 0.354
## GroupMED21:Session2:Timepoint8 9.135e-01 3.999e-01 4.972e+04 2.284
## GroupMED1:Session3:Timepoint8 9.294e-01 4.217e-01 4.972e+04 2.204
## GroupMED21:Session3:Timepoint8 1.225e+00 3.964e-01 4.972e+04 3.089
## GroupMED1:Session2:Timepoint9 1.154e-01 4.217e-01 4.972e+04 0.274
## GroupMED21:Session2:Timepoint9 8.176e-01 3.999e-01 4.972e+04 2.045
## GroupMED1:Session3:Timepoint9 1.164e-01 4.217e-01 4.972e+04 0.276
## GroupMED21:Session3:Timepoint9 8.727e-01 3.964e-01 4.972e+04 2.201
## GroupMED1:Session2:Timepoint10 -4.685e-01 4.217e-01 4.972e+04 -1.111
## GroupMED21:Session2:Timepoint10 2.471e-01 3.999e-01 4.972e+04 0.618
## GroupMED1:Session3:Timepoint10 1.335e+00 4.217e-01 4.972e+04 3.165
## GroupMED21:Session3:Timepoint10 1.598e+00 3.964e-01 4.972e+04 4.030
## GroupMED1:Session2:RoIC 3.254e-01 4.452e-01 4.972e+04 0.731
## GroupMED21:Session2:RoIC 7.180e-01 4.219e-01 4.972e+04 1.702
## GroupMED1:Session3:RoIC -9.102e-02 4.449e-01 4.972e+04 -0.205
## GroupMED21:Session3:RoIC 3.238e-01 4.180e-01 4.972e+04 0.775
## GroupMED1:Session2:RoIP 2.341e-01 4.458e-01 4.972e+04 0.525
## GroupMED21:Session2:RoIP 8.020e-01 4.225e-01 4.972e+04 1.898
## GroupMED1:Session3:RoIP 1.684e-02 4.456e-01 4.972e+04 0.038
## GroupMED21:Session3:RoIP -4.962e-02 4.186e-01 4.972e+04 -0.119
## GroupMED1:Timepoint2:RoIC 4.571e-03 4.450e-01 4.972e+04 0.010
## GroupMED21:Timepoint2:RoIC 5.887e-02 4.182e-01 4.972e+04 0.141
## GroupMED1:Timepoint3:RoIC 3.003e-01 4.450e-01 4.972e+04 0.675
## GroupMED21:Timepoint3:RoIC 4.008e-01 4.182e-01 4.972e+04 0.958
## GroupMED1:Timepoint4:RoIC 1.228e-01 4.450e-01 4.972e+04 0.276
## GroupMED21:Timepoint4:RoIC 9.273e-02 4.182e-01 4.972e+04 0.222
## GroupMED1:Timepoint5:RoIC -2.185e-01 4.450e-01 4.972e+04 -0.491
## GroupMED21:Timepoint5:RoIC -1.999e-01 4.182e-01 4.972e+04 -0.478
## GroupMED1:Timepoint6:RoIC 4.979e-02 4.450e-01 4.972e+04 0.112
## GroupMED21:Timepoint6:RoIC -1.516e-02 4.182e-01 4.972e+04 -0.036
## GroupMED1:Timepoint7:RoIC 1.556e-01 4.450e-01 4.972e+04 0.350
## GroupMED21:Timepoint7:RoIC 7.360e-01 4.182e-01 4.972e+04 1.760
## GroupMED1:Timepoint8:RoIC 2.383e-01 4.450e-01 4.972e+04 0.536
## GroupMED21:Timepoint8:RoIC 7.570e-01 4.182e-01 4.972e+04 1.810
## GroupMED1:Timepoint9:RoIC 1.408e-01 4.450e-01 4.972e+04 0.316
## GroupMED21:Timepoint9:RoIC 3.430e-01 4.182e-01 4.972e+04 0.820
## GroupMED1:Timepoint10:RoIC -5.585e-01 4.450e-01 4.972e+04 -1.255
## GroupMED21:Timepoint10:RoIC 6.729e-01 4.182e-01 4.972e+04 1.609
## GroupMED1:Timepoint2:RoIP 1.067e-01 4.466e-01 4.972e+04 0.239
## GroupMED21:Timepoint2:RoIP -8.823e-02 4.194e-01 4.972e+04 -0.210
## GroupMED1:Timepoint3:RoIP 8.834e-01 4.466e-01 4.972e+04 1.978
## GroupMED21:Timepoint3:RoIP 1.736e-01 4.194e-01 4.972e+04 0.414
## GroupMED1:Timepoint4:RoIP 1.109e+00 4.466e-01 4.972e+04 2.484
## GroupMED21:Timepoint4:RoIP -1.210e-01 4.194e-01 4.972e+04 -0.289
## GroupMED1:Timepoint5:RoIP 8.121e-01 4.466e-01 4.972e+04 1.818
## GroupMED21:Timepoint5:RoIP -5.637e-01 4.194e-01 4.972e+04 -1.344
## GroupMED1:Timepoint6:RoIP 7.848e-01 4.466e-01 4.972e+04 1.757
## GroupMED21:Timepoint6:RoIP -2.793e-02 4.194e-01 4.972e+04 -0.067
## GroupMED1:Timepoint7:RoIP 1.839e+00 4.466e-01 4.972e+04 4.118
## GroupMED21:Timepoint7:RoIP 1.314e+00 4.194e-01 4.972e+04 3.134
## GroupMED1:Timepoint8:RoIP 1.948e+00 4.466e-01 4.972e+04 4.362
## GroupMED21:Timepoint8:RoIP 7.544e-01 4.194e-01 4.972e+04 1.799
## GroupMED1:Timepoint9:RoIP 1.749e+00 4.466e-01 4.972e+04 3.917
## GroupMED21:Timepoint9:RoIP 7.623e-01 4.194e-01 4.972e+04 1.818
## GroupMED1:Timepoint10:RoIP 1.072e+00 4.466e-01 4.972e+04 2.401
## GroupMED21:Timepoint10:RoIP 1.173e+00 4.194e-01 4.972e+04 2.797
## Session2:Timepoint2:RoIC 1.977e-01 4.583e-01 4.972e+04 0.431
## Session3:Timepoint2:RoIC 1.372e-01 4.579e-01 4.972e+04 0.300
## Session2:Timepoint3:RoIC 1.564e-01 4.583e-01 4.972e+04 0.341
## Session3:Timepoint3:RoIC 4.645e-01 4.579e-01 4.972e+04 1.014
## Session2:Timepoint4:RoIC 3.704e-01 4.583e-01 4.972e+04 0.808
## Session3:Timepoint4:RoIC 3.930e-01 4.579e-01 4.972e+04 0.858
## Session2:Timepoint5:RoIC -1.372e-03 4.583e-01 4.972e+04 -0.003
## Session3:Timepoint5:RoIC -1.135e-01 4.579e-01 4.972e+04 -0.248
## Session2:Timepoint6:RoIC 2.490e-02 4.583e-01 4.972e+04 0.054
## Session3:Timepoint6:RoIC 2.079e-01 4.579e-01 4.972e+04 0.454
## Session2:Timepoint7:RoIC 6.662e-01 4.583e-01 4.972e+04 1.454
## Session3:Timepoint7:RoIC 7.628e-01 4.579e-01 4.972e+04 1.666
## Session2:Timepoint8:RoIC 8.602e-01 4.583e-01 4.972e+04 1.877
## Session3:Timepoint8:RoIC 5.213e-01 4.579e-01 4.972e+04 1.139
## Session2:Timepoint9:RoIC 8.541e-01 4.583e-01 4.972e+04 1.863
## Session3:Timepoint9:RoIC 1.859e-01 4.579e-01 4.972e+04 0.406
## Session2:Timepoint10:RoIC 4.952e-01 4.583e-01 4.972e+04 1.080
## Session3:Timepoint10:RoIC 2.971e-01 4.579e-01 4.972e+04 0.649
## Session2:Timepoint2:RoIP 2.283e-01 4.592e-01 4.972e+04 0.497
## Session3:Timepoint2:RoIP 4.178e-01 4.588e-01 4.972e+04 0.911
## Session2:Timepoint3:RoIP 4.483e-02 4.592e-01 4.972e+04 0.098
## Session3:Timepoint3:RoIP 5.435e-01 4.588e-01 4.972e+04 1.185
## Session2:Timepoint4:RoIP 3.937e-01 4.592e-01 4.972e+04 0.857
## Session3:Timepoint4:RoIP 8.133e-01 4.588e-01 4.972e+04 1.773
## Session2:Timepoint5:RoIP -1.912e-01 4.592e-01 4.972e+04 -0.416
## Session3:Timepoint5:RoIP 9.130e-02 4.588e-01 4.972e+04 0.199
## Session2:Timepoint6:RoIP 1.982e-01 4.592e-01 4.972e+04 0.432
## Session3:Timepoint6:RoIP 4.716e-01 4.588e-01 4.972e+04 1.028
## Session2:Timepoint7:RoIP 1.028e+00 4.592e-01 4.972e+04 2.238
## Session3:Timepoint7:RoIP 1.251e+00 4.588e-01 4.972e+04 2.726
## Session2:Timepoint8:RoIP 1.015e+00 4.592e-01 4.972e+04 2.211
## Session3:Timepoint8:RoIP 7.636e-01 4.588e-01 4.972e+04 1.664
## Session2:Timepoint9:RoIP 1.284e+00 4.592e-01 4.972e+04 2.796
## Session3:Timepoint9:RoIP 3.578e-01 4.588e-01 4.972e+04 0.780
## Session2:Timepoint10:RoIP 9.435e-01 4.592e-01 4.972e+04 2.055
## Session3:Timepoint10:RoIP 3.263e-01 4.588e-01 4.972e+04 0.711
## GroupMED1:Session2:Timepoint2:RoIC -1.445e-01 6.296e-01 4.972e+04 -0.230
## GroupMED21:Session2:Timepoint2:RoIC -1.733e-01 5.967e-01 4.972e+04 -0.290
## GroupMED1:Session3:Timepoint2:RoIC -3.799e-01 6.293e-01 4.972e+04 -0.604
## GroupMED21:Session3:Timepoint2:RoIC 1.421e-01 5.911e-01 4.972e+04 0.240
## GroupMED1:Session2:Timepoint3:RoIC -3.805e-01 6.296e-01 4.972e+04 -0.604
## GroupMED21:Session2:Timepoint3:RoIC -5.291e-01 5.967e-01 4.972e+04 -0.887
## GroupMED1:Session3:Timepoint3:RoIC -8.567e-01 6.293e-01 4.972e+04 -1.361
## GroupMED21:Session3:Timepoint3:RoIC -4.558e-01 5.911e-01 4.972e+04 -0.771
## GroupMED1:Session2:Timepoint4:RoIC -4.561e-01 6.296e-01 4.972e+04 -0.724
## GroupMED21:Session2:Timepoint4:RoIC -9.029e-01 5.967e-01 4.972e+04 -1.513
## GroupMED1:Session3:Timepoint4:RoIC -4.223e-01 6.293e-01 4.972e+04 -0.671
## GroupMED21:Session3:Timepoint4:RoIC -1.720e-01 5.911e-01 4.972e+04 -0.291
## GroupMED1:Session2:Timepoint5:RoIC -1.180e-01 6.296e-01 4.972e+04 -0.187
## GroupMED21:Session2:Timepoint5:RoIC -9.952e-02 5.967e-01 4.972e+04 -0.167
## GroupMED1:Session3:Timepoint5:RoIC 2.367e-01 6.293e-01 4.972e+04 0.376
## GroupMED21:Session3:Timepoint5:RoIC 1.060e-01 5.911e-01 4.972e+04 0.179
## GroupMED1:Session2:Timepoint6:RoIC -1.868e-01 6.296e-01 4.972e+04 -0.297
## GroupMED21:Session2:Timepoint6:RoIC -3.610e-01 5.967e-01 4.972e+04 -0.605
## GroupMED1:Session3:Timepoint6:RoIC -3.797e-01 6.293e-01 4.972e+04 -0.603
## GroupMED21:Session3:Timepoint6:RoIC -2.173e-01 5.911e-01 4.972e+04 -0.368
## GroupMED1:Session2:Timepoint7:RoIC -6.021e-01 6.296e-01 4.972e+04 -0.956
## GroupMED21:Session2:Timepoint7:RoIC -1.077e+00 5.967e-01 4.972e+04 -1.806
## GroupMED1:Session3:Timepoint7:RoIC -3.650e-01 6.293e-01 4.972e+04 -0.580
## GroupMED21:Session3:Timepoint7:RoIC -1.052e+00 5.911e-01 4.972e+04 -1.779
## GroupMED1:Session2:Timepoint8:RoIC -8.040e-01 6.296e-01 4.972e+04 -1.277
## GroupMED21:Session2:Timepoint8:RoIC -1.392e+00 5.967e-01 4.972e+04 -2.333
## GroupMED1:Session3:Timepoint8:RoIC -7.003e-02 6.293e-01 4.972e+04 -0.111
## GroupMED21:Session3:Timepoint8:RoIC -1.003e+00 5.911e-01 4.972e+04 -1.697
## GroupMED1:Session2:Timepoint9:RoIC -4.418e-01 6.296e-01 4.972e+04 -0.702
## GroupMED21:Session2:Timepoint9:RoIC -9.843e-01 5.967e-01 4.972e+04 -1.650
## GroupMED1:Session3:Timepoint9:RoIC 1.512e-02 6.293e-01 4.972e+04 0.024
## GroupMED21:Session3:Timepoint9:RoIC 1.057e-01 5.911e-01 4.972e+04 0.179
## GroupMED1:Session2:Timepoint10:RoIC -3.739e-02 6.296e-01 4.972e+04 -0.059
## GroupMED21:Session2:Timepoint10:RoIC -7.431e-01 5.967e-01 4.972e+04 -1.245
## GroupMED1:Session3:Timepoint10:RoIC 2.072e-01 6.293e-01 4.972e+04 0.329
## GroupMED21:Session3:Timepoint10:RoIC -8.719e-01 5.911e-01 4.972e+04 -1.475
## GroupMED1:Session2:Timepoint2:RoIP 5.201e-01 6.305e-01 4.972e+04 0.825
## GroupMED21:Session2:Timepoint2:RoIP 8.825e-02 5.975e-01 4.972e+04 0.148
## GroupMED1:Session3:Timepoint2:RoIP -4.239e-01 6.302e-01 4.972e+04 -0.673
## GroupMED21:Session3:Timepoint2:RoIP 5.988e-01 5.919e-01 4.972e+04 1.012
## GroupMED1:Session2:Timepoint3:RoIP 1.957e-01 6.305e-01 4.972e+04 0.310
## GroupMED21:Session2:Timepoint3:RoIP 4.210e-02 5.975e-01 4.972e+04 0.070
## GroupMED1:Session3:Timepoint3:RoIP -9.732e-01 6.302e-01 4.972e+04 -1.544
## GroupMED21:Session3:Timepoint3:RoIP 4.569e-01 5.919e-01 4.972e+04 0.772
## GroupMED1:Session2:Timepoint4:RoIP -4.777e-01 6.305e-01 4.972e+04 -0.758
## GroupMED21:Session2:Timepoint4:RoIP -9.447e-01 5.975e-01 4.972e+04 -1.581
## GroupMED1:Session3:Timepoint4:RoIP -1.229e+00 6.302e-01 4.972e+04 -1.950
## GroupMED21:Session3:Timepoint4:RoIP 6.242e-01 5.919e-01 4.972e+04 1.055
## GroupMED1:Session2:Timepoint5:RoIP 3.040e-02 6.305e-01 4.972e+04 0.048
## GroupMED21:Session2:Timepoint5:RoIP -6.989e-02 5.975e-01 4.972e+04 -0.117
## GroupMED1:Session3:Timepoint5:RoIP -1.187e-02 6.302e-01 4.972e+04 -0.019
## GroupMED21:Session3:Timepoint5:RoIP 9.470e-01 5.919e-01 4.972e+04 1.600
## GroupMED1:Session2:Timepoint6:RoIP 6.668e-01 6.305e-01 4.972e+04 1.058
## GroupMED21:Session2:Timepoint6:RoIP -4.695e-01 5.975e-01 4.972e+04 -0.786
## GroupMED1:Session3:Timepoint6:RoIP -3.547e-01 6.302e-01 4.972e+04 -0.563
## GroupMED21:Session3:Timepoint6:RoIP 5.212e-01 5.919e-01 4.972e+04 0.880
## GroupMED1:Session2:Timepoint7:RoIP -6.638e-01 6.305e-01 4.972e+04 -1.053
## GroupMED21:Session2:Timepoint7:RoIP -1.842e+00 5.975e-01 4.972e+04 -3.083
## GroupMED1:Session3:Timepoint7:RoIP -1.480e+00 6.302e-01 4.972e+04 -2.348
## GroupMED21:Session3:Timepoint7:RoIP -7.564e-01 5.919e-01 4.972e+04 -1.278
## GroupMED1:Session2:Timepoint8:RoIP -6.396e-01 6.305e-01 4.972e+04 -1.014
## GroupMED21:Session2:Timepoint8:RoIP -1.367e+00 5.975e-01 4.972e+04 -2.287
## GroupMED1:Session3:Timepoint8:RoIP -1.833e+00 6.302e-01 4.972e+04 -2.909
## GroupMED21:Session3:Timepoint8:RoIP 3.882e-01 5.919e-01 4.972e+04 0.656
## GroupMED1:Session2:Timepoint9:RoIP -4.050e-01 6.305e-01 4.972e+04 -0.642
## GroupMED21:Session2:Timepoint9:RoIP -1.583e+00 5.975e-01 4.972e+04 -2.650
## GroupMED1:Session3:Timepoint9:RoIP -9.752e-01 6.302e-01 4.972e+04 -1.548
## GroupMED21:Session3:Timepoint9:RoIP 5.886e-01 5.919e-01 4.972e+04 0.994
## GroupMED1:Session2:Timepoint10:RoIP -6.417e-02 6.305e-01 4.972e+04 -0.102
## GroupMED21:Session2:Timepoint10:RoIP -1.743e+00 5.975e-01 4.972e+04 -2.917
## GroupMED1:Session3:Timepoint10:RoIP -1.514e-01 6.302e-01 4.972e+04 -0.240
## GroupMED21:Session3:Timepoint10:RoIP -2.312e-01 5.919e-01 4.972e+04 -0.391
## Pr(>|t|)
## (Intercept) < 2e-16 ***
## GroupMED1 0.111774
## GroupMED21 0.032855 *
## Session2 0.002754 **
## Session3 0.002842 **
## Timepoint2 0.428377
## Timepoint3 0.164865
## Timepoint4 0.013610 *
## Timepoint5 0.909535
## Timepoint6 0.284918
## Timepoint7 0.016915 *
## Timepoint8 0.519952
## Timepoint9 0.077605 .
## Timepoint10 0.030433 *
## RoIC 1.41e-08 ***
## RoIP < 2e-16 ***
## GroupMED1:Session2 0.771058
## GroupMED21:Session2 0.226747
## GroupMED1:Session3 0.003200 **
## GroupMED21:Session3 0.135618
## GroupMED1:Timepoint2 0.584214
## GroupMED21:Timepoint2 0.703633
## GroupMED1:Timepoint3 0.420165
## GroupMED21:Timepoint3 0.376790
## GroupMED1:Timepoint4 0.308994
## GroupMED21:Timepoint4 0.854993
## GroupMED1:Timepoint5 0.397871
## GroupMED21:Timepoint5 0.535041
## GroupMED1:Timepoint6 0.599943
## GroupMED21:Timepoint6 0.876001
## GroupMED1:Timepoint7 0.521917
## GroupMED21:Timepoint7 0.201108
## GroupMED1:Timepoint8 0.811845
## GroupMED21:Timepoint8 0.359404
## GroupMED1:Timepoint9 0.947401
## GroupMED21:Timepoint9 0.160610
## GroupMED1:Timepoint10 0.646167
## GroupMED21:Timepoint10 0.249695
## Session2:Timepoint2 0.152906
## Session3:Timepoint2 0.218786
## Session2:Timepoint3 0.682431
## Session3:Timepoint3 0.057213 .
## Session2:Timepoint4 0.015251 *
## Session3:Timepoint4 0.002392 **
## Session2:Timepoint5 0.290050
## Session3:Timepoint5 0.170605
## Session2:Timepoint6 0.078588 .
## Session3:Timepoint6 0.002438 **
## Session2:Timepoint7 0.015318 *
## Session3:Timepoint7 0.000151 ***
## Session2:Timepoint8 0.051958 .
## Session3:Timepoint8 0.454723
## Session2:Timepoint9 0.005468 **
## Session3:Timepoint9 0.498635
## Session2:Timepoint10 0.119327
## Session3:Timepoint10 0.001132 **
## GroupMED1:RoIC 0.051657 .
## GroupMED21:RoIC 0.975456
## GroupMED1:RoIP 0.621019
## GroupMED21:RoIP 0.941303
## Session2:RoIC 0.045285 *
## Session3:RoIC 0.195177
## Session2:RoIP 0.087564 .
## Session3:RoIP 0.316873
## Timepoint2:RoIC 0.852137
## Timepoint3:RoIC 0.718544
## Timepoint4:RoIC 0.715071
## Timepoint5:RoIC 0.074555 .
## Timepoint6:RoIC 0.237382
## Timepoint7:RoIC 0.606080
## Timepoint8:RoIC 0.983487
## Timepoint9:RoIC 0.969251
## Timepoint10:RoIC 0.655281
## Timepoint2:RoIP 0.829892
## Timepoint3:RoIP 0.544569
## Timepoint4:RoIP 0.933232
## Timepoint5:RoIP 0.007227 **
## Timepoint6:RoIP 0.067764 .
## Timepoint7:RoIP 0.640614
## Timepoint8:RoIP 0.532735
## Timepoint9:RoIP 0.899052
## Timepoint10:RoIP 0.268259
## GroupMED1:Session2:Timepoint2 0.936385
## GroupMED21:Session2:Timepoint2 0.406847
## GroupMED1:Session3:Timepoint2 0.274005
## GroupMED21:Session3:Timepoint2 0.116852
## GroupMED1:Session2:Timepoint3 0.241802
## GroupMED21:Session2:Timepoint3 0.405341
## GroupMED1:Session3:Timepoint3 0.215928
## GroupMED21:Session3:Timepoint3 0.139373
## GroupMED1:Session2:Timepoint4 0.134608
## GroupMED21:Session2:Timepoint4 0.001935 **
## GroupMED1:Session3:Timepoint4 0.001967 **
## GroupMED21:Session3:Timepoint4 0.005693 **
## GroupMED1:Session2:Timepoint5 0.754853
## GroupMED21:Session2:Timepoint5 0.033917 *
## GroupMED1:Session3:Timepoint5 0.164837
## GroupMED21:Session3:Timepoint5 0.009426 **
## GroupMED1:Session2:Timepoint6 0.358047
## GroupMED21:Session2:Timepoint6 0.047773 *
## GroupMED1:Session3:Timepoint6 0.002154 **
## GroupMED21:Session3:Timepoint6 0.000175 ***
## GroupMED1:Session2:Timepoint7 0.067998 .
## GroupMED21:Session2:Timepoint7 0.002773 **
## GroupMED1:Session3:Timepoint7 0.000348 ***
## GroupMED21:Session3:Timepoint7 1.79e-06 ***
## GroupMED1:Session2:Timepoint8 0.723263
## GroupMED21:Session2:Timepoint8 0.022353 *
## GroupMED1:Session3:Timepoint8 0.027541 *
## GroupMED21:Session3:Timepoint8 0.002011 **
## GroupMED1:Session2:Timepoint9 0.784343
## GroupMED21:Session2:Timepoint9 0.040908 *
## GroupMED1:Session3:Timepoint9 0.782582
## GroupMED21:Session3:Timepoint9 0.027715 *
## GroupMED1:Session2:Timepoint10 0.266602
## GroupMED21:Session2:Timepoint10 0.536607
## GroupMED1:Session3:Timepoint10 0.001551 **
## GroupMED21:Session3:Timepoint10 5.58e-05 ***
## GroupMED1:Session2:RoIC 0.464849
## GroupMED21:Session2:RoIC 0.088782 .
## GroupMED1:Session3:RoIC 0.837911
## GroupMED21:Session3:RoIC 0.438519
## GroupMED1:Session2:RoIP 0.599449
## GroupMED21:Session2:RoIP 0.057663 .
## GroupMED1:Session3:RoIP 0.969857
## GroupMED21:Session3:RoIP 0.905639
## GroupMED1:Timepoint2:RoIC 0.991805
## GroupMED21:Timepoint2:RoIC 0.888050
## GroupMED1:Timepoint3:RoIC 0.499793
## GroupMED21:Timepoint3:RoIC 0.337860
## GroupMED1:Timepoint4:RoIC 0.782527
## GroupMED21:Timepoint4:RoIC 0.824517
## GroupMED1:Timepoint5:RoIC 0.623368
## GroupMED21:Timepoint5:RoIC 0.632591
## GroupMED1:Timepoint6:RoIC 0.910918
## GroupMED21:Timepoint6:RoIC 0.971080
## GroupMED1:Timepoint7:RoIC 0.726674
## GroupMED21:Timepoint7:RoIC 0.078440 .
## GroupMED1:Timepoint8:RoIC 0.592274
## GroupMED21:Timepoint8:RoIC 0.070271 .
## GroupMED1:Timepoint9:RoIC 0.751709
## GroupMED21:Timepoint9:RoIC 0.412140
## GroupMED1:Timepoint10:RoIC 0.209434
## GroupMED21:Timepoint10:RoIC 0.107641
## GroupMED1:Timepoint2:RoIP 0.811268
## GroupMED21:Timepoint2:RoIP 0.833366
## GroupMED1:Timepoint3:RoIP 0.047946 *
## GroupMED21:Timepoint3:RoIP 0.678996
## GroupMED1:Timepoint4:RoIP 0.013005 *
## GroupMED21:Timepoint4:RoIP 0.772953
## GroupMED1:Timepoint5:RoIP 0.069014 .
## GroupMED21:Timepoint5:RoIP 0.178936
## GroupMED1:Timepoint6:RoIP 0.078910 .
## GroupMED21:Timepoint6:RoIP 0.946908
## GroupMED1:Timepoint7:RoIP 3.83e-05 ***
## GroupMED21:Timepoint7:RoIP 0.001726 **
## GroupMED1:Timepoint8:RoIP 1.29e-05 ***
## GroupMED21:Timepoint8:RoIP 0.072073 .
## GroupMED1:Timepoint9:RoIP 8.98e-05 ***
## GroupMED21:Timepoint9:RoIP 0.069118 .
## GroupMED1:Timepoint10:RoIP 0.016365 *
## GroupMED21:Timepoint10:RoIP 0.005164 **
## Session2:Timepoint2:RoIC 0.666150
## Session3:Timepoint2:RoIC 0.764397
## Session2:Timepoint3:RoIC 0.732899
## Session3:Timepoint3:RoIC 0.310410
## Session2:Timepoint4:RoIC 0.418982
## Session3:Timepoint4:RoIC 0.390756
## Session2:Timepoint5:RoIC 0.997611
## Session3:Timepoint5:RoIC 0.804206
## Session2:Timepoint6:RoIC 0.956664
## Session3:Timepoint6:RoIC 0.649806
## Session2:Timepoint7:RoIC 0.146069
## Session3:Timepoint7:RoIC 0.095731 .
## Session2:Timepoint8:RoIC 0.060553 .
## Session3:Timepoint8:RoIC 0.254876
## Session2:Timepoint9:RoIC 0.062398 .
## Session3:Timepoint9:RoIC 0.684704
## Session2:Timepoint10:RoIC 0.279931
## Session3:Timepoint10:RoIC 0.516371
## Session2:Timepoint2:RoIP 0.619046
## Session3:Timepoint2:RoIP 0.362462
## Session2:Timepoint3:RoIP 0.922233
## Session3:Timepoint3:RoIP 0.236114
## Session2:Timepoint4:RoIP 0.391200
## Session3:Timepoint4:RoIP 0.076262 .
## Session2:Timepoint5:RoIP 0.677144
## Session3:Timepoint5:RoIP 0.842256
## Session2:Timepoint6:RoIP 0.666060
## Session3:Timepoint6:RoIP 0.303940
## Session2:Timepoint7:RoIP 0.025241 *
## Session3:Timepoint7:RoIP 0.006405 **
## Session2:Timepoint8:RoIP 0.027051 *
## Session3:Timepoint8:RoIP 0.096035 .
## Session2:Timepoint9:RoIP 0.005183 **
## Session3:Timepoint9:RoIP 0.435382
## Session2:Timepoint10:RoIP 0.039924 *
## Session3:Timepoint10:RoIP 0.476935
## GroupMED1:Session2:Timepoint2:RoIC 0.818447
## GroupMED21:Session2:Timepoint2:RoIC 0.771457
## GroupMED1:Session3:Timepoint2:RoIC 0.545991
## GroupMED21:Session3:Timepoint2:RoIC 0.810041
## GroupMED1:Session2:Timepoint3:RoIC 0.545580
## GroupMED21:Session2:Timepoint3:RoIC 0.375185
## GroupMED1:Session3:Timepoint3:RoIC 0.173397
## GroupMED21:Session3:Timepoint3:RoIC 0.440613
## GroupMED1:Session2:Timepoint4:RoIC 0.468824
## GroupMED21:Session2:Timepoint4:RoIC 0.130229
## GroupMED1:Session3:Timepoint4:RoIC 0.502105
## GroupMED21:Session3:Timepoint4:RoIC 0.771097
## GroupMED1:Session2:Timepoint5:RoIC 0.851359
## GroupMED21:Session2:Timepoint5:RoIC 0.867532
## GroupMED1:Session3:Timepoint5:RoIC 0.706752
## GroupMED21:Session3:Timepoint5:RoIC 0.857635
## GroupMED1:Session2:Timepoint6:RoIC 0.766701
## GroupMED21:Session2:Timepoint6:RoIC 0.545131
## GroupMED1:Session3:Timepoint6:RoIC 0.546222
## GroupMED21:Session3:Timepoint6:RoIC 0.713177
## GroupMED1:Session2:Timepoint7:RoIC 0.338902
## GroupMED21:Session2:Timepoint7:RoIC 0.070959 .
## GroupMED1:Session3:Timepoint7:RoIC 0.561914
## GroupMED21:Session3:Timepoint7:RoIC 0.075200 .
## GroupMED1:Session2:Timepoint8:RoIC 0.201579
## GroupMED21:Session2:Timepoint8:RoIC 0.019659 *
## GroupMED1:Session3:Timepoint8:RoIC 0.911381
## GroupMED21:Session3:Timepoint8:RoIC 0.089685 .
## GroupMED1:Session2:Timepoint9:RoIC 0.482811
## GroupMED21:Session2:Timepoint9:RoIC 0.099030 .
## GroupMED1:Session3:Timepoint9:RoIC 0.980832
## GroupMED21:Session3:Timepoint9:RoIC 0.858058
## GroupMED1:Session2:Timepoint10:RoIC 0.952646
## GroupMED21:Session2:Timepoint10:RoIC 0.212985
## GroupMED1:Session3:Timepoint10:RoIC 0.741933
## GroupMED21:Session3:Timepoint10:RoIC 0.140204
## GroupMED1:Session2:Timepoint2:RoIP 0.409424
## GroupMED21:Session2:Timepoint2:RoIP 0.882585
## GroupMED1:Session3:Timepoint2:RoIP 0.501107
## GroupMED21:Session3:Timepoint2:RoIP 0.311748
## GroupMED1:Session2:Timepoint3:RoIP 0.756214
## GroupMED21:Session2:Timepoint3:RoIP 0.943833
## GroupMED1:Session3:Timepoint3:RoIP 0.122497
## GroupMED21:Session3:Timepoint3:RoIP 0.440187
## GroupMED1:Session2:Timepoint4:RoIP 0.448681
## GroupMED21:Session2:Timepoint4:RoIP 0.113869
## GroupMED1:Session3:Timepoint4:RoIP 0.051228 .
## GroupMED21:Session3:Timepoint4:RoIP 0.291614
## GroupMED1:Session2:Timepoint5:RoIP 0.961543
## GroupMED21:Session2:Timepoint5:RoIP 0.906889
## GroupMED1:Session3:Timepoint5:RoIP 0.984967
## GroupMED21:Session3:Timepoint5:RoIP 0.109648
## GroupMED1:Session2:Timepoint6:RoIP 0.290223
## GroupMED21:Session2:Timepoint6:RoIP 0.432028
## GroupMED1:Session3:Timepoint6:RoIP 0.573501
## GroupMED21:Session3:Timepoint6:RoIP 0.378594
## GroupMED1:Session2:Timepoint7:RoIP 0.292449
## GroupMED21:Session2:Timepoint7:RoIP 0.002050 **
## GroupMED1:Session3:Timepoint7:RoIP 0.018856 *
## GroupMED21:Session3:Timepoint7:RoIP 0.201284
## GroupMED1:Session2:Timepoint8:RoIP 0.310351
## GroupMED21:Session2:Timepoint8:RoIP 0.022195 *
## GroupMED1:Session3:Timepoint8:RoIP 0.003627 **
## GroupMED21:Session3:Timepoint8:RoIP 0.511963
## GroupMED1:Session2:Timepoint9:RoIP 0.520652
## GroupMED21:Session2:Timepoint9:RoIP 0.008058 **
## GroupMED1:Session3:Timepoint9:RoIP 0.121738
## GroupMED21:Session3:Timepoint9:RoIP 0.320049
## GroupMED1:Session2:Timepoint10:RoIP 0.918928
## GroupMED21:Session2:Timepoint10:RoIP 0.003534 **
## GroupMED1:Session3:Timepoint10:RoIP 0.810104
## GroupMED21:Session3:Timepoint10:RoIP 0.696058
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1##
## Correlation matrix not shown by default, as p = 270 > 12.
## Use print(x, correlation=TRUE) or
## vcov(x) if you need it## optimizer (nloptwrap) convergence code: 0 (OK)
## boundary (singular) fit: see ?isSingular#Model 5 Effects. This plot looks more realistic to how the data looks like to show Group differences at the sessions, timepoints and RoI.
ae.m.T5 <- allEffects(m.T5)
ae.m.df.T5 <- as.data.frame(ae.m.T5[[1]])
#Ordering RoI
ae.m.df.T5$RoI <- factor(ae.m.df.T5$RoI, levels=c('F', 'C', 'P'))
#Ordering Timepoint for plotting sick.
ae.m.df.T5$Timepoint <- as.character(ae.m.df.T5$Timepoint)
ae.m.df.T5$Timepoint <- as.numeric(ae.m.df.T5$Timepoint)
ae.T5 <- ggplot(ae.m.df.T5, aes(x=Timepoint, y=fit, ymin=lower, ymax=upper))+
geom_point(aes(color = Group))+
geom_path(aes(x=Timepoint, y=fit, color=Group))+
geom_ribbon(aes(ymin=lower, ymax=upper, fill=Group), alpha=0.2)+
ylab("Mean Theta Power (n.u.)")+
scale_x_continuous(name="Timpoint", breaks=c(1,2,3,4,5,6,7,8,9,10) , limits=c(1,10))+
ggtitle("muTheta Group x Session x Timepoint x RoI")+
theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank(),panel.background = element_blank(), axis.line = element_line(colour = "black"))+
facet_grid(RoI ~ Session)
plot(ae.T5)
#Modelfit checking No.2.
AIC(m.T0, m.T1, m.T2, m.T2.5, m.T3, m.T4, m.T5)## df AIC
## m.T0 56 213031.6
## m.T1 57 212806.7
## m.T2 61 209397.3
## m.T2.5 58 212073.7
## m.T3 65 208354.3
## m.T4 76 207423.0
## m.T5 273 212613.6####Posthocs? I don’t think we need it as we should rely on plots.###
#Method 1: To get summary of posthocs, simple but not the most correct since they are only conparing against Control. #summary(glht(m.ThetaRanSS1, linfct = mcp(Group = “Tukey”)), test = adjusted(“holm”)) #lsmeans(m.ThetaRanSS1, pairwise~Group, adjust=“tukey”) #emmeans(m.ThetaRanSS1, pairwise~Group)
#summary(glht(m.ThetaRanSS2, linfct = mcp(Group = “Tukey”)), test = adjusted(“holm”)) #lsmeans(m.ThetaRanSS2, pairwise~GroupSession, adjust=“tukey”) #emmeans(m.ThetaRanSS2, pairwise~GroupSession)
#lsmeans(m.ThetaRanSS3, pairwise~GroupSessionTimepoint, adjust=“tukey”) #emmeans(m.ThetaRanSS3, pairwise~GroupSessionTimepoint)
#lsmeans(m.ThetaRanSS4, pairwise~GroupSessionTimepointRoI, adjust=“tukey”) #emmeans(m.ThetaRanSS3, pairwise~GroupSession*Timepoint)
#Method 2: To get summary of posthocs - This is the more correct way. #emmeans(MODEL_NAME, pairwise ~ size) #Overall emmeans -> Hasty and not correct #Create new interaction matrices #EG. emm_NAME <- emmeans(MODEL_NAME, pairwise ~ size | type) #The last term is the conditioning on type.