SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
comparing results by cuffdiff, edgeR, DESeq PFS Bioinformatics 5 03-12-2014 03:01 AM
cummeRbund andplotting for selected treatments among multiple treatments ataheri RNA Sequencing 1 04-10-2013 06:41 PM
A few questions regarding RNA-seq results Kocur RNA Sequencing 0 03-26-2013 12:56 AM
Several questions about our RNA-seq results sazz RNA Sequencing 3 01-29-2013 12:25 PM
edgeR glmFit results capricy Bioinformatics 1 11-19-2012 02:43 PM

Reply
 
Thread Tools
Old 08-30-2013, 07:47 AM   #1
lran2008
Member
 
Location: quebec

Join Date: Apr 2013
Posts: 35
Default EdgeR DE results of two treatments questions

Hi all,

I am using edger to find DE genes among two time points (P7, P28) compared with control. In my experiment, 6 cows were subjected to treatA and 5 cows were subjected to treatB. I have the miRNA seq data for control, P7 and P28 of each cow.

Now My Question is that:
The first time I got the DE genes of treatA (P28 vs control, P7 vs control) by just importing the seq data of the 6 cows in treatA to EdgeR. (I perform the same work for treat B separately.)

The second time I imported all the seq data of the two treatments together into Edger and I still could get the DE genes for treat A and treat B. (In this way, I can compare the effects of the two treatments)

HOWEVER, the results from the two approaches for the same treatment are different. I understand that the normalization factor changed when the data of the two treatments were put together.

BUT which one is more reliable? Should I analyze the DE genes of two treatments separately or together?
lran2008 is offline   Reply With Quote
Old 08-30-2013, 11:43 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

Well, the proper comparison would be to compare one of the treatments at the various timepoints to an appropriate sham treatment at the same time points, but that train has pretty much left the station (i.e., it's too late to complain about that since you already did the experiment).

In real biological situations, there are typically changes overtime both due to a treatment and in the absence of one (biology is messy). In the first case, you're picking up changes that simply happen randomly over time and those that happen due to the treatment. In the second comparison, you somewhat account for random changes over time, but not completely because treatment B can cause its own host of time-dependent changes. Lacking any information about the actual treatments, I would predict that the second comparison is more reliable. In the future, though, you should talk to someone who's done animal experiments like this prior to doing them so you can setup the proper controls...your life will then be much easier.
dpryan is offline   Reply With Quote
Old 09-03-2013, 05:08 PM   #3
mikep
Member
 
Location: Singapore

Join Date: Feb 2011
Posts: 45
Default

if you use the GLM functions in edgeR then you can set your experiment up as the 2-way repeated measures anova that it actually is.

Section 3.3 of the edgeR users guide covers this.
mikep is offline   Reply With Quote
Old 09-05-2013, 04:38 AM   #4
lran2008
Member
 
Location: quebec

Join Date: Apr 2013
Posts: 35
Default

Quote:
Originally Posted by mikep View Post
if you use the GLM functions in edgeR then you can set your experiment up as the 2-way repeated measures anova that it actually is.

Section 3.3 of the edgeR users guide covers this.
Thanks. This is exactly what I need.
lran2008 is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 09:25 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO