bootstrap replicate data

24 views
Skip to first unread message

Kathryn Greiner-Ferris

unread,
Aug 12, 2025, 1:41:56 PMAug 12
to Stacks
I would like to calculate the confidence intervals for the pairwise fst values calculated by STACKS after bootstrapping. Is there a way to extract the  values for each pop pair bootstrap replicate? I was hoping that this could be found in populations.bsarchive copy.tsv file, but the content does not change even after bootstrapping, so I assume it is just the the original per-locus data from your initial sampling. Will using --log-fst-comp give me the data I am looking for? 

If so.. I got an error when I would did add --log-fst-comp to the end of my code: 

populations: 
unrecognized option '--log-fst-comp'
populations 2.62
Usage:
populations -P dir [-O dir] [-M popmap] (filters) [--fstats] [-k [--sigma=150000] [--bootstrap [-N 100]]] (output formats)
populations -V vcf -O dir [-M popmap] (filters) [--fstats] [-k [--sigma=150000] [--bootstrap [-N 100]]] (output formats)


This is how I used it: sbatch ./scripts3/whitelist.slurm -i ./optimize_sync_new/refmap -o ./optimize_sync_new/refmap/output_bootstrap -M ./info/synsc_popmap.tsv -W ./info/synsc_WL.tsv -X ' --fstats --fst-correction --bootstrap-fst --bootstrap-reps 1000 --log-fst-comp'





Reply all
Reply to author
Forward
0 new messages