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

loi de fisher sous R




  1. #1
    goz

    loi de fisher sous R

    Bonjour,
    je souhaite récuperer la valeur d'une loi de fisher sous R, voici la commade que j'ecris :
    f=qf(0.95,length(type1[3])-1,length(type2[3])-1)
    (type1 et type2 sont des matrices)
    et voici la réponse de R :
    production de NaN in: qf(p, df1, df2, log)

    j'ai essayé en remplacant les longueurs des vecteur par des nombres ca marche.

    -----


  2. Publicité
  3. #2
    goz

    Re : loi de fisher sous R

    au cas ou, voici le code de la fonction ...

    moy<-function(a,b){
    for (i in 3:dim(a)[2]) {
    x=var(a[i])
    y=var(a[i])
    t=length(a[i])-1
    y=length(b[i])-1
    f=qf(0.95,t,y)
    f
    if (x>y) { fc=x/y} else {fc=y/x}
    if (fc<f) { result[i]<-t.test(a[i],b[i],var.equal=TRUE)$p.value} else { result[i]<-t.test(a[i],b[i],var.equal=FALSE)$p.value}
    }
    }

  4. #3
    HigginsVincent

    Re : loi de fisher sous R

    Mais quelle est la question ???


  5. #4
    kinette

    Re : loi de fisher sous R

    Bonjour,
    Je n'ai pas non plus compris la question... de plus beaucoup de fonctions de base sont déjà présentes par défaut dans R...

    K
    Nomina si nescis, perit et cognito rerum.

  6. #5
    goz

    Re : loi de fisher sous R

    c'est ici qu'il y a un problème :
    production de NaN in: qf(p, df1, df2, log)
    c'est ici le problème, du coup cette ligne ne renvoie aucune valeure, ce qui pose problème pour la comparaison
    if (fc<f)
    le but est donc de savoir comment faire pour que le qf(...) renvoie une valeure
    en tout cas merci de votre interêt

  7. A voir en vidéo sur Futura
  8. #6
    goz

    Re : loi de fisher sous R

    résolution du problème, en fait l'utilisation d'une autre version de R m'a fait changer qques trucs de syntaxe et m'a permis de débloquer la situation. En fait le length de a[i] renvoie 1, il fallait faire length a[,i]
    voici le code final, si jamais cela peut servir a qqun

    function(a,b){
    for (i in 3:dim(a)[2]) {
    x=var(a[,i])
    y=var(b[,i])
    t=length(a[,i])-1
    y=length(b[,i])-1
    f=qf(0.95,t,y)
    f
    if (x>y) { fc=x/y} else {fc=y/x}
    if (fc<f) { result[i]<-t.test(a[,i],b[,i],var.equal=TRUE)$p.value} else { result[i]<-t.test(a[,i],b[,i],var.equal=FALSE)$p.value}
    }
    return (result)
    }

  9. #7
    kinette

    Re : loi de fisher sous R

    Bonjour,
    Explique clairement:
    - ce que tu as au départ: valeurs expérimentales, organisées comment...
    - ce que tu souhaites obtenir

    C'est très difficile à partir d'un code sans plus d'explication de voir ce que tu souhaites faire...

    K
    Nomina si nescis, perit et cognito rerum.

  10. Publicité
  11. #8
    goz

    Re : loi de fisher sous R

    a et b sont des matrices de réels (valeurs expérimentales, spécimens en lignes, parramètres en colonnes), il s'agit en fait d'un test de comparaison de moyenne entre les valeurs de a et b (les colonnes correspondent au mêmes paramètes dans ces deux matrices)
    pour cela, je fait d'abord un test sur les variances (j'avais des problèmes avec le var.test, c'est pour cela que je l'ai écrit) puis ensuite le test sur les moyennes selon les résultats du test de variances.

    sinon, je viens de penser à un truc, monter que les variances ne sont pas significativement identiques revient il à montrer qu'elles not significativement différentes ??
    voici les commentaires en plus

    moy<-function(a,b){
    for (i in 2:dim(a)[2]) { // la première colonne correspond au nom des echantillons
    x=var(a[,i]) variance sur les données d'un paramètre
    y=var(b[,i])
    t=length(a[,i])-1
    y=length(b[,i])-1
    f=qf(0.95,t,y) calcul de variance théorique admissible pour le paramètre
    f
    if (x>y) { fc=x/y} else {fc=y/x} fc: variance calculée
    if (fc<f) { result[i]<-t.test(a[,i],b[,i],var.equal=TRUE)$p.value} else { result[i]<-t.test(a[,i],b[,i],var.equal=FALSE)$p.value}
    } // si les variances sont identiques student, sinon welsh, et je stoque les résultats dans le vecteur result
    return (result)
    }

Sur le même thème :

Discussions similaires

  1. Karl fisher
    Par luce dans le forum Chimie
    Réponses: 5
    Dernier message: 25/02/2009, 18h36
  2. Fisher
    Par Salfania dans le forum Chimie
    Réponses: 4
    Dernier message: 20/10/2007, 12h53
  3. Matrices de Fisher
    Par physeb dans le forum Mathématiques du supérieur
    Réponses: 4
    Dernier message: 12/04/2007, 17h29
  4. synthèse de kiliani fisher
    Par casar900 dans le forum Chimie
    Réponses: 5
    Dernier message: 03/03/2005, 21h53
  5. A propos de Karl-Fisher...
    Par grego7 dans le forum Chimie
    Réponses: 6
    Dernier message: 17/02/2005, 12h39