[R-es] Como calcular la media de una gran cantidad de flow.frame en un ciclo "for" ...
Oscar Perpiñan Lamigueiro
oscar.perpinan en upm.es
Vie Dic 23 11:50:55 CET 2011
Hola Eric,
Dado que los flowFrame son clases definidas en el paquete de
Bioconductor flowCore, es probable que obtengas ayuda más fácilmente si
escribes a la lista de correo de Bioconductor:
http://www.bioconductor.org/help/mailing-list/#bioconductor.
Saludos.
Oscar.
El Thu, 22
Dec 2011 13:57:04 -0300 Eric <ericconchamunoz en gmail.com> escribió:
> Queridos amigos tengo el siguiente problema:
>
> Tengo alrededor de 100 archivos del tipo flowFrame (datos de
> citometria) que he importado a R. Estos archivos tienen 6 columnas y
> 10.000 (diez mil) filas. Tengo que calcular el promedio de la quinta
> columna y graficarlo contra el promedio de la segunda columna. La
> segunda columna tiene el mismo valor para las 10.000 filas en cada
> archivo asi es que sabemos el promedio de antemano, pero quiero
> incluirlo en el calculo de todos modos. Para calcular las medias he
> creado el siguiente nano-algoritmo que incluye un ciclo "for" :
>
> lista <- ls(pattern="dat12") # para usar solo los archivos que son
> del mes 12 que son alrededor de 100
> medias <- c(1:length(lista))
> for (i in length(lista)){
> medias[i] <- mean(exprs(lista[i][,5]))
> }
>
> pero obtengo el siguiente error:
>
> Error in exprs(lista[i][, 5]) :
> error in evaluating the argument 'object' in selecting a method for
> function 'exprs': Error in lista[i][, 5] : incorrect number of
> dimensions
>
> parece que usar *lista[i][,5]* como argumento de la funcion exprs es
> un error , porque lista[i] no se comporta como un nombre de archivo
> como yo crei que haria y esto es lo que no se como hacer, es decir,
> como se hace para poder utilizar los nombres contenidos en el vector
> "lista" y construir el argumento de la funcion "exprs"
> correctamente ... o como se puede hacer de otro modo esto mismo??
>
> me explico bien ??
>
> entre parentesis es la primera vez que trato de hacer un ciclo en R y
> quiza hasta eso este mal ...
>
> Saludos y muchas gracias por su ayuda,
>
> Eric.
>
>
>
>
>
>
>
Más información sobre la lista de distribución R-help-es