Go Back   SEQanswers > Applications Forums > RNA Sequencing

Similar Threads
Thread Thread Starter Forum Replies Last Post
EdgeR or GLM model design question younko Bioinformatics 10 08-31-2014 11:22 PM
edgeR with uneven groups sindrle Bioinformatics 6 10-22-2013 03:08 PM
Experimental design with edgeR and DESeq packages yvan.wenger Bioinformatics 0 11-16-2012 02:53 AM
seeking pool-seq experimental design advice bluesquid Genomic Resequencing 0 08-07-2012 04:15 AM
EdgeR design-matrix design extended.wobble RNA Sequencing 3 07-11-2011 06:58 AM

Thread Tools
Old 08-17-2015, 01:00 PM   #1
Junior Member
Location: Penn

Join Date: Jun 2013
Posts: 4
Default edgeR design: pool groups


I am using edgeR post kallisto to look for DGE.

I have RNAseq samples from 4 different tissue types and my design currently looks like this:

> head(design)
(Intercept)$samples[, 1]CA Endothelium
1 1 1
2 1 1
3 1 1
4 1 1
5 1 1
6 1 1$samples[, 1]CA Smooth Muscle$samples[, 1]Coronary Artery
1 0 0
2 0 0
3 0 0
4 0 0
5 0 0
6 0 0

[1] 0 1 1 1
attr(,"contrasts")$`$samples[, 1]`
[1] "contr.treatment"

So my four tissue types are Aorta (intercept), CA Endothelium, CA Smooth Muscle and Coronary Artery.

As I understand I can compare gene expression between any tissue using:
lrt.ECvsAo<-glmLRT(, contrast=c(1,-1,0,0))

If I use this contrast
lrt.ECvsAll<-glmLRT(, contrast=c(1,-1,1,1))
am I asking for DGE between the second column (CA Endothelium) and EITHER other tissue?

I'd like to find genes that are differentially expressed between CA Endothelium and ALL Aorta, CA Smooth Muscle and Coronary Artery.

Is this doable?

Many thanks
sylvian is offline   Reply With Quote
Old 08-17-2015, 04:17 PM   #2
Location: Berkeley, cA

Join Date: Feb 2010
Posts: 40

If you're quantifying with kallisto you may want to take a look at sleuth
lpachter is offline   Reply With Quote
Old 11-11-2015, 07:39 AM   #3
Junior Member
Location: United States-NJ

Join Date: Jul 2015
Posts: 6

Hi all,
I am using edgeR flow to get gene lists of multifactorial experiment. I have 6 conditions with 3 samples per each. The problem is some of the gene lists I got end at FDR much less than what I entered and so I did not get the whole gene list. For example, the upregulated genes were completely missing in some lists. I tried increasing the number of genes the program gives using the
< gene_table<-topTags(lrt, n=40000)
but did not help. Then I tried very restreignent FDR as follow
< DE<-tp$table[tp$table$FDR<0.000001,]$GeneID
But again did not help. The last gene in one of my lists was with
Symbol Chr logFC logCPM LR PValue FDR
Hpx 7 -7.186931303 13.83923493 57.6010426 3.21E-014 3.23E-014

Any help would be appreciated!
TRA is offline   Reply With Quote


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:56 PM.

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