Error after running simulations

[charlottebourbon]

Hello,
I have one population and 6 scenarii for determining bottlnecks and effective population size variations. After checking the data and the prior settings, I ran the simulations but after fews runs, a problem occured and the program stopped. The run log finished by "coal_pop : population size <1 (0) ". I don't understand the error. Could you help me ?

There is the run log :

debut
headerRF.txt found, choosing it over header.txt
simulating data sets with all summary statistics
I have read RNGs' states from file ./RNG_state_0000.bin
debut de l'action r
./test2.DIYABC.snp.bin
insert title here
Fichier ./test2.DIYABC.snp : SNP
lecture du fichier binaire
smaf=0.05
THIS->MAF=0.05
sexratio=0.714286
nloc=42427
nsample=1
type 0 catexist =1
type 1 catexist =0
type 2 catexist =0
type 3 catexist =0
type 4 catexist =0
type 0
36
this->maf=0.05
dans libin
this->locus[0].haplosnp[0][0]=1
sexratio=0.714286
fin de la lecture du fichier binaire
avant missingdata
apres missingdata
nparamtot=13
nstat=4
----------------------------------------apres readHeaderDebut
6 scenario(s)
----------------------------------------apres readHeaderScenarios
historical parameters priors (13,9)
----------------------------------------apres readHeaderHistParam
fichier SNP dans redheaderLoci
42427 G1 from 1
prem = 0
groupe[1].nloc=42427
avant recopie des locus
dataobs.locus.size()=42427
dataobs.nloc = 42427
----------------------------------------apres readHeaderLoci
----------------------------------------apres readHeaderGroupPrior
----------------------------------------apres buildSuperScen
debut de readHeaderAllStat
nsamp=1 nstat=4 nparamut=0
./headerRFRF.txt
nstat=4

this->entete=
scenario Ne N1 N2 N3 N4 N5 N6 N7 N8 N9 t1 t2 t3 HWm_1 HWv_1 FST1m_1 FST1v_1
----------------------------------------apres readHeaderAllStat
----------------------------------------apres buildMutParam
fin de readheader nloc retenus=42427
Renaming headerRFRF.txt to headerRF.txt
calstatobs (./statobsRF.txt)
stat_obs[0]=0.39952
stat_obs[1]=0.0133501
stat_obs[2]=0
stat_obs[3]=0

verification de l'integrite de la table de reference
fichier./reftableRF.bin
nrec = 2390
nrecscen[0] = 395
nrecscen[1] = 441
nrecscen[2] = 368
nrecscen[3] = 396
nrecscen[4] = 411
nrecscen[5] = 379
nparam[0] = 1
nparam[1] = 3
nparam[2] = 3
nparam[3] = 5
nparam[4] = 5
nparam[5] = 7
nstat = 4
6 scenario(s)
40
80
120
160
200
240
280
320
360
400
440
480
520
560
600
640
680
720
760
800
840
880
920
960
1000
1040
1080
1120
1160
1200
1240
1280
1320
1360
1400
1440
1480
1520
1560
1600
1640
1680
1720
1760
1800
1840
1880
1920
1960
2000
2040
2080
2120
2160
2200
2240
2280
2320
2360
2400
fichier reftable OK

dans rt.readheader this->filerefscen = ./reftabscen.txt
nparam[0] = 1
nparam[1] = 3
nparam[2] = 3
nparam[3] = 5
nparam[4] = 5
nparam[5] = 7
debut de sethistparamname
nparamut=0 nscenarios=6
avant la boucle des scenarios nscenarios =6
scenario 0 scenario[i].nparam=1 nparamvar=1
scenario 1 scenario[i].nparam=3 nparamvar=3
scenario 2 scenario[i].nparam=3 nparamvar=3
scenario 3 scenario[i].nparam=5 nparamvar=5
scenario 4 scenario[i].nparam=5 nparamvar=5
scenario 5 scenario[i].nparam=7 nparamvar=7
fin de sethistparamname
apres readheader k=0

nloc = 42427 nsample = 1 k=0
DEBUT nrecneeded=300000 rt.nrec=2390 rt.nstat=4 nscenarios=6

Ouverture du fichier
./reftableRF.bin

nparammax=16
n simu data = 2430
n simu data = 2470
n simu data = 2510
n simu data = 2550
n simu data = 2590
coal_pop : population size <1 (0)

[gdurif]

Hi, I am really sorry for the very late reply.

Due to work overload, the development, support and maintenance for the project diyabcGUI is currently paused. Feel free to contact us (by submitting an issue) if you want to take over the development for this project. Pull request are also welcome to fix bug or implement missing functionalities.

I have updated the README file to announce this. This is the official announcement to clarify a situation that has been going on for quite some time.

You can still use the command line version of the software which is still maintained/supported. Check the documentation page here.

All my apologies for the trouble,
Best regards