Discussion :

[karaudrey88]

Bonjour
J'aimerais comparer les colonnes de deux matrices :

matrice1
matrice2

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

var.test(matrice1$V1,matrice2$V1)

Mais mes deux matrices étant très grandes, j'aimerais automatiser le processus.
J'ai déjà essayé le script suivant mais cela ne fonctionne pas :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
> m<-read.table ("essai.txt")
> p<-read.table ("essai3.txt")
> for (i in 0:1000){
+ var.test(m$Vi,p$Vi)}
Error in var.test.default(m$Vi, p$Vi) : not enough 'x' observations

Est-ce que quelqu'un saurait comment faire?

Merci.

[hoccha]

Salut,

Je ne sais pas si ceci vous convient:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
m <- matrix(runif(1000),ncol=10)
p <- matrix(runif(1000),ncol=10)
 
res <- numeric(ncol(m))
for(i in 1:ncol(m))
	res[i] <- var.test(m[,i]  , p[,i])$p.value  # ici je prends les p-value
 
print(res)

Bon week-end

[karaudrey88]

Votre réponse me convient tout à fait.
Merci de votre aide.
Bon week-end.