Je tente de modéliser une variable ordinale à 3 modalités par 2 prédicteurs nominaux. J'utilise la fonction clm() du package ordinal:
Code : Tout sélectionner
fit.clm<-clm(RindFleshColorContrast~Variety+Repetition)
Je voudrais savoir si certains d'entre vous ont déjà eu cette erreur:
Code : Tout sélectionner
Warning message:
(1) Hessian is numerically singular: parameters are not uniquely determined
In addition: Absolute convergence criterion was met, but relative criterion was not met
Je pense que cette erreur est dûe aux effectifs des modalités de la variable à expliquer pour chaque modalité du prédicteur "variety":
Code : Tout sélectionner
> table(RindFleshColorContrast,Variety)
Variety
RindFleshColorContrast V1 V10 V11 V12 V13 V14 V15 V16 V17 V18 V19 V2 V20 V21 V22 V23 V24 V25 V26 V27 V28 V29 V3 V30
2 0 1 0 3 1 0 0 0 4 0 1 0 1 4 13 2 0 3 1 4 6 0 1 0
3 1 6 0 3 4 6 10 3 19 2 5 3 8 20 8 10 1 3 4 0 3 0 3 6
4 21 3 12 3 7 5 2 9 1 1 5 19 3 0 3 0 11 0 0 0 2 11 2 2
Variety
RindFleshColorContrast V31 V32 V33 V34 V35 V36 V37 V38 V39 V4 V40 V41 V42 V43 V44 V45 V5 V6 V7 V8 V9
2 0 0 1 0 1 2 0 0 0 0 0 0 0 3 0 7 2 1 0 0 1
3 7 7 4 9 3 9 11 8 0 18 3 12 0 16 6 5 10 2 7 3 3
4 0 4 3 2 4 1 1 4 12 6 9 4 20 5 10 0 0 3 0 2 4
La modalité 2 est beaucoup moins représentée que les deux autres.
De plus, les effectifs des modalités de la variable à expliquer ne sont pas équilibrés:
Code : Tout sélectionner
> table(RindFleshColorContrast)
RindFleshColorContrast
2 3 4
63 271 216
Je ne sais pas comment m'y prendre pour régler ce problème... Faudrait-il rassembler les modalités?
Si quelqu'un a une idée je suis preneuse.
Merci!
Chloé