schraiber / selection Goto Github PK

Dear Schraiber
I am interest in your paper, and tried the asip data, while I cannot get your result, I post all my files and command here, please tell me is there something wrong about it?
sample file
0 10 -0.078 -0.078
1 22 -0.051 -0.051
15 20 -0.014 -0.014
12 20 -0.011 -0.011
15 36 -0.004 -0.004
18 38 -0.002 -0.002
popfile
16000 0 -0.0390625
110000 -20.956 -0.2109375
110000 0 -0.307109375
195019 0 -0.414019531
107519 0 -0.540960938
82404 0 -0.691695313
124078 0 -0.870679688
171778 0 -1.083210938
181878 0 -1.335566406
112777 0 -1.63521875
45266 0 -1.991035156
20816 0 -2.41353125
16110 0 -2.915207031
23901 0 -3.510910156
53320 0 -4.21825
105982 0 -5.05815625
138531 0 -6.055472656
127587 0 -7.239695313
98370 0 -8.645859375
74423 0 -10.31555469
63543 0 -12.29817578
66523 0 -14.65235547
83610 0 -17.44775391
112343 0 -20.76702734
142428 0 -24.70838672
161440 0 -Inf
./sr -D asip.input -G 8 -N 16000 -n 500000 -d 0.001 -F 20 -f 1000 -s 100 -P horse.all.pop -e 8067 -a -o asip
Thankyou so much!
Best
xue

Hi Joshua,

I have my files formatted to run for allele age + selection estimation, but I'm thinking I must have gone wrong somewhere as the program keeps throwing an error "Times are not in decreasing order!"

The command I run is: /ohta2/julia.kreiner/software/selection/sr -G 1 -N 109868 -D test_derivedcounts.txt -P dom_popsizehistory_deorder.txt -o $path/output/${lname} -a, where I've kept years and pop size unscaled.

Here's the head of my popfile and count input file, and I've also attached a full version for inpsection/testing.

test_derivedcounts.txt
13 16 0 0
0 0 -89 -89
2 2 -90 -90
4 4 -91 -91
4 4 -97 -97
2 2 -98 -98
2 2 -99 -99
2 2 -100 -100
2 2 -410 -410
0 0 -418 -418

dom_popsizehistory_deorder.txt
109868.20575394991 0 0
35014.67057441156 0 -89
35014.67057441156 0 -90
35014.67057441156 0 -91
35150.23096022847 0 -97
35150.23096022847 0 -98
35150.23096022847 0 -99
35150.23096022847 0 -100
58640.33186376603 0 -410
58640.33186376603 0 -418

Thanks for your help!
Julia