clmm不计算随机效应

时间:2019-11-07 15:22:15

标签: r statistics mixed-models

对我的研究参与者(n = N)排名4张图片,产生4 * N个案例。 现在我想为参与者和图片进行控制。我包括了joinerID(here = PerNO)和pictureID,但仅适用于PicID。我做错了什么?

非常感谢您提出任何想法。

在此处显示带有参与者ID的输出,之后您将找到具有PictureID的有效输出

> summary(var.with.random)
Cumulative Link Mixed Model fitted with the Laplace approximation

formula: Int ~ PicSex + Posture + (1 | PerNO)
data:    x

link  threshold nobs logLik  AIC     niter    max.grad cond.H
logit flexible  668  -917.62 1847.24 423(331) 3.80e-06 NaN   

Random effects:
Groups Name        Variance Std.Dev. 
PerNO  (Intercept) 5.38e-07 0.0007335
Number of groups:  PerNO 167 

Coefficients:
     Estimate Std. Error z value Pr(>|z|)
PicSex1   0.4107         NA      NA       NA
Posture  -0.4035         NA      NA       NA

Threshold coefficients:

 Estimate Std. Error z value
1|2 -1.12189         NA      NA
2|3 -0.01118         NA      NA
3|4  1.11189         NA      NA 
```

**Output with PictureID:**
```
> var.with.random <- clmm(Int~PicSex + Posture +(1|PicID), data = x)
> summary(var.with.random)
Cumulative Link Mixed Model fitted with the Laplace approximation

formula: Int ~ PicSex + Posture + (1 | PicID)
data:    x

link  threshold nobs logLik  AIC     niter    max.grad cond.H 
logit flexible  668  -896.18 1804.35 317(954) 7.20e-04 1.2e+02

Random effects:
Groups Name        Variance Std.Dev.
PicID  (Intercept) 0.3262   0.5711  
Number of groups:  PicID 16 

Coefficients:
     Estimate Std. Error z value Pr(>|z|)
PicSex1   0.4428     0.3195   1.386    0.166
Posture  -0.4528     0.3198  -1.416    0.157

Threshold coefficients:
 Estimate Std. Error z value
1|2 -1.21490    0.28418  -4.275
2|3 -0.01622    0.27953  -0.058
3|4  1.18597    0.28276   4.194

0 个答案:

没有答案