J'ai un problème avec la 2l.norm
méthode d'imputation multiniveau dans mice
.
Malheureusement, je ne peux pas publier d'exemple reproductible en raison de la taille de mes données - lorsque je réduis la taille, le problème disparaît.
Pour une variable particulière, mice
génère les erreurs et avertissements suivants:
Error in chol.default(inv.sigma2[class] * X.SS[[class]] + inv.psi) :
the leading minor of order 1 is not positive definite
In addition: Warning messages:
1: In rgamma(n.class, n.g/2 + 1/(2 * theta), scale = 2 * theta/(ss * :
NAs produced
2: In rgamma(1, n.class/(2 * theta) + 1, scale = 2 * theta * H/n.class) :
NAs produced
3: In rgamma(1, n.class/2 - 1, scale = 2/(n.class * (sigma2.0/H - log(sigma2.0) + :
NAs produced
Si j'utilise les 2l.pan
, norm
ou les pmm
méthodes, le problème ne se produit pas.
La variable a la distribution suivante:
Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
50.0 117.0 136.0 136.7 155.0 249.0 3124
De plus, les tailles de classe ont la distribution suivante:
Min. 1st Qu. Median Mean 3rd Qu. Max.
3.00 50.00 80.00 88.52 111.00 350.00
la source