bootstrapping help
0
0
Entering edit mode
3.6 years ago

Hi,

I have a script to do bootstrapping, and I've successfully bootstrapped a vcf file of a single chromosome into 50,000 bp chunks, and resampled until the size of the original vcf file was reached. However, because when you bootstrap you end up with duplicates within each replicate, I am receiving the following error when inputting the data into a program (SMC++ but I don't think it's important here): "WARNING Multiple entries found at 498,975 positions; skipped all but the first".

Is there a way to assign new positions to every region, so that even the duplicated regions are no longer read that way into the program?

Thank you!

genome bootstrap vcf bootstrapping SNP • 777 views
ADD COMMENT

Login before adding your answer.

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