hola!
estoy trabajando con la funcion cor.test, y quiero calcular un coeficiente de correlacion para muchos archivos.
En la funcion creo dos directorios, y en cada uno de ellos tengo 200 archivos, a partir de los cuales debo calcular el coeficiente.
Mi problema es que si leo los archivos (dentro de un for)y despues intento calcularlo, el resultado que obtengo solo me lo calcula para el primer fichero.
Por ejemplo:
(estas sentencias estan dentro de una funcion, q recorre con un for los archivos en nam1 y nam2, que son 200 archivos con columnas de 30 datos cada uno de ellos).
write.table(x1$V2,file=nam1,append=F,quote=T,sep="",row.names=F,col.names=F) write.table(y1$V2,file=nam2,append=F,quote=T,sep="",row.names=F,col.names=F)
a1=read.table(nam1,sep="")
a2=read.table(nam2,sep="")
m1<-cor.test(a1$V1,a2$V1,method="spearman")
pval=m1$estimate
p=rbind(m1$estimate)
print(p)
el resultado que obtengo es:
rho
[1,] 0.2282342
lo que inteno es obtener el rho para los 200 archivos.
Espero haberme explicado con claridad, si alguien sabe como puedo conseguirlo?
muchas gracias
La instrucción for entiende que el cuerpo del bucle es la instrucción siguiente, si escribes:
for (i in 1:10)
print(i);
print(i^2);
Únicamente el primer print se ejecuta dentro del for, si deseas ejecutar un grupo de instrucciones debes encerrarlas entre llaves. Para el caso anterior:
for (i in 1:10)
{
print(i);
print(i^2);
}
Saludos.
for (i in 1:10)
print(i);
print(i^2);
Únicamente el primer print se ejecuta dentro del for, si deseas ejecutar un grupo de instrucciones debes encerrarlas entre llaves. Para el caso anterior:
for (i in 1:10)
{
print(i);
print(i^2);
}
Saludos.