EMG of motor learning analysis
Russell Chan
2023-04-26
#Load packages
## Warning: package 'readxl' was built under R version 4.2.2## Warning: package 'haven' was built under R version 4.2.2## Warning: package 'tidyr' was built under R version 4.2.2## Warning: package 'tidyverse' was built under R version 4.2.2## Warning: package 'ggplot2' was built under R version 4.2.2## Warning: package 'readr' was built under R version 4.2.2## Warning: package 'purrr' was built under R version 4.2.2## Warning: package 'dplyr' was built under R version 4.2.2## Warning: package 'stringr' was built under R version 4.2.2## Warning: package 'forcats' was built under R version 4.2.2## Warning: package 'lubridate' was built under R version 4.2.2## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr 1.1.0 ✔ purrr 1.0.1
## ✔ forcats 1.0.0 ✔ readr 2.1.4
## ✔ ggplot2 3.4.1 ✔ stringr 1.5.0
## ✔ lubridate 1.9.2 ✔ tibble 3.1.8
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag() masks stats::lag()
## ℹ Use the ]8;;http://conflicted.r-lib.org/conflicted package]8;; to force all conflicts to become errors## Warning: package 'lme4' was built under R version 4.2.2## Loading required package: Matrix## Warning: package 'Matrix' was built under R version 4.2.2##
## Attaching package: 'Matrix'
##
## The following objects are masked from 'package:tidyr':
##
## expand, pack, unpack## Warning: package 'effects' was built under R version 4.2.2## Loading required package: carData## Warning: package 'carData' was built under R version 4.2.2## lattice theme set by effectsTheme()
## See ?effectsTheme for details.## Warning: package 'car' was built under R version 4.2.2##
## Attaching package: 'car'
##
## The following object is masked from 'package:dplyr':
##
## recode
##
## The following object is masked from 'package:purrr':
##
## some## Warning: package 'knitr' was built under R version 4.2.2## Warning: package 'reshape2' was built under R version 4.2.2##
## Attaching package: 'reshape2'
##
## The following object is masked from 'package:tidyr':
##
## smiths## Warning: package 'DHARMa' was built under R version 4.2.2## This is DHARMa 0.4.6. For overview type '?DHARMa'. For recent changes, type news(package = 'DHARMa')## Warning: package 'phia' was built under R version 4.2.2## Warning: package 'lsmeans' was built under R version 4.2.2## Loading required package: emmeans## Warning: package 'emmeans' was built under R version 4.2.2## The 'lsmeans' package is now basically a front end for 'emmeans'.
## Users are encouraged to switch the rest of the way.
## See help('transition') for more information, including how to
## convert old 'lsmeans' objects and scripts to work with 'emmeans'.## Warning: package 'multcomp' was built under R version 4.2.2## Loading required package: mvtnorm
## Loading required package: survival
## Loading required package: TH.data## Warning: package 'TH.data' was built under R version 4.2.2## Loading required package: MASS
##
## Attaching package: 'MASS'
##
## The following object is masked from 'package:dplyr':
##
## select
##
##
## Attaching package: 'TH.data'
##
## The following object is masked from 'package:MASS':
##
## geyser
##
##
## Attaching package: 'nlme'
##
## The following object is masked from 'package:lme4':
##
## lmList
##
## The following object is masked from 'package:dplyr':
##
## collapse## Warning: package 'gplots' was built under R version 4.2.2##
## Attaching package: 'gplots'
##
## The following object is masked from 'package:stats':
##
## lowess## Warning: package 'htmltools' was built under R version 4.2.2## Warning: package 'Hmisc' was built under R version 4.2.2## Loading required package: Formula## Warning: package 'Formula' was built under R version 4.2.2##
## Attaching package: 'Hmisc'
##
## The following objects are masked from 'package:dplyr':
##
## src, summarize
##
## The following objects are masked from 'package:base':
##
## format.pval, units#Load dataframes
Behav <- read_excel("TIMEERROR.xlsx")
EMG <- read_excel("COACTIVATION TABLE 3.xlsx")#Factors for EMG & Behav.
#Factors for Behav
Behav$PARTICIPANT <- factor(Behav$PARTICIPANT)
Behav$BLOCK <- factor(Behav$BLOCK)
Behav$CONDITION <- factor(Behav$CONDITION)
#Factors for EMG
EMG$PARTICIPANT <- factor(EMG$PARTICIPANT)
EMG$BLOCK <- factor(EMG$BLOCK)
EMG$PHASE <- factor(EMG$PHASE)
EMG$CONDITION <- factor(EMG$CONDITION)Model Behav
#RT
m.RT <- lmer(TIME ~ BLOCK * CONDITION + (1|PARTICIPANT), data = Behav)
Anova(m.RT)## Analysis of Deviance Table (Type II Wald chisquare tests)
##
## Response: TIME
## Chisq Df Pr(>Chisq)
## BLOCK 143.7178 9 <2e-16 ***
## CONDITION 0.1392 1 0.7090
## BLOCK:CONDITION 4.8662 9 0.8458
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1summary(m.RT)## Linear mixed model fit by REML ['lmerMod']
## Formula: TIME ~ BLOCK * CONDITION + (1 | PARTICIPANT)
## Data: Behav
##
## REML criterion at convergence: 3083.8
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -4.2352 -0.4917 -0.0123 0.4429 5.5761
##
## Random effects:
## Groups Name Variance Std.Dev.
## PARTICIPANT (Intercept) 7783 88.22
## Residual 2028 45.03
## Number of obs: 300, groups: PARTICIPANT, 30
##
## Fixed effects:
## Estimate Std. Error t value
## (Intercept) 320.88 25.57 12.547
## BLOCK2 -48.24 16.44 -2.934
## BLOCK3 -52.00 16.44 -3.162
## BLOCK4 -59.36 16.44 -3.610
## BLOCK5 -65.48 16.44 -3.982
## BLOCK6 -48.92 16.44 -2.975
## BLOCK7 -83.96 16.44 -5.106
## BLOCK8 -85.76 16.44 -5.215
## BLOCK9 -102.76 16.44 -6.249
## BLOCK10 -102.52 16.44 -6.234
## CONDITION2 12.72 36.17 0.352
## BLOCK2:CONDITION2 24.20 23.26 1.041
## BLOCK3:CONDITION2 12.32 23.26 0.530
## BLOCK4:CONDITION2 0.96 23.26 0.041
## BLOCK5:CONDITION2 -2.84 23.26 -0.122
## BLOCK6:CONDITION2 -18.76 23.26 -0.807
## BLOCK7:CONDITION2 -9.96 23.26 -0.428
## BLOCK8:CONDITION2 -5.28 23.26 -0.227
## BLOCK9:CONDITION2 3.08 23.26 0.132
## BLOCK10:CONDITION2 -9.16 23.26 -0.394##
## Correlation matrix not shown by default, as p = 20 > 12.
## Use print(x, correlation=TRUE) or
## vcov(x) if you need it#Errors
m.error <- lmer(ERRORS ~ BLOCK * CONDITION + (1|PARTICIPANT), data = Behav)
Anova(m.error)## Analysis of Deviance Table (Type II Wald chisquare tests)
##
## Response: ERRORS
## Chisq Df Pr(>Chisq)
## BLOCK 247.1471 9 <2e-16 ***
## CONDITION 0.5117 1 0.4744
## BLOCK:CONDITION 5.2798 9 0.8093
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1summary(m.error)## Linear mixed model fit by REML ['lmerMod']
## Formula: ERRORS ~ BLOCK * CONDITION + (1 | PARTICIPANT)
## Data: Behav
##
## REML criterion at convergence: 2926.2
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -3.0335 -0.5440 -0.0325 0.5404 3.6496
##
## Random effects:
## Groups Name Variance Std.Dev.
## PARTICIPANT (Intercept) 2664 51.61
## Residual 1220 34.93
## Number of obs: 300, groups: PARTICIPANT, 30
##
## Fixed effects:
## Estimate Std. Error t value
## (Intercept) 239.000 16.091 14.853
## BLOCK2 -51.267 12.754 -4.020
## BLOCK3 -65.600 12.754 -5.144
## BLOCK4 -71.600 12.754 -5.614
## BLOCK5 -74.200 12.754 -5.818
## BLOCK6 -97.933 12.754 -7.679
## BLOCK7 -97.267 12.754 -7.626
## BLOCK8 -107.400 12.754 -8.421
## BLOCK9 -108.867 12.754 -8.536
## BLOCK10 -113.467 12.754 -8.897
## CONDITION2 -26.333 22.756 -1.157
## BLOCK2:CONDITION2 23.867 18.037 1.323
## BLOCK3:CONDITION2 12.867 18.037 0.713
## BLOCK4:CONDITION2 19.600 18.037 1.087
## BLOCK5:CONDITION2 3.333 18.037 0.185
## BLOCK6:CONDITION2 3.133 18.037 0.174
## BLOCK7:CONDITION2 1.133 18.037 0.063
## BLOCK8:CONDITION2 18.267 18.037 1.013
## BLOCK9:CONDITION2 24.933 18.037 1.382
## BLOCK10:CONDITION2 18.333 18.037 1.016##
## Correlation matrix not shown by default, as p = 20 > 12.
## Use print(x, correlation=TRUE) or
## vcov(x) if you need it#Check if errors is affecting RT, using simple LM not really needed as below is a mixed model.
m.errorRT <- lm(TIME ~ ERRORS, data = Behav)
Anova(m.errorRT)## Anova Table (Type II tests)
##
## Response: TIME
## Sum Sq Df F value Pr(>F)
## ERRORS 112897 1 11.418 0.0008244 ***
## Residuals 2946633 298
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1summary(m.errorRT)##
## Call:
## lm(formula = TIME ~ ERRORS, data = Behav)
##
## Residuals:
## Min 1Q Median 3Q Max
## -184.86 -66.60 -8.38 56.36 452.27
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 218.67161 14.06782 15.544 < 2e-16 ***
## ERRORS 0.28300 0.08375 3.379 0.000824 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 99.44 on 298 degrees of freedom
## Multiple R-squared: 0.0369, Adjusted R-squared: 0.03367
## F-statistic: 11.42 on 1 and 298 DF, p-value: 0.0008244plot(m.errorRT)
#RT * Error model
m.RT.ER <- lmer(TIME ~ BLOCK * CONDITION * ERRORS + (1|PARTICIPANT), data = Behav)
Anova(m.RT.ER)## Analysis of Deviance Table (Type II Wald chisquare tests)
##
## Response: TIME
## Chisq Df Pr(>Chisq)
## BLOCK 53.1673 9 2.717e-08 ***
## CONDITION 0.2282 1 0.632872
## ERRORS 11.0919 1 0.000867 ***
## BLOCK:CONDITION 2.3138 9 0.985480
## BLOCK:ERRORS 1.4964 9 0.997173
## CONDITION:ERRORS 0.9029 1 0.342009
## BLOCK:CONDITION:ERRORS 7.0006 9 0.637062
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1summary(m.RT.ER)## Linear mixed model fit by REML ['lmerMod']
## Formula: TIME ~ BLOCK * CONDITION * ERRORS + (1 | PARTICIPANT)
## Data: Behav
##
## REML criterion at convergence: 3086.3
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -4.2203 -0.4645 -0.0056 0.4070 5.1291
##
## Random effects:
## Groups Name Variance Std.Dev.
## PARTICIPANT (Intercept) 7950 89.16
## Residual 2017 44.91
## Number of obs: 300, groups: PARTICIPANT, 30
##
## Fixed effects:
## Estimate Std. Error t value
## (Intercept) 254.41268 49.37244 5.153
## BLOCK2 -4.59144 50.59842 -0.091
## BLOCK3 10.28548 52.39018 0.196
## BLOCK4 -1.10420 48.98079 -0.023
## BLOCK5 -34.23324 49.05395 -0.698
## BLOCK6 -9.59480 46.04440 -0.208
## BLOCK7 -55.39388 49.83086 -1.112
## BLOCK8 -78.18337 57.94078 -1.349
## BLOCK9 -105.66953 51.44698 -2.054
## BLOCK10 -80.78868 54.07831 -1.494
## CONDITION2 8.85264 63.66348 0.139
## ERRORS 0.27811 0.17619 1.578
## BLOCK2:CONDITION2 -7.03963 65.99368 -0.107
## BLOCK3:CONDITION2 -51.73719 69.46090 -0.745
## BLOCK4:CONDITION2 -77.59512 68.41105 -1.134
## BLOCK5:CONDITION2 -15.47391 68.80140 -0.225
## BLOCK6:CONDITION2 -27.41962 66.59408 -0.412
## BLOCK7:CONDITION2 23.25460 68.22655 0.341
## BLOCK8:CONDITION2 18.16966 72.87285 0.249
## BLOCK9:CONDITION2 36.49737 68.36537 0.534
## BLOCK10:CONDITION2 38.80660 72.57065 0.535
## BLOCK2:ERRORS -0.15656 0.22481 -0.696
## BLOCK3:ERRORS -0.25399 0.24708 -1.028
## BLOCK4:ERRORS -0.22905 0.22618 -1.013
## BLOCK5:ERRORS -0.06439 0.22991 -0.280
## BLOCK6:ERRORS -0.08570 0.21425 -0.400
## BLOCK7:ERRORS -0.01069 0.25972 -0.041
## BLOCK8:ERRORS 0.16939 0.35828 0.473
## BLOCK9:ERRORS 0.25502 0.29512 0.864
## BLOCK10:ERRORS 0.07826 0.33196 0.236
## CONDITION2:ERRORS 0.05262 0.22819 0.231
## BLOCK2:CONDITION2:ERRORS 0.13849 0.29928 0.463
## BLOCK3:CONDITION2:ERRORS 0.37411 0.34528 1.084
## BLOCK4:CONDITION2:ERRORS 0.46244 0.34223 1.351
## BLOCK5:CONDITION2:ERRORS 0.09841 0.36769 0.268
## BLOCK6:CONDITION2:ERRORS 0.09153 0.39747 0.230
## BLOCK7:CONDITION2:ERRORS -0.24663 0.41053 -0.601
## BLOCK8:CONDITION2:ERRORS -0.18192 0.45413 -0.401
## BLOCK9:CONDITION2:ERRORS -0.27637 0.40321 -0.685
## BLOCK10:CONDITION2:ERRORS -0.40357 0.47185 -0.855##
## Correlation matrix not shown by default, as p = 40 > 12.
## Use print(x, correlation=TRUE) or
## vcov(x) if you need it#I only see a Block * Time and BLock * Errors interaction. and later then I run the 3-way I see a Time and Errors interaction. I don't see any significant interactions by condition.
#Here I will just perform the emmeans on the 1st model
emmeans(m.RT, pairwise ~ CONDITION | BLOCK)## $emmeans
## BLOCK = 1:
## CONDITION emmean SE df lower.CL upper.CL
## 1 321 25.6 42 269 372
## 2 334 25.6 42 282 385
##
## BLOCK = 2:
## CONDITION emmean SE df lower.CL upper.CL
## 1 273 25.6 42 221 324
## 2 310 25.6 42 258 361
##
## BLOCK = 3:
## CONDITION emmean SE df lower.CL upper.CL
## 1 269 25.6 42 217 320
## 2 294 25.6 42 242 346
##
## BLOCK = 4:
## CONDITION emmean SE df lower.CL upper.CL
## 1 262 25.6 42 210 313
## 2 275 25.6 42 224 327
##
## BLOCK = 5:
## CONDITION emmean SE df lower.CL upper.CL
## 1 255 25.6 42 204 307
## 2 265 25.6 42 214 317
##
## BLOCK = 6:
## CONDITION emmean SE df lower.CL upper.CL
## 1 272 25.6 42 220 324
## 2 266 25.6 42 214 318
##
## BLOCK = 7:
## CONDITION emmean SE df lower.CL upper.CL
## 1 237 25.6 42 185 289
## 2 240 25.6 42 188 291
##
## BLOCK = 8:
## CONDITION emmean SE df lower.CL upper.CL
## 1 235 25.6 42 184 287
## 2 243 25.6 42 191 294
##
## BLOCK = 9:
## CONDITION emmean SE df lower.CL upper.CL
## 1 218 25.6 42 167 270
## 2 234 25.6 42 182 286
##
## BLOCK = 10:
## CONDITION emmean SE df lower.CL upper.CL
## 1 218 25.6 42 167 270
## 2 222 25.6 42 170 274
##
## Degrees-of-freedom method: kenward-roger
## Confidence level used: 0.95
##
## $contrasts
## BLOCK = 1:
## contrast estimate SE df t.ratio p.value
## CONDITION1 - CONDITION2 -12.72 36.2 42 -0.352 0.7268
##
## BLOCK = 2:
## contrast estimate SE df t.ratio p.value
## CONDITION1 - CONDITION2 -36.92 36.2 42 -1.021 0.3132
##
## BLOCK = 3:
## contrast estimate SE df t.ratio p.value
## CONDITION1 - CONDITION2 -25.04 36.2 42 -0.692 0.4925
##
## BLOCK = 4:
## contrast estimate SE df t.ratio p.value
## CONDITION1 - CONDITION2 -13.68 36.2 42 -0.378 0.7072
##
## BLOCK = 5:
## contrast estimate SE df t.ratio p.value
## CONDITION1 - CONDITION2 -9.88 36.2 42 -0.273 0.7861
##
## BLOCK = 6:
## contrast estimate SE df t.ratio p.value
## CONDITION1 - CONDITION2 6.04 36.2 42 0.167 0.8682
##
## BLOCK = 7:
## contrast estimate SE df t.ratio p.value
## CONDITION1 - CONDITION2 -2.76 36.2 42 -0.076 0.9395
##
## BLOCK = 8:
## contrast estimate SE df t.ratio p.value
## CONDITION1 - CONDITION2 -7.44 36.2 42 -0.206 0.8380
##
## BLOCK = 9:
## contrast estimate SE df t.ratio p.value
## CONDITION1 - CONDITION2 -15.80 36.2 42 -0.437 0.6645
##
## BLOCK = 10:
## contrast estimate SE df t.ratio p.value
## CONDITION1 - CONDITION2 -3.56 36.2 42 -0.098 0.9221
##
## Degrees-of-freedom method: kenward-roger#You can see of course results corraborates that we don't have a condition and Block intereaction
#then if you just check on the block level.
emmeans(m.RT, pairwise ~ BLOCK)## NOTE: Results may be misleading due to involvement in interactions## $emmeans
## BLOCK emmean SE df lower.CL upper.CL
## 1 327 18.1 42 291 364
## 2 291 18.1 42 255 328
## 3 281 18.1 42 245 318
## 4 268 18.1 42 232 305
## 5 260 18.1 42 224 297
## 6 269 18.1 42 232 305
## 7 238 18.1 42 202 275
## 8 239 18.1 42 202 275
## 9 226 18.1 42 190 263
## 10 220 18.1 42 184 257
##
## Results are averaged over the levels of: CONDITION
## Degrees-of-freedom method: kenward-roger
## Confidence level used: 0.95
##
## $contrasts
## contrast estimate SE df t.ratio p.value
## BLOCK1 - BLOCK2 36.14 11.6 252 3.108 0.0636
## BLOCK1 - BLOCK3 45.84 11.6 252 3.942 0.0041
## BLOCK1 - BLOCK4 58.88 11.6 252 5.064 <.0001
## BLOCK1 - BLOCK5 66.90 11.6 252 5.753 <.0001
## BLOCK1 - BLOCK6 58.30 11.6 252 5.014 <.0001
## BLOCK1 - BLOCK7 88.94 11.6 252 7.649 <.0001
## BLOCK1 - BLOCK8 88.40 11.6 252 7.603 <.0001
## BLOCK1 - BLOCK9 101.22 11.6 252 8.705 <.0001
## BLOCK1 - BLOCK10 107.10 11.6 252 9.211 <.0001
## BLOCK2 - BLOCK3 9.70 11.6 252 0.834 0.9980
## BLOCK2 - BLOCK4 22.74 11.6 252 1.956 0.6308
## BLOCK2 - BLOCK5 30.76 11.6 252 2.645 0.2024
## BLOCK2 - BLOCK6 22.16 11.6 252 1.906 0.6652
## BLOCK2 - BLOCK7 52.80 11.6 252 4.541 0.0004
## BLOCK2 - BLOCK8 52.26 11.6 252 4.494 0.0004
## BLOCK2 - BLOCK9 65.08 11.6 252 5.597 <.0001
## BLOCK2 - BLOCK10 70.96 11.6 252 6.103 <.0001
## BLOCK3 - BLOCK4 13.04 11.6 252 1.121 0.9821
## BLOCK3 - BLOCK5 21.06 11.6 252 1.811 0.7276
## BLOCK3 - BLOCK6 12.46 11.6 252 1.072 0.9870
## BLOCK3 - BLOCK7 43.10 11.6 252 3.707 0.0096
## BLOCK3 - BLOCK8 42.56 11.6 252 3.660 0.0112
## BLOCK3 - BLOCK9 55.38 11.6 252 4.763 0.0001
## BLOCK3 - BLOCK10 61.26 11.6 252 5.268 <.0001
## BLOCK4 - BLOCK5 8.02 11.6 252 0.690 0.9996
## BLOCK4 - BLOCK6 -0.58 11.6 252 -0.050 1.0000
## BLOCK4 - BLOCK7 30.06 11.6 252 2.585 0.2302
## BLOCK4 - BLOCK8 29.52 11.6 252 2.539 0.2532
## BLOCK4 - BLOCK9 42.34 11.6 252 3.641 0.0120
## BLOCK4 - BLOCK10 48.22 11.6 252 4.147 0.0019
## BLOCK5 - BLOCK6 -8.60 11.6 252 -0.740 0.9992
## BLOCK5 - BLOCK7 22.04 11.6 252 1.895 0.6722
## BLOCK5 - BLOCK8 21.50 11.6 252 1.849 0.7031
## BLOCK5 - BLOCK9 34.32 11.6 252 2.952 0.0972
## BLOCK5 - BLOCK10 40.20 11.6 252 3.457 0.0221
## BLOCK6 - BLOCK7 30.64 11.6 252 2.635 0.2070
## BLOCK6 - BLOCK8 30.10 11.6 252 2.589 0.2285
## BLOCK6 - BLOCK9 42.92 11.6 252 3.691 0.0101
## BLOCK6 - BLOCK10 48.80 11.6 252 4.197 0.0015
## BLOCK7 - BLOCK8 -0.54 11.6 252 -0.046 1.0000
## BLOCK7 - BLOCK9 12.28 11.6 252 1.056 0.9883
## BLOCK7 - BLOCK10 18.16 11.6 252 1.562 0.8651
## BLOCK8 - BLOCK9 12.82 11.6 252 1.103 0.9841
## BLOCK8 - BLOCK10 18.70 11.6 252 1.608 0.8432
## BLOCK9 - BLOCK10 5.88 11.6 252 0.506 1.0000
##
## Results are averaged over the levels of: CONDITION
## Degrees-of-freedom method: kenward-roger
## P value adjustment: tukey method for comparing a family of 10 estimates#You do see that there are significant differences between the blocks, in that block 10 is faster than block 1.Effects and Plots
#Effects just on the RT model
ae.m.RT<-allEffects(m.RT)
ae.m.RT.df<-as.data.frame(ae.m.RT[[1]])
plot(ae.m.RT)
#The plot
plot.m.RT <- ggplot(ae.m.RT.df, aes(x=BLOCK, y=fit, group=CONDITION))+
geom_ribbon(aes(ymin=lower, ymax=upper, fill=CONDITION), alpha = 0.2) +
geom_line(aes(color=CONDITION), linewidth = 1) +
geom_point(aes(color=CONDITION, shape=CONDITION), size = 3.5)+
ylab("Response Time (s)")+
xlab("Blocks")+
ggtitle("Response time between groups")+
theme_classic()
#facet_wrap(~BLOCK * CONDITION) not needed here
print (plot.m.RT)
EMG
#Coactivation
m.coact <- lmer(COACTIVATION.INDEX ~ BLOCK * CONDITION * PHASE + (1|PARTICIPANT), data = EMG)
Anova(m.coact)## Analysis of Deviance Table (Type II Wald chisquare tests)
##
## Response: COACTIVATION.INDEX
## Chisq Df Pr(>Chisq)
## BLOCK 14.6702 9 0.1004
## CONDITION 0.9161 1 0.3385
## PHASE 70.4279 9 1.255e-11 ***
## BLOCK:CONDITION 14.5032 9 0.1055
## BLOCK:PHASE 72.5833 81 0.7366
## CONDITION:PHASE 8.6894 9 0.4664
## BLOCK:CONDITION:PHASE 73.1454 81 0.7209
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1summary(m.coact)## Linear mixed model fit by REML ['lmerMod']
## Formula: COACTIVATION.INDEX ~ BLOCK * CONDITION * PHASE + (1 | PARTICIPANT)
## Data: EMG
##
## REML criterion at convergence: 133291.6
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -4.762 -0.008 0.004 0.028 104.872
##
## Random effects:
## Groups Name Variance Std.Dev.
## PARTICIPANT (Intercept) 4.441 2.107
## Residual 610.271 24.704
## Number of obs: 14510, groups: PARTICIPANT, 30
##
## Fixed effects:
## Estimate Std. Error t value
## (Intercept) 1.013e+02 2.923e+00 34.656
## BLOCK2 -1.324e+00 4.075e+00 -0.325
## BLOCK3 1.287e-01 4.061e+00 0.032
## BLOCK4 2.212e-04 4.048e+00 0.000
## BLOCK5 2.130e-02 4.048e+00 0.005
## BLOCK6 -5.560e-02 4.061e+00 -0.014
## BLOCK7 1.296e-01 4.048e+00 0.032
## BLOCK8 2.893e-02 4.075e+00 0.007
## BLOCK9 -4.285e-02 4.048e+00 -0.011
## BLOCK10 -4.690e-02 4.048e+00 -0.012
## CONDITION1 -1.313e-01 4.147e+00 -0.032
## PHASE20 1.224e-01 4.089e+00 0.030
## PHASE30 -2.582e+00 4.048e+00 -0.638
## PHASE40 -2.818e-02 4.048e+00 -0.007
## PHASE50 -2.682e+00 4.075e+00 -0.658
## PHASE60 -4.090e-02 4.061e+00 -0.010
## PHASE70 -2.575e+00 4.075e+00 -0.632
## PHASE80 2.106e-01 4.061e+00 0.052
## PHASE90 -2.645e+00 4.075e+00 -0.649
## PHASE100 -2.566e+00 4.312e+00 -0.595
## BLOCK2:CONDITION1 1.141e+00 5.763e+00 0.198
## BLOCK3:CONDITION1 6.243e-01 5.753e+00 0.109
## BLOCK4:CONDITION1 -3.968e-01 5.734e+00 -0.069
## BLOCK5:CONDITION1 3.054e+00 5.734e+00 0.533
## BLOCK6:CONDITION1 -4.002e-01 5.763e+00 -0.069
## BLOCK7:CONDITION1 1.041e+00 5.744e+00 0.181
## BLOCK8:CONDITION1 6.881e-01 5.754e+00 0.120
## BLOCK9:CONDITION1 3.550e+01 5.754e+00 6.170
## BLOCK10:CONDITION1 -6.719e-01 5.734e+00 -0.117
## BLOCK2:PHASE20 1.143e+00 5.783e+00 0.198
## BLOCK3:PHASE20 -8.556e-02 5.754e+00 -0.015
## BLOCK4:PHASE20 -2.209e-01 5.744e+00 -0.038
## BLOCK5:PHASE20 -1.547e-01 5.744e+00 -0.027
## BLOCK6:PHASE20 7.978e-02 5.754e+00 0.014
## BLOCK7:PHASE20 -2.896e-01 5.744e+00 -0.050
## BLOCK8:PHASE20 -1.845e-01 5.783e+00 -0.032
## BLOCK9:PHASE20 -7.681e-02 5.744e+00 -0.013
## BLOCK10:PHASE20 -3.063e-02 5.744e+00 -0.005
## BLOCK2:PHASE30 1.516e+00 5.774e+00 0.263
## BLOCK3:PHASE30 -1.021e-01 5.744e+00 -0.018
## BLOCK4:PHASE30 8.491e-03 5.715e+00 0.001
## BLOCK5:PHASE30 6.137e-02 5.715e+00 0.011
## BLOCK6:PHASE30 1.029e-01 5.724e+00 0.018
## BLOCK7:PHASE30 -1.724e-01 5.715e+00 -0.030
## BLOCK8:PHASE30 -2.754e-02 5.754e+00 -0.005
## BLOCK9:PHASE30 5.534e-02 5.715e+00 0.010
## BLOCK10:PHASE30 7.281e-02 5.715e+00 0.013
## BLOCK2:PHASE40 1.575e+00 5.754e+00 0.274
## BLOCK3:PHASE40 -1.028e-01 5.744e+00 -0.018
## BLOCK4:PHASE40 5.188e-02 5.715e+00 0.009
## BLOCK5:PHASE40 2.157e-02 5.715e+00 0.004
## BLOCK6:PHASE40 7.237e-02 5.724e+00 0.013
## BLOCK7:PHASE40 -3.242e-01 5.715e+00 -0.057
## BLOCK8:PHASE40 -6.692e-02 5.754e+00 -0.012
## BLOCK9:PHASE40 4.439e-02 5.715e+00 0.008
## BLOCK10:PHASE40 1.290e-01 5.715e+00 0.023
## BLOCK2:PHASE50 1.221e+00 5.773e+00 0.211
## BLOCK3:PHASE50 -2.256e-02 5.744e+00 -0.004
## BLOCK4:PHASE50 5.689e-02 5.734e+00 0.010
## BLOCK5:PHASE50 2.118e-01 5.734e+00 0.037
## BLOCK6:PHASE50 1.495e-01 5.744e+00 0.026
## BLOCK7:PHASE50 -1.025e-01 5.734e+00 -0.018
## BLOCK8:PHASE50 -4.778e-02 5.783e+00 -0.008
## BLOCK9:PHASE50 1.206e-01 5.734e+00 0.021
## BLOCK10:PHASE50 1.870e-01 5.744e+00 0.033
## BLOCK2:PHASE60 1.400e+00 5.763e+00 0.243
## BLOCK3:PHASE60 -1.241e+00 5.734e+00 -0.216
## BLOCK4:PHASE60 6.348e-03 5.734e+00 0.001
## BLOCK5:PHASE60 1.740e-01 5.724e+00 0.030
## BLOCK6:PHASE60 -6.633e-03 5.734e+00 -0.001
## BLOCK7:PHASE60 -1.375e-01 5.724e+00 -0.024
## BLOCK8:PHASE60 1.809e-01 5.784e+00 0.031
## BLOCK9:PHASE60 9.197e-02 5.724e+00 0.016
## BLOCK10:PHASE60 1.908e-01 5.724e+00 0.033
## BLOCK2:PHASE70 1.318e+00 5.783e+00 0.228
## BLOCK3:PHASE70 -2.017e-01 5.763e+00 -0.035
## BLOCK4:PHASE70 -2.411e-02 5.744e+00 -0.004
## BLOCK5:PHASE70 -1.858e-01 5.734e+00 -0.032
## BLOCK6:PHASE70 1.223e-01 5.753e+00 0.021
## BLOCK7:PHASE70 -3.501e-01 5.754e+00 -0.061
## BLOCK8:PHASE70 -2.006e-02 5.783e+00 -0.003
## BLOCK9:PHASE70 -3.060e-03 5.754e+00 -0.001
## BLOCK10:PHASE70 4.818e-02 5.744e+00 0.008
## BLOCK2:PHASE80 1.170e+00 5.773e+00 0.203
## BLOCK3:PHASE80 -3.920e-01 5.734e+00 -0.068
## BLOCK4:PHASE80 -1.222e-01 5.724e+00 -0.021
## BLOCK5:PHASE80 -2.231e-01 5.724e+00 -0.039
## BLOCK6:PHASE80 -1.218e-01 5.734e+00 -0.021
## BLOCK7:PHASE80 -3.510e-01 5.734e+00 -0.061
## BLOCK8:PHASE80 -2.066e-01 5.773e+00 -0.036
## BLOCK9:PHASE80 -1.899e-01 5.744e+00 -0.033
## BLOCK10:PHASE80 -1.825e-01 5.724e+00 -0.032
## BLOCK2:PHASE90 1.404e+00 5.783e+00 0.243
## BLOCK3:PHASE90 -6.145e-02 5.744e+00 -0.011
## BLOCK4:PHASE90 1.486e-01 5.744e+00 0.026
## BLOCK5:PHASE90 6.287e-02 5.734e+00 0.011
## BLOCK6:PHASE90 1.547e-01 5.744e+00 0.027
## BLOCK7:PHASE90 -4.203e-02 5.734e+00 -0.007
## BLOCK8:PHASE90 -2.283e-02 5.783e+00 -0.004
## BLOCK9:PHASE90 1.401e-01 5.734e+00 0.024
## BLOCK10:PHASE90 9.656e-02 5.744e+00 0.017
## BLOCK2:PHASE100 1.349e+00 6.122e+00 0.220
## BLOCK3:PHASE100 -6.985e-02 6.083e+00 -0.011
## BLOCK4:PHASE100 -6.823e-02 6.074e+00 -0.011
## BLOCK5:PHASE100 3.519e-02 6.074e+00 0.006
## BLOCK6:PHASE100 7.238e-02 6.083e+00 0.012
## BLOCK7:PHASE100 -1.836e-01 6.074e+00 -0.030
## BLOCK8:PHASE100 8.337e-02 6.121e+00 0.014
## BLOCK9:PHASE100 8.290e-02 6.089e+00 0.014
## BLOCK10:PHASE100 6.892e-02 6.074e+00 0.011
## CONDITION1:PHASE20 -7.775e-02 5.773e+00 -0.013
## CONDITION1:PHASE30 2.582e-01 5.744e+00 0.045
## CONDITION1:PHASE40 2.673e-02 5.744e+00 0.005
## CONDITION1:PHASE50 3.930e-01 5.773e+00 0.068
## CONDITION1:PHASE60 1.658e-01 5.763e+00 0.029
## CONDITION1:PHASE70 -3.562e-04 5.773e+00 0.000
## CONDITION1:PHASE80 -2.727e-01 5.763e+00 -0.047
## CONDITION1:PHASE90 2.529e-01 5.793e+00 0.044
## CONDITION1:PHASE100 -1.134e-01 6.107e+00 -0.019
## BLOCK2:CONDITION1:PHASE20 -1.586e+00 8.157e+00 -0.194
## BLOCK3:CONDITION1:PHASE20 -1.141e+00 8.130e+00 -0.140
## BLOCK4:CONDITION1:PHASE20 9.170e-01 8.117e+00 0.113
## BLOCK5:CONDITION1:PHASE20 -3.333e+00 8.123e+00 -0.410
## BLOCK6:CONDITION1:PHASE20 6.669e-01 8.151e+00 0.082
## BLOCK7:CONDITION1:PHASE20 -7.893e-01 8.130e+00 -0.097
## BLOCK8:CONDITION1:PHASE20 5.912e-01 8.151e+00 0.073
## BLOCK9:CONDITION1:PHASE20 -3.053e+01 8.130e+00 -3.756
## BLOCK10:CONDITION1:PHASE20 2.071e+00 8.117e+00 0.255
## BLOCK2:CONDITION1:PHASE30 -1.571e+00 8.151e+00 -0.193
## BLOCK3:CONDITION1:PHASE30 -8.171e-01 8.137e+00 -0.100
## BLOCK4:CONDITION1:PHASE30 2.815e-01 8.096e+00 0.035
## BLOCK5:CONDITION1:PHASE30 -3.396e+00 8.102e+00 -0.419
## BLOCK6:CONDITION1:PHASE30 2.120e-01 8.130e+00 0.026
## BLOCK7:CONDITION1:PHASE30 -8.896e-01 8.116e+00 -0.110
## BLOCK8:CONDITION1:PHASE30 -3.036e+00 8.123e+00 -0.374
## BLOCK9:CONDITION1:PHASE30 -3.723e+01 8.109e+00 -4.591
## BLOCK10:CONDITION1:PHASE30 1.514e+00 8.096e+00 0.187
## BLOCK2:CONDITION1:PHASE40 -1.186e-02 8.137e+00 -0.001
## BLOCK3:CONDITION1:PHASE40 -7.341e-01 8.123e+00 -0.090
## BLOCK4:CONDITION1:PHASE40 3.593e-01 8.102e+00 0.044
## BLOCK5:CONDITION1:PHASE40 -1.491e+00 8.102e+00 -0.184
## BLOCK6:CONDITION1:PHASE40 4.804e-01 8.123e+00 0.059
## BLOCK7:CONDITION1:PHASE40 -2.077e+00 8.123e+00 -0.256
## BLOCK8:CONDITION1:PHASE40 1.455e+01 8.130e+00 1.790
## BLOCK9:CONDITION1:PHASE40 -3.021e+01 8.109e+00 -3.725
## BLOCK10:CONDITION1:PHASE40 -1.237e-01 8.096e+00 -0.015
## BLOCK2:CONDITION1:PHASE50 -1.440e+00 8.164e+00 -0.176
## BLOCK3:CONDITION1:PHASE50 -1.346e+00 8.137e+00 -0.165
## BLOCK4:CONDITION1:PHASE50 2.542e-01 8.123e+00 0.031
## BLOCK5:CONDITION1:PHASE50 -3.531e+00 8.130e+00 -0.434
## BLOCK6:CONDITION1:PHASE50 1.265e-01 8.144e+00 0.016
## BLOCK7:CONDITION1:PHASE50 -1.142e+00 8.130e+00 -0.140
## BLOCK8:CONDITION1:PHASE50 -4.450e-01 8.151e+00 -0.055
## BLOCK9:CONDITION1:PHASE50 -3.632e+01 8.130e+00 -4.468
## BLOCK10:CONDITION1:PHASE50 3.265e+00 8.123e+00 0.402
## BLOCK2:CONDITION1:PHASE60 -1.680e+00 8.157e+00 -0.206
## BLOCK3:CONDITION1:PHASE60 7.170e-01 8.123e+00 0.088
## BLOCK4:CONDITION1:PHASE60 6.228e-01 8.116e+00 0.077
## BLOCK5:CONDITION1:PHASE60 -3.348e+00 8.116e+00 -0.413
## BLOCK6:CONDITION1:PHASE60 2.880e-01 8.137e+00 0.035
## BLOCK7:CONDITION1:PHASE60 -1.546e+00 8.130e+00 -0.190
## BLOCK8:CONDITION1:PHASE60 1.402e+00 8.151e+00 0.172
## BLOCK9:CONDITION1:PHASE60 -3.534e+01 8.123e+00 -4.350
## BLOCK10:CONDITION1:PHASE60 4.316e+00 8.116e+00 0.532
## BLOCK2:CONDITION1:PHASE70 -2.134e+00 8.178e+00 -0.261
## BLOCK3:CONDITION1:PHASE70 -3.171e-01 8.144e+00 -0.039
## BLOCK4:CONDITION1:PHASE70 5.840e-01 8.123e+00 0.072
## BLOCK5:CONDITION1:PHASE70 -2.780e+00 8.130e+00 -0.342
## BLOCK6:CONDITION1:PHASE70 5.234e-01 8.150e+00 0.064
## BLOCK7:CONDITION1:PHASE70 -7.969e-01 8.144e+00 -0.098
## BLOCK8:CONDITION1:PHASE70 -6.945e-01 8.151e+00 -0.085
## BLOCK9:CONDITION1:PHASE70 -3.335e+01 8.151e+00 -4.091
## BLOCK10:CONDITION1:PHASE70 9.537e-01 8.123e+00 0.117
## BLOCK2:CONDITION1:PHASE80 -3.059e-01 8.179e+00 -0.037
## BLOCK3:CONDITION1:PHASE80 -1.618e-01 8.130e+00 -0.020
## BLOCK4:CONDITION1:PHASE80 -8.319e-01 8.123e+00 -0.102
## BLOCK5:CONDITION1:PHASE80 -2.669e+00 8.116e+00 -0.329
## BLOCK6:CONDITION1:PHASE80 -6.322e-01 8.144e+00 -0.078
## BLOCK7:CONDITION1:PHASE80 2.143e+00 8.130e+00 0.264
## BLOCK8:CONDITION1:PHASE80 1.946e+00 8.151e+00 0.239
## BLOCK9:CONDITION1:PHASE80 -2.881e+01 8.144e+00 -3.538
## BLOCK10:CONDITION1:PHASE80 1.636e+00 8.109e+00 0.202
## BLOCK2:CONDITION1:PHASE90 -1.240e+00 8.193e+00 -0.151
## BLOCK3:CONDITION1:PHASE90 -7.183e-01 8.151e+00 -0.088
## BLOCK4:CONDITION1:PHASE90 1.080e-01 8.151e+00 0.013
## BLOCK5:CONDITION1:PHASE90 -3.207e+00 8.145e+00 -0.394
## BLOCK6:CONDITION1:PHASE90 2.377e-01 8.165e+00 0.029
## BLOCK7:CONDITION1:PHASE90 -1.178e+00 8.151e+00 -0.145
## BLOCK8:CONDITION1:PHASE90 8.358e-01 8.172e+00 0.102
## BLOCK9:CONDITION1:PHASE90 -3.697e+01 8.145e+00 -4.540
## BLOCK10:CONDITION1:PHASE90 7.312e-01 8.138e+00 0.090
## BLOCK2:CONDITION1:PHASE100 -4.926e-01 8.647e+00 -0.057
## BLOCK3:CONDITION1:PHASE100 -4.821e-01 8.620e+00 -0.056
## BLOCK4:CONDITION1:PHASE100 6.053e-01 8.597e+00 0.070
## BLOCK5:CONDITION1:PHASE100 -2.908e+00 8.607e+00 -0.338
## BLOCK6:CONDITION1:PHASE100 6.513e-01 8.626e+00 0.076
## BLOCK7:CONDITION1:PHASE100 -7.604e-01 8.635e+00 -0.088
## BLOCK8:CONDITION1:PHASE100 2.345e-01 8.630e+00 0.027
## BLOCK9:CONDITION1:PHASE100 -3.596e+01 8.620e+00 -4.172
## BLOCK10:CONDITION1:PHASE100 1.710e+00 8.597e+00 0.199##
## Correlation matrix not shown by default, as p = 200 > 12.
## Use print(x, correlation=TRUE) or
## vcov(x) if you need it#Hi there are errors in your COACTIVATION file. Basically at the 100% line of every participant if you check your original file, there is a Warning. So that made the Participant into 514 levels which is not true. I suggest you check your original dataframes and how you found those interactions because I didn't find any after I filtered out the warning lines.
#If you fix it then maybe you can continue the same way as above and try and reach the plots as well as above on your own.