0

我正在使用 anova() 函数比较 R 中的两个多级模型。一个模型包含一个控制变量，另一个模型包含一个实验变量。当我比较这两者时，我得到一个奇怪的结果，其中卡方为 0，p 值为 1。我将其解释为模型没有显着差异，但这对于我的数据和其他分析没有意义已经完成了这个实验变量。有人可以帮我理解这个输出吗？

为了解释变量，block_order（控制）是问题的平衡。这是一个有5个级别的因素。team_num 是 2 级随机效果；这是他们所属的参与者的团队。cent_team_wm_agg 是球队保持健康体重的愿望。它是一个连续变量。exer_vig 是连续因变量，它是人们锻炼的频率。

这是让我感到困惑的模型比较输出：

anova(m2_ev_full_team, m1_ev_control_block_team) 

refitting model(s) with ML (instead of REML)
Data: clean_data_0_nona
Models:
m2_ev_full_team: exer_vig ~ 1 + cent_team_wm_agg + (1 | team_num)
m1_ev_control_block_team: exer_vig ~ 1 + block_order + (1 | team_num)
                         Df    AIC    BIC  logLik deviance Chisq Chi Df Pr(>Chisq)
m2_ev_full_team           4 523.75 536.27 -257.88   515.75                        
m1_ev_control_block_team  8 533.96 559.00 -258.98   517.96     0      4          1

如果这有帮助，这里是模型本身。这是带有实验变量的那个：

summary(m2_ev_full_team <- lmer(exer_vig ~ 1 + cent_team_wm_agg + (1 |team_num), data = clean_data_0_nona))
Linear mixed model fit by REML. t-tests use Satterthwaite's method ['lmerModLmerTest']
Formula: exer_vig ~ 1 + cent_team_wm_agg + (1 | team_num)
   Data: clean_data_0_nona

REML criterion at convergence: 519.7

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-1.7585 -0.5819 -0.2432  0.5531  2.5569 

Random effects:
 Groups   Name        Variance Std.Dev.
 team_num (Intercept) 0.1004   0.3168  
 Residual             1.1628   1.0783  
Number of obs: 169, groups:  team_num, 58

Fixed effects:
                 Estimate Std. Error       df t value Pr(>|t|)    
(Intercept)       2.65955    0.09478 42.39962  28.061   <2e-16 ***
cent_team_wm_agg  0.73291    0.23572 64.27148   3.109   0.0028 ** 
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr)
cnt_tm_wm_g -0.004

和一个控制：

summary(m1_ev_control_block_team <- lmer(exer_vig ~ 1 + block_order + (1  |team_num), data = clean_data_0_nona))
Linear mixed model fit by REML. t-tests use Satterthwaite's method ['lmerModLmerTest']
Formula: exer_vig ~ 1 + block_order + (1 | team_num)
   Data: clean_data_0_nona

REML criterion at convergence: 525.1

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-1.6796 -0.6597 -0.1625  0.5291  2.0941 

Random effects:
 Groups   Name        Variance Std.Dev.
 team_num (Intercept) 0.2499   0.4999  
 Residual             1.1003   1.0490  
Number of obs: 169, groups:  team_num, 58

Fixed effects:
                                Estimate Std. Error       df t value Pr(>|t|)    
(Intercept)                       3.0874     0.2513 155.4960  12.284   <2e-16 ***
block_orderBlock2|Block4|Block3  -0.2568     0.3057 154.8652  -0.840   0.4020    
block_orderBlock3|Block2|Block4  -0.3036     0.3438 160.8279  -0.883   0.3785    
block_orderBlock3|Block4|Block2  -0.6204     0.3225 161.5186  -1.924   0.0561 .  
block_orderBlock4|Block2|Block3  -0.4215     0.3081 151.2908  -1.368   0.1733    
block_orderBlock4|Block3|Block2  -0.7306     0.3178 156.5548  -2.299   0.0228 *  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Correlation of Fixed Effects:
            (Intr) b_B2|B b_B3|B2 b_B3|B4 b_B4|B2
bl_B2|B4|B3 -0.757                               
bl_B3|B2|B4 -0.687  0.557                        
bl_B3|B4|B2 -0.733  0.585  0.543                 
bl_B4|B2|B3 -0.741  0.601  0.545   0.577         
bl_B4|B3|B2 -0.734  0.586  0.535   0.561   0.575

编辑：如果我不得不猜测，我认为这是因为控制模型比实验模型具有更多的自由度，但这就是我能想到的。我已经尝试在翻转模型顺序的情况下运行 anova，但它并没有改变任何东西。如果是这种情况，我不知道为什么 dfs 的数量会影响比较哪个更好。

谢谢！

4

0 回答 0