Des données sont toujours disponible à cette adresse:
ftp://ftp.cirad.fr/pub/group-r/groupe-r ... metrie.txt
J'utilise les contrasts "contr.sum"
Sur ce jeu de données on peut regarder le modèle suivant :
lm1<-lm(log(mstoter)~log(dbm)+variete)
Je n’étudie pas ce modèle en raison du déséquilibre du nombre de mesures mais à titre d’exemple cela permet de voir le problème rencontré avec la fonction dummy.coef().
> coef(lm1)
(Intercept) log(dbm) variete1 variete2 variete3 variete4 variete5
-2.93441516 2.68139874 0.15415229 -0.15055904 0.12163052 -0.12665775 -0.01721356
>dummy.coef(lm1)
Full coefficients are
(Intercept): -2.934415
log(dbm): 0
variete: jose irwi cogs kens tomm kent
0.15415229 -0.15055904 0.12163052 -0.12665775 -0.01721356 0.01864755
Lorsque je lançait un modèle similaire sur l’ancienne version de R la valeur de log(dbm) apparaissait dans les dummy.coef() comme dans la fonction coef().
Je ne sais pas si cela a un rapport avec la nouvelle version de R mais à part ça je n’ai rien changé ni à mes analyses ni à mes données.
Est ce qu’il existe une autre fonction que dummy.coef() pour récupérer tous les coefficients ?
Je vous remercie de votre aide.
Bien Cordialement
Claire