Hi all,
I am trying to use fastsimcoal2 (fsc28) to estimate the population history and gene flow, but here is something wrong. I will show my data and work step, hoping someone can help me.
*****here is detail******
I have two populations. SNP called and filtered with gatk. easySFS used to calculate the folded SFS with command, --proj 24,28 --order Ax,Ex --dtype int -a -f, *_jointMAFpop1_0.obs file directly used for fastsimcoal2.
****here is command for fastsimcoal2******
# 如果文件是Windows换行,先转一下
#dos2unix AxEx_AncientANDRecentGeneFlow.tpl AxEx_AncientANDRecentGeneFlow.est AxEx_AncientANDRecentGeneFlow_jointMAFpop1_0.obs 2>/dev/null
ulimit -s 1024000
#fsc28模型测试
${fsc28} \
-t ${tpl} \
-e ${est} \
-m -M \
-n 10 \
-L 2 \
-c 1 \
--removeZeroSFS
*_jointMAFpop1_0.obs, .tpl and .est file is uploaded.
Thanks for your help !
kami