Question

我需要帮助解释实验性2x2阶乘设计的平均边际效应。我已经应用了一个lmer模型，该模型呈现出不重要的主效应，但是却具有显着的交互作用。

因此，我想研究一下边际效应，以了解是什么驱动了这种互动效应。

我已经使用了R的margins包。我设法计算了平均边际效应。

library(lme4)
library(margins)

model_cre <- lmer(AVG_CREATIVITY~IV_PTP*IV_PR + (1|Studentno), sl_data)
summary(model_cre)

margin_cre <- margins(model_cre, at=list(IV_PTP=0:1, IV_PR=0:1))
summary(margin_cre)

Linear mixed model fit by REML ['lmerMod']
Formula: AVG_CREATIVITY ~ IV_PTP * IV_PR + (1 | Studentno)
  Data: sl_data

REML criterion at convergence: 3345

Scaled residuals: 
   Min      1Q  Median      3Q     Max 
-2.6399 -0.6239 -0.0822  0.4710  4.2569 

Random effects:
Groups    Name        Variance Std.Dev.
Studentno (Intercept) 0.09803  0.3131  
Residual              0.65872  0.8116  
Number of obs: 1326, groups:  Studentno, 178

Fixed effects:
             Estimate Std. Error t value
(Intercept)   2.948545   0.069900  42.183
IV_PTP       -0.093376   0.099162  -0.942
IV_PR         0.008009   0.096209   0.083
IV_PTP:IV_PR  0.282937   0.135182   2.093

Correlation of Fixed Effects:
           (Intr) IV_PTP IV_PR 
IV_PTP      -0.705              
IV_PR       -0.727  0.512       
IV_PTP:IV_P  0.517 -0.734 -0.712

Average Marginal Effects
factor IV_PTP  IV_PR     AME     SE       z      p   lower  upper
 IV_PR 0.0000 0.0000  0.0080 0.0962  0.0833 0.9337 -0.1806 0.1966
 IV_PR 0.0000 1.0000  0.0080 0.0962  0.0833 0.9337 -0.1806 0.1966
 IV_PR 1.0000 0.0000  0.2909 0.0950  3.0638 0.0022  0.1048 0.4771
 IV_PR 1.0000 1.0000  0.2909 0.0950  3.0638 0.0022  0.1048 0.4771
IV_PTP 0.0000 0.0000 -0.0934 0.0992 -0.9417 0.3464 -0.2877 0.1010
IV_PTP 0.0000 1.0000  0.1896 0.0919  2.0632 0.0391  0.0095 0.3696
IV_PTP 1.0000 0.0000 -0.0934 0.0992 -0.9417 0.3464 -0.2877 0.1010
IV_PTP 1.0000 1.0000  0.1896 0.0919  2.0633 0.0391  0.0095 0.3696

我对R非常陌生（一周后），我想我能从中得出的结论是：

factor IV_PTP  IV_PR     AME     SE       z      p   lower  upper
 IV_PR 1.0000 0.0000  0.2909 0.0950  3.0638 0.0022**  0.1048 0.4771
 IV_PR 1.0000 1.0000  0.2909 0.0950  3.0638 0.0022**  0.1048 0.4771

表示当IV_PTP和IV_PR = 1时-与仅IV_PTP = 1时相比，这会产生更高的平均创造力。

但是，我对“平均边际效应”摘要显示8行这一事实感到困惑。 2x2设计的结果是在4个条件下得出的，而不是8个。我想知道为什么最后4行的结果不同，即使它应该表示IV_PTP和IV_PR的相同组合。

也有重要结果。

factor IV_PTP  IV_PR     AME     SE       z      p   lower  upper 
IV_PTP 0.0000 1.0000  0.1896 0.0919  2.0632 0.0391*  0.0095 0.3696
IV_PTP 1.0000 1.0000  0.1896 0.0919  2.0633 0.0391*  0.0095 0.3696

那是否还表明与仅IV_PR = 1时相比，IV_PTP和IV_PR都为1时，创造力会更高？

希望您能理解我的问题-对R来说是非常陌生的-让我知道是否以及如何可以提供有关数据的其他信息。

R：lmer模型的平均边际效应的解释

0 个答案: