Problème d'estimation des coefficients clm()

Postez ici vos questions, réponses, commentaires ou suggestions - Les sujets seront ultérieurement répartis dans les archives par les modérateurs

Modérateur : Groupe des modérateurs

Chloé Magnin
Messages : 16
Enregistré le : 24 Avr 2017, 06:45

Problème d'estimation des coefficients clm()

Messagepar Chloé Magnin » 08 Juil 2017, 13:46

Bonjour à tous !

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é

Florent Aubry
Messages : 324
Enregistré le : 25 Juin 2010, 10:21

Re: Problème d'estimation des coefficients clm()

Messagepar Florent Aubry » 10 Juil 2017, 08:41

A vue de nez, RindFleshColorContrast et Variety ne sont pas indépendants entre eux. ce qui expliquerait le message "Hessian is numerically singular: parameters are not uniquely determined".


Retourner vers « Questions en cours »

Qui est en ligne

Utilisateurs parcourant ce forum : Google [Bot] et 1 invité