Hi. I've used Affy 133 Plus 2 Arrays on a series of samples and normalised the data using rma. Can anyone tell me how to run eBayes to compare two subsets of samples (with and without monosomy of chromosome 10) using only probesets mapped to chromosome 10? Thanks. Dr. Martin G. McCabe Cancer Research UK Clinical Research Training Fellow Cambridge University Department of Pathology Addenbrooke's Hospital Cambridge

Hi Martin, Martin McCabe wrote: > Hi. I've used Affy 133 Plus 2 Arrays on a series of samples and > normalised the data using rma. Can anyone tell me how to run eBayes to > compare two subsets of samples (with and without monosomy of chromosome > 10) using only probesets mapped to chromosome 10? I will assume that 'run eBayes' means you want to fit a linear model to your data using the limma package. The eBayes step assumes that you are using all the genes on your chip, so you should fit the model up to the eBayes step and then filter on chromosome. To filter on chromosome you can do something like [Some code to get to this point] fit2 <- eBayes(fit) library(hgu133plus2.db) chr10 <- mget("10", revmap(hgu133plus2CHR))[[1]] just.chr10 <- fit2[chr10,] Best, Jim > > Thanks. > > Dr. Martin G. McCabe > > Cancer Research UK Clinical Research Training Fellow > Cambridge University Department of Pathology > Addenbrooke's Hospital > Cambridge > > _______________________________________________ > 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 -- James W. MacDonald, M.S. Biostatistician Hildebrandt Lab 8220D MSRB III 1150 W. Medical Center Drive Ann Arbor MI 48109-0646 734-936-8662