Hello,
I am trying to run RAiSD on simulated data. My input file:
sfscode2ms 24 2
701933490 701933490 701933490
//
segsites: 84
positions: 0.011958 0.018982 0.025858 0.082044 0.085840 0.096328 0.098668 0.106622 0.125632 0.125632 0.137494 0.137494 0.138510 0.168142 0.168142 0.171742 0.176526 0.176526 0.189940 0.190036 0.190036 0.204400 0.204400 0.244744 0.260109 0.267261 0.273693 0.278781 0.278781 0.288465 0.300219 0.309801 0.310421 0.324129 0.358413 0.374283 0.425157 0.444129 0.463365 0.474057 0.530187 0.543843 0.550053 0.556365 0.556365 0.575541 0.575541 0.608667 0.617751 0.617751 0.623763 0.682605 0.694169 0.695707 0.712821 0.712821 0.721119 0.753524 0.787634 0.790106 0.790358 0.796694 0.812612 0.816926 0.818768 0.818768 0.823580 0.836424 0.857144 0.857144 0.869936 0.869936 0.884714 0.895322 0.925208 0.929114 0.935684 0.959486 0.964562 0.964562 0.965168 0.965168 0.972530 0.972644
010011100000000100110010000010010000000000010001100000100100000000000010100000101000
000000000101001001110010000010010000000000010000000000100000000110001010100000101000
000000100000000100010000111000100001000000010010100100100000100000000010101101101011
000000000000000000010000111000100000000000010000000000100000000010001010100000101000
000001100000000100010000111000100000000100010001100000100000000010001010100000101000
000001100000000100010000111000100000000100010010000010100000000010001010100000101000
000001100000000100010000111000100000000100010001100000100000000010001010100000101000
000001100000000100010000111000100000000100010001100000100000000010001010100000101000
000000000000000000010000000000000010000100010001100000100000000010001010100000101000
000100000101001001010010000010001000011000010010000011100001000010001010100000101000
000000100000000100010000111000100001000000010001100000100001000010001010100000101000
000001100000000000010000111000100000000101010000000000100000000010001010100000101000
000000001010010010101100000100010000000000001100000010010000000000000101100000010100
000000001010110010001100000100001000000010001100000010010100000000000001100000010100
000000001010010010101100000100001000000000001100000010011000000000000101100000010100
000000001010010010101101000100010000000000001100000010010000000000000101100000010100
001000000000010010101100000100010000000000001100000010010000000000000001101101010101
000100001010010010001100000100010100000010001100001010010000000000000001101101010101
000000001010010010001100000100001000100000001100000010010000001000000101100000010100
100000010000000000101100000100001000000000001100000010010000001000100101101111010101
000000001010010010101100000100001000000000001100010000010010000001000101110000010100
000000001010010010101100000100010100000010101100010010010000000000000001101000010100
000000001010010010101100000100001000000000001100000010010000000000010101100000010100
000000001010010010101100000101001000000010001100010010010000010000000001101101010101
//
segsites: 86
positions: 0.023872 0.044224 0.044620 0.072280 0.083470 0.083470 0.089938 0.089938 0.092692 0.111400 0.120796 0.131428 0.134266 0.138838 0.138838 0.159436 0.159436 0.166940 0.170206 0.176602 0.178666 0.194926 0.212568 0.235180 0.242508 0.299673 0.305631 0.306947 0.325563 0.336723 0.343749 0.345113 0.345819 0.348981 0.354795 0.354795 0.363247 0.375117 0.387897 0.406881 0.429075 0.439677 0.458025 0.458025 0.462255 0.481113 0.487425 0.506165 0.541941 0.570917 0.583817 0.627855 0.633297 0.633297 0.634113 0.645747 0.680435 0.692799 0.698219 0.701031 0.707139 0.715697 0.720473 0.734169 0.736353 0.757294 0.760242 0.768112 0.768112 0.776294 0.800666 0.818856 0.826316 0.826316 0.836252 0.842078 0.866162 0.908414 0.908414 0.918956 0.924590 0.939092 0.939092 0.940856 0.964070 0.980398
00000101000001000000000000010000000111000100100000000000000000100000101010000010010001
00000101100001000000100100010000000100000000100000000000000000100100101010000010000000
01000101000001001000010010000000000100000001000000000110001000010010100010010010110000
00000101000001000000100001000000000100000100100000000000100000100000101010000010010000
00000101000001001000000000000000000100000001000000010110000000000000100010000010010000
00000101000001001100001000000000000100000001000010000000000000100000100010000010010000
00000101000011000000100000100110110000000000000000000000010000100000101010000010000000
00000101010001000000100000000000000101000100100000000000000000000000101010000010010001
00000101000001000000000000101110100100000000100000000000000010000000000010000010010000
00000101000001000000100000010000000101000100100000000000000000100000101010000011010001
00000101000001001000000000100110100000000000010000000000000000100000000010000010000000
00000101000001001010000000000000000100000001000010000100000001000000100010000010010100
00001010000100110000000000000000000000000000000000000001000000000000000001000100001000
01011010000000110000000000000000001000000010000010001000000000000001000001000100001000
00001010000000100000000000000001001000000010000000101000000000000001000001000100001000
00001010000000110000000000000000001000000010000000001010001100000001000001000100001000
01011010000000110000000000000000001000000010000010001000000000000001000001000100001000
01001010000100110001000000000000001000100010000010001000001100000001000001000100001000
01001010000100100000000000000000001000000010000010001000000000001001000101100100001010
00101010000000110000000000000000001000001010001010001000000000000001000001000100001000
00001010001000110000000000000000001000000010000010001010000000000001000001000100001000
00001010000000110000000000000000001000010010000001001010001000010001010001000100001000
00001010000000110000000000000000001000000010001010001000000000000001000001000100001000
11001010000100110000000000000000001000000010000110001000000000000001000001001100001000
My command to run RAiSD:
RAiSD -I testneutral_2nm_mig_set9_S75.ms -n testneutral_2nm_mig_set9_75.out -L 500000 -R -w 40
I then get the following error:
RAiSD: sources/RAiSD.c:256: main: Assertion `(uint64_t)RSDDataset->setSNPs==RSDDataset->preLoadedsetSNPs' failed
I am using version 2.9. I am unsure what this error means - other ms files that I ran with similar number of SNPs worked fine. I would appreciate any help with this error.
Thank you!
Ashleigh
Hello,
I am trying to run RAiSD on simulated data. My input file:
sfscode2ms 24 2
701933490 701933490 701933490
//
segsites: 84
positions: 0.011958 0.018982 0.025858 0.082044 0.085840 0.096328 0.098668 0.106622 0.125632 0.125632 0.137494 0.137494 0.138510 0.168142 0.168142 0.171742 0.176526 0.176526 0.189940 0.190036 0.190036 0.204400 0.204400 0.244744 0.260109 0.267261 0.273693 0.278781 0.278781 0.288465 0.300219 0.309801 0.310421 0.324129 0.358413 0.374283 0.425157 0.444129 0.463365 0.474057 0.530187 0.543843 0.550053 0.556365 0.556365 0.575541 0.575541 0.608667 0.617751 0.617751 0.623763 0.682605 0.694169 0.695707 0.712821 0.712821 0.721119 0.753524 0.787634 0.790106 0.790358 0.796694 0.812612 0.816926 0.818768 0.818768 0.823580 0.836424 0.857144 0.857144 0.869936 0.869936 0.884714 0.895322 0.925208 0.929114 0.935684 0.959486 0.964562 0.964562 0.965168 0.965168 0.972530 0.972644
010011100000000100110010000010010000000000010001100000100100000000000010100000101000
000000000101001001110010000010010000000000010000000000100000000110001010100000101000
000000100000000100010000111000100001000000010010100100100000100000000010101101101011
000000000000000000010000111000100000000000010000000000100000000010001010100000101000
000001100000000100010000111000100000000100010001100000100000000010001010100000101000
000001100000000100010000111000100000000100010010000010100000000010001010100000101000
000001100000000100010000111000100000000100010001100000100000000010001010100000101000
000001100000000100010000111000100000000100010001100000100000000010001010100000101000
000000000000000000010000000000000010000100010001100000100000000010001010100000101000
000100000101001001010010000010001000011000010010000011100001000010001010100000101000
000000100000000100010000111000100001000000010001100000100001000010001010100000101000
000001100000000000010000111000100000000101010000000000100000000010001010100000101000
000000001010010010101100000100010000000000001100000010010000000000000101100000010100
000000001010110010001100000100001000000010001100000010010100000000000001100000010100
000000001010010010101100000100001000000000001100000010011000000000000101100000010100
000000001010010010101101000100010000000000001100000010010000000000000101100000010100
001000000000010010101100000100010000000000001100000010010000000000000001101101010101
000100001010010010001100000100010100000010001100001010010000000000000001101101010101
000000001010010010001100000100001000100000001100000010010000001000000101100000010100
100000010000000000101100000100001000000000001100000010010000001000100101101111010101
000000001010010010101100000100001000000000001100010000010010000001000101110000010100
000000001010010010101100000100010100000010101100010010010000000000000001101000010100
000000001010010010101100000100001000000000001100000010010000000000010101100000010100
000000001010010010101100000101001000000010001100010010010000010000000001101101010101
//
segsites: 86
positions: 0.023872 0.044224 0.044620 0.072280 0.083470 0.083470 0.089938 0.089938 0.092692 0.111400 0.120796 0.131428 0.134266 0.138838 0.138838 0.159436 0.159436 0.166940 0.170206 0.176602 0.178666 0.194926 0.212568 0.235180 0.242508 0.299673 0.305631 0.306947 0.325563 0.336723 0.343749 0.345113 0.345819 0.348981 0.354795 0.354795 0.363247 0.375117 0.387897 0.406881 0.429075 0.439677 0.458025 0.458025 0.462255 0.481113 0.487425 0.506165 0.541941 0.570917 0.583817 0.627855 0.633297 0.633297 0.634113 0.645747 0.680435 0.692799 0.698219 0.701031 0.707139 0.715697 0.720473 0.734169 0.736353 0.757294 0.760242 0.768112 0.768112 0.776294 0.800666 0.818856 0.826316 0.826316 0.836252 0.842078 0.866162 0.908414 0.908414 0.918956 0.924590 0.939092 0.939092 0.940856 0.964070 0.980398
00000101000001000000000000010000000111000100100000000000000000100000101010000010010001
00000101100001000000100100010000000100000000100000000000000000100100101010000010000000
01000101000001001000010010000000000100000001000000000110001000010010100010010010110000
00000101000001000000100001000000000100000100100000000000100000100000101010000010010000
00000101000001001000000000000000000100000001000000010110000000000000100010000010010000
00000101000001001100001000000000000100000001000010000000000000100000100010000010010000
00000101000011000000100000100110110000000000000000000000010000100000101010000010000000
00000101010001000000100000000000000101000100100000000000000000000000101010000010010001
00000101000001000000000000101110100100000000100000000000000010000000000010000010010000
00000101000001000000100000010000000101000100100000000000000000100000101010000011010001
00000101000001001000000000100110100000000000010000000000000000100000000010000010000000
00000101000001001010000000000000000100000001000010000100000001000000100010000010010100
00001010000100110000000000000000000000000000000000000001000000000000000001000100001000
01011010000000110000000000000000001000000010000010001000000000000001000001000100001000
00001010000000100000000000000001001000000010000000101000000000000001000001000100001000
00001010000000110000000000000000001000000010000000001010001100000001000001000100001000
01011010000000110000000000000000001000000010000010001000000000000001000001000100001000
01001010000100110001000000000000001000100010000010001000001100000001000001000100001000
01001010000100100000000000000000001000000010000010001000000000001001000101100100001010
00101010000000110000000000000000001000001010001010001000000000000001000001000100001000
00001010001000110000000000000000001000000010000010001010000000000001000001000100001000
00001010000000110000000000000000001000010010000001001010001000010001010001000100001000
00001010000000110000000000000000001000000010001010001000000000000001000001000100001000
11001010000100110000000000000000001000000010000110001000000000000001000001001100001000
My command to run RAiSD:
RAiSD -I testneutral_2nm_mig_set9_S75.ms -n testneutral_2nm_mig_set9_75.out -L 500000 -R -w 40
I then get the following error:
RAiSD: sources/RAiSD.c:256: main: Assertion `(uint64_t)RSDDataset->setSNPs==RSDDataset->preLoadedsetSNPs' failed
I am using version 2.9. I am unsure what this error means - other ms files that I ran with similar number of SNPs worked fine. I would appreciate any help with this error.
Thank you!
Ashleigh