clmm не вычисляет случайный эффект - PullRequest
Новая
       20

clmm не вычисляет случайный эффект

1 голос
/ 07 ноября 2019

для моих участников исследования (n = N) ранжировались 4 рисунка, что привело к 4 * N случаям. Теперь я хочу контролировать участников и фотографии. Я включил идентификатор участника (здесь = PerNO) и pictureID, но он работает только для PicID. Что я делаю не так?

Большое спасибо за любую идею.

Здесь вывод с memberID, после чего вы найдете рабочий вывод с 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
...