likelihoods always 0.000

225 views
Skip to first unread message

Jeremy Davis

unread,
May 11, 2022, 11:59:52 AM5/11/22
to fastsimcoal2
Hello,
I am trying to get several models working with some different parameters to compare most likely scenarios. After some troubleshooting I have gotten Fsc to run, however my results for every iteration so far have had a likelihood of "0.000" and is the same no matter what parameters are attempted. 

My question is, why is this the likelihood I am getting for every iteration, and where am I going wrong?

I have 64 diploid individuals in one population, and got my .obs file from ANGSD. I have attached the output file as well as my .obs file. The output file also includes my command line, but here it is again:

fsc2705 -t model0.tpl -n 1000 -m -e model0.est -M -L 40 -q 

TPL:
//Number of population samples (demes)
1
//Population effective sizes (number of genes)
NNA
//Sample sizes
128
//Growth rates : negative growth implies population expansion
0
//Number of migration matrices : 0 implies no migration between demes
0
//historical event: time, source, sink, migrants, new size, new growth rate, migr. matrix  
2 historical event
$TENDBOT$ 0 0 1 $NBOT$ 0 0 absoluteResize
$TBOT$ 0 0 1 $NEU$ 0 0 absoluteResize
//Number of independent loci [chromosome]
1 0
//Per chromosome: Number of linkage blocks
1
//per Block: data type, num loci, rec. rate and mut rate + optional parameters
FREQ 1 0 3.5e-9

EST:
// Priors and rules file  
// *********************  
[PARAMETERS]  
//#isInt? #name #dist.#min #max  
//all Ns are in number of haploid individuals  
1 $NEU$ unif 100000 100000000 output  
1 $NBOT$ unif 100 50000 output
1 $NNA$ unif 100000 10000000 output
1 $TBOT$ unif 100 300 output
1 $TENDBOT$ unif 50 $TBOT$ output paramInRange
model0_MAFpop0.obs
model0.out

Jack Harper

unread,
Jan 7, 2025, 1:10:10 PMJan 7
to fastsimcoal2
Hi Jeremy,

Did you manage to find a solution to this? I am currently struggling with the same issue.

Thanks!
Jack

Laurent Excoffier

unread,
Jan 8, 2025, 4:48:32 AMJan 8
to fastsimcoal2
Hi,
it seems as if your obs file is not properly formatted (Missing header). Have a look at the manual on how to do that
best
Laurent

Reply all
Reply to author
Forward
0 new messages