problème de niveaux sur une ANOVA

Version imprimable

Bonjour à tous !

Je cherche à déterminer quels facteurs influencent le plus une teneur en champignon (lnEF1a_Fgraminearum, qui est en log) après une récolte ... qui pourrait s'expliquer par plusieurs variables (toutes qualitatives). Je décide donc de faire une ANOVA

Voici un bout de mon script :
Code:

1 2 3 4 5 6 res.lm <- lm( lnEF1a_Fgraminearum ~ -1 + V3 + V4 + V21 + V23 + V16 + V17 + V21 + V35 , data = data) anova(res.lm) summary(res.lm) residuals(res.lm) hist(residuals(res.lm)) curve(dnorm(x,0,sd=sd(residuals(res.lm))), add=T, lwd=2, col="red")
J'ai quelques soucis : lorsque je veux rajouter une variable explicative, j'ai le message d'erreur suivant :
Erreur dans `contrasts<-`(`*tmp*`, value = contr.funs[1 + isOF[nn]]) :
les contrastes ne peuvent être appliqués qu'aux facteurs ayant au moins deux niveaux

Je ne sais pas pourquoi R refuse que je rajoute certaines variables, avez vous une idée ?

Merci

11/04/2013, 12h08
Sengar

Citation:

Envoyé par Fnouch

Bonjour à tous !
Erreur dans `contrasts<-`(`*tmp*`, value = contr.funs[1 + isOF[nn]]) :
les contrastes ne peuvent être appliqués qu'aux facteurs ayant au moins deux niveaux

Je ne sais pas pourquoi R refuse que je rajoute certaines variables, avez vous une idée ?

Hello,

Vu l'erreur je pense qu'avant toute chose il serait bon que tu postes le type de la nouvelle variable que tu ajoutes, un table() de cette dernière, et un levels().

Ma variable V18 par exemple, quand je la mets, l'anova bug !

Code:

1
2
3
4
5
6
> table(data$V18)
 
    BROYES BROYES SPE     BRULES    ENFOUIS   EXPORTES NON BROYES 
       551         33          5         29         68         85 
> levels(data$V18)
[1] "BROYES"     "BROYES SPE" "BRULES"     "ENFOUIS"    "EXPORTES"   "NON BROYES"

Citation:

Envoyé par Fnouch

Erreur dans `contrasts<-`(`*tmp*`, value = contr.funs[1 + isOF[nn]]) :
les contrastes ne peuvent être appliqués qu'aux facteurs ayant au moins deux niveaux

:ccool:
Code:

1 2 3 4 5 6 > table(data$V18) BROYES BROYES SPE BRULES ENFOUIS EXPORTES NON BROYES 551 33 5 29 68 85 > levels(data$V18) [1] "BROYES" "BROYES SPE" "BRULES" "ENFOUIS" "EXPORTES" "NON BROYES"
:aie:

Le agatha christie du R :mouarf:

Allez jtente mon hercule poirot:
Au moins une de tes autres variables dans le lm a des valeurs manquantes sur tous les facteurs de V18 sauf un.

Edit: pour vérifier si c'est ca le problème tu fais:

Code:

table(data$V18[!is.na(data$V1) & !is.na(data$V2) & ....])

...etc pour toutes les variables explicatives.