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