Emotion Paper Models with Word Attack, Feb 26, 2020
Sage Pickren
Attaching Data
Checking Packages
library("lavaan")
library("effects")
library("psych")
library("lme4")
library("lmerTest")
library("sjstats")
library("sjmisc")
library("tidyverse")
library("tinytex")Centering Variables
Corrected_EmotionPaper_Data_Quartiles$ID<- as.factor(Corrected_EmotionPaper_Data_Quartiles$ID)
Corrected_EmotionPaper_Data_Quartiles$zArousal<-scale(Corrected_EmotionPaper_Data_Quartiles$Arousal_Mean_Sum, scale=TRUE, center=TRUE)
Corrected_EmotionPaper_Data_Quartiles$zValence<-scale(Corrected_EmotionPaper_Data_Quartiles$Valence_Mean_Sum, scale=TRUE, center=TRUE)
Corrected_EmotionPaper_Data_Quartiles$zWordConcreteness<-scale(Corrected_EmotionPaper_Data_Quartiles$Word_concreteness_z_score, scale=TRUE, center=TRUE)
Corrected_EmotionPaper_Data_Quartiles$zNarrativity<-scale(Corrected_EmotionPaper_Data_Quartiles$Narrativity, scale=TRUE, center=TRUE)
Corrected_EmotionPaper_Data_Quartiles$zDeepCohesion<-scale(Corrected_EmotionPaper_Data_Quartiles$Deep_cohesion_z_score, scale=TRUE, center=TRUE)
Corrected_EmotionPaper_Data_Quartiles$zRefCohesion<-scale(Corrected_EmotionPaper_Data_Quartiles$Referential_cohesion_z_score, scale=TRUE, center=TRUE)
Corrected_EmotionPaper_Data_Quartiles$zSyntacticSimplicity<-scale(Corrected_EmotionPaper_Data_Quartiles$Syntactic_simplicity_z_score, scale=TRUE, center=TRUE)
Corrected_EmotionPaper_Data_Quartiles$zWordAttack<-scale(Corrected_EmotionPaper_Data_Quartiles$WJ.III.Word.Attack.Standard.Score, scale=TRUE, center=TRUE)
Corrected_EmotionPaper_Data_Quartiles$zLWID<-scale(Corrected_EmotionPaper_Data_Quartiles$WJ.III.Letter.Word.Identification.Standard.Score, scale=TRUE, center=TRUE)
BASE Models with LWID
M3 = glmer(GatesScore ~ zLWID + zNarrativity + zWordConcreteness + zRefCohesion + zDeepCohesion
+ zSyntacticSimplicity
+ (1|ID),data = Corrected_EmotionPaper_Data_Quartiles, family = binomial)
summary(M3)## Generalized linear mixed model fit by maximum likelihood (Laplace
## Approximation) [glmerMod]
## Family: binomial ( logit )
## Formula:
## GatesScore ~ zLWID + zNarrativity + zWordConcreteness + zRefCohesion +
## zDeepCohesion + zSyntacticSimplicity + (1 | ID)
## Data: Corrected_EmotionPaper_Data_Quartiles
##
## AIC BIC logLik deviance df.resid
## 7752.7 7808.5 -3868.3 7736.7 7912
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -6.3715 -0.5808 0.3277 0.5367 2.8730
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 0.8209 0.906
## Number of obs: 7920, groups: ID, 165
##
## Fixed effects:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 1.20728 0.07812 15.454 < 2e-16 ***
## zLWID 0.92881 0.07736 12.006 < 2e-16 ***
## zNarrativity 0.20829 0.04232 4.921 8.59e-07 ***
## zWordConcreteness -0.40273 0.04621 -8.716 < 2e-16 ***
## zRefCohesion -0.19281 0.04697 -4.105 4.04e-05 ***
## zDeepCohesion 0.34958 0.03082 11.344 < 2e-16 ***
## zSyntacticSimplicity -0.40721 0.05013 -8.122 4.57e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) zLWID zNrrtv zWrdCn zRfChs zDpChs
## zLWID 0.063
## zNarrativty 0.028 0.017
## zWrdCncrtns -0.051 -0.027 0.391
## zRefCohesin -0.027 -0.014 -0.589 0.050
## zDeepCohesn 0.061 0.036 0.028 -0.237 -0.193
## zSyntctcSmp -0.051 -0.025 -0.102 0.617 0.580 -0.166M4 = glmer(GatesScore ~ zLWID + zNarrativity + zWordConcreteness + zRefCohesion + zDeepCohesion
+ (1|ID),data = Corrected_EmotionPaper_Data_Quartiles, family = binomial)
summary(M4)## Generalized linear mixed model fit by maximum likelihood (Laplace
## Approximation) [glmerMod]
## Family: binomial ( logit )
## Formula:
## GatesScore ~ zLWID + zNarrativity + zWordConcreteness + zRefCohesion +
## zDeepCohesion + (1 | ID)
## Data: Corrected_EmotionPaper_Data_Quartiles
##
## AIC BIC logLik deviance df.resid
## 7817.7 7866.6 -3901.9 7803.7 7913
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -5.2700 -0.6103 0.3241 0.5413 3.0411
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 0.8068 0.8982
## Number of obs: 7920, groups: ID, 165
##
## Fixed effects:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 1.18914 0.07741 15.361 < 2e-16 ***
## zLWID 0.92170 0.07670 12.017 < 2e-16 ***
## zNarrativity 0.17934 0.04196 4.274 1.92e-05 ***
## zWordConcreteness -0.17888 0.03702 -4.831 1.36e-06 ***
## zRefCohesion 0.03009 0.03869 0.778 0.437
## zDeepCohesion 0.31316 0.02993 10.463 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) zLWID zNrrtv zWrdCn zRfChs
## zLWID 0.064
## zNarrativty 0.018 0.016
## zWrdCncrtns -0.025 -0.015 0.587
## zRefCohesin 0.007 0.001 -0.653 -0.481
## zDeepCohesn 0.049 0.034 -0.033 -0.160 -0.093Arousal and Valence while controlling for LWID
E1 = glmer(GatesScore ~ zLWID + zNarrativity + zWordConcreteness + zRefCohesion + zDeepCohesion
+ zSyntacticSimplicity +zArousal
+ (1|ID),data = Corrected_EmotionPaper_Data_Quartiles, family = binomial)
summary(E1)## Generalized linear mixed model fit by maximum likelihood (Laplace
## Approximation) [glmerMod]
## Family: binomial ( logit )
## Formula:
## GatesScore ~ zLWID + zNarrativity + zWordConcreteness + zRefCohesion +
## zDeepCohesion + zSyntacticSimplicity + zArousal + (1 | ID)
## Data: Corrected_EmotionPaper_Data_Quartiles
##
## AIC BIC logLik deviance df.resid
## 7724.2 7787.0 -3853.1 7706.2 7911
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -6.1816 -0.5751 0.3278 0.5283 2.7302
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 0.8286 0.9103
## Number of obs: 7920, groups: ID, 165
##
## Fixed effects:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 1.21372 0.07848 15.465 < 2e-16 ***
## zLWID 0.93300 0.07771 12.006 < 2e-16 ***
## zNarrativity 0.24795 0.04340 5.713 1.11e-08 ***
## zWordConcreteness -0.44783 0.04653 -9.625 < 2e-16 ***
## zRefCohesion -0.08364 0.05066 -1.651 0.0987 .
## zDeepCohesion 0.27021 0.03363 8.035 9.36e-16 ***
## zSyntacticSimplicity -0.32751 0.05141 -6.371 1.88e-10 ***
## zArousal 0.22823 0.04116 5.545 2.95e-08 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) zLWID zNrrtv zWrdCn zRfChs zDpChs zSyntS
## zLWID 0.063
## zNarrativty 0.036 0.020
## zWrdCncrtns -0.054 -0.031 0.335
## zRefCohesin -0.014 -0.005 -0.476 -0.031
## zDeepCohesn 0.043 0.025 -0.036 -0.121 -0.314
## zSyntctcSmp -0.041 -0.019 -0.065 0.528 0.617 -0.244
## zArousal 0.030 0.019 0.177 -0.198 0.370 -0.411 0.248E2 = glmer(GatesScore ~ zLWID + zNarrativity + zWordConcreteness + zRefCohesion + zDeepCohesion
+zArousal
+ (1|ID),data = Corrected_EmotionPaper_Data_Quartiles, family = binomial)
summary(E2)## Generalized linear mixed model fit by maximum likelihood (Laplace
## Approximation) [glmerMod]
## Family: binomial ( logit )
## Formula:
## GatesScore ~ zLWID + zNarrativity + zWordConcreteness + zRefCohesion +
## zDeepCohesion + zArousal + (1 | ID)
## Data: Corrected_EmotionPaper_Data_Quartiles
##
## AIC BIC logLik deviance df.resid
## 7763.3 7819.1 -3873.7 7747.3 7912
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -5.6037 -0.5923 0.3240 0.5359 2.7239
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 0.8203 0.9057
## Number of obs: 7920, groups: ID, 165
##
## Fixed effects:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 1.20171 0.07805 15.396 < 2e-16 ***
## zLWID 0.92892 0.07731 12.015 < 2e-16 ***
## zNarrativity 0.23248 0.04309 5.395 6.84e-08 ***
## zWordConcreteness -0.29557 0.03985 -7.416 1.20e-13 ***
## zRefCohesion 0.11636 0.04016 2.897 0.00376 **
## zDeepCohesion 0.22217 0.03214 6.912 4.77e-12 ***
## zArousal 0.29096 0.03896 7.469 8.09e-14 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) zLWID zNrrtv zWrdCn zRfChs zDpChs
## zLWID 0.064
## zNarrativty 0.028 0.019
## zWrdCncrtns -0.038 -0.024 0.451
## zRefCohesin 0.016 0.008 -0.564 -0.534
## zDeepCohesn 0.031 0.022 -0.074 0.014 -0.194
## zArousal 0.040 0.026 0.181 -0.394 0.282 -0.367E3 = glmer(GatesScore ~ zLWID + zNarrativity + zWordConcreteness + zRefCohesion + zDeepCohesion
+ zSyntacticSimplicity +zValence
+ (1|ID),data = Corrected_EmotionPaper_Data_Quartiles, family = binomial)## Warning in checkConv(attr(opt, "derivs"), opt$par, ctrl =
## control$checkConv, : Model failed to converge with max|grad| = 0.00172486
## (tol = 0.001, component 1)summary(E3)## Generalized linear mixed model fit by maximum likelihood (Laplace
## Approximation) [glmerMod]
## Family: binomial ( logit )
## Formula:
## GatesScore ~ zLWID + zNarrativity + zWordConcreteness + zRefCohesion +
## zDeepCohesion + zSyntacticSimplicity + zValence + (1 | ID)
## Data: Corrected_EmotionPaper_Data_Quartiles
##
## AIC BIC logLik deviance df.resid
## 7745.5 7808.3 -3863.8 7727.5 7911
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -6.4455 -0.5829 0.3238 0.5367 3.1464
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 0.8241 0.9078
## Number of obs: 7920, groups: ID, 165
##
## Fixed effects:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 1.20782 0.07826 15.434 < 2e-16 ***
## zLWID 0.93054 0.07750 12.006 < 2e-16 ***
## zNarrativity 0.16455 0.04456 3.693 0.000222 ***
## zWordConcreteness -0.46398 0.05048 -9.191 < 2e-16 ***
## zRefCohesion -0.19564 0.04686 -4.175 2.98e-05 ***
## zDeepCohesion 0.41973 0.03849 10.905 < 2e-16 ***
## zSyntacticSimplicity -0.54553 0.06787 -8.038 9.13e-16 ***
## zValence -0.16037 0.05264 -3.047 0.002314 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) zLWID zNrrtv zWrdCn zRfChs zDpChs zSyntS
## zLWID 0.064
## zNarrativty 0.022 0.012
## zWrdCncrtns -0.050 -0.030 0.474
## zRefCohesin -0.026 -0.014 -0.549 0.050
## zDeepCohesn 0.053 0.037 -0.174 -0.410 -0.168
## zSyntctcSmp -0.043 -0.027 0.151 0.688 0.439 -0.496
## zValence -0.008 -0.013 0.323 0.403 0.022 -0.596 0.675
## convergence code: 0
## Model failed to converge with max|grad| = 0.00172486 (tol = 0.001, component 1)Interactions with LWID and narrativity
M5 = glmer(GatesScore ~ zNarrativity*zLWID + zWordConcreteness + zRefCohesion + zDeepCohesion
+ zSyntacticSimplicity
+ (1|ID),data = Corrected_EmotionPaper_Data_Quartiles, family = binomial)
summary(M5)## Generalized linear mixed model fit by maximum likelihood (Laplace
## Approximation) [glmerMod]
## Family: binomial ( logit )
## Formula:
## GatesScore ~ zNarrativity * zLWID + zWordConcreteness + zRefCohesion +
## zDeepCohesion + zSyntacticSimplicity + (1 | ID)
## Data: Corrected_EmotionPaper_Data_Quartiles
##
## AIC BIC logLik deviance df.resid
## 7754.6 7817.4 -3868.3 7736.6 7911
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -6.3991 -0.5831 0.3278 0.5371 2.8476
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 0.8209 0.906
## Number of obs: 7920, groups: ID, 165
##
## Fixed effects:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 1.208129 0.078201 15.449 < 2e-16 ***
## zNarrativity 0.210309 0.043075 4.882 1.05e-06 ***
## zLWID 0.929374 0.077402 12.007 < 2e-16 ***
## zWordConcreteness -0.403525 0.046339 -8.708 < 2e-16 ***
## zRefCohesion -0.193673 0.047111 -4.111 3.94e-05 ***
## zDeepCohesion 0.349853 0.030839 11.344 < 2e-16 ***
## zSyntacticSimplicity -0.408349 0.050387 -8.104 5.31e-16 ***
## zNarrativity:zLWID 0.007283 0.029859 0.244 0.807
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) zNrrtv zLWID zWrdCn zRfChs zDpChs zSyntS
## zNarrativty 0.036
## zLWID 0.065 0.023
## zWrdCncrtns -0.054 0.370 -0.029
## zRefCohesin -0.030 -0.591 -0.016 0.056
## zDeepCohesn 0.063 0.035 0.038 -0.239 -0.196
## zSyntctcSmp -0.055 -0.118 -0.028 0.619 0.583 -0.170
## zNrrtv:LWID 0.044 0.187 0.031 -0.071 -0.074 0.037 -0.093M6 = glmer(GatesScore ~ zNarrativity*zLWID + zWordConcreteness + zRefCohesion + zDeepCohesion
+ (1|ID),data = Corrected_EmotionPaper_Data_Quartiles, family = binomial)
summary(M6)## Generalized linear mixed model fit by maximum likelihood (Laplace
## Approximation) [glmerMod]
## Family: binomial ( logit )
## Formula:
## GatesScore ~ zNarrativity * zLWID + zWordConcreteness + zRefCohesion +
## zDeepCohesion + (1 | ID)
## Data: Corrected_EmotionPaper_Data_Quartiles
##
## AIC BIC logLik deviance df.resid
## 7819.6 7875.5 -3901.8 7803.6 7912
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -5.2893 -0.6048 0.3225 0.5405 3.0779
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 0.8067 0.8982
## Number of obs: 7920, groups: ID, 165
##
## Fixed effects:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 1.188116 0.077459 15.339 < 2e-16 ***
## zNarrativity 0.176867 0.042586 4.153 3.28e-05 ***
## zLWID 0.921109 0.076710 12.008 < 2e-16 ***
## zWordConcreteness -0.178537 0.037021 -4.823 1.42e-06 ***
## zRefCohesion 0.030217 0.038696 0.781 0.435
## zDeepCohesion 0.313340 0.029944 10.464 < 2e-16 ***
## zNarrativity:zLWID -0.009773 0.029535 -0.331 0.741
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) zNrrtv zLWID zWrdCn zRfChs zDpChs
## zNarrativty 0.023
## zLWID 0.064 0.019
## zWrdCncrtns -0.026 0.574 -0.015
## zRefCohesin 0.007 -0.645 0.000 -0.481
## zDeepCohesn 0.048 -0.036 0.034 -0.160 -0.093
## zNrrtv:LWID 0.038 0.169 0.022 -0.027 -0.009 -0.019Interactions with word attack and Word concreteness
M7 = glmer(GatesScore ~ zNarrativity+ zLWID*zWordConcreteness + zRefCohesion + zDeepCohesion
+ zSyntacticSimplicity
+ (1|ID),data = Corrected_EmotionPaper_Data_Quartiles, family = binomial)
summary(M7)## Generalized linear mixed model fit by maximum likelihood (Laplace
## Approximation) [glmerMod]
## Family: binomial ( logit )
## Formula:
## GatesScore ~ zNarrativity + zLWID * zWordConcreteness + zRefCohesion +
## zDeepCohesion + zSyntacticSimplicity + (1 | ID)
## Data: Corrected_EmotionPaper_Data_Quartiles
##
## AIC BIC logLik deviance df.resid
## 7753.4 7816.1 -3867.7 7735.4 7911
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -6.5047 -0.5897 0.3285 0.5325 2.8378
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 0.8212 0.9062
## Number of obs: 7920, groups: ID, 165
##
## Fixed effects:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 1.21009 0.07819 15.476 < 2e-16 ***
## zNarrativity 0.21032 0.04234 4.967 6.79e-07 ***
## zLWID 0.93103 0.07741 12.027 < 2e-16 ***
## zWordConcreteness -0.41155 0.04671 -8.810 < 2e-16 ***
## zRefCohesion -0.19462 0.04699 -4.142 3.45e-05 ***
## zDeepCohesion 0.35049 0.03087 11.355 < 2e-16 ***
## zSyntacticSimplicity -0.41006 0.05035 -8.145 3.81e-16 ***
## zLWID:zWordConcreteness -0.03411 0.02938 -1.161 0.246
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) zNrrtv zLWID zWrdCn zRfChs zDpChs zSyntS
## zNarrativty 0.029
## zLWID 0.065 0.019
## zWrdCncrtns -0.059 0.376 -0.033
## zRefCohesin -0.028 -0.588 -0.015 0.054
## zDeepCohesn 0.062 0.030 0.037 -0.239 -0.197
## zSyntctcSmp -0.053 -0.106 -0.027 0.614 0.582 -0.171
## zLWID:zWrdC -0.033 -0.043 -0.030 0.159 0.032 -0.026 0.049M8 = glmer(GatesScore ~ zNarrativity + zLWID*zWordConcreteness + zRefCohesion + zDeepCohesion
+ (1|ID),data = Corrected_EmotionPaper_Data_Quartiles, family = binomial)
summary(M8)## Generalized linear mixed model fit by maximum likelihood (Laplace
## Approximation) [glmerMod]
## Family: binomial ( logit )
## Formula:
## GatesScore ~ zNarrativity + zLWID * zWordConcreteness + zRefCohesion +
## zDeepCohesion + (1 | ID)
## Data: Corrected_EmotionPaper_Data_Quartiles
##
## AIC BIC logLik deviance df.resid
## 7818.9 7874.7 -3901.4 7802.9 7912
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -5.2421 -0.6064 0.3206 0.5422 3.0103
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 0.8071 0.8984
## Number of obs: 7920, groups: ID, 165
##
## Fixed effects:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 1.19132 0.07747 15.377 < 2e-16 ***
## zNarrativity 0.18047 0.04196 4.301 1.70e-05 ***
## zLWID 0.92327 0.07674 12.032 < 2e-16 ***
## zWordConcreteness -0.18642 0.03786 -4.923 8.51e-07 ***
## zRefCohesion 0.02999 0.03863 0.776 0.437
## zDeepCohesion 0.31284 0.02995 10.447 < 2e-16 ***
## zLWID:zWordConcreteness -0.02769 0.02963 -0.934 0.350
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) zNrrtv zLWID zWrdCn zRfChs zDpChs
## zNarrativty 0.019
## zLWID 0.065 0.017
## zWrdCncrtns -0.035 0.564 -0.022
## zRefCohesin 0.007 -0.652 0.000 -0.471
## zDeepCohesn 0.048 -0.033 0.034 -0.153 -0.094
## zLWID:zWrdC -0.032 -0.030 -0.026 0.209 0.002 0.011Interaction with Word attack and cohesions
M9 = glmer(GatesScore ~ zNarrativity + zWordConcreteness + zRefCohesion*zLWID + zDeepCohesion
+ zSyntacticSimplicity
+ (1|ID),data = Corrected_EmotionPaper_Data_Quartiles, family = binomial)
summary(M9)## Generalized linear mixed model fit by maximum likelihood (Laplace
## Approximation) [glmerMod]
## Family: binomial ( logit )
## Formula: GatesScore ~ zNarrativity + zWordConcreteness + zRefCohesion *
## zLWID + zDeepCohesion + zSyntacticSimplicity + (1 | ID)
## Data: Corrected_EmotionPaper_Data_Quartiles
##
## AIC BIC logLik deviance df.resid
## 7749.0 7811.8 -3865.5 7731.0 7911
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -6.2926 -0.5828 0.3294 0.5315 2.9871
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 0.8216 0.9064
## Number of obs: 7920, groups: ID, 165
##
## Fixed effects:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 1.20378 0.07814 15.405 < 2e-16 ***
## zNarrativity 0.21055 0.04229 4.979 6.38e-07 ***
## zWordConcreteness -0.39807 0.04605 -8.645 < 2e-16 ***
## zRefCohesion -0.20732 0.04692 -4.418 9.95e-06 ***
## zLWID 0.92805 0.07741 11.989 < 2e-16 ***
## zDeepCohesion 0.34631 0.03071 11.276 < 2e-16 ***
## zSyntacticSimplicity -0.40489 0.04996 -8.104 5.33e-16 ***
## zRefCohesion:zLWID -0.06982 0.02894 -2.413 0.0158 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) zNrrtv zWrdCn zRfChs zLWID zDpChs zSyntS
## zNarrativty 0.028
## zWrdCncrtns -0.051 0.394
## zRefCohesin -0.031 -0.582 0.041
## zLWID 0.063 0.018 -0.026 -0.018
## zDeepCohesn 0.062 0.031 -0.236 -0.194 0.035
## zSyntctcSmp -0.051 -0.100 0.616 0.570 -0.025 -0.167
## zRfChs:LWID 0.014 -0.019 -0.039 0.118 -0.007 0.037 -0.019M10 = glmer(GatesScore ~ zNarrativity + zWordConcreteness + zRefCohesion*zLWID + zDeepCohesion
+ (1|ID),data = Corrected_EmotionPaper_Data_Quartiles, family = binomial)
summary(M10)## Generalized linear mixed model fit by maximum likelihood (Laplace
## Approximation) [glmerMod]
## Family: binomial ( logit )
## Formula: GatesScore ~ zNarrativity + zWordConcreteness + zRefCohesion *
## zLWID + zDeepCohesion + (1 | ID)
## Data: Corrected_EmotionPaper_Data_Quartiles
##
## AIC BIC logLik deviance df.resid
## 7813.7 7869.5 -3898.9 7797.7 7912
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -5.0249 -0.5948 0.3238 0.5411 3.1640
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 0.8074 0.8985
## Number of obs: 7920, groups: ID, 165
##
## Fixed effects:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 1.18567 0.07742 15.314 < 2e-16 ***
## zNarrativity 0.18250 0.04193 4.352 1.35e-05 ***
## zWordConcreteness -0.17603 0.03695 -4.764 1.90e-06 ***
## zRefCohesion 0.01150 0.03911 0.294 0.7686
## zLWID 0.92093 0.07673 12.002 < 2e-16 ***
## zDeepCohesion 0.30989 0.02982 10.391 < 2e-16 ***
## zRefCohesion:zLWID -0.07206 0.02910 -2.476 0.0133 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) zNrrtv zWrdCn zRfChs zLWID zDpChs
## zNarrativty 0.018
## zWrdCncrtns -0.026 0.588
## zRefCohesin 0.002 -0.642 -0.479
## zLWID 0.063 0.017 -0.014 -0.006
## zDeepCohesn 0.049 -0.030 -0.158 -0.095 0.033
## zRfChs:LWID 0.012 -0.026 -0.028 0.178 -0.007 0.037M11 = glmer(GatesScore ~ zNarrativity + zWordConcreteness + zRefCohesion + zDeepCohesion*zLWID
+ zSyntacticSimplicity
+ (1|ID),data = Corrected_EmotionPaper_Data_Quartiles, family = binomial)
summary(M11)## Generalized linear mixed model fit by maximum likelihood (Laplace
## Approximation) [glmerMod]
## Family: binomial ( logit )
## Formula: GatesScore ~ zNarrativity + zWordConcreteness + zRefCohesion +
## zDeepCohesion * zLWID + zSyntacticSimplicity + (1 | ID)
## Data: Corrected_EmotionPaper_Data_Quartiles
##
## AIC BIC logLik deviance df.resid
## 7752.7 7815.5 -3867.4 7734.7 7911
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -6.5697 -0.5871 0.3267 0.5366 2.7336
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 0.8214 0.9063
## Number of obs: 7920, groups: ID, 165
##
## Fixed effects:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 1.21257 0.07827 15.492 < 2e-16 ***
## zNarrativity 0.20651 0.04247 4.863 1.16e-06 ***
## zWordConcreteness -0.40832 0.04648 -8.785 < 2e-16 ***
## zRefCohesion -0.19258 0.04718 -4.082 4.47e-05 ***
## zDeepCohesion 0.36114 0.03188 11.329 < 2e-16 ***
## zLWID 0.93267 0.07746 12.041 < 2e-16 ***
## zSyntacticSimplicity -0.41241 0.05026 -8.206 2.29e-16 ***
## zDeepCohesion:zLWID 0.04225 0.02978 1.419 0.156
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) zNrrtv zWrdCn zRfChs zDpChs zLWID zSyntS
## zNarrativty 0.026
## zWrdCncrtns -0.056 0.390
## zRefCohesin -0.027 -0.592 0.050
## zDeepCohesn 0.078 0.012 -0.259 -0.179
## zLWID 0.066 0.016 -0.031 -0.013 0.049
## zSyntctcSmp -0.055 -0.106 0.620 0.579 -0.183 -0.029
## zDpChs:LWID 0.050 -0.031 -0.089 0.006 0.254 0.041 -0.075M12 = glmer(GatesScore ~ zNarrativity + zWordConcreteness + zRefCohesion + zDeepCohesion*zLWID
+ (1|ID),data = Corrected_EmotionPaper_Data_Quartiles, family = binomial)
summary(M12)## Generalized linear mixed model fit by maximum likelihood (Laplace
## Approximation) [glmerMod]
## Family: binomial ( logit )
## Formula: GatesScore ~ zNarrativity + zWordConcreteness + zRefCohesion +
## zDeepCohesion * zLWID + (1 | ID)
## Data: Corrected_EmotionPaper_Data_Quartiles
##
## AIC BIC logLik deviance df.resid
## 7819.1 7875.0 -3901.6 7803.1 7912
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -5.5575 -0.6049 0.3217 0.5402 2.9580
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 0.8071 0.8984
## Number of obs: 7920, groups: ID, 165
##
## Fixed effects:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 1.19171 0.07751 15.374 < 2e-16 ***
## zNarrativity 0.17700 0.04213 4.202 2.65e-05 ***
## zWordConcreteness -0.17989 0.03707 -4.852 1.22e-06 ***
## zRefCohesion 0.03260 0.03890 0.838 0.402
## zDeepCohesion 0.31862 0.03071 10.375 < 2e-16 ***
## zLWID 0.92338 0.07675 12.031 < 2e-16 ***
## zDeepCohesion:zLWID 0.02297 0.02956 0.777 0.437
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) zNrrtv zWrdCn zRfChs zDpChs zLWID
## zNarrativty 0.014
## zWrdCncrtns -0.027 0.589
## zRefCohesin 0.011 -0.656 -0.483
## zDeepCohesn 0.060 -0.053 -0.168 -0.065
## zLWID 0.065 0.013 -0.016 0.004 0.042
## zDpChs:LWID 0.044 -0.073 -0.037 0.085 0.225 0.031plot(effect("zRefCohesion:zLWID", M9), multiline=TRUE)
Interaction with LWID and syntactic simplicity
M13 = glmer(GatesScore ~ zNarrativity + zWordConcreteness + zRefCohesion + zDeepCohesion + zSyntacticSimplicity*zLWID
+ (1|ID),data = Corrected_EmotionPaper_Data_Quartiles, family = binomial)
summary(M13)## Generalized linear mixed model fit by maximum likelihood (Laplace
## Approximation) [glmerMod]
## Family: binomial ( logit )
## Formula: GatesScore ~ zNarrativity + zWordConcreteness + zRefCohesion +
## zDeepCohesion + zSyntacticSimplicity * zLWID + (1 | ID)
## Data: Corrected_EmotionPaper_Data_Quartiles
##
## AIC BIC logLik deviance df.resid
## 7751.9 7814.7 -3866.9 7733.9 7911
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -6.2420 -0.5841 0.3278 0.5311 2.8739
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 0.8212 0.9062
## Number of obs: 7920, groups: ID, 165
##
## Fixed effects:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 1.20563 0.07813 15.431 < 2e-16 ***
## zNarrativity 0.20722 0.04228 4.902 9.50e-07 ***
## zWordConcreteness -0.39968 0.04607 -8.675 < 2e-16 ***
## zRefCohesion -0.19019 0.04692 -4.053 5.05e-05 ***
## zDeepCohesion 0.34633 0.03079 11.247 < 2e-16 ***
## zSyntacticSimplicity -0.39363 0.05075 -7.757 8.72e-15 ***
## zLWID 0.92826 0.07739 11.995 < 2e-16 ***
## zSyntacticSimplicity:zLWID 0.04949 0.02935 1.686 0.0918 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) zNrrtv zWrdCn zRfChs zDpChs zSyntS zLWID
## zNarrativty 0.028
## zWrdCncrtns -0.051 0.392
## zRefCohesin -0.027 -0.585 0.047
## zDeepCohesn 0.061 0.026 -0.239 -0.198
## zSyntctcSmp -0.048 -0.103 0.609 0.579 -0.178
## zLWID 0.063 0.017 -0.027 -0.013 0.035 -0.022
## zSyntS:LWID -0.009 -0.014 0.037 0.034 -0.061 0.154 0.003Interaction with Word attack and emotion
M14 = glmer(GatesScore ~ zNarrativity + zWordConcreteness + zRefCohesion + zDeepCohesion + zSyntacticSimplicity + zLWID*zArousal
+ (1|ID),data = Corrected_EmotionPaper_Data_Quartiles, family = binomial)
summary(M14)## Generalized linear mixed model fit by maximum likelihood (Laplace
## Approximation) [glmerMod]
## Family: binomial ( logit )
## Formula: GatesScore ~ zNarrativity + zWordConcreteness + zRefCohesion +
## zDeepCohesion + zSyntacticSimplicity + zLWID * zArousal +
## (1 | ID)
## Data: Corrected_EmotionPaper_Data_Quartiles
##
## AIC BIC logLik deviance df.resid
## 7715.9 7785.6 -3847.9 7695.9 7910
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -6.4256 -0.5705 0.3203 0.5325 2.8358
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 0.8305 0.9113
## Number of obs: 7920, groups: ID, 165
##
## Fixed effects:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 1.21893 0.07862 15.504 < 2e-16 ***
## zNarrativity 0.25267 0.04358 5.798 6.72e-09 ***
## zWordConcreteness -0.44868 0.04655 -9.638 < 2e-16 ***
## zRefCohesion -0.08233 0.05057 -1.628 0.1035
## zDeepCohesion 0.26774 0.03368 7.949 1.88e-15 ***
## zSyntacticSimplicity -0.32844 0.05095 -6.447 1.14e-10 ***
## zLWID 0.93888 0.07786 12.059 < 2e-16 ***
## zArousal 0.25422 0.04212 6.035 1.59e-09 ***
## zLWID:zArousal 0.09425 0.02910 3.239 0.0012 **
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) zNrrtv zWrdCn zRfChs zDpChs zSyntS zLWID zArosl
## zNarrativty 0.037
## zWrdCncrtns -0.055 0.333
## zRefCohesin -0.014 -0.480 -0.032
## zDeepCohesn 0.042 -0.035 -0.117 -0.311
## zSyntctcSmp -0.042 -0.072 0.532 0.613 -0.239
## zLWID 0.065 0.022 -0.031 -0.004 0.024 -0.020
## zArousal 0.043 0.177 -0.205 0.362 -0.415 0.232 0.032
## zLWID:zArsl 0.027 0.036 -0.013 0.008 -0.022 -0.009 0.037 0.190M15 = glmer(GatesScore ~ zNarrativity + zWordConcreteness + zRefCohesion + zDeepCohesion + zSyntacticSimplicity + zLWID*zValence
+ (1|ID),data = Corrected_EmotionPaper_Data_Quartiles, family = binomial)
summary(M15)## Generalized linear mixed model fit by maximum likelihood (Laplace
## Approximation) [glmerMod]
## Family: binomial ( logit )
## Formula: GatesScore ~ zNarrativity + zWordConcreteness + zRefCohesion +
## zDeepCohesion + zSyntacticSimplicity + zLWID * zValence +
## (1 | ID)
## Data: Corrected_EmotionPaper_Data_Quartiles
##
## AIC BIC logLik deviance df.resid
## 7745.9 7815.6 -3862.9 7725.9 7910
##
## Scaled residuals:
## Min 1Q Median 3Q Max
## -6.6288 -0.5879 0.3256 0.5342 3.2778
##
## Random effects:
## Groups Name Variance Std.Dev.
## ID (Intercept) 0.8245 0.908
## Number of obs: 7920, groups: ID, 165
##
## Fixed effects:
## Estimate Std. Error z value Pr(>|z|)
## (Intercept) 1.21100 0.07834 15.459 < 2e-16 ***
## zNarrativity 0.16619 0.04467 3.720 0.000199 ***
## zWordConcreteness -0.46877 0.05071 -9.245 < 2e-16 ***
## zRefCohesion -0.19768 0.04702 -4.204 2.62e-05 ***
## zDeepCohesion 0.42470 0.03884 10.934 < 2e-16 ***
## zSyntacticSimplicity -0.55103 0.06797 -8.107 5.20e-16 ***
## zLWID 0.93327 0.07756 12.033 < 2e-16 ***
## zValence -0.15387 0.05296 -2.906 0.003667 **
## zLWID:zValence 0.03830 0.02974 1.288 0.197819
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Correlation of Fixed Effects:
## (Intr) zNrrtv zWrdCn zRfChs zDpChs zSyntS zLWID zValnc
## zNarrativty 0.023
## zWrdCncrtns -0.053 0.467
## zRefCohesin -0.028 -0.552 0.056
## zDeepCohesn 0.057 -0.169 -0.417 -0.166
## zSyntctcSmp -0.045 0.144 0.691 0.441 -0.500
## zLWID 0.065 0.014 -0.033 -0.015 0.041 -0.030
## zValence -0.003 0.322 0.392 0.019 -0.583 0.664 -0.008
## zLWID:zVlnc 0.034 0.029 -0.077 -0.032 0.102 -0.064 0.033 0.091plot(effect("zLWID:zArousal", M14), multiline=TRUE)