Répondre à la discussion
Affichage des résultats 1 à 5 sur 5

[R] Problème sur une ANOVA



  1. #1
    petitelectron

    Question [R] Problème sur une ANOVA


    ------

    Bonjour à tous !

    je vous explique mon sujet : 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 ! (je ne sais pas si c'est la meilleure idée mais bon..)

    Voici un bout de mon script :
    Code:
    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

    1) Je ne sais pas pourquoi R refuse que je rajoute certaines variables, avez vous une idée ?
    2) Lorsque je fais mon histogramme de la distribution des résidus, la fonction curve() que j'applique renvois une courbe qui ne suit pas du tout une loi normale
    residu anova.jpeg

    j'ai du rater quelque chose dans ma ligne curve()


    ...

    Merci !

    -----

  2. Publicité
  3. #2
    toothpick-charlie

    Re : [R] Problème sur une ANOVA

    le message d'erreur laisse penser qu'une variable du modèle est constante.

    et quant à la densité gaussienne, il faut la multiplier par le nombre de résidus (qui est le nombre d'observations).

  4. #3
    petitelectron

    Re : [R] Problème sur une ANOVA

    Euu merci ca me donne une bien plus jolie répartition (même si je suis pas sûre du tout que ce soit "autorisé")
    mais ca ne change rien à ma curve !!! :'(

  5. #4
    toothpick-charlie

    Re : [R] Problème sur une ANOVA

    il ne faut pas utiliser la fonction curve.

    ou sinon, tu peux faire

    hist(residuals(res.lm),freq=F)
    curve(dnorm(x,0,sd=sd(residual s(res.lm))), add=T, lwd=2, col="red")

  6. #5
    petitelectron

    Re : [R] Problème sur une ANOVA

    Ça ne marche pas non plus ...

    J'avais l'habitude d'utiliser curve() et tout allait bien ..

    je vais essayer avec line() ou abline()

  7. A voir en vidéo sur Futura

Sur le même thème :

Discussions similaires

  1. gros problème ANOVA deux facteurs
    Par xpressmusic dans le forum Mathématiques du supérieur
    Réponses: 0
    Dernier message: 08/04/2012, 16h17
  2. problème ANOVA sous R
    Par holaquimica dans le forum Mathématiques du supérieur
    Réponses: 4
    Dernier message: 29/06/2011, 14h21
  3. Problème pour réaliser une ANOVA
    Par Foufy68 dans le forum Mathématiques du supérieur
    Réponses: 14
    Dernier message: 10/06/2009, 10h44
  4. Statistique sur R (anova)
    Par Kisssmo dans le forum Mathématiques du supérieur
    Réponses: 1
    Dernier message: 27/05/2009, 16h22
  5. question sur le nombre minimum d'observations à avoir pour pouvoir faire une ANOVA?
    Par aliciavasseur dans le forum Mathématiques du supérieur
    Réponses: 2
    Dernier message: 22/05/2009, 11h11
Découvrez nos comparatifs produits sur l'informatique et les technologies.