Model 1
library(MASS)
AM.lm1 = lm(AMefficiency ~ AMcorrect * AMavRT * AMincorrect * AMsdRT * AMmdRT *
AMresponses, data = AM, na.action = na.omit)
summary(AM.lm1)$coefficients
## Estimate Std. Error t value
## (Intercept) 2.878e+02 1.080e+02 2.663453
## AMcorrect 4.190e+01 1.163e+00 36.014848
## AMavRT -1.281e-01 6.598e-02 -1.941793
## AMsdRT -2.869e-02 3.567e-02 -0.804402
## AMmdRT -4.090e-02 5.398e-02 -0.757652
## AMcorrect:AMavRT -3.848e-03 7.416e-04 -5.189143
## AMcorrect:AMincorrect -1.448e+00 4.663e-01 -3.105346
## AMcorrect:AMsdRT -7.682e-05 5.208e-04 -0.147502
## AMavRT:AMsdRT 1.160e-05 1.045e-05 1.110681
## AMcorrect:AMmdRT -1.725e-03 6.571e-04 -2.624583
## AMavRT:AMmdRT 2.072e-05 1.097e-05 1.887880
## AMsdRT:AMmdRT 1.422e-05 2.006e-05 0.708775
## AMcorrect:AMavRT:AMincorrect 6.132e-04 2.767e-04 2.216352
## AMcorrect:AMavRT:AMsdRT -3.000e-08 9.479e-08 -0.316452
## AMcorrect:AMincorrect:AMsdRT 9.640e-05 1.575e-04 0.612064
## AMcorrect:AMavRT:AMmdRT 3.283e-07 1.157e-07 2.836585
## AMcorrect:AMincorrect:AMmdRT 2.336e-04 2.299e-04 1.016032
## AMcorrect:AMsdRT:AMmdRT 8.279e-08 2.063e-07 0.401377
## AMavRT:AMsdRT:AMmdRT -2.900e-09 1.438e-09 -2.016328
## AMcorrect:AMavRT:AMincorrect:AMsdRT -4.257e-08 4.651e-08 -0.915403
## AMcorrect:AMavRT:AMincorrect:AMmdRT -1.022e-07 4.680e-08 -2.184967
## AMcorrect:AMavRT:AMsdRT:AMmdRT -9.485e-14 1.401e-11 -0.006769
## AMcorrect:AMincorrect:AMsdRT:AMmdRT -6.804e-08 8.854e-08 -0.768410
## AMcorrect:AMavRT:AMincorrect:AMsdRT:AMmdRT 1.255e-11 6.488e-12 1.934619
## Pr(>|t|)
## (Intercept) 8.040e-03
## AMcorrect 6.443e-129
## AMavRT 5.285e-02
## AMsdRT 4.216e-01
## AMmdRT 4.491e-01
## AMcorrect:AMavRT 3.333e-07
## AMcorrect:AMincorrect 2.033e-03
## AMcorrect:AMsdRT 8.828e-01
## AMavRT:AMsdRT 2.674e-01
## AMcorrect:AMmdRT 9.000e-03
## AMavRT:AMmdRT 5.975e-02
## AMsdRT:AMmdRT 4.789e-01
## AMcorrect:AMavRT:AMincorrect 2.722e-02
## AMcorrect:AMavRT:AMsdRT 7.518e-01
## AMcorrect:AMincorrect:AMsdRT 5.408e-01
## AMcorrect:AMavRT:AMmdRT 4.787e-03
## AMcorrect:AMincorrect:AMmdRT 3.102e-01
## AMcorrect:AMsdRT:AMmdRT 6.884e-01
## AMavRT:AMsdRT:AMmdRT 4.442e-02
## AMcorrect:AMavRT:AMincorrect:AMsdRT 3.605e-01
## AMcorrect:AMavRT:AMincorrect:AMmdRT 2.946e-02
## AMcorrect:AMavRT:AMsdRT:AMmdRT 9.946e-01
## AMcorrect:AMincorrect:AMsdRT:AMmdRT 4.427e-01
## AMcorrect:AMavRT:AMincorrect:AMsdRT:AMmdRT 5.373e-02
summary(AM.lm1)$r.squared
## [1] 0.9969
# step <- stepAIC(AM.lm1, direction='both')
Model 2
AM.lm2 = lm(AMefficiency ~ AMcorrect * AMavRT * AMincorrect * AMsdRT * AMmdRT,
data = AM, na.action = na.omit)
summary(AM.lm2)$coefficients
## Estimate Std. Error t value
## (Intercept) 2.878e+02 1.080e+02 2.663453
## AMcorrect 4.190e+01 1.163e+00 36.014848
## AMavRT -1.281e-01 6.598e-02 -1.941793
## AMsdRT -2.869e-02 3.567e-02 -0.804402
## AMmdRT -4.090e-02 5.398e-02 -0.757652
## AMcorrect:AMavRT -3.848e-03 7.416e-04 -5.189143
## AMcorrect:AMincorrect -1.448e+00 4.663e-01 -3.105346
## AMcorrect:AMsdRT -7.682e-05 5.208e-04 -0.147502
## AMavRT:AMsdRT 1.160e-05 1.045e-05 1.110681
## AMcorrect:AMmdRT -1.725e-03 6.571e-04 -2.624583
## AMavRT:AMmdRT 2.072e-05 1.097e-05 1.887880
## AMsdRT:AMmdRT 1.422e-05 2.006e-05 0.708775
## AMcorrect:AMavRT:AMincorrect 6.132e-04 2.767e-04 2.216352
## AMcorrect:AMavRT:AMsdRT -3.000e-08 9.479e-08 -0.316452
## AMcorrect:AMincorrect:AMsdRT 9.640e-05 1.575e-04 0.612064
## AMcorrect:AMavRT:AMmdRT 3.283e-07 1.157e-07 2.836585
## AMcorrect:AMincorrect:AMmdRT 2.336e-04 2.299e-04 1.016032
## AMcorrect:AMsdRT:AMmdRT 8.279e-08 2.063e-07 0.401377
## AMavRT:AMsdRT:AMmdRT -2.900e-09 1.438e-09 -2.016328
## AMcorrect:AMavRT:AMincorrect:AMsdRT -4.257e-08 4.651e-08 -0.915403
## AMcorrect:AMavRT:AMincorrect:AMmdRT -1.022e-07 4.680e-08 -2.184967
## AMcorrect:AMavRT:AMsdRT:AMmdRT -9.485e-14 1.401e-11 -0.006769
## AMcorrect:AMincorrect:AMsdRT:AMmdRT -6.804e-08 8.854e-08 -0.768410
## AMcorrect:AMavRT:AMincorrect:AMsdRT:AMmdRT 1.255e-11 6.488e-12 1.934619
## Pr(>|t|)
## (Intercept) 8.040e-03
## AMcorrect 6.443e-129
## AMavRT 5.285e-02
## AMsdRT 4.216e-01
## AMmdRT 4.491e-01
## AMcorrect:AMavRT 3.333e-07
## AMcorrect:AMincorrect 2.033e-03
## AMcorrect:AMsdRT 8.828e-01
## AMavRT:AMsdRT 2.674e-01
## AMcorrect:AMmdRT 9.000e-03
## AMavRT:AMmdRT 5.975e-02
## AMsdRT:AMmdRT 4.789e-01
## AMcorrect:AMavRT:AMincorrect 2.722e-02
## AMcorrect:AMavRT:AMsdRT 7.518e-01
## AMcorrect:AMincorrect:AMsdRT 5.408e-01
## AMcorrect:AMavRT:AMmdRT 4.787e-03
## AMcorrect:AMincorrect:AMmdRT 3.102e-01
## AMcorrect:AMsdRT:AMmdRT 6.884e-01
## AMavRT:AMsdRT:AMmdRT 4.442e-02
## AMcorrect:AMavRT:AMincorrect:AMsdRT 3.605e-01
## AMcorrect:AMavRT:AMincorrect:AMmdRT 2.946e-02
## AMcorrect:AMavRT:AMsdRT:AMmdRT 9.946e-01
## AMcorrect:AMincorrect:AMsdRT:AMmdRT 4.427e-01
## AMcorrect:AMavRT:AMincorrect:AMsdRT:AMmdRT 5.373e-02
summary(AM.lm2)$r.squared
## [1] 0.9969
# step <- stepAIC(AM.lm2, direction='both')
Model 3
AM.lm3 = lm(AMefficiency ~ AMcorrect * AMavRT * AMincorrect, data = AM, na.action = na.omit)
summary(AM.lm3)$coefficients
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -2.121e+02 5.553e+01 -3.819 1.540e-04
## AMcorrect 4.275e+01 8.755e-01 48.831 7.453e-177
## AMavRT 5.089e-02 1.143e-02 4.451 1.094e-05
## AMcorrect:AMavRT -4.022e-03 1.902e-04 -21.141 2.101e-68
## AMcorrect:AMincorrect 4.105e-01 2.362e-01 1.738 8.300e-02
## AMcorrect:AMavRT:AMincorrect -1.059e-04 4.886e-05 -2.167 3.077e-02
summary(AM.lm3)$r.squared
## [1] 0.9645
# step <- stepAIC(AM.lm3, direction='both')
Model 4
AM.lm4 = lm(AMefficiency ~ AMcorrect * AMavRT, data = AM, na.action = na.omit)
summary(AM.lm4)$coefficients
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.191e+02 1.328e+01 -8.966 9.586e-18
## AMcorrect 4.256e+01 8.683e-01 49.012 6.858e-178
## AMavRT 2.689e-02 2.829e-03 9.503 1.461e-19
## AMcorrect:AMavRT -3.969e-03 1.893e-04 -20.967 1.038e-67
summary(AM.lm4)$r.squared
## [1] 0.9641
# step <- stepAIC(AM.lm4, direction='both')
Model 5
AM.lm5 = lm(AMefficiency ~ AMcorrect * AMavRT * AMmdRT, data = AM, na.action = na.omit)
summary(AM.lm5)$coefficients
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -1.245e+02 1.431e+01 -8.701 7.297e-17
## AMcorrect 4.604e+01 9.701e-01 47.458 5.923e-172
## AMavRT 2.288e-02 4.766e-03 4.801 2.189e-06
## AMmdRT 3.376e-02 7.386e-03 4.571 6.379e-06
## AMcorrect:AMavRT -4.328e-03 3.200e-04 -13.525 6.441e-35
## AMcorrect:AMmdRT -3.131e-03 4.895e-04 -6.395 4.234e-10
## AMavRT:AMmdRT -5.400e-06 5.683e-07 -9.502 1.522e-19
## AMcorrect:AMavRT:AMmdRT 5.058e-07 4.226e-08 11.970 1.177e-28
summary(AM.lm5)$r.squared
## [1] 0.9903
# step <- stepAIC(AM.lm5, direction='both')
Model 6
AM.lm6 = lm(AMefficiency ~ AMcorrect * AMavRT * AMsdRT, data = AM, na.action = na.omit)
summary(AM.lm6)$coefficients
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -3.372e+01 1.145e+01 -2.9443 3.415e-03
## AMcorrect 3.973e+01 7.483e-01 53.0972 1.089e-189
## AMavRT 1.176e-02 3.557e-03 3.3076 1.021e-03
## AMsdRT -7.131e-03 3.471e-03 -2.0544 4.055e-02
## AMcorrect:AMavRT -3.582e-03 2.315e-04 -15.4742 3.686e-43
## AMcorrect:AMsdRT -2.948e-04 2.369e-04 -1.2442 2.141e-01
## AMavRT:AMsdRT 2.533e-07 4.129e-07 0.6134 5.400e-01
## AMcorrect:AMavRT:AMsdRT 1.131e-07 2.914e-08 3.8831 1.195e-04
summary(AM.lm6)$r.squared
## [1] 0.9917
# step <- stepAIC(AM.lm6, direction='both')
Model 7
AM.lm7 = lm(AMefficiency ~ AMcorrect * AMsdRT * AMmdRT, data = AM, na.action = na.omit)
summary(AM.lm7)$coefficients
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) -9.479e+01 1.839e+01 -5.155 3.896e-07
## AMcorrect 4.384e+01 1.232e+00 35.580 1.610e-129
## AMsdRT 1.116e-02 4.571e-03 2.442 1.500e-02
## AMmdRT 3.121e-02 6.067e-03 5.144 4.121e-07
## AMcorrect:AMsdRT -2.361e-03 3.206e-04 -7.363 9.391e-13
## AMcorrect:AMmdRT -4.985e-03 4.076e-04 -12.230 1.114e-29
## AMsdRT:AMmdRT -4.045e-06 9.433e-07 -4.289 2.226e-05
## AMcorrect:AMsdRT:AMmdRT 4.747e-07 6.931e-08 6.849 2.623e-11
summary(AM.lm7)$r.squared
## [1] 0.9808
# step <- stepAIC(AM.lm5, direction='both')