Hi everyone:
Now I want to normalize ChIP-seq bed file reads to certain number by randomly choose certain number of reads from the original bed file.
Dose anyone have such kind of scripts to do this? I tried to do the perl script by ran() operator. It will work for small bed file. But it will calculate forever if the reads if around 38 million. So is there any way I can unsort the @array element? So that I can use head to form a certain number of reads by randomly selection. Or Does anyone have better idea about that?
Thank you
Lijing
Now I want to normalize ChIP-seq bed file reads to certain number by randomly choose certain number of reads from the original bed file.
Dose anyone have such kind of scripts to do this? I tried to do the perl script by ran() operator. It will work for small bed file. But it will calculate forever if the reads if around 38 million. So is there any way I can unsort the @array element? So that I can use head to form a certain number of reads by randomly selection. Or Does anyone have better idea about that?
Thank you
Lijing
Comment