ngsLD sometimes hangs indefinitely
James-S-Santangelo opened this issue · 8 comments
Hello!
I've come across an interesting issue where ngsLD
will sometimes hang on "Reading data from file..." and never proceed any further, even after approximately 48 hours.
I noticed this problem while estimate pairwise LD in parallel across 16 separate chromosomes. Some chromosomes would finish within a couple hours, while others would just hang and eventually get terminated by the cluster due to timing out.
I'm estimating pairwise LD from genotype likelihoods that I calculated in ANGSD
using the following command:
angsd -GL 1 \
-out {0}/angsd/{{wildcards.chrom}}/{{wildcards.chrom}}_genolike_allSamples \
-nThreads {{threads}} \
-doGlf 2 \
-doMajorMinor 1 \
-SNP_pval 1e-6 \
-doMaf 1 \
-doCounts 1 \
-setMinDepthInd 3 \
-setMaxDepth 4500 \
-baq 2 \
-ref {1} \
-minInd 96 \
-minQ 20 \
-minMapQ 30 \
-doSaf 1 \
-anc {1} \
-r {{wildcards.chrom}} \
-bam {{input.bams}} 2> {{log}}
This command succeeds without issue for all chromosomes and generates all output files, including the BEAGLE format genotype likelihoods used as input to ngsLD
. Here are the first 10 lines of one of these files:
marker allele1 allele2 Ind0 Ind0 Ind0 Ind1 Ind1 Ind1 Ind2 Ind2 Ind2 Ind3 Ind3 Ind3 Ind4 Ind4 Ind4 Ind5 Ind5 Ind5 Ind6 Ind6 Ind6 Ind7 Ind7 Ind7 Ind8 Ind8 Ind8 Ind9 Ind9 Ind9 Ind10 Ind10 Ind10 Ind11 Ind11 Ind11 Ind12 Ind12 Ind12 Ind13 Ind13 Ind13 Ind14 Ind14 Ind14 Ind15 Ind15 Ind15 Ind16 Ind16 Ind16 Ind17 Ind17 Ind17 Ind18 Ind18 Ind18 Ind19 Ind19 Ind19 Ind20 Ind20 Ind20 Ind21 Ind21 Ind21 Ind22 Ind22 Ind22 Ind23 Ind23 Ind23 Ind24 Ind24 Ind24 Ind25 Ind25 Ind25 Ind26 Ind26 Ind26 Ind27 Ind27 Ind27 Ind28 Ind28 Ind28 Ind29 Ind29 Ind29 Ind30 Ind30 Ind30 Ind31 Ind31 Ind31 Ind32 Ind32 Ind32 Ind33 Ind33 Ind33 Ind34 Ind34 Ind34 Ind35 Ind35 Ind35 Ind36 Ind36 Ind36 Ind37 Ind37 Ind37 Ind38 Ind38 Ind38 Ind39 Ind39 Ind39 Ind40 Ind40 Ind40 Ind41 Ind41 Ind41 Ind42 Ind42 Ind42 Ind43 Ind43 Ind43 Ind44 Ind44 Ind44 Ind45 Ind45 Ind45 Ind46 Ind46 Ind46 Ind47 Ind47 Ind47 Ind48 Ind48 Ind48 Ind49 Ind49 Ind49 Ind50 Ind50 Ind50 Ind51 Ind51 Ind51 Ind52 Ind52 Ind52 Ind53 Ind53 Ind53 Ind54 Ind54 Ind54 Ind55 Ind55 Ind55 Ind56 Ind56 Ind56 Ind57 Ind57 Ind57 Ind58 Ind58 Ind58 Ind59 Ind59 Ind59 Ind60 Ind60 Ind60 Ind61 Ind61 Ind61 Ind62 Ind62 Ind62 Ind63 Ind63 Ind63 Ind64 Ind64 Ind64 Ind65 Ind65 Ind65 Ind66 Ind66 Ind66 Ind67 Ind67 Ind67 Ind68 Ind68 Ind68 Ind69 Ind69 Ind69 Ind70 Ind70 Ind70 Ind71 Ind71 Ind71 Ind72 Ind72 Ind72 Ind73 Ind73 Ind73 Ind74 Ind74 Ind74 Ind75 Ind75 Ind75 Ind76 Ind76 Ind76 Ind77 Ind77 Ind77 Ind78 Ind78 Ind78 Ind79 Ind79 Ind79 Ind80 Ind80 Ind80 Ind81 Ind81 Ind81 Ind82 Ind82 Ind82 Ind83 Ind83 Ind83 Ind84 Ind84 Ind84 Ind85 Ind85 Ind85 Ind86 Ind86 Ind86 Ind87 Ind87 Ind87 Ind88 Ind88 Ind88 Ind89 Ind89 Ind89 Ind90 Ind90 Ind90 Ind91 Ind91 Ind91 Ind92 Ind92 Ind92 Ind93 Ind93 Ind93 Ind94 Ind94 Ind94 Ind95 Ind95 Ind95 Ind96 Ind96 Ind96 Ind97 Ind97 Ind97 Ind98 Ind98 Ind98 Ind99 Ind99 Ind99 Ind100 Ind100 Ind100 Ind101 Ind101 Ind101 Ind102 Ind102 Ind102 Ind103 Ind103 Ind103 Ind104 Ind104 Ind104 Ind105 Ind105 Ind105 Ind106 Ind106 Ind106 Ind107 Ind107 Ind107 Ind108 Ind108 Ind108 Ind109 Ind109 Ind109 Ind110 Ind110 Ind110 Ind111 Ind111 Ind111 Ind112 Ind112 Ind112 Ind113 Ind113 Ind113 Ind114 Ind114 Ind114 Ind115 Ind115 Ind115 Ind116 Ind116 Ind116 Ind117 Ind117 Ind117 Ind118 Ind118 Ind118 Ind119 Ind119 Ind119
CM019102.1_30951 2 0 0.984616 0.015384 0.000000 0.998051 0.001949 0.000000 0.999992 0.000008 0.000000 0.999512 0.000488 0.000000 0.996109 0.003891 0.000000 0.333333 0.333333 0.333333 0.333333 0.333333 0.333333 0.969698 0.030302 0.000000 0.000041 0.999959 0.000000 0.999024 0.000976 0.000000 0.999024 0.000976 0.000000 0.996109 0.003891 0.000000 0.998051 0.001949 0.000000 0.333333 0.333333 0.333333 0.399098 0.600902 0.000000 0.999878 0.000122 0.000000 0.941178 0.058822 0.000000 0.010308 0.989692 0.000000 0.999512 0.000488 0.000000 0.992249 0.007751 0.000000 0.996109 0.003891 0.000000 0.092094 0.907906 0.000000 0.001593 0.998406 0.000000 0.333333 0.333333 0.333333 0.333333 0.333333 0.333333 0.333333 0.333333 0.333333 0.142436 0.857564 0.000000 0.999756 0.000244 0.000000 0.999969 0.000031 0.000000 0.024890 0.975110 0.000000 0.999024 0.000976 0.000000 0.992249 0.007751 0.000000 0.333333 0.333333 0.333333 0.449468 0.550532 0.000000 0.000020 0.999980 0.000000 0.076885 0.923115 0.000000 0.888891 0.111109 0.000000 0.006310 0.993690 0.000000 0.996109 0.003891 0.000000 0.998051 0.001949 0.000000 0.024890 0.975110 0.000000 0.999756 0.000244 0.000000 0.333333 0.333333 0.333333 0.333333 0.333333 0.333333 0.000001 0.999999 0.000000 0.765539 0.234461 0.000000 0.999939 0.000061 0.000000 0.000020 0.999980 0.000000 0.000073 0.999927 0.000000 0.998051 0.001949 0.000000 0.941178 0.058822 0.000000 0.941178 0.058822 0.000000 0.888891 0.111109 0.000000 0.999985 0.000015 0.000000 0.888891 0.111109 0.000000 0.999024 0.000976 0.000000 0.969698 0.030302 0.000000 0.984616 0.015384 0.000000 0.333333 0.333333 0.333333 0.024890 0.975110 0.000000 0.333333 0.333333 0.333333 0.169534 0.830466 0.000000 0.984616 0.015384 0.000000 0.941178 0.058822 0.000000 0.169534 0.830466 0.000000 0.969698 0.030302 0.000000 0.999512 0.000488 0.000000 0.999024 0.000976 0.000000 0.888891 0.111109 0.000000 0.992249 0.007751 0.000000 0.333333 0.333333 0.333333 0.984616 0.015384 0.000000 0.333333 0.333333 0.333333 0.992249 0.007751 0.000000 0.941178 0.058822 0.000000 0.999756 0.000244 0.000000 0.048566 0.951434 0.000000 0.941178 0.058822 0.000000 0.984616 0.015384 0.000000 0.333333 0.333333 0.333333 0.761488 0.238512 0.000000 0.984616 0.015384 0.000000 0.333333 0.333333 0.333333 0.999024 0.000976 0.000000 0.984616 0.015384 0.000000 0.333333 0.333333 0.333333 0.888891 0.111109 0.000000 0.333333 0.333333 0.333333 0.289899 0.710101 0.000000 0.888891 0.111109 0.000000 0.004109 0.995891 0.000000 0.000013 0.999987 0.000000 0.984616 0.015384 0.000000 0.992249 0.007751 0.000000 0.996109 0.003891 0.000000 0.969698 0.030302 0.000000 0.888891 0.111109 0.000000 0.449468 0.550532 0.000000 0.000081 0.999919 0.000000 0.000003 0.999997 0.000000 0.992249 0.007751 0.000000 0.000000 1.000000 0.000000 0.941178 0.058822 0.000000 0.999939 0.000061 0.000000 0.999878 0.000122 0.000000 0.941178 0.058822 0.000000 0.992249 0.007751 0.000000 0.048566 0.951434 0.000000 0.999996 0.000004 0.000000 0.074818 0.925182 0.000000 0.999512 0.000488 0.000000 0.000000 1.000000 0.000000 0.992249 0.007751 0.000000 0.998051 0.001949 0.000000 0.996109 0.003891 0.000000 0.999024 0.000976 0.000000 0.992249 0.007751 0.000000 0.996109 0.003891 0.000000 0.984616 0.015384 0.000000 0.000010 0.999990 0.000000
CM019102.1_30958 0 2 0.000000 1.000000 0.000000 0.999024 0.000976 0.000000 0.981211 0.018789 0.000000 0.245248 0.754752 0.000000 0.998051 0.001949 0.000000 0.333333 0.333333 0.333333 0.333333 0.333333 0.333333 0.000020 0.999980 0.000000 0.999512 0.000488 0.000000 0.962152 0.037848 0.000000 0.000000 1.000000 0.000000 0.996109 0.003891 0.000000 0.000000 1.000000 0.000000 0.969698 0.030302 0.000000 0.000041 0.999959 0.000000 0.999985 0.000015 0.000000 0.941178 0.058822 0.000000 1.000000 0.000000 0.000000 0.620161 0.379839 0.000000 0.024890 0.975110 0.000000 0.444653 0.555347 0.000000 0.000497 0.999503 0.000000 0.001593 0.998405 0.000001 0.001593 0.998406 0.000000 0.333333 0.333333 0.333333 0.333333 0.333333 0.333333 0.999880 0.000120 0.000000 0.000000 1.000000 0.000000 0.999992 0.000008 0.000000 0.998051 0.001949 0.000000 0.864413 0.135587 0.000000 0.992249 0.007751 0.000000 0.333333 0.333333 0.333333 0.001301 0.998699 0.000000 0.449468 0.550532 0.000000 0.142775 0.857225 0.000000 0.003181 0.996819 0.000000 0.001593 0.998403 0.000004 0.000000 1.000000 0.000000 0.996109 0.003891 0.000000 0.999512 0.000488 0.000000 0.999996 0.000004 0.000000 0.333333 0.333333 0.333333 0.333333 0.333333 0.333333 0.010308 0.989692 0.000000 0.001014 0.998986 0.000000 0.999999 0.000001 0.000000 0.999878 0.000122 0.000000 0.999512 0.000488 0.000000 0.999756 0.000244 0.000000 0.941178 0.058822 0.000000 0.557681 0.442319 0.000000 0.888891 0.111109 0.000000 0.999969 0.000031 0.000000 0.969698 0.030302 0.000000 0.927171 0.072829 0.000000 0.999024 0.000976 0.000000 0.999939 0.000061 0.000000 0.333333 0.333333 0.333333 0.998051 0.001949 0.000000 0.333333 0.333333 0.333333 0.777677 0.222323 0.000000 0.048566 0.951434 0.000000 0.941178 0.058822 0.000000 0.999024 0.000976 0.000000 0.999024 0.000976 0.000000 0.000000 1.000000 0.000000 0.999024 0.000976 0.000000 0.001593 0.998406 0.000000 0.286207 0.713793 0.000000 0.333333 0.333333 0.333333 0.996109 0.003891 0.000000 0.333333 0.333333 0.333333 0.992249 0.007751 0.000000 0.984616 0.015384 0.000000 0.999878 0.000122 0.000000 0.000000 1.000000 0.000000 0.003181 0.996819 0.000000 0.000010 0.999990 0.000000 0.333333 0.333333 0.333333 0.999512 0.000488 0.000000 0.006342 0.993658 0.000000 0.888891 0.111109 0.000000 0.998051 0.001949 0.000000 0.996109 0.003891 0.000000 0.969698 0.030302 0.000000 0.000005 0.999995 0.000000 0.333333 0.333333 0.333333 0.999992 0.000008 0.000000 0.888891 0.111109 0.000000 0.000025 0.999975 0.000000 0.998051 0.001949 0.000000 0.167221 0.832779 0.000000 0.999024 0.000976 0.000000 0.092626 0.907374 0.000000 0.000174 0.999826 0.000000 0.941178 0.058822 0.000000 0.999999 0.000001 0.000000 0.999985 0.000015 0.000000 0.999024 0.000976 0.000000 0.996109 0.003891 0.000000 0.999878 0.000122 0.000000 0.984616 0.015384 0.000000 0.867191 0.132809 0.000000 0.001301 0.998699 0.000000 0.006342 0.993658 0.000000 0.996109 0.003891 0.000000 0.761488 0.238512 0.000000 0.000000 1.000000 0.000000 0.999939 0.000061 0.000000 0.001585 0.998415 0.000000 0.999512 0.000488 0.000000 0.286207 0.713793 0.000000 0.999512 0.000488 0.000000 0.008237 0.991763 0.000000 0.000662 0.999338 0.000000 0.996109 0.003891 0.000000 0.169534 0.830466 0.000000 0.996109 0.003891 0.000000 0.169534 0.830466 0.000000
CM019102.1_30963 1 3 0.000000 1.000000 0.000000 0.000020 0.999980 0.000000 0.000016 0.999984 0.000000 0.999996 0.000004 0.000000 0.028342 0.971658 0.000000 0.333333 0.333333 0.333333 0.333333 0.333333 0.333333 0.092626 0.907374 0.000000 0.449468 0.550532 0.000000 0.999024 0.000976 0.000000 0.999985 0.000015 0.000000 0.996109 0.003891 0.000000 0.998051 0.001949 0.000000 0.941178 0.058822 0.000000 0.000010 0.999990 0.000000 0.079841 0.920159 0.000000 0.941178 0.058822 0.000000 0.972459 0.027541 0.000000 0.999512 0.000488 0.000000 0.024890 0.975110 0.000000 0.941178 0.058822 0.000000 0.992249 0.007751 0.000000 0.969698 0.030302 0.000000 0.003181 0.996819 0.000000 0.333333 0.333333 0.333333 0.333333 0.333333 0.333333 0.981106 0.018894 0.000000 0.000115 0.999885 0.000000 0.000178 0.999822 0.000000 0.000107 0.999893 0.000000 0.999756 0.000244 0.000000 0.615213 0.384787 0.000000 0.333333 0.333333 0.333333 0.000016 0.999984 0.000000 0.000081 0.999919 0.000000 0.000270 0.999730 0.000000 0.444653 0.555347 0.000000 0.888891 0.111109 0.000000 0.999024 0.000976 0.000000 0.000034 0.999966 0.000000 0.000000 1.000000 0.000000 0.140893 0.859107 0.000000 0.333333 0.333333 0.333333 0.333333 0.333333 0.333333 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.308111 0.691889 0.000000 0.002545 0.997455 0.000000 0.119771 0.880229 0.000000 0.006342 0.993658 0.000000 0.992249 0.007751 0.000000 0.941178 0.058822 0.000000 0.004615 0.995385 0.000000 0.969698 0.030302 0.000000 0.999024 0.000976 0.000000 0.031769 0.968231 0.000000 0.000004 0.999996 0.000000 0.333333 0.333333 0.333333 0.998051 0.001949 0.000000 0.333333 0.333333 0.333333 0.000012 0.999988 0.000000 0.660786 0.339214 0.000000 0.012603 0.987397 0.000000 0.000026 0.999974 0.000000 0.000129 0.999871 0.000000 0.653200 0.346800 0.000000 0.999024 0.000976 0.000000 0.888891 0.111109 0.000000 0.014671 0.985329 0.000000 0.333333 0.333333 0.333333 0.048566 0.951434 0.000000 0.333333 0.333333 0.333333 0.998051 0.001949 0.000000 0.984616 0.015384 0.000000 0.000000 1.000000 0.000000 0.999024 0.000976 0.000000 0.969698 0.030302 0.000000 0.048566 0.951434 0.000000 0.333333 0.333333 0.333333 0.998051 0.001949 0.000000 0.015770 0.984230 0.000000 0.888891 0.111109 0.000000 0.000786 0.999214 0.000000 0.984616 0.015384 0.000000 0.969698 0.030302 0.000000 0.992249 0.007751 0.000000 0.333333 0.333333 0.333333 0.000002 0.999998 0.000000 0.941178 0.058822 0.000000 0.000059 0.999941 0.000000 0.000161 0.999839 0.000000 0.000000 1.000000 0.000000 0.028342 0.971658 0.000000 0.000098 0.999902 0.000000 0.123066 0.876934 0.000000 0.003181 0.996819 0.000000 0.000000 1.000000 0.000000 0.999969 0.000031 0.000000 0.169534 0.830466 0.000000 0.169534 0.830466 0.000000 0.000000 1.000000 0.000000 0.048566 0.951434 0.000000 0.932763 0.067237 0.000000 0.000000 1.000000 0.000000 0.888891 0.111109 0.000000 0.996109 0.003891 0.000000 0.999024 0.000976 0.000000 0.000108 0.999892 0.000000 0.765539 0.234461 0.000000 0.764131 0.235869 0.000000 0.101860 0.898140 0.000000 0.996109 0.003891 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.998051 0.001949 0.000000 0.999024 0.000976 0.000000 0.000000 1.000000 0.000000 0.984616 0.015384 0.000000 0.999512 0.000488 0.000000
CM019102.1_30964 3 0 0.999756 0.000244 0.000000 0.169534 0.830466 0.000000 0.997611 0.002389 0.000000 0.999996 0.000004 0.000000 0.999512 0.000488 0.000000 0.333333 0.333333 0.333333 0.333333 0.333333 0.333333 0.169534 0.830466 0.000000 0.999878 0.000122 0.000000 0.999878 0.000122 0.000000 0.999985 0.000015 0.000000 0.998051 0.001949 0.000000 0.999512 0.000488 0.000000 0.992249 0.007751 0.000000 0.999998 0.000002 0.000000 1.000000 0.000000 0.000000 0.941178 0.058822 0.000000 1.000000 0.000000 0.000000 0.999878 0.000122 0.000000 0.992249 0.007751 0.000000 0.996109 0.003891 0.000000 0.999024 0.000976 0.000000 0.969698 0.030302 0.000000 0.941178 0.058822 0.000000 0.333333 0.333333 0.333333 0.333333 0.333333 0.333333 1.000000 0.000000 0.000000 0.999939 0.000061 0.000000 1.000000 0.000000 0.000000 0.092626 0.907374 0.000000 0.999756 0.000244 0.000000 0.999024 0.000976 0.000000 0.333333 0.333333 0.333333 0.140893 0.859107 0.000000 0.020125 0.979875 0.000000 0.999701 0.000299 0.000000 0.992249 0.007751 0.000000 0.333333 0.333333 0.333333 0.999024 0.000976 0.000000 0.449468 0.550532 0.000000 0.999878 0.000122 0.000000 0.999996 0.000004 0.000000 0.333333 0.333333 0.333333 0.333333 0.333333 0.333333 0.000091 0.999909 0.000000 0.270465 0.729535 0.000000 0.076885 0.923115 0.000000 0.999878 0.000122 0.000000 0.999756 0.000244 0.000000 0.999756 0.000244 0.000000 0.969698 0.030302 0.000000 0.998051 0.001949 0.000000 0.941178 0.058822 0.000000 0.157587 0.842413 0.000000 0.969698 0.030302 0.000000 0.999512 0.000488 0.000000 0.289899 0.710101 0.000000 0.010308 0.989692 0.000000 0.333333 0.333333 0.333333 0.999024 0.000976 0.000000 0.333333 0.333333 0.333333 0.000003 0.999997 0.000000 0.998051 0.001949 0.000000 0.984616 0.015384 0.000000 0.999756 0.000244 0.000000 0.999756 0.000244 0.000000 0.867191 0.132809 0.000000 0.999939 0.000061 0.000000 0.888891 0.111109 0.000000 0.996109 0.003891 0.000000 0.333333 0.333333 0.333333 0.996109 0.003891 0.000000 0.333333 0.333333 0.333333 0.999024 0.000976 0.000000 0.984616 0.015384 0.000000 0.999878 0.000122 0.000000 0.999512 0.000488 0.000000 0.969698 0.030302 0.000000 0.048566 0.951434 0.000000 0.333333 0.333333 0.333333 0.999512 0.000488 0.000000 0.060202 0.939798 0.000000 0.888891 0.111109 0.000000 0.864413 0.135587 0.000000 0.996109 0.003891 0.000000 0.992249 0.007751 0.000000 0.996109 0.003891 0.000000 0.333333 0.333333 0.333333 0.928866 0.071134 0.000000 0.941178 0.058822 0.000000 0.146738 0.853262 0.000000 0.998051 0.001949 0.000000 0.984616 0.015384 0.000000 0.289899 0.710101 0.000000 0.000163 0.999837 0.000000 0.286207 0.713793 0.000000 0.012603 0.987397 0.000000 0.188440 0.811560 0.000000 0.999992 0.000008 0.000000 0.999024 0.000976 0.000000 0.620161 0.379839 0.000000 0.009795 0.990205 0.000000 0.996109 0.003891 0.000000 0.963118 0.036882 0.000000 0.014840 0.985160 0.000000 0.969698 0.030302 0.000000 0.998051 0.001949 0.000000 0.999024 0.000976 0.000000 1.000000 0.000000 0.000000 0.999969 0.000031 0.000000 0.999939 0.000061 0.000000 0.999878 0.000122 0.000000 0.996109 0.003891 0.000000 0.672647 0.327353 0.000000 0.999969 0.000031 0.000000 0.999939 0.000061 0.000000 0.999024 0.000976 0.000000 0.000041 0.999959 0.000000 0.996109 0.003891 0.000000 0.999756 0.000244 0.000000
CM019102.1_30969 1 3 0.765539 0.234461 0.000000 0.000002 0.999998 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.333333 0.333333 0.333333 0.000001 0.999982 0.000017 0.999939 0.000061 0.000000 0.028960 0.971040 0.000000 0.867191 0.132809 0.000000 0.014840 0.985160 0.000000 0.867191 0.132809 0.000000 0.005115 0.994885 0.000000 0.615213 0.384787 0.000000 1.000000 0.000000 0.000000 0.999925 0.000075 0.000000 0.003181 0.996819 0.000000 0.000000 1.000000 0.000000 0.000004 0.999996 0.000000 0.012603 0.987397 0.000000 0.996109 0.003891 0.000000 0.449468 0.550532 0.000000 0.941178 0.058822 0.000000 0.024890 0.975110 0.000000 0.333333 0.333333 0.333333 0.000000 0.999999 0.000001 0.419912 0.580088 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000041 0.999959 0.000000 0.000000 1.000000 0.000000 0.092626 0.907374 0.000000 0.333333 0.333333 0.333333 0.963118 0.036882 0.000000 0.928866 0.071134 0.000000 0.011908 0.988092 0.000000 0.999878 0.000122 0.000000 0.003181 0.996819 0.000000 0.000000 1.000000 0.000000 0.000024 0.999976 0.000000 0.000033 0.999967 0.000000 0.000050 0.999950 0.000000 0.333333 0.333333 0.333333 0.888891 0.111109 0.000000 0.004185 0.995815 0.000000 0.055784 0.944216 0.000000 0.294053 0.705947 0.000000 0.045531 0.954469 0.000000 0.000607 0.999393 0.000000 0.000001 0.999999 0.000000 0.000020 0.999980 0.000000 0.000000 1.000000 0.000000 0.000000 0.999997 0.000003 0.000000 1.000000 0.000000 0.941178 0.058822 0.000000 0.009795 0.990205 0.000000 0.996109 0.003891 0.000000 0.000000 1.000000 0.000000 0.333333 0.333333 0.333333 0.928307 0.071693 0.000000 0.333333 0.333333 0.333333 0.990515 0.009485 0.000000 0.999939 0.000061 0.000000 0.000003 0.999997 0.000000 0.999939 0.000061 0.000000 0.000000 1.000000 0.000000 0.999939 0.000061 0.000000 0.010171 0.989829 0.000000 0.000000 0.998406 0.001593 0.000254 0.999746 0.000000 0.333333 0.333333 0.333333 0.992249 0.007751 0.000000 0.333333 0.333333 0.333333 0.001945 0.998055 0.000000 0.000000 0.993658 0.006342 0.094645 0.905355 0.000000 0.998051 0.001949 0.000000 0.000000 1.000000 0.000000 0.996109 0.003891 0.000000 0.333333 0.333333 0.333333 0.620161 0.379839 0.000000 0.765539 0.234461 0.000000 0.888891 0.111109 0.000000 0.999969 0.000031 0.000000 0.000010 0.999990 0.000000 0.000016 0.999984 0.000000 0.999024 0.000976 0.000000 0.333333 0.333333 0.333333 0.000003 0.999997 0.000000 0.000000 0.999995 0.000005 0.000000 1.000000 0.000000 0.024890 0.975110 0.000000 0.000161 0.999839 0.000000 0.010308 0.989692 0.000000 0.963118 0.036882 0.000000 0.101860 0.898140 0.000000 0.969698 0.030302 0.000000 0.000000 1.000000 0.000000 0.000001 0.999999 0.000000 0.999878 0.000122 0.000000 0.000326 0.999674 0.000000 0.000326 0.999674 0.000000 0.999024 0.000976 0.000000 0.000004 0.999996 0.000000 0.640032 0.359968 0.000000 0.969698 0.030302 0.000000 0.000000 1.000000 0.000000 0.001945 0.998055 0.000000 0.999998 0.000002 0.000000 0.999756 0.000244 0.000000 0.990515 0.009485 0.000000 0.003836 0.996164 0.000000 0.001297 0.998703 0.000000 0.998051 0.001949 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.449468 0.550532 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.996109 0.003891 0.000000
CM019102.1_30981 2 0 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.997742 0.002258 0.000000 1.000000 0.000000 0.000000 0.928866 0.071134 0.000000 0.984616 0.015384 0.000000 0.992249 0.007751 0.000000 0.000003 0.999997 0.000000 0.026658 0.973342 0.000000 0.006476 0.993524 0.000000 1.000000 0.000000 0.000000 0.246871 0.753129 0.000000 1.000000 0.000000 0.000000 0.012866 0.987134 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.000017 0.999983 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.998792 0.001208 0.000000 0.999024 0.000976 0.000000 0.914174 0.085826 0.000000 0.999512 0.000488 0.000000 0.999024 0.000976 0.000000 0.984616 0.015384 0.000000 0.999756 0.000244 0.000000 0.000002 0.999998 0.000000 1.000000 0.000000 0.000000 0.994164 0.005836 0.000000 0.571173 0.428827 0.000000 0.999963 0.000037 0.000000 0.007556 0.992444 0.000000 0.888891 0.111109 0.000000 0.999927 0.000073 0.000000 1.000000 0.000000 0.000000 0.999419 0.000581 0.000000 1.000000 0.000000 0.000000 0.984616 0.015384 0.000000 0.999969 0.000031 0.000000 0.000013 0.999987 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.333333 0.333333 0.333333 0.969698 0.030302 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.999998 0.000002 0.000000 1.000000 0.000000 0.000000 0.999756 0.000244 0.000000 1.000000 0.000000 0.000000 0.996109 0.003891 0.000000 0.001494 0.998506 0.000000 0.934301 0.065699 0.000000 0.999997 0.000003 0.000000 0.992249 0.007751 0.000000 1.000000 0.000000 0.000000 0.333333 0.333333 0.333333 1.000000 0.000000 0.000000 0.999985 0.000015 0.000000 0.999985 0.000015 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.999981 0.000019 0.000000 0.996109 0.003891 0.000000 0.999996 0.000004 0.000000 0.333333 0.333333 0.333333 0.005181 0.994819 0.000000 0.333333 0.333333 0.333333 0.001757 0.998243 0.000000 0.992249 0.007751 0.000000 1.000000 0.000000 0.000000 0.620161 0.379839 0.000000 0.999999 0.000001 0.000000 0.999969 0.000031 0.000000 0.996109 0.003891 0.000000 0.000338 0.999662 0.000000 1.000000 0.000000 0.000000 0.000001 0.999999 0.000000 0.997611 0.002389 0.000000 0.999024 0.000976 0.000000 0.999969 0.000031 0.000000 0.999939 0.000061 0.000000 0.333333 0.333333 0.333333 1.000000 0.000000 0.000000 0.999992 0.000008 0.000000 1.000000 0.000000 0.000000 0.999939 0.000061 0.000000 0.999969 0.000031 0.000000 1.000000 0.000000 0.000000 0.999999 0.000001 0.000000 0.999756 0.000244 0.000000 0.999756 0.000244 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.992448 0.007552 0.000000 1.000000 0.000000 0.000000 0.995215 0.004785 0.000000 0.000002 0.999998 0.000000 1.000000 0.000000 0.000000 0.999756 0.000244 0.000000 0.999878 0.000122 0.000000 0.999261 0.000739 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.891133 0.108867 0.000000 0.999998 0.000002 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.999850 0.000150 0.000000 1.000000 0.000000 0.000000 0.999024 0.000976 0.000000 0.955158 0.044842 0.000000
CM019102.1_30983 2 0 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.002598 0.997402 0.000000 0.984616 0.015384 0.000000 0.996109 0.003891 0.000000 1.000000 0.000000 0.000000 0.158516 0.841484 0.000000 0.999985 0.000015 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.999969 0.000031 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.992249 0.007751 0.000000 1.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.999999 0.000001 0.000000 0.999512 0.000488 0.000000 1.000000 0.000000 0.000000 0.999512 0.000488 0.000000 0.998051 0.001949 0.000000 0.984616 0.015384 0.000000 0.999939 0.000061 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.999969 0.000031 0.000000 0.941178 0.058822 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.005181 0.994819 0.000000 0.984616 0.015384 0.000000 0.999992 0.000008 0.000000 0.999998 0.000002 0.000000 1.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.333333 0.333333 0.333333 0.969698 0.030302 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.765539 0.234461 0.000000 1.000000 0.000000 0.000000 0.992249 0.007751 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.996109 0.003891 0.000000 1.000000 0.000000 0.000000 0.333333 0.333333 0.333333 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 1.000000 0.000000 1.000000 0.000000 0.000000 0.999963 0.000037 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.999024 0.000976 0.000000 1.000000 0.000000 0.000000 0.333333 0.333333 0.333333 0.999998 0.000002 0.000000 0.333333 0.333333 0.333333 0.999980 0.000020 0.000000 0.996109 0.003891 0.000000 1.000000 0.000000 0.000000 0.999878 0.000122 0.000000 0.993959 0.006041 0.000000 0.999985 0.000015 0.000000 0.996109 0.003891 0.000000 0.999999 0.000001 0.000000 1.000000 0.000000 0.000000 0.969698 0.030302 0.000000 1.000000 0.000000 0.000000 0.998051 0.001949 0.000000 0.999878 0.000122 0.000000 0.999969 0.000031 0.000000 0.333333 0.333333 0.333333 1.000000 0.000000 0.000000 0.999998 0.000002 0.000000 0.999313 0.000687 0.000000 0.999985 0.000015 0.000000 0.999969 0.000031 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.999512 0.000488 0.000000 0.999512 0.000488 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.999999 0.000001 0.000000 1.000000 0.000000 0.000000 0.999996 0.000004 0.000000 0.414286 0.585714 0.000000 1.000000 0.000000 0.000000 0.999756 0.000244 0.000000 0.999985 0.000015 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.999969 0.000031 0.000000 0.999992 0.000008 0.000000 0.000681 0.999319 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.999512 0.000488 0.000000 0.999995 0.000005 0.000000
CM019102.1_30984 1 2 1.000000 0.000000 0.000000 0.000000 1.000000 0.000000 1.000000 0.000000 0.000000 0.999713 0.000287 0.000000 0.999985 0.000015 0.000000 0.984616 0.015384 0.000000 0.996109 0.003891 0.000000 1.000000 0.000000 0.000000 0.000000 1.000000 0.000000 0.002598 0.997402 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.999992 0.000008 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.941178 0.058822 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.990515 0.009485 0.000000 0.999024 0.000976 0.000000 0.996463 0.003537 0.000000 0.101860 0.898140 0.000000 0.000322 0.999678 0.000000 0.969698 0.030302 0.000000 0.999939 0.000061 0.000000 0.994593 0.005407 0.000000 1.000000 0.000000 0.000000 0.999428 0.000572 0.000000 0.399779 0.600221 0.000000 0.999998 0.000002 0.000000 0.999969 0.000031 0.000000 0.941178 0.058822 0.000000 0.000000 1.000000 0.000000 0.000253 0.999747 0.000000 0.000000 1.000000 0.000000 0.000128 0.999872 0.000000 0.969698 0.030302 0.000000 0.999985 0.000015 0.000000 0.999998 0.000002 0.000000 0.001523 0.998477 0.000000 1.000000 0.000000 0.000000 0.333333 0.333333 0.333333 0.984616 0.015384 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.881719 0.118281 0.000000 0.955158 0.044842 0.000000 0.942591 0.057409 0.000000 0.000025 0.999975 0.000000 0.249851 0.750149 0.000000 0.999969 0.000031 0.000000 1.000000 0.000000 0.000000 0.992249 0.007751 0.000000 0.000000 1.000000 0.000000 0.727041 0.272959 0.000000 0.493040 0.506960 0.000000 0.984616 0.015384 0.000000 0.000013 0.999987 0.000000 0.333333 0.333333 0.333333 0.977062 0.022938 0.000000 0.064072 0.935928 0.000000 0.002598 0.997402 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.999024 0.000976 0.000000 1.000000 0.000000 0.000000 0.333333 0.333333 0.333333 0.000008 0.999992 0.000000 0.333333 0.333333 0.333333 1.000000 0.000000 0.000000 0.998051 0.001949 0.000000 1.000000 0.000000 0.000000 0.000326 0.999674 0.000000 0.188440 0.811560 0.000000 0.999985 0.000015 0.000000 0.000161 0.999839 0.000000 0.999996 0.000004 0.000000 0.249851 0.750149 0.000000 0.969698 0.030302 0.000000 0.727041 0.272959 0.000000 0.999024 0.000976 0.000000 0.990231 0.009769 0.000000 0.999939 0.000061 0.000000 0.333333 0.333333 0.333333 0.777190 0.222810 0.000000 0.371164 0.628836 0.000000 0.998083 0.001917 0.000000 0.001301 0.998699 0.000000 0.999939 0.000061 0.000000 0.000773 0.999227 0.000000 1.000000 0.000000 0.000000 0.999512 0.000488 0.000000 0.000041 0.999959 0.000000 0.000400 0.999600 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.999756 0.000244 0.000000 0.999985 0.000015 0.000000 0.999998 0.000002 0.000000 0.000000 1.000000 0.000000 0.118420 0.881580 0.000000 0.000009 0.999991 0.000000 0.914174 0.085826 0.000000 0.999969 0.000031 0.000000 0.867191 0.132809 0.000000 0.000000 1.000000 0.000000 1.000000 0.000000 0.000000 0.399779 0.600221 0.000000 0.963560 0.036440 0.000000 0.999512 0.000488 0.000000 0.991127 0.008873 0.000000
CM019102.1_30985 3 1 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.928866 0.071134 0.000000 0.984616 0.015384 0.000000 0.048566 0.951434 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.867191 0.132809 0.000000 1.000000 0.000000 0.000000 0.000000 1.000000 0.000000 1.000000 0.000000 0.000000 0.999992 0.000008 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.024890 0.975110 0.000000 1.000000 0.000000 0.000000 0.999995 0.000005 0.000000 0.999395 0.000605 0.000000 0.999024 0.000976 0.000000 1.000000 0.000000 0.000000 0.999512 0.000488 0.000000 0.999024 0.000976 0.000000 0.984616 0.015384 0.000000 0.999939 0.000061 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.999701 0.000299 0.000000 0.999999 0.000001 0.000000 0.928866 0.071134 0.000000 0.047944 0.952056 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.984616 0.015384 0.000000 0.999996 0.000004 0.000000 0.995235 0.004765 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.333333 0.333333 0.333333 0.984616 0.015384 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.998804 0.001196 0.000000 1.000000 0.000000 0.000000 0.999939 0.000061 0.000000 1.000000 0.000000 0.000000 0.996109 0.003891 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.992249 0.007751 0.000000 0.399779 0.600221 0.000000 0.333333 0.333333 0.333333 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.999996 0.000004 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.999963 0.000037 0.000000 0.999024 0.000976 0.000000 0.997611 0.002389 0.000000 0.333333 0.333333 0.333333 0.999999 0.000001 0.000000 0.333333 0.333333 0.333333 1.000000 0.000000 0.000000 0.999024 0.000976 0.000000 1.000000 0.000000 0.000000 0.999969 0.000031 0.000000 1.000000 0.000000 0.000000 0.999939 0.000061 0.000000 0.996109 0.003891 0.000000 0.000178 0.999822 0.000000 1.000000 0.000000 0.000000 0.969698 0.030302 0.000000 0.999963 0.000037 0.000000 0.998051 0.001949 0.000000 0.999969 0.000031 0.000000 0.999969 0.000031 0.000000 0.333333 0.333333 0.333333 1.000000 0.000000 0.000000 0.999998 0.000002 0.000000 1.000000 0.000000 0.000000 0.999992 0.000008 0.000000 0.999969 0.000031 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.999756 0.000244 0.000000 0.999756 0.000244 0.000000 1.000000 0.000000 0.000000 0.999981 0.000019 0.000000 1.000000 0.000000 0.000000 0.076885 0.923115 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.000163 0.999837 0.000000 0.010308 0.989692 0.000000 0.999999 0.000001 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.942420 0.057580 0.000000 0.999996 0.000004 0.000000 1.000000 0.000000 0.000000 1.000000 0.000000 0.000000 0.999701 0.000299 0.000000 1.000000 0.000000 0.000000 0.000501 0.999499 0.000000 1.000000 0.000000 0.000000
Here is the ngsLD
command that I've been running:
( NUM_SITES=$(cat {input.pos} | wc -l) &&
ngsLD --geno {input.gls} \
--pos {input.pos} \
--n_ind 120 \
--n_sites $NUM_SITES \
--probs \
--n_threads {threads} \
--max_kb_dist 25 | gzip --best > {output} ) 2> {log}
where the first few lines of the position file looks like:
CM019102.1 30951
CM019102.1 30958
CM019102.1 30963
CM019102.1 30964
CM019102.1 30969
CM019102.1 30981
CM019102.1 30983
CM019102.1 30984
CM019102.1 30985
CM019102.1 30989
For reference, here is what the log file looks like for a successful run:
==> Input Arguments:
geno: ../results/population_structure/angsd_gl/full/CM019101.1/CM019101.1_genolike_allSamples.beagle.gz
probs: true
log_scale: false
n_ind: 120
n_sites: 2487093
pos: ../results/population_structure/angsd_gl/full/CM019101.1/CM019101.1_angsdGL.pos (WITHOUT header)
max_kb_dist (kb): 25
max_snp_dist: 0
min_maf: 0.001000
ignore_miss_data: false
call_geno: false
N_thresh: 0.000000
call_thresh: 0.000000
rnd_sample: 1.000000
seed: 1609882280
extend_out: false
out: (null) (WITHOUT header)
n_threads: 16
verbose: 1
version: 1.1.1 (Dec 22 2020 @ 00:24:22)
==> GZIP input file (not BINARY)
> Reading data from file...
> Header found! Skipping line...
==> Calculating MAF for all sites...
==> Getting sites coordinates
==> Launching threads...
==> Waiting for all threads to finish...
==> Freeing memory...
Done!
And here is what the log file looks like for a failed run that hangs indefinitely:
==> Input Arguments:
geno: ../results/population_structure/angsd/CM019105.1/CM019105.1_genolike_allSamples.beagle.gz
probs: true
log_scale: false
n_ind: 120
n_sites: 2036161
pos: ../results/population_structure/angsd/CM019105.1/CM019105.1_angsdGL.pos (WITHOUT header)
max_kb_dist (kb): 25
max_snp_dist: 0
min_maf: 0.001000
ignore_miss_data: false
call_geno: false
N_thresh: 0.000000
call_thresh: 0.000000
rnd_sample: 1.000000
seed: 1610049659
extend_out: false
out: (null) (WITHOUT header)
n_threads: 16
verbose: 1
version: 1.1.1 (Dec 22 2020 @ 00:24:22)
==> GZIP input file (not BINARY)
> Reading data from file...
Finally, I should mention that the commands are being run inside of a singularity container (see HERE) into which I've installed ngsLD
.
Also, it doesn't appear to be a problem with particular chromosomes. If I restart the failed runs, they will eventually succeed so I'm not sure what is causing the occasional hanging.
Thanks!
Hi James,
and sorry for the late reply.
This is indeed strange, specially since you say that these failed runs eventually finish, if you launch them again, right?
That makes me think that it might not be ngsLD
related.
Can you check if you have enough memory to launch the 16 chromosome separately? If you don't, you might be using SWAP memory that is considerably slower.
This is happening to me as well, in a slightly different way. I have over 100 datasets with the same number of SNPs in each (from different species and different chromosomes -- I use the same size of subsample of SNPs from each), and all of them run nicely except for one, which hangs indefinitely. In my case, it is one particular dataset that hangs every time.
I know this might be difficult to troubleshoot, I could potentially send some data if it would help.
Thanks,
-Teresa
It's possible that my issue is caused by an input beagle file that is composed of mostly missing data (the vast majority of genotype likelihoods were 0.33333). Could it be getting stuck because of that?
Is it stuck (as in, no output being generated), or is it just slow?
It seems to get stuck -- I haven't let it run for more than 20 hours before killing it, but this is on a dataset of a size that typically takes about 10 minutes to run, so I don't think it's doing anything.
I was able to get it to complete when I add the options --min_maf 0 and --ignore_miss_data. Obviously this particular dataset is pretty worthless, it's mostly missing, but hopefully it's helpful to reveal a particular circumstance that can cause ngsLD to hang up and never finish!