Statistique sur R (anova)
Répondre à la discussion
Affichage des résultats 1 à 2 sur 2

Statistique sur R (anova)



  1. #1
    invite8b118ae7

    Exclamation Statistique sur R (anova)


    ------

    Bonjour,


    je suis étudiant et actuellement en stage... et en plein désespoir..

    J'ai enfin fini d'entrer mes données sur excel... mais probleme: mes tests ne fonctionne pas! et j'ai un probleme pour la rédaction de l'un d'eux...


    Déjà le choix du test:

    je compare 2 espèces qui ont subis 6 traitements différents et j'ai relevé différents caractères. (oui c'est de la bio.. je sais pas si je peux poster ici.. mais je tente quand même ^^ )

    sur excel mes colones sont:
    - espèce: AT ou AJ (initial de leur nom scientifique..)
    - traitement: le n° code du traitement
    - nom de l'individu
    - masse
    - taille
    - longévité
    - fécondité
    etc


    Donc je pensais faire une ANOVA de tout ça sur R mais il faut vraiment que je tienne compte de l'espèce et du traitement à chaque fois... Je n'arrive pas à réussir mon test de normalité et surtout d'homoscédasticité!!

    En fait, j'avais ma matrice sur R. J'ai fait une ANOVA pour la longévité en tapant ceci sur R (ici en plus clair):

    > Longévité=aov(longévité~traite ment*espèce)

    Ensuite:
    > summary(longévité)
    Df Sum Sq Mean Sq F value Pr(>F)
    traitement 3 7162 2387 2.9114 0.03654 *
    espèce 1 16824 16824 20.5174 1.219e-05 ***
    traitement:espèce 3 1290 430 0.5243 0.66623
    Residuals 146 119719 820


    Donc là 1er problème! p<<0,05!

    Au cas où j'essaie mon test de Levene:
    levene.test(Longévité,interact ion(traitement*espèce))

    là ça marche pas du tout! Comment je peux faire pour que mon test de Levene fonctionne? Comment je peux lier le traitement et l'espèce pour mon test?

    Voilà j'espère que vous pourrez m'aider
    Merci d'avance

    -----

  2. #2
    invite8b118ae7

    Exclamation Re : Statistique sur R (anova)

    dsl petit probleme dans mes explications^^

    et le menu editer n'a pas voulu fonctionner....

    voilà le "vrai" message


    Bonjour,


    je suis étudiant et actuellement en stage... et en plein désespoir..

    J'ai enfin fini d'entrer mes données sur excel... mais probleme: mes tests ne fonctionne pas! et j'ai un probleme pour la rédaction de l'un d'eux...


    Déjà le choix du test:

    je compare 2 espèces qui ont subis 6 traitements différents et j'ai relevé différents caractères. (oui c'est de la bio.. je sais pas si je peux poster ici.. mais je tente quand même ^^ )

    sur excel mes colones sont:
    - espèce: AT ou AJ (initial de leur nom scientifique..)
    - traitement: le n° code du traitement
    - nom de l'individu
    - masse
    - taille
    - longévité
    - fécondité
    etc

    (PS quand je n'ai pas de données pour certains individus (soit en longévité soit en fécondité par exemple) je mets des NA... donc plein de NA dans mon tableau)
    Donc je pensais faire une ANOVA de tout ça sur R mais il faut vraiment que je tienne compte de l'espèce et du traitement à chaque fois... Je n'arrive pas à réussir mon test de normalité et surtout d'homoscédasticité!!

    En fait, j'avais ma matrice sur R. J'ai fait une ANOVA pour la longévité en tapant ceci sur R (ici en plus clair):

    > Longévité=aov(longévité~traite ment*espèce)

    Ensuite:
    > summary(Longévité)
    Df Sum Sq Mean Sq F value Pr(>F)
    traitement 3 7162 2387 2.9114 0.03654 *
    espèce 1 16824 16824 20.5174 1.219e-05 ***
    traitement:espèce 3 1290 430 0.5243 0.66623
    Residuals 146 119719 820

    Je fais mon test de Shapiro:
    shapiro.test(Longévité)---> fonctionne pas
    shapiro.test(Longévité$residua ls)
    Shapiro-Wilk normality test

    Longévité$residuals
    W = 0.969, p-value = 0.001510

    Donc là 1er problème! p<<0,05!

    Au cas où j'essaie mon test de Levene:
    levene.test(long,interaction(t raitement*espèce))

    J'ai tenté tout betement: levene.test(Longévité) et ça me donne quelque chose mais je ne sais pas si ça prend en compte l'interaction espèce/traitement pour vérifier l'homoscédasticité...

    Avec levene.test(Longévité$residual s)
    là ça marche pas du tout! Comment je peux faire pour que mon test de Levene fonctionne? Comment je peux lier le traitement et l'espèce pour mon test?

    Voilà j'espère que vous pourrez m'aider
    Merci d'avance

Discussions similaires

  1. Statistique sur R
    Par invite8b118ae7 dans le forum Mathématiques du supérieur
    Réponses: 2
    Dernier message: 29/05/2009, 13h40
  2. question sur le nombre minimum d'observations à avoir pour pouvoir faire une ANOVA?
    Par invite7323c4b6 dans le forum Mathématiques du supérieur
    Réponses: 2
    Dernier message: 22/05/2009, 10h11
  3. ANOVA sur SPSS
    Par invited7ea5151 dans le forum Mathématiques du supérieur
    Réponses: 4
    Dernier message: 11/06/2008, 11h33
  4. SPSS - anova 2x2
    Par invite3c9f3cce dans le forum Mathématiques du supérieur
    Réponses: 0
    Dernier message: 12/01/2008, 00h25
  5. [statistiques] p-value ANOVA
    Par inviteb1ba0069 dans le forum Biologie
    Réponses: 1
    Dernier message: 27/06/2007, 14h02