Go Back   SEQanswers > Applications Forums > RNA Sequencing

Similar Threads
Thread Thread Starter Forum Replies Last Post
DESeq: defining conditions mouchkam Bioinformatics 2 04-10-2012 05:45 PM
Is more than two conditions possible in DESEQ? greener RNA Sequencing 5 05-09-2011 04:10 PM
RNASeq experiments with 2 conditions and more than 3 replicates Jouneau Luc RNA Sequencing 1 04-05-2011 04:00 AM
PCR conditions for PE adapters Luke Illumina/Solexa 3 03-02-2011 09:48 AM

Thread Tools
Old 05-19-2014, 11:15 AM   #21
Junior Member
Location: Maine, USA

Join Date: Jul 2013
Posts: 7

Continuing on my post above, I am now trying to generate the fold change between treatments corrected by the effect of time. Using the nomenclature from my post above this is the fold change I would like to generate : (Drug.2h/Drug.0h)/(Placebo.2h/Placebo.0h).

If you will remember, DEU was tested for using the following formulae :

formulaFullModel = ˜ sample + exon + treatment:exon + time:exon + treatment:time:exon
formulaReducedModel = ˜ sample + exon + treatment:exon + time:exon

I am now trying to generate the fold changes that occur using the testForDEU() function. However, this command only takes one variable and my model uses two.

In my scenario (see post above for full explanation) I am thinking that the model that needs to be fit in testForDEU() is '~ sample + time*treatment* exon. With the denominator being 'Placebo'

My questions are: is this the right model to use to obtain the desired fold changes?
can this model be used in testForDEU() with some small modifications?
GreboGuru is offline   Reply With Quote
Old 03-09-2015, 10:25 AM   #22
Junior Member
Location: USA, virginia

Join Date: Feb 2014
Posts: 1

I'm using DEXSeq for the analysis of 16 human tissues. I'm trying to identify exons that are differentially expressed in one tissue against all other tissues. I ran the whole pipeline of DEXSeq on all the sample for all tissues. What I understand is when I have exon with differential usage in the results table, that means it is differentially expressed in at least one tissue but I can't determine which one. I also can't determine if it is differentially expressed on more than one tissue, is that right? is there anyway to know such information?
emanbadr is offline   Reply With Quote
Old 04-25-2019, 09:00 PM   #23
Junior Member
Location: india

Join Date: Apr 2019
Posts: 1


I want to conduct DEU analysis for three different tissues of different conditions namely F, LR and LR.
Here I used the code in DEXseq package of R to make sample table and design. The input data is counts.txt file from FeatureCounts Here I used the code

sampleTable = data.frame(
  row.names = c("F1","F2","LR1","LR2","LV1","LV2"),
  condition = factor(rep(c("F","LR","LV"),each=2)),
  libType = c("paired-end","paired-end","paired-end", "paired-end","paired-end", "paired-end"))


dxd <- DEXSeqDataSet(cc,design=~sample + exon + condition:exon, featureID=as.factor(exon_ids),
but I only get 'exon fold change' in between F and LR, F and LV but I could not get the exon fold change between LV and LR. kindly suggest the solution to get comparision in between all three conditions.

Thank you,

aniket 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 07:14 PM.

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