Base Model
M1 = glmer(GatesScore ~ zNarrativity + zWordConcreteness + zSyntacticSimplicity + zRefCohesion + zDeepCohesion + zArousal
+ (1|ID),data = Emotion_Data, family = binomial)
summary(M1)
## Generalized linear mixed model fit by maximum likelihood (Laplace
## Approximation) [glmerMod]
## Family: binomial ( logit )
## Formula:
## GatesScore ~ zNarrativity + zWordConcreteness + zSyntacticSimplicity +
## zRefCohesion + zDeepCohesion + zArousal + (1 | ID)
## Data: Emotion_Data
##
## AIC BIC logLik deviance df.resid
## 7828.7 7884.5 -3906.4 7812.7 7912
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -7.3094 -0.5803 0.3305 0.5283 2.7691
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 1.776 1.333
## Number of obs: 7920, groups: ID, 165
##
## Fixed effects:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 1.22822 0.10927 11.240 < 2e-16 ***
## zNarrativity 0.24729 0.04328 5.713 1.11e-08 ***
## zWordConcreteness -0.44924 0.04642 -9.677 < 2e-16 ***
## zSyntacticSimplicity -0.32845 0.05127 -6.406 1.50e-10 ***
## zRefCohesion -0.08313 0.05052 -1.646 0.0999 .
## zDeepCohesion 0.27155 0.03354 8.097 5.65e-16 ***
## zArousal 0.22800 0.04107 5.552 2.83e-08 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) zNrrtv zWrdCn zSyntS zRfChs zDpChs
## zNarrativty 0.025
## zWrdCncrtns -0.039 0.336
## zSyntctcSmp -0.029 -0.064 0.528
## zRefCohesin -0.010 -0.476 -0.031 0.617
## zDeepCohesn 0.031 -0.036 -0.120 -0.243 -0.315
## zArousal 0.021 0.177 -0.199 0.248 0.370 -0.411
Base Model with word attack main effect in thirds
M2 = glmer(GatesScore ~ WordAttack.Thirds + zNarrativity + zWordConcreteness + zSyntacticSimplicity + zRefCohesion + zDeepCohesion + zArousal
+ (1|ID),data = Emotion_Data, family = binomial)
summary(M2)
## Generalized linear mixed model fit by maximum likelihood (Laplace
## Approximation) [glmerMod]
## Family: binomial ( logit )
## Formula: GatesScore ~ WordAttack.Thirds + zNarrativity + zWordConcreteness +
## zSyntacticSimplicity + zRefCohesion + zDeepCohesion + zArousal +
## (1 | ID)
## Data: Emotion_Data
##
## AIC BIC logLik deviance df.resid
## 7719.4 7789.1 -3849.7 7699.4 7862
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -7.2161 -0.5821 0.3246 0.5253 2.7620
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 1.048 1.024
## Number of obs: 7872, groups: ID, 164
##
## Fixed effects:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 1.77001 0.14744 12.005 < 2e-16 ***
## WordAttack.Thirdslow -1.76262 0.20939 -8.418 < 2e-16 ***
## WordAttack.Thirdsmedium -0.00792 0.21045 -0.038 0.9700
## zNarrativity 0.24901 0.04349 5.726 1.03e-08 ***
## zWordConcreteness -0.45241 0.04663 -9.702 < 2e-16 ***
## zSyntacticSimplicity -0.33429 0.05153 -6.487 8.73e-11 ***
## zRefCohesion -0.08438 0.05075 -1.663 0.0964 .
## zDeepCohesion 0.27026 0.03370 8.019 1.06e-15 ***
## zArousal 0.22938 0.04126 5.560 2.70e-08 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) WrdAttck.Thrdsl WrdAttck.Thrdsm zNrrtv zWrdCn zSyntS
## WrdAttck.Thrdsl -0.704
## WrdAttck.Thrdsm -0.692 0.491
## zNarrativty 0.025 -0.014 0.001
## zWrdCncrtns -0.039 0.022 0.001 0.334
## zSyntctcSmp -0.029 0.014 0.001 -0.065 0.529
## zRefCohesin -0.009 0.003 0.000 -0.476 -0.030 0.617
## zDeepCohesn 0.032 -0.018 -0.001 -0.036 -0.121 -0.245
## zArousal 0.022 -0.014 0.000 0.177 -0.198 0.248
## zRfChs zDpChs
## WrdAttck.Thrdsl
## WrdAttck.Thrdsm
## zNarrativty
## zWrdCncrtns
## zSyntctcSmp
## zRefCohesin
## zDeepCohesn -0.314
## zArousal 0.370 -0.411
#Interactions with Word Attack in Thirds
M3 = glmer(GatesScore ~ WordAttack.Thirds*zNarrativity + zWordConcreteness + zSyntacticSimplicity + zRefCohesion + zDeepCohesion + zArousal
+ (1|ID),data = Emotion_Data, family = binomial)
summary(M3)
## Generalized linear mixed model fit by maximum likelihood (Laplace
## Approximation) [glmerMod]
## Family: binomial ( logit )
## Formula: GatesScore ~ WordAttack.Thirds * zNarrativity + zWordConcreteness +
## zSyntacticSimplicity + zRefCohesion + zDeepCohesion + zArousal +
## (1 | ID)
## Data: Emotion_Data
##
## AIC BIC logLik deviance df.resid
## 7716.5 7800.1 -3846.2 7692.5 7860
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -7.9657 -0.5891 0.3257 0.5268 2.9152
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 1.049 1.024
## Number of obs: 7872, groups: ID, 164
##
## Fixed effects:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 1.774971 0.147717 12.016 < 2e-16 ***
## WordAttack.Thirdslow -1.769155 0.209623 -8.440 < 2e-16 ***
## WordAttack.Thirdsmedium 0.008761 0.211070 0.042 0.9669
## zNarrativity 0.270654 0.063689 4.250 2.14e-05 ***
## zWordConcreteness -0.461349 0.046968 -9.823 < 2e-16 ***
## zSyntacticSimplicity -0.343136 0.052053 -6.592 4.34e-11 ***
## zRefCohesion -0.090698 0.050966 -1.780 0.0751 .
## zDeepCohesion 0.272344 0.033691 8.084 6.29e-16 ***
## zArousal 0.234080 0.041320 5.665 1.47e-08 ***
## WordAttack.Thirdslow:zNarrativity -0.101636 0.071208 -1.427 0.1535
## WordAttack.Thirdsmedium:zNarrativity 0.087744 0.078477 1.118 0.2635
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) WrdAttck.Thrdsl WrdAttck.Thrdsm zNrrtv zWrdCn zSyntS
## WrdAttck.Thrdsl -0.704
## WrdAttck.Thrdsm -0.692 0.491
## zNarrativty 0.058 -0.041 -0.036
## zWrdCncrtns -0.042 0.025 -0.001 0.197
## zSyntctcSmp -0.032 0.017 0.000 -0.074 0.533
## zRefCohesin -0.012 0.006 -0.001 -0.345 -0.021 0.623
## zDeepCohesn 0.033 -0.019 -0.001 -0.012 -0.128 -0.256
## zArousal 0.024 -0.016 0.001 0.139 -0.200 0.243
## WrdAttck.Thrdsl:N -0.052 0.040 0.033 -0.675 0.062 0.059
## WrdAttck.Thrdsm:N -0.042 0.030 0.070 -0.595 -0.015 -0.005
## zRfChs zDpChs zArosl WrdAttck.Thrdsl:N
## WrdAttck.Thrdsl
## WrdAttck.Thrdsm
## zNarrativty
## zWrdCncrtns
## zSyntctcSmp
## zRefCohesin
## zDeepCohesn -0.320
## zArousal 0.364 -0.406
## WrdAttck.Thrdsl:N 0.043 -0.023 -0.036
## WrdAttck.Thrdsm:N -0.003 0.002 0.011 0.532
M4 = glmer(GatesScore ~ zNarrativity + WordAttack.Thirds*zWordConcreteness + zSyntacticSimplicity + zRefCohesion + zDeepCohesion + zArousal
+ (1|ID),data = Emotion_Data, family = binomial)
## Warning in checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
## Model failed to converge with max|grad| = 0.00394573 (tol = 0.001, component 1)
summary(M4)
## Generalized linear mixed model fit by maximum likelihood (Laplace
## Approximation) [glmerMod]
## Family: binomial ( logit )
## Formula: GatesScore ~ zNarrativity + WordAttack.Thirds * zWordConcreteness +
## zSyntacticSimplicity + zRefCohesion + zDeepCohesion + zArousal +
## (1 | ID)
## Data: Emotion_Data
##
## AIC BIC logLik deviance df.resid
## 7722.8 7806.5 -3849.4 7698.8 7860
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -7.2329 -0.5824 0.3259 0.5256 2.8247
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 1.049 1.024
## Number of obs: 7872, groups: ID, 164
##
## Fixed effects:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 1.77331 0.14759 12.015 < 2e-16
## zNarrativity 0.25018 0.04352 5.748 9.02e-09
## WordAttack.Thirdslow -1.76616 0.20951 -8.430 < 2e-16
## WordAttack.Thirdsmedium -0.01082 0.21065 -0.051 0.9590
## zWordConcreteness -0.47994 0.06363 -7.542 4.62e-14
## zSyntacticSimplicity -0.33558 0.05165 -6.497 8.18e-11
## zRefCohesion -0.08533 0.05077 -1.681 0.0928
## zDeepCohesion 0.27067 0.03373 8.025 1.01e-15
## zArousal 0.22928 0.04126 5.558 2.74e-08
## WordAttack.Thirdslow:zWordConcreteness 0.04911 0.06837 0.718 0.4725
## WordAttack.Thirdsmedium:zWordConcreteness 0.02458 0.07353 0.334 0.7382
##
## (Intercept) ***
## zNarrativity ***
## WordAttack.Thirdslow ***
## WordAttack.Thirdsmedium
## zWordConcreteness ***
## zSyntacticSimplicity ***
## zRefCohesion .
## zDeepCohesion ***
## zArousal ***
## WordAttack.Thirdslow:zWordConcreteness
## WordAttack.Thirdsmedium:zWordConcreteness
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) zNrrtv WrdAttck.Thrdsl WrdAttck.Thrdsm zWrdCn zSyntS
## zNarrativty 0.026
## WrdAttck.Thrdsl -0.704 -0.015
## WrdAttck.Thrdsm -0.693 0.001 0.491
## zWrdCncrtns -0.055 0.227 0.035 0.024
## zSyntctcSmp -0.029 -0.067 0.014 0.000 0.398
## zRefCohesin -0.010 -0.476 0.004 -0.001 -0.013 0.619
## zDeepCohesn 0.032 -0.035 -0.018 -0.001 -0.096 -0.247
## zArousal 0.022 0.178 -0.014 0.000 -0.144 0.248
## WrdAttck.Thrdsl:WC 0.034 0.037 -0.026 -0.022 -0.620 -0.030
## WrdAttck.Thrdsm:WC 0.029 -0.001 -0.021 -0.039 -0.569 0.005
## zRfChs zDpChs zArosl WrdAttck.Thrdsl:WC
## zNarrativty
## WrdAttck.Thrdsl
## WrdAttck.Thrdsm
## zWrdCncrtns
## zSyntctcSmp
## zRefCohesin
## zDeepCohesn -0.316
## zArousal 0.369 -0.411
## WrdAttck.Thrdsl:WC -0.023 0.016 -0.003
## WrdAttck.Thrdsm:WC 0.004 -0.002 0.000 0.531
## convergence code: 0
## Model failed to converge with max|grad| = 0.00394573 (tol = 0.001, component 1)
M5 = glmer(GatesScore ~ zNarrativity + zWordConcreteness + WordAttack.Thirds*zSyntacticSimplicity + zRefCohesion + zDeepCohesion + zArousal
+ (1|ID),data = Emotion_Data, family = binomial)
## Warning in checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
## Model failed to converge with max|grad| = 0.00220419 (tol = 0.001, component 1)
summary(M5)
## Generalized linear mixed model fit by maximum likelihood (Laplace
## Approximation) [glmerMod]
## Family: binomial ( logit )
## Formula: GatesScore ~ zNarrativity + zWordConcreteness + WordAttack.Thirds *
## zSyntacticSimplicity + zRefCohesion + zDeepCohesion + zArousal +
## (1 | ID)
## Data: Emotion_Data
##
## AIC BIC logLik deviance df.resid
## 7720.4 7804.1 -3848.2 7696.4 7860
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -7.3306 -0.5789 0.3245 0.5282 2.8027
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 1.049 1.024
## Number of obs: 7872, groups: ID, 164
##
## Fixed effects:
## Estimate Std. Error z value
## (Intercept) 1.767371 0.147493 11.983
## zNarrativity 0.248697 0.043488 5.719
## zWordConcreteness -0.452181 0.046639 -9.695
## WordAttack.Thirdslow -1.759602 0.209477 -8.400
## WordAttack.Thirdsmedium -0.002055 0.210567 -0.010
## zSyntacticSimplicity -0.264578 0.067303 -3.931
## zRefCohesion -0.083475 0.050764 -1.644
## zDeepCohesion 0.269533 0.033740 7.989
## zArousal 0.229207 0.041243 5.557
## WordAttack.Thirdslow:zSyntacticSimplicity -0.083652 0.067761 -1.235
## WordAttack.Thirdsmedium:zSyntacticSimplicity -0.120141 0.071878 -1.671
## Pr(>|z|)
## (Intercept) < 2e-16 ***
## zNarrativity 1.07e-08 ***
## zWordConcreteness < 2e-16 ***
## WordAttack.Thirdslow < 2e-16 ***
## WordAttack.Thirdsmedium 0.9922
## zSyntacticSimplicity 8.45e-05 ***
## zRefCohesion 0.1001
## zDeepCohesion 1.36e-15 ***
## zArousal 2.74e-08 ***
## WordAttack.Thirdslow:zSyntacticSimplicity 0.2170
## WordAttack.Thirdsmedium:zSyntacticSimplicity 0.0946 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) zNrrtv zWrdCn WrdAttck.Thrdsl WrdAttck.Thrdsm zSyntS
## zNarrativty 0.025
## zWrdCncrtns -0.039 0.334
## WrdAttck.Thrdsl -0.704 -0.014 0.022
## WrdAttck.Thrdsm -0.692 0.001 0.000 0.491
## zSyntctcSmp -0.023 -0.056 0.412 0.012 0.003
## zRefCohesin -0.009 -0.476 -0.030 0.004 0.000 0.486
## zDeepCohesn 0.032 -0.036 -0.122 -0.018 -0.001 -0.207
## zArousal 0.022 0.177 -0.199 -0.014 0.000 0.185
## WrdAttck.Thrdsl:SS 0.007 0.012 -0.024 -0.007 -0.003 -0.586
## WrdAttck.Thrdsm:SS 0.003 0.002 0.006 -0.003 -0.017 -0.536
## zRfChs zDpChs zArosl WrdAttck.Thrdsl:SS
## zNarrativty
## zWrdCncrtns
## WrdAttck.Thrdsl
## WrdAttck.Thrdsm
## zSyntctcSmp
## zRefCohesin
## zDeepCohesn -0.316
## zArousal 0.369 -0.410
## WrdAttck.Thrdsl:SS -0.026 0.045 0.010
## WrdAttck.Thrdsm:SS -0.010 0.003 0.001 0.531
## convergence code: 0
## Model failed to converge with max|grad| = 0.00220419 (tol = 0.001, component 1)
M6 = glmer(GatesScore ~ zNarrativity + zWordConcreteness + zSyntacticSimplicity + WordAttack.Thirds*zRefCohesion + zDeepCohesion + zArousal
+ (1|ID),data = Emotion_Data, family = binomial)
## Warning in checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
## Model failed to converge with max|grad| = 0.00114104 (tol = 0.001, component 1)
summary(M6)
## Generalized linear mixed model fit by maximum likelihood (Laplace
## Approximation) [glmerMod]
## Family: binomial ( logit )
## Formula:
## GatesScore ~ zNarrativity + zWordConcreteness + zSyntacticSimplicity +
## WordAttack.Thirds * zRefCohesion + zDeepCohesion + zArousal + (1 | ID)
## Data: Emotion_Data
##
## AIC BIC logLik deviance df.resid
## 7716.3 7800.0 -3846.2 7692.3 7860
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -7.4426 -0.5811 0.3204 0.5328 2.7635
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 1.05 1.025
## Number of obs: 7872, groups: ID, 164
##
## Fixed effects:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 1.764134 0.147504 11.960 < 2e-16 ***
## zNarrativity 0.247673 0.043529 5.690 1.27e-08 ***
## zWordConcreteness -0.452404 0.046686 -9.690 < 2e-16 ***
## zSyntacticSimplicity -0.333413 0.051555 -6.467 9.99e-11 ***
## WordAttack.Thirdslow -1.756029 0.209527 -8.381 < 2e-16 ***
## WordAttack.Thirdsmedium 0.008921 0.210750 0.042 0.96624
## zRefCohesion -0.180367 0.065420 -2.757 0.00583 **
## zDeepCohesion 0.271155 0.033761 8.032 9.62e-16 ***
## zArousal 0.229430 0.041274 5.559 2.72e-08 ***
## WordAttack.Thirdslow:zRefCohesion 0.100999 0.067382 1.499 0.13390
## WordAttack.Thirdsmedium:zRefCohesion 0.194012 0.072661 2.670 0.00758 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) zNrrtv zWrdCn zSyntS WrdAttck.Thrdsl WrdAttck.Thrdsm
## zNarrativty 0.025
## zWrdCncrtns -0.039 0.335
## zSyntctcSmp -0.029 -0.066 0.529
## WrdAttck.Thrdsl -0.703 -0.014 0.022 0.014
## WrdAttck.Thrdsm -0.692 -0.001 0.000 0.001 0.490
## zRefCohesin -0.005 -0.363 -0.035 0.471 0.000 -0.006
## zDeepCohesn 0.032 -0.037 -0.123 -0.244 -0.018 0.000
## zArousal 0.022 0.177 -0.198 0.248 -0.014 0.000
## WrdAttck.Thrdsl:RC -0.010 0.001 0.028 0.010 0.009 0.006
## WrdAttck.Thrdsm:RC -0.008 -0.024 -0.002 0.012 0.006 0.031
## zRfChs zDpChs zArosl WrdAttck.Thrdsl:RC
## zNarrativty
## zWrdCncrtns
## zSyntctcSmp
## WrdAttck.Thrdsl
## WrdAttck.Thrdsm
## zRefCohesin
## zDeepCohesn -0.239
## zArousal 0.290 -0.410
## WrdAttck.Thrdsl:RC -0.573 -0.026 -0.009
## WrdAttck.Thrdsm:RC -0.520 0.022 0.001 0.516
## convergence code: 0
## Model failed to converge with max|grad| = 0.00114104 (tol = 0.001, component 1)
plot(effect("WordAttack.Thirds:zRefCohesion", M6), multiline=TRUE)
M7 = glmer(GatesScore ~ zNarrativity + zWordConcreteness + zSyntacticSimplicity + zRefCohesion + WordAttack.Thirds*zDeepCohesion + zArousal
+ (1|ID),data = Emotion_Data, family = binomial)
## Warning in checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
## Model failed to converge with max|grad| = 0.00119292 (tol = 0.001, component 1)
summary(M7)
## Generalized linear mixed model fit by maximum likelihood (Laplace
## Approximation) [glmerMod]
## Family: binomial ( logit )
## Formula:
## GatesScore ~ zNarrativity + zWordConcreteness + zSyntacticSimplicity +
## zRefCohesion + WordAttack.Thirds * zDeepCohesion + zArousal + (1 | ID)
## Data: Emotion_Data
##
## AIC BIC logLik deviance df.resid
## 7721.6 7805.3 -3848.8 7697.6 7860
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -7.4004 -0.5820 0.3247 0.5266 2.6932
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 1.048 1.024
## Number of obs: 7872, groups: ID, 164
##
## Fixed effects:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 1.770779 0.147605 11.997 < 2e-16 ***
## zNarrativity 0.247878 0.043593 5.686 1.30e-08 ***
## zWordConcreteness -0.456411 0.046826 -9.747 < 2e-16 ***
## zSyntacticSimplicity -0.338289 0.051641 -6.551 5.72e-11 ***
## zRefCohesion -0.084099 0.050910 -1.652 0.0986 .
## WordAttack.Thirdslow -1.763537 0.209500 -8.418 < 2e-16 ***
## WordAttack.Thirdsmedium 0.001315 0.210863 0.006 0.9950
## zDeepCohesion 0.272455 0.056102 4.856 1.20e-06 ***
## zArousal 0.229131 0.041232 5.557 2.74e-08 ***
## WordAttack.Thirdslow:zDeepCohesion -0.040668 0.069868 -0.582 0.5605
## WordAttack.Thirdsmedium:zDeepCohesion 0.053798 0.076317 0.705 0.4809
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) zNrrtv zWrdCn zSyntS zRfChs WrdAttck.Thrdsl
## zNarrativty 0.024
## zWrdCncrtns -0.041 0.334
## zSyntctcSmp -0.031 -0.068 0.533
## zRefCohesin -0.009 -0.479 -0.029 0.616
## WrdAttck.Thrdsl -0.704 -0.013 0.024 0.016 0.003
## WrdAttck.Thrdsm -0.692 0.000 0.000 0.000 0.000 0.491
## zDeepCohesn 0.060 -0.029 -0.101 -0.174 -0.188 -0.042
## zArousal 0.022 0.178 -0.196 0.248 0.369 -0.014
## WrdAttck.Thrdsl:DC -0.048 0.014 0.056 0.052 -0.001 0.036
## WrdAttck.Thrdsm:DC -0.041 -0.009 -0.012 -0.006 0.007 0.031
## WrdAttck.Thrdsm zDpChs zArosl WrdAttck.Thrdsl:DC
## zNarrativty
## zWrdCncrtns
## zSyntctcSmp
## zRefCohesin
## WrdAttck.Thrdsl
## WrdAttck.Thrdsm
## zDeepCohesn -0.039
## zArousal 0.000 -0.249
## WrdAttck.Thrdsl:DC 0.031 -0.731 0.004
## WrdAttck.Thrdsm:DC 0.063 -0.662 -0.001 0.531
## convergence code: 0
## Model failed to converge with max|grad| = 0.00119292 (tol = 0.001, component 1)
M8 = glmer(GatesScore ~ zNarrativity + zWordConcreteness + zSyntacticSimplicity + zRefCohesion + zDeepCohesion + WordAttack.Thirds*zArousal
+ (1|ID),data = Emotion_Data, family = binomial)
## Warning in checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
## Model failed to converge with max|grad| = 0.00252078 (tol = 0.001, component 1)
summary(M8)
## Generalized linear mixed model fit by maximum likelihood (Laplace
## Approximation) [glmerMod]
## Family: binomial ( logit )
## Formula:
## GatesScore ~ zNarrativity + zWordConcreteness + zSyntacticSimplicity +
## zRefCohesion + zDeepCohesion + WordAttack.Thirds * zArousal + (1 | ID)
## Data: Emotion_Data
##
## AIC BIC logLik deviance df.resid
## 7713.2 7796.9 -3844.6 7689.2 7860
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -7.0459 -0.5754 0.3267 0.5278 2.6498
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 1.053 1.026
## Number of obs: 7872, groups: ID, 164
##
## Fixed effects:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 1.78303 0.14788 12.057 < 2e-16 ***
## zNarrativity 0.25092 0.04358 5.758 8.51e-09 ***
## zWordConcreteness -0.45388 0.04666 -9.727 < 2e-16 ***
## zSyntacticSimplicity -0.33464 0.05136 -6.516 7.23e-11 ***
## zRefCohesion -0.08332 0.05073 -1.642 0.10051
## zDeepCohesion 0.26869 0.03375 7.962 1.69e-15 ***
## WordAttack.Thirdslow -1.77591 0.20989 -8.461 < 2e-16 ***
## WordAttack.Thirdsmedium -0.02232 0.21093 -0.106 0.91573
## zArousal 0.37138 0.06077 6.112 9.86e-10 ***
## WordAttack.Thirdslow:zArousal -0.19571 0.06864 -2.851 0.00436 **
## WordAttack.Thirdsmedium:zArousal -0.20510 0.07441 -2.756 0.00585 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) zNrrtv zWrdCn zSyntS zRfChs zDpChs WrdAttck.Thrdsl
## zNarrativty 0.026
## zWrdCncrtns -0.040 0.334
## zSyntctcSmp -0.029 -0.068 0.530
## zRefCohesin -0.009 -0.478 -0.031 0.616
## zDeepCohesn 0.030 -0.035 -0.119 -0.242 -0.313
## WrdAttck.Thrdsl -0.704 -0.015 0.022 0.014 0.003 -0.017
## WrdAttck.Thrdsm -0.693 0.000 0.001 0.001 -0.001 -0.001 0.492
## zArousal 0.051 0.132 -0.152 0.163 0.258 -0.296 -0.035
## WrdAttck.Thrdsl:A -0.036 -0.026 0.015 0.005 -0.009 0.022 0.026
## WrdAttck.Thrdsm:A -0.033 0.000 0.021 0.001 -0.007 0.011 0.023
## WrdAttck.Thrdsm zArosl WrdAttck.Thrdsl:A
## zNarrativty
## zWrdCncrtns
## zSyntctcSmp
## zRefCohesin
## zDeepCohesn
## WrdAttck.Thrdsl
## WrdAttck.Thrdsm
## zArousal -0.027
## WrdAttck.Thrdsl:A 0.023 -0.668
## WrdAttck.Thrdsm:A 0.024 -0.608 0.529
## convergence code: 0
## Model failed to converge with max|grad| = 0.00252078 (tol = 0.001, component 1)
plot(effect("WordAttack.Thirds:zArousal", M8), multiline=TRUE)
##Interactions with Word Attack and Arousal in Quarters
M9 = glmer(GatesScore ~ zNarrativity + zWordConcreteness + zSyntacticSimplicity + zRefCohesion + zDeepCohesion + WordAttack.Quarters*zArousal
+ (1|ID),data = Emotion_Data, family = binomial)
## Warning in checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
## Model failed to converge with max|grad| = 0.00694542 (tol = 0.001, component 1)
summary(M9)
## Generalized linear mixed model fit by maximum likelihood (Laplace
## Approximation) [glmerMod]
## Family: binomial ( logit )
## Formula:
## GatesScore ~ zNarrativity + zWordConcreteness + zSyntacticSimplicity +
## zRefCohesion + zDeepCohesion + WordAttack.Quarters * zArousal +
## (1 | ID)
## Data: Emotion_Data
##
## AIC BIC logLik deviance df.resid
## 7710.9 7808.5 -3841.5 7682.9 7858
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -6.8066 -0.5675 0.3227 0.5325 2.7321
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 1.036 1.018
## Number of obs: 7872, groups: ID, 164
##
## Fixed effects:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 1.98666 0.17215 11.540 < 2e-16
## zNarrativity 0.25305 0.04365 5.797 6.76e-09
## zWordConcreteness -0.45455 0.04668 -9.738 < 2e-16
## zSyntacticSimplicity -0.33561 0.05119 -6.556 5.54e-11
## zRefCohesion -0.08316 0.05071 -1.640 0.10104
## zDeepCohesion 0.26822 0.03378 7.941 2.00e-15
## WordAttack.Quarterslow -2.20568 0.24394 -9.042 < 2e-16
## WordAttack.Quartersmedium_high -0.27321 0.24503 -1.115 0.26484
## WordAttack.Quartersmedium_low -0.74726 0.23781 -3.142 0.00168
## zArousal 0.31378 0.07026 4.466 7.98e-06
## WordAttack.Quarterslow:zArousal -0.18585 0.08195 -2.268 0.02334
## WordAttack.Quartersmedium_high:zArousal 0.08625 0.08783 0.982 0.32607
## WordAttack.Quartersmedium_low:zArousal -0.15529 0.08307 -1.869 0.06158
##
## (Intercept) ***
## zNarrativity ***
## zWordConcreteness ***
## zSyntacticSimplicity ***
## zRefCohesion
## zDeepCohesion ***
## WordAttack.Quarterslow ***
## WordAttack.Quartersmedium_high
## WordAttack.Quartersmedium_low **
## zArousal ***
## WordAttack.Quarterslow:zArousal *
## WordAttack.Quartersmedium_high:zArousal
## WordAttack.Quartersmedium_low:zArousal .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation matrix not shown by default, as p = 13 > 12.
## Use print(x, correlation=TRUE) or
## vcov(x) if you need it
## convergence code: 0
## Model failed to converge with max|grad| = 0.00694542 (tol = 0.001, component 1)
plot(effect("WordAttack.Quarters:zArousal", M9), multiline=TRUE)