Hi all,
I'm having trouble restarting an IMa2 run from the .mcf file (using the -r3 option). I am using this command line (which is identical to the command line in the original analysis, except for a different seed, shorter burnin and use of -r3 and -f followed by the name of the .mcf file):
./IMa2 -b1000 -hfg -hn120 -ha0.999 -hb0.3 -i PwMtEste.txt -l6.0 -q50.0 -m0.5 -t10.0 -p23567 -u1 -s1454584 -o PWMtEste.out -y1 -r3 -f PWMtEste.out.mcf
And after the program reads the loci info, and is about to start the Markov chain, I get this error message:
Assertion failed: (pdgoldsum == C[ci]->allpcalc.pdg), function changet_RY1, file update_t_RY.cpp, line 442.
Abort trap: 6
Anyone else having this error? Is this a memory error or something? The .mcf file is big, about 200 MB...
Any advice will be highly appreciated, the computer crashed when IMa2 had been running for about 400 hours and I wouldn't want to start over...
Thanks,