MDS plot for miRNA data
0
0
Entering edit mode
6.1 years ago
Adrian Pelin ★ 2.5k

Hello,

I have sequenced human miRNA (total RNA enriched for 15-25nt) in different conditions. I mapped my reads with a minimum overlap of 15 + perfect match to the human genome, and extracted read counts using bedtools of miRBase miRNA annotations.

I now want to generally compare my different conditions and see which ones are more similar to each other at the small RNA level. I remember it is possible to do that with cummeRBound package in R, but it require cufflinks output, and it is my understanding that cufflinks is not necessarly a good option for miRNA.

Anyone have any idea how to do this?

Also, I read in some papers that people add a pseudocount of 1 to all read counts so as not to encounter division by 0 errors.

Thanks,

Adrian

RNA-Seq smallRNA MDS • 1.6k views
ADD COMMENT
1
Entering edit mode

Bioconductor's RNAseq workflow may be of help. There are sections on PCA and MDS, and the data expected are just read counts, so your bedtools counts should suffice.

ADD REPLY
0
Entering edit mode

Nice, I see it. Thanks, this is the answer.

ADD REPLY

Login before adding your answer.

Traffic: 2770 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