How can you check if the downsampler in indeed random and efficient
1
1
Entering edit mode
5.6 years ago
shruti ▴ 10

Hi, I am using the Picard downsampler for downsampling around 50 million reads from the entire human transciptome for my experiment. I want to know if the downsampler is indeed random and not creating bias. Is there anyway I can check that ? 

-Shruti

RNA-Seq • 1.4k views
ADD COMMENT
0
Entering edit mode

Downsampling reads from a bam file using DownsampleSam?

ADD REPLY
0
Entering edit mode

Yes, i've used tophat for generating the BAM files from which I have used DownsampleSam from Picard.

ADD REPLY
2
Entering edit mode
5.6 years ago

One way is to validate the results by generating different random sets and compare them wrt your goal. For e.g. If you have annotations, you can do quantification of genes, and check if they correlate between different random sets to get an idea.

Otherwise, if you go into details, it uses java.util.Random, which depends on ["Linear congruential generator"] algorithm to generate random numbers.

ADD COMMENT
0
Entering edit mode

Thank you @Goutham ! 

ADD REPLY

Login before adding your answer.

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