Hi,

I was processing some data and wanted to consolidate expression values by gene, instead of by probe id as given the expression dataset. The hgu133a chips has 22283 probe ids corresponding to 12704 genes. Is there any function in the package that will automatically consolidate the expression value genewise? In other words (from the code given below), how can I get a 12704 X 38 matrix from the probe level 22283 X 38 matrix?
> dim(emat)  # expression data for 22283 probes and 38 samples
[1] 22283    38

> probes <- rownames(emat)
> gs <- unlist(mget(probes,hgu133aSYMBOL))
> length(probes)
[1] 22283

> length(unique(gs))
[1] 12704
many thanks!


      
	[[alternative HTML version deleted]]

