对我的研究参与者(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