SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Exome data: how to filter by cases and controls? Sini Bioinformatics 2 02-17-2013 10:47 PM
Tag reads from cases and controls meher Bioinformatics 0 10-25-2012 11:48 AM
edgeR Puva Bioinformatics 2 05-19-2011 10:04 AM

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

Join Date: Apr 2013
Posts: 35
Default edgeR for paired cases

Hi all,

I am using EdgeR for DE analysis. I would appreciate if somebody could help me check my experiment design.

I have 6 cows (2,63,66,67,73,74) for a treatment, each with 3 time points(m14,p7,p28). I want to compare between p28 vs m14, p7 vs m14 and p28 vs p7.


cow <- factor(c(2,2,2,63,63,63,66,66,66,67,67,67,73,73,73,74,74,74))
timep <- factor(c("m14","p7","p28","m14","p7","p28","m14","p7","p28","m14","p7","p28","m14","p7","p28","m14","p7","p28"))
data.frame(Sample=colnames(y),cow,timep)
design<-model.matrix(~cow+timep)

This follows the EdgerR user guide (31 March 2013 revised), case study 4.4 on page 53.

Thanks.
lran2008 is offline   Reply With Quote
Old 08-27-2013, 02:39 PM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

The short answer is that you want to simply use constrasts, such as accepted by glmLRT. Here's one thread on the bioconductor email list (the edgeR authors are pretty active there) that gives some suggestions to get you started.
dpryan is offline   Reply With Quote
Reply

Tags
edger, paired

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 08:35 AM.


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