Question: Comparing Gene Expression From Strand-Specific Library And Unstranded Library In Rnaseq
gravatar for camelbbs
6.4 years ago by
camelbbs660 wrote:

As I know, the strand-specific and unstranded library have different counting method, so does that make sense to compare the gene expression (rpkm or rpm) from the strand-specific library and from unstranded library? If yes. which normalization method is better for this. Thanks.

rnaseq strand • 1.8k views
ADD COMMENTlink modified 6.4 years ago by Nicolas Rosewick8.3k • written 6.4 years ago by camelbbs660
gravatar for Nicolas Rosewick
6.4 years ago by
Belgium, Brussels
Nicolas Rosewick8.3k wrote:

I think you could use a multifactor design with DESeq (R - bioconductor) by specifying a design matrix like that :

Sample    condition    library
A        untreated    unstranded
B       untreated    unstranded
C       treated     unstranded
D       treated     unstranded
E       untreated    stranded
F        untreated    stranded
G        treated     stranded
H        treated     stranded

and with your count matrix :

cdsFull = newCountDataSet( countTable, designMatrix )

after that follow DESeq manual at "Multi Factor Designs" section.

I think edgeR has also this type of feature.

ADD COMMENTlink written 6.4 years ago by Nicolas Rosewick8.3k
Please log in to add an answer.


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