1
$\begingroup$

I am using this R package "sjPlot" to plot the forest-plot of standardized beta values. My model includes interaction terms

m1.3.5=lmer(Q3~1+Q4+gender+gcOUT*gender+maxXcorr*gender+(1+Q4+gcOUT*gender+maxXcorr*gender|sub))

The "gender" variable is binary, and the others are continuous variables.

The results of Anova (car package)

Anova(m1.3.5)
Analysis of Deviance Table (Type II Wald chisquare tests)

Response: Q3
                   Chisq Df Pr(>Chisq)    
Q4              228.6415  1  < 2.2e-16 ***
gender            0.4101  1   0.521927    
gcOUT             7.8720  1   0.005021 ** 
maxXcorr          6.2570  1   0.012370 *  
gender:gcOUT      2.3550  1   0.124884    
gender:maxXcorr   0.5016  1   0.478812  

Here is the plot:

plot_model(m1.3.5, type = "std",se=TRUE,line.size=1.5,dot.size=4,show.values = TRUE)

enter image description here

The plot and the Anova() results seems to be quite different. For example, the gcOUT effect is significant in Anova. However, the SE of the plot of that variable crossed 0, and it does not have *, suggesting it's not significant.

May I ask why the results of these 2 were inconsistent?

Note that if without interaction variables, the results of these two are the same.

Thank you all for your great help in advance :)

$\endgroup$

0

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Browse other questions tagged or ask your own question.