DRIMSeq with technical replicates
1
0
Entering edit mode
@poecilepal-22574
Last seen 10 months 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 • 468 views
0
Entering edit mode
0
Entering edit mode

Hello Kevin,

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

Best regards, Poecile palustris

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.

0
Entering edit mode
@kevin
Last seen 5 hours 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

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