edgeR - expressen between genes
1
0
Entering edit mode
@sridhara-gupta-kunjeti-4449
Last seen 9.7 years ago
United States
Hi All, I am using RNA-seq to study the expression levels of gene globally and as well as group of genes. Plotsmear have given me a meaningful information. I believe the estimated norm.factor is directly included in the model to test DE, but the values themselves are not modified. My question is if I want to check compare the expression between genes (within a group) what values do I need to consider for this analysis. Any help or suggestions will be highly appreciated. Thanks in advance! Sridhara -- Sridhara G Kunjeti PhD Candidate University of Delaware Department of Plant and Soil Science email- sridhara@udel.edu Ph: 832-566-0011 [[alternative HTML version deleted]]
• 1.2k views
ADD COMMENT
0
Entering edit mode
@sridhara-gupta-kunjeti-4449
Last seen 9.7 years ago
United States
Opps I am sorry there is a typo in the subject - expression Thanks, Sridhara On Fri, Jan 28, 2011 at 6:37 PM, Sridhara Gupta Kunjeti <sridhara@udel.edu>wrote: > Hi All, > I am using RNA-seq to study the expression levels of gene globally and as > well as group of genes. Plotsmear have given me a meaningful information. I > believe the estimated norm.factor is directly included in the model to test > DE, but the values themselves are not modified. My question is if I want to > check compare the expression between genes (within a group) what values do I > need to consider for this analysis. Any help or suggestions will be highly > appreciated. > > Thanks in advance! > Sridhara > > > > -- > Sridhara G Kunjeti > PhD Candidate > University of Delaware > Department of Plant and Soil Science > email- sridhara@udel.edu > Ph: 832-566-0011 > -- Sridhara G Kunjeti PhD Candidate University of Delaware Department of Plant and Soil Science email- sridhara@udel.edu Ph: 832-566-0011 [[alternative HTML version deleted]]
ADD COMMENT
0
Entering edit mode
Hi Sridhara. >> I believe the estimated norm.factor is directly included in the model to test >> DE, but the values themselves are not modified. This is correct: the normalization factors are built in as "offsets" in the statistical model, so the count data does not get transformed. >> My question is if I want to >> check compare the expression between genes (within a group) what values do I >> need to consider for this analysis. Comparing expression b/w genes (as opposed to across conditions for a given gene) is not a standard use case within edgeR (yet) and there are a number of issues that should be dealt with if you want to do this well. I'm sure others on the list have more experience and can comment more on this, but it seems you need to get a very good handle on gene length (can you always trust the annotation?) and other biases that (e.g. GC content, positional biases) that may introduce differences in read counts that are different from gene to gene. That said, edgeR has fairly general routines to analyze count data, if you can sort out these other factors. Cheers, Mark On 2011-01-29, at 10:37 AM, Sridhara Gupta Kunjeti wrote: > Opps I am sorry there is a typo in the subject - expression > Thanks, > Sridhara > > > > On Fri, Jan 28, 2011 at 6:37 PM, Sridhara Gupta Kunjeti > <sridhara at="" udel.edu="">wrote: > >> Hi All, >> I am using RNA-seq to study the expression levels of gene globally and as >> well as group of genes. Plotsmear have given me a meaningful information. I >> believe the estimated norm.factor is directly included in the model to test >> DE, but the values themselves are not modified. My question is if I want to >> check compare the expression between genes (within a group) what values do I >> need to consider for this analysis. Any help or suggestions will be highly >> appreciated. >> >> Thanks in advance! >> Sridhara >> >> >> >> -- >> Sridhara G Kunjeti >> PhD Candidate >> University of Delaware >> Department of Plant and Soil Science >> email- sridhara at udel.edu >> Ph: 832-566-0011 >> > > > > -- > Sridhara G Kunjeti > PhD Candidate > University of Delaware > Department of Plant and Soil Science > email- sridhara at udel.edu > Ph: 832-566-0011 > > [[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 ------------------------------ Mark Robinson, PhD (Melb) Epigenetics Laboratory, Garvan Bioinformatics Division, WEHI e: mrobinson at wehi.edu.au e: m.robinson at garvan.org.au p: +61 (0)3 9345 2628 f: +61 (0)3 9347 0852 ------------------------------ ______________________________________________________________________ The information in this email is confidential and intend...{{dropped:6}}
ADD REPLY

Login before adding your answer.

Traffic: 587 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