options on outputing gene expression based on gene level instead of probe set level for justRMA
3
0
Entering edit mode
Jack Luo ▴ 440
@jack-luo-4241
Last seen 9.6 years ago
Hi, I am trying to RMA normalize a bunch of cel files and output the gene expression matrix to gene level (instead of probe set level), is there any option in justRMA that can do this (or any other function that can do this?) I searched the options on justRMA, does not seem to have this option. Thanks, -Jack [[alternative HTML version deleted]]
probe probe • 1.1k views
ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 1 hour ago
United States
Hi Jack, On 3/6/2013 11:48 AM, Jack Luo wrote: > Hi, > > I am trying to RMA normalize a bunch of cel files and output the gene A bunch of cel files of what Array type? > expression matrix to gene level (instead of probe set level), is there any > option in justRMA that can do this (or any other function that can do > this?) I searched the options on justRMA, does not seem to have this > option. > > Thanks, > > -Jack > > [[alternative HTML version deleted]] > > _______________________________________________ > Bioconductor mailing list > Bioconductor at r-project.org > 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 University of Washington Environmental and Occupational Health Sciences 4225 Roosevelt Way NE, # 100 Seattle WA 98105-6099
ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 1 hour ago
United States
Please don't take things off list. On 3/6/2013 7:59 PM, Jack Luo wrote: > Gene 1.0 ST. At this point I guess I need to know what you mean by gene level. The affy package by definition only summarizes at the transcript level. But that doesn't mean that you only have one probeset per gene, as there may still be duplicates (splice variants or whatnot). If you just want a single expression value per gene, you can either use the genefilter package to subset by Entrez ID, or you could alternatively use one of the MBNI re-mapped cdf packages. Best, Jim > > On Wed, Mar 6, 2013 at 2:53 PM, James W. MacDonald <jmacdon at="" uw.edu=""> <mailto:jmacdon at="" uw.edu="">> wrote: > > Hi Jack, > > > On 3/6/2013 11:48 AM, Jack Luo wrote: > > Hi, > > I am trying to RMA normalize a bunch of cel files and output > the gene > > > A bunch of cel files of what Array type? > > > expression matrix to gene level (instead of probe set level), > is there any > option in justRMA that can do this (or any other function that > can do > this?) I searched the options on justRMA, does not seem to > have this > option. > > Thanks, > > -Jack > > [[alternative HTML version deleted]] > > _______________________________________________ > Bioconductor mailing list > Bioconductor at r-project.org <mailto:bioconductor at="" r-project.org=""> > 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 > University of Washington > Environmental and Occupational Health Sciences > 4225 Roosevelt Way NE, # 100 > Seattle WA 98105-6099 > > -- James W. MacDonald, M.S. Biostatistician University of Washington Environmental and Occupational Health Sciences 4225 Roosevelt Way NE, # 100 Seattle WA 98105-6099
ADD COMMENT
0
Entering edit mode
Thanks, Jim. That's very helpful. -Jack On Thu, Mar 7, 2013 at 8:57 AM, James W. MacDonald <jmacdon@uw.edu> wrote: > Please don't take things off list. > > On 3/6/2013 7:59 PM, Jack Luo wrote: > >> Gene 1.0 ST. >> > > At this point I guess I need to know what you mean by gene level. The affy > package by definition only summarizes at the transcript level. But that > doesn't mean that you only have one probeset per gene, as there may still > be duplicates (splice variants or whatnot). If you just want a single > expression value per gene, you can either use the genefilter package to > subset by Entrez ID, or you could alternatively use one of the MBNI > re-mapped cdf packages. > > Best, > > Jim > > > >> On Wed, Mar 6, 2013 at 2:53 PM, James W. MacDonald <jmacdon@uw.edu<mailto:>> jmacdon@uw.edu>> wrote: >> >> Hi Jack, >> >> >> On 3/6/2013 11:48 AM, Jack Luo wrote: >> >> Hi, >> >> I am trying to RMA normalize a bunch of cel files and output >> the gene >> >> >> A bunch of cel files of what Array type? >> >> >> expression matrix to gene level (instead of probe set level), >> is there any >> option in justRMA that can do this (or any other function that >> can do >> this?) I searched the options on justRMA, does not seem to >> have this >> option. >> >> Thanks, >> >> -Jack >> >> [[alternative HTML version deleted]] >> >> ______________________________**_________________ >> Bioconductor mailing list >> Bioconductor@r-project.org <mailto:bioconductor@r-**project.org<bioconductor@r-project.org> >> > >> >> https://stat.ethz.ch/mailman/**listinfo/bioconductor<https: stat.ethz.ch="" mailman="" listinfo="" bioconductor=""> >> Search the archives: >> http://news.gmane.org/gmane.**science.biology.informatics.** >> conductor<http: news.gmane.org="" gmane.science.biology.informatics.c="" onductor=""> >> >> >> -- James W. MacDonald, M.S. >> Biostatistician >> University of Washington >> Environmental and Occupational Health Sciences >> 4225 Roosevelt Way NE, # 100 >> Seattle WA 98105-6099 >> >> >> > -- > James W. MacDonald, M.S. > Biostatistician > University of Washington > Environmental and Occupational Health Sciences > 4225 Roosevelt Way NE, # 100 > Seattle WA 98105-6099 > > [[alternative HTML version deleted]]
ADD REPLY
0
Entering edit mode
@peter-langfelder-4469
Last seen 23 days ago
United States
On Wed, Mar 6, 2013 at 8:48 AM, Jack Luo <jluo.rhelp at="" gmail.com=""> wrote: > Hi, > > I am trying to RMA normalize a bunch of cel files and output the gene > expression matrix to gene level (instead of probe set level), is there any > option in justRMA that can do this (or any other function that can do > this?) I searched the options on justRMA, does not seem to have this > option. You can output probe-set level data, then use the function collapseRows in the WGCNA package (CRAN) to summarize probesets for each gene to a single value. Best, Peter
ADD COMMENT
0
Entering edit mode

I am new to bioinformatics, and I am faced with this exact problem of summarizing probe-set level data to a single value for each gene.

Your WGCNA suggestion sounds simple enough, but I wanted to check if there is a newer suggested approach since your answer is over 6 years old?

James MacDonald suggests using either genefilter or MBNI re-mapped cdf packages, but these solutions are less clear to me. Do you think these alternatives have potential advantages?

Thank you, James

ADD REPLY

Login before adding your answer.

Traffic: 685 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6