[BioC] merging columns of large data sets
Martin Morgan
mtmorgan at fhcrc.org
Sat Sep 26 22:11:32 CEST 2009
Amit Kumar wrote:
> Hi! All
>
> I am trying to merge very large data sets. Here fullset1 contains 13
> ExpressionSets. Each data set has columns containing Cy5 and Cy3
> intensities which I need to merge.
>
> Here I am trying to merge columns (containing Cy5 intensity)
> 2,6,10,14....till end for all 13
> data sets in fullset1. But I am only getting 2nd column from each of
> the 13 data sets.
>
> Rchan1 = sapply(1:length(fullset1), function(i) exprs(fullset1[[i]])[ ,2])
> Rch1 = as.list(Rchan1)
> red1 = do.call(cbind,Rch1)
looks like you want columns in
seq(2, ncol(exprs), 4)
so, assuming fullset1 is a list,
red1 = lapply(fullset1,
function(eset) exprs(est)[,seq(2, ncol(eset), 4)])
red = do.call(cbind, red1)
No need to cross-post (to the R and Bioc mailing lists)
Martin
>
> Likewise, here I need to merge columns (containing Cy3 intensity)
> 3,7,11,15...till end for all 13
> data sets in fullset1. And I am getting 3rd column from each of the 13
> data sets.
>
> Gchan1 = sapply(1:length(fullset1), function(i) exprs(fullset1[[i]])[ ,3])
> Gch1 = as.list(Gchan1)
> green1 = do.call(cbind,Gch1)
>
> I am stuck! Please Help! If its not possible this way please suggest
> me alternative way to do it.
>
> Cheers!
> Amit
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
--
Martin Morgan
Computational Biology / Fred Hutchinson Cancer Research Center
1100 Fairview Ave. N.
PO Box 19024 Seattle, WA 98109
Location: Arnold Building M1 B861
Phone: (206) 667-2793
More information about the Bioconductor
mailing list