Question: Normalized transcript/exon counts for network analysis
gravatar for maya.kappil
4 months ago by
maya.kappil10 wrote:


I would like to generate a cosplicing network similar to how it's described here ( My understanding is that I would need normalized exon/transcript counts prior to generating the network. For DEXSeq, I was planning to input the normalized exon counts from the DEXSEqDataSet generated as follows:

dxd <- DEXSeqDataSet(,, design=~sample + exon + group:exon, featureID=counts(d)$featureid, groupID=counts(d)$geneid)

system.time({ dxd <- estimateSizeFactors(dxd) dxd <- estimateDispersions(dxd, quiet=TRUE) dxd <- nbinomLRT(dxd, reduced=~sample + exon) })


I was wondering whether a similar option is available to output normalized transcript levels from DRIMSeq?



dexseq sva wgcna drimseq • 158 views
ADD COMMENTlink modified 4 months ago by Gosia Nowicka40 • written 4 months ago by maya.kappil10
Answer: Normalized transcript/exon counts for network analysis
gravatar for Gosia Nowicka
4 months ago by
Gosia Nowicka40 wrote:

Hello Maya,

no, there is no way to get normalized data from the DRIMSeq workflow. This is because DRIMSeq does not need to use normalized data when comparing proportions.

What you could do is to input your transcript counts into DEXSeq in the same way you did with exon counts.

ADD COMMENTlink written 4 months ago by Gosia Nowicka40

Ah, ok - that makes sense! Thanks!!

ADD REPLYlink written 3 months ago by maya.kappil10
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 264 users visited in the last hour