DRIMSeq with technical replicates
1
0
Entering edit mode
@poecilepal-22574
Last seen 3.1 years ago
Russia

Good afternoon,

I use DRIMSeq-workflow from this page http://bioconductor.org/packages/release/workflows/vignettes/rnaseqDTU/inst/doc/rnaseqDTU.html#statistical-analysis-of-differential-transcript-usage to explore differential transcript usage. Could you please point me out, where and how should I collapse technical replicates in DRIMSeq? For example, for DEG search in DESeq2 I use collapseReplicates function before running DESeq.

Best regards, Poecile palustris

DRIMSeq R DTU technical replicates deseq2 • 1.4k views
ADD COMMENT
0
Entering edit mode
ADD REPLY
0
Entering edit mode

Hello Kevin,

Yes, thanks, I have duplicated this post on another forum just in case.

Best regards, Poecile palustris

ADD REPLY
0
Entering edit mode

Oh, I read about disadvantages of cross-hosting, I apologize for my ignorance. It seems that post https://www.biostars.org/p/450220/ has already been deleted.

ADD REPLY
0
Entering edit mode
Kevin Blighe ★ 4.0k
@kevin
Last seen 6 weeks ago
Republic of Ireland

It seems that the DRIMSeq authors have not yet checked - you may consider contacting them. Also, you could still use the DESeq2 functionality here and then export the 'collapsed' counts (via counts(dds, normalised = FALSE)) for use with DRIMSeq. You could also apply your own function manually for collapsing your replicates.

Kevin

ADD COMMENT
0
Entering edit mode

Thank you very much for your response! I used DESeq2 functionality: firstlycollapseReplicates function from DESeq2, then divided each value by the number of replicates.

Best regards, Poecile palustris

ADD REPLY

Login before adding your answer.

Traffic: 1255 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6