ANGSD/ngsPSMC

Error: -> Problem opening file: '-dospline'

Opened this issue · 2 comments

I tested out running ngsPSMC using the instructions outlined in the tutorial.

I first ran this script on a single .bam file to generate an input file

./angsd -i S13_2949_trim_cleaned.bam -ref canFam3_withY.fa -out S13_2949_psmcinput -dopsmc 1 -gl 1 -minq 30 -minmapq 30

Looking at the bottom of my .arg file it appears it worked and generated the following output files

abcPSMC.cpp:
-writePSMC 1
1: (still beta, not really working)
->"S13_2949_psmcinput.arg"
->"S13_2949_psmcinput.psmc.gz"
->"S13_2949_psmcinput.psmc.pos.gz"
->"S13_2949_psmcinput.psmc.idx"

    [ALL done] cpu-time used =  4683.03 sec
    [ALL done] walltime used =  4776.00 sec

I then tried running ngsPSMC and am getting the following error:

./ngsPSMC S13_2949_psmcinput.psmc.idx -p "1*4+25*2+1*4+1*6" -dospline 0 -nthreads 8 -nIter 20 -init 1  -theta 0.000233095 -rho 0.005357 > S13_2949.psmc.ml.txt
-> we are in file: main_psmc.cpp function: main_psmc line:187
    -> Version of fname: '/home/sophiepq/GrayFoxWGS/GrayFox1/angsd/ngsPSMC/S13_2949_psmcinput.psmc.idx' is:1
    -> Assuming .psmc.gz file: /home/sophiepq/GrayFoxWGS/GrayFox1/angsd/ngsPSMC/S13_2949_psmcinput.psmc.gz
    -> Assuming .psmc.pos.gz: /home/sophiepq/GrayFoxWGS/GrayFox1/angsd/ngsPSMC/S13_2949_psmcinput.psmc.pos.gz
    -> Problem opening file: '-dospline'

I don't really know what this error means or how to address it. I believe -dospline 0 or -dospline 1 are requirements to run program. I'm not sure whether there is another input file it's looking for that I am missing here.

Hi @squisquater I am having the exact same issue – were you able to resolve it somehow?

Alternatively @ANGSD if you have any ideas why this might be happening your guidance would be much appreciated