Parametric bootstrap

135 views
Skip to first unread message

Vicente Castillo

unread,
Apr 25, 2024, 4:05:11 PM4/25/24
to fastsimcoal2
Hello all,

I have some problems understanding the parameric bootstrap. 
As far I know, a way to obtain a CI for parameter estimates is to simulate a dataset with the ML values of the selected model.

But there are 2 parameters in fastsimcoal that i dont really understand.

If my dataset contains 33,000 no linked SNPs,
the following lines are ok? 

//Number of independent loci [chromosome]
33000 0

//per Block: data type, num loci, rec. rate and mut rate + optional parameters
FREQ 1 0 1.5e-9 OUTEXP

thanks




Chenfei Zheng

unread,
Apr 26, 2024, 3:40:28 AM4/26/24
to fastsimcoal2
Hi Vicent,
I think this number of independent loci [chromosome] should be 1 when you set the data type as "FREQ"
Hope it will be helpful

yy Wu

unread,
Jan 10, 2025, 3:49:43 AMJan 10
to fastsimcoal2
I wonder if it's normal to have only one build file *_boot_1.simparam when running bootstrap? This is my command fsc28 -i *_boot.par -n1 -j -m -s0 -I -x -q, I don't know what the problem is, any suggestions are appreciated!
Reply all
Reply to author
Forward
0 new messages