I am using -Sc to simulate selection, it works for this parameter combination:
msms 120 1000 -I 1 120 -t 18000 -N 100000 -r 15600 -en 0.0005 1 0.24 -SI .005 1 0.001 -Sc 0.001 1 1600 800 0 -Sp 0.5 -Smark -SFC -threads 30
But it will not work when change the time from 0.001 to 0.( actually less than 0.001, eg 0.0001)
msms 120 1000 -I 1 120 -t 18000 -N 100000 -r 15600 -en 0.0005 1 0.24 -SI .005 1 0.001 -Sc 0 1 1600 800 0 -Sp 0.5 -Smark -SFC -threads 30
java.lang.ArrayIndexOutOfBoundsException
at at.mabs.model.selection.SuperFrequencyTrace.getFrequencys(SuperFrequencyTrace.java:142)
at at.mabs.model.selection.SelectionData.getFrequency(SelectionData.java:268)
at at.mabs.model.selection.SelectionData.coalescentCumulantIntegration(SelectionData.java:296)
at at.mabs.coalescent.CoalescentEventCalculator.nextCoalescentSelection(CoalescentEventCalculator.java:361)
I am using -Sc to simulate selection, it works for this parameter combination: msms 120 1000 -I 1 120 -t 18000 -N 100000 -r 15600 -en 0.0005 1 0.24 -SI .005 1 0.001 -Sc 0.001 1 1600 800 0 -Sp 0.5 -Smark -SFC -threads 30
But it will not work when change the time from 0.001 to 0.( actually less than 0.001, eg 0.0001) msms 120 1000 -I 1 120 -t 18000 -N 100000 -r 15600 -en 0.0005 1 0.24 -SI .005 1 0.001 -Sc 0 1 1600 800 0 -Sp 0.5 -Smark -SFC -threads 30
java.lang.ArrayIndexOutOfBoundsException at at.mabs.model.selection.SuperFrequencyTrace.getFrequencys(SuperFrequencyTrace.java:142) at at.mabs.model.selection.SelectionData.getFrequency(SelectionData.java:268) at at.mabs.model.selection.SelectionData.coalescentCumulantIntegration(SelectionData.java:296) at at.mabs.coalescent.CoalescentEventCalculator.nextCoalescentSelection(CoalescentEventCalculator.java:361)
Look forward to your reply !
Best