loi de fisher sous R

[invite003c405b]

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.
-----

[invite003c405b]

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}
}
}

[invitec5eb4b89]

Mais quelle est la question ???

[kinette]

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

[A voir en vidéo sur Futura]

[invite003c405b]

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

[invite003c405b]

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)
}

[kinette]

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

[invite003c405b]

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)
}