|
Commit message:
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bayescan/ commit 33637968e1e32c02d7765a6701e930a0ea0dd903 |
|
added:
bayescan.xml test-data/pilot.txt test-data/result.out test-data/test_binary_AFLP.txt test-data/verif.txt |
| b |
| diff -r 000000000000 -r 72f3a333f155 bayescan.xml --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/bayescan.xml Fri Mar 17 15:20:38 2017 -0400 |
| [ |
| @@ -0,0 +1,147 @@ +<tool id="BayeScan" name="BayeScan" version="2.1"> + <description>Detecting natural selection from population-based genetic data</description> + <requirements> + <requirement type="package" version="2.0.1">bayescan</requirement> + </requirements> + <command> + + <![CDATA[ + mkdir 'output_dir'; + + bayescan2 '$input' + -od output_dir + #if '$loci_file.loci' == "1" + -d '$loci_file.input_loci' + #end if + '$snp_genotypes_matrix' + '$fstats' + '$pilot_runs' + '$allele_frequency' + -o bayescan + -n '$sample_size' + -thin '$thinning_interval' + -nbp '$num_pilot_runs' + -pilot '$length_pilot_run' + -burn '$burn' + -pr_odds '$prior_odds' + -lb_fis '$lower_prior' + -hb_fis '$higher_prior' + -aflp_pc '$threshold' > '$output' + ]]> + + </command> + + <inputs> + + <param name="input" type="data" format="tabular,txt" label="Input genotype data file" help="must be space/tab delimitted plain text file" /> + <conditional name="loci_file"> + <param name="loci" type="select" label="Discard loci?" help="" > + <option value="0">No</option> + <option value="1">Yes</option> + </param> + <when value="0"></when> + <when value="1"> + <param name="input_loci" type="data" format="tabular,txt" label="Discard loci file" help="Optional input file containing list of loci to discard" /> + </when> + </conditional> + <param name="snp_genotypes_matrix" type="boolean" checked="false" truevalue="-fstat" falsevalue="" label="SNP genotypes matrix data" help="Use SNP genotypes matrix"/> + + <param name="fstats" type="boolean" checked="false" truevalue="-snp" falsevalue="" label="Only estimate F-stats (no selection)"/> + + <param name="sample_size" type="integer" value="5000" label="Number of outputted iterations"/> + <param name="thinning_interval" type="integer" value="10" label="Thinning interval size"/> + <param name="num_pilot_runs" type="integer" value="20" label="Number of pilot runs" /> + <param name="length_pilot_run" type="integer" value="5000" label="Length of pilot runs" /> + <param name="burn" type="integer" value="50000" label="Additional burn-in length" help="Nlength of an interval between two draws in MCMC"/> + + <param name="prior_odds" type="integer" value="10" size="3" label="Prior odds for the neutral model"/> + <param name="lower_prior" type="float" value="0.0" label="Lower bound for uniform prior on Fis (dominant data)"/> + <param name="higher_prior" type="float" value="1.0" label="Higher bound for uniform prior on Fis (dominant data)"/> + <param name="threshold" type="float" value="0.1" label="Threshold for the recessive genotype as a fraction of maximum band intensity"/> + + <param name="pilot_runs" type="boolean" checked="false" truevalue="-out_pilot" falsevalue="" label="Optional output file for pilot runs"/> + <param name="allele_frequency" type="boolean" checked="false" truevalue="-out_freq" falsevalue="" label="Optional output file for allele frequencies"/> + + + </inputs> + + <outputs> + <data format="txt" name="output"> + <discover_datasets pattern="__designation_and_ext__" directory="output_dir" visible="true" /> + </data> + </outputs> + + <tests> + <test> + <param name="input" value="test_binary_AFLP.txt" /> + <param name="loci" value="0" /> + <param name="snp_genotypes_matrix" value="true" /> + <param name="fstats" value="true"/> + + <param name="sample_size" value="5000" /> + <param name="thinning_interval" value="10" /> + <param name="num_pilot_runs" value="20" /> + <param name="length_pilot_run" value="5000" /> + <param name="burn" value="50000" /> + + <param name="prior_odds" value="10"/> + <param name="lower_prior" value="0.0" /> + <param name="higher_prior" value="1.0"/> + <param name="threshold" value="0.1"/> + + <param name="pilot_runs" value="true"/> + <param name="allele_frequency" value="true" /> + + <output name="output" file="result.out" ftype="txt"> + <discovered_dataset designation="bayescan" ftype="sel"> + <assert_contents> + <has_text text="logL Fis1 Fis2 Fis3 Fis4 Fis5 Fis6 Fis7 Fis8 Fis9 Fis10 Fst1 Fst2 Fst3 Fst4 Fst5 Fst6 Fst7 Fst8 Fst9 Fst10" /> + </assert_contents> + </discovered_dataset> + <discovered_dataset designation="bayescan_Verif" ftype="txt" value="verif.txt" /> + <discovered_dataset designation="bayescan_AccRte" ftype="txt"> + <assert_contents> + <has_line_matching expression="alpha beta ances freq a_p.*" /> + </assert_contents> + </discovered_dataset> + <discovered_dataset designation="bayescan_prop" ftype="txt" value="pilot.txt" /> + <discovered_dataset designation="bayescan_freq" ftype="txt"> + <assert_contents> + <has_text text="locus1 locus2 locus3 locus4 locus5 locus6 locus7 locus8 locus9 locus10" /> + </assert_contents> + </discovered_dataset> + </output> + </test> + </tests> + + <help><![CDATA[ + +**What it does** + +This program, BayeScan aims at identifying candidate loci under natural selection from genetic data, using differences in allele frequencies between populations. + +`BayeScan`_ is based on the multinomial-Dirichlet model. + +One of the simplest possible scenarios covered consists of an island model in which subpopulation allele frequencies are correlated through a common migrant gene pool from which they differ in varying degrees. The difference in allele frequency between this common gene pool and each subpopulation is measured by a subpopulation specific FST coefficient. + +Therefore, this formulation can consider realistic ecological scenarios where the effective size and the immigration rate may differ among subpopulations. + +.. _Bayescan: http://cmpg.unibe.ch/software/BayeScan/ + +.. class:: infomark + +**Input file** + +BayeScan uses its own input file formats, which depend on the type of data used. All input files are simply in text format. + +Read the `manual`_ to create the input file. + +.. _manual: http://cmpg.unibe.ch/software/BayeScan/files/BayeScan2.1_manual.pdf + +]]> + + </help> + <citations> + <citation type="doi">10.1534/genetics.108.092221</citation> + </citations> +</tool> |
| b |
| diff -r 000000000000 -r 72f3a333f155 test-data/pilot.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/pilot.txt Fri Mar 17 15:20:38 2017 -0400 |
| b |
| @@ -0,0 +1,106 @@ + +Variance of alpha proposal: 0.0261 +Variance of beta proposal: 0.583 +Variance of a proposal: 0.00522 +Variance of allele freq proposal: 7.67 +Variance of ancestral allele freq proposal: 0.293 +Mean and variance alpha_1 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_2 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_3 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_4 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_5 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_6 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_7 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_8 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_9 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_10 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_11 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_12 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_13 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_14 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_15 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_16 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_17 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_18 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_19 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_20 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_21 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_22 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_23 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_24 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_25 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_26 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_27 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_28 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_29 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_30 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_31 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_32 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_33 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_34 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_35 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_36 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_37 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_38 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_39 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_40 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_41 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_42 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_43 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_44 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_45 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_46 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_47 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_48 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_49 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_50 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_51 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_52 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_53 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_54 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_55 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_56 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_57 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_58 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_59 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_60 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_61 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_62 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_63 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_64 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_65 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_66 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_67 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_68 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_69 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_70 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_71 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_72 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_73 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_74 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_75 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_76 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_77 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_78 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_79 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_80 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_81 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_82 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_83 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_84 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_85 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_86 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_87 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_88 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_89 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_90 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_91 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_92 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_93 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_94 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_95 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_96 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_97 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_98 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_99 RJ proposal: 0.00 , 0.00 +Mean and variance alpha_100 RJ proposal: 0.00 , 0.00 |
| b |
| diff -r 000000000000 -r 72f3a333f155 test-data/result.out --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/result.out Fri Mar 17 15:20:38 2017 -0400 |
| b |
| @@ -0,0 +1,4 @@ +Pilot runs... +1% 2% 3% 4% 5% 6% 7% 8% 9% 10% 11% 12% 13% 14% 15% 16% 17% 18% 19% 20% 21% 22% 23% 24% 25% 26% 27% 28% 29% 30% 31% 32% 33% 34% 35% 36% 37% 38% 39% 40% 41% 42% 43% 44% 45% 46% 47% 48% 49% 50% 51% 52% 53% 54% 55% 56% 57% 58% 59% 60% 61% 62% 63% 64% 65% 66% 67% 68% 69% 70% 71% 72% 73% 74% 75% 76% 77% 78% 79% 80% 81% 82% 83% 84% 85% 86% 87% 88% 89% 90% 91% 92% 93% 94% 95% 96% 97% 98% 99% +Calculation... +1% 2% 3% 4% 5% 6% 7% 8% 9% 10% 11% 12% 13% 14% 15% 16% 17% 18% 19% 20% 21% 22% 23% 24% 25% 26% 27% 28% 29% 30% 31% 32% 33% 34% 35% 36% 37% 38% 39% 40% 41% 42% 43% 44% 45% 46% 47% 48% 49% 50% 51% 52% 53% 54% 55% 56% 57% 58% 59% 60% 61% 62% 63% 64% 65% 66% 67% 68% 69% 70% 71% 72% 73% 74% 75% 76% 77% 78% 79% 80% 81% 82% 83% 84% 85% 86% 87% 88% 89% 90% 91% 92% 93% 94% 95% 96% 97% 98% 99% |
| b |
| diff -r 000000000000 -r 72f3a333f155 test-data/test_binary_AFLP.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test_binary_AFLP.txt Fri Mar 17 15:20:38 2017 -0400 |
| [ |
| b'@@ -0,0 +1,1024 @@\n+[loci]=100\n+\n+[populations]=10\n+\n+[pop]=1\n+ 1 30 26\n+ 2 30 2\n+ 3 30 29\n+ 4 30 0\n+ 5 30 30\n+ 6 30 20\n+ 7 30 0\n+ 8 30 25\n+ 9 30 29\n+ 10 30 30\n+ 11 30 0\n+ 12 30 0\n+ 13 30 30\n+ 14 30 30\n+ 15 30 26\n+ 16 30 0\n+ 17 30 30\n+ 18 30 29\n+ 19 30 27\n+ 20 30 30\n+ 21 30 0\n+ 22 30 0\n+ 23 30 30\n+ 24 30 1\n+ 25 30 28\n+ 26 30 0\n+ 27 30 29\n+ 28 30 0\n+ 29 30 21\n+ 30 30 27\n+ 31 30 29\n+ 32 30 0\n+ 33 30 1\n+ 34 30 2\n+ 35 30 13\n+ 36 30 23\n+ 37 30 6\n+ 38 30 28\n+ 39 30 30\n+ 40 30 30\n+ 41 30 5\n+ 42 30 0\n+ 43 30 0\n+ 44 30 29\n+ 45 30 25\n+ 46 30 1\n+ 47 30 22\n+ 48 30 21\n+ 49 30 24\n+ 50 30 17\n+ 51 30 0\n+ 52 30 13\n+ 53 30 29\n+ 54 30 30\n+ 55 30 13\n+ 56 30 27\n+ 57 30 30\n+ 58 30 0\n+ 59 30 0\n+ 60 30 0\n+ 61 30 1\n+ 62 30 29\n+ 63 30 27\n+ 64 30 18\n+ 65 30 29\n+ 66 30 30\n+ 67 30 29\n+ 68 30 30\n+ 69 30 25\n+ 70 30 0\n+ 71 30 9\n+ 72 30 16\n+ 73 30 25\n+ 74 30 0\n+ 75 30 30\n+ 76 30 30\n+ 77 30 25\n+ 78 30 5\n+ 79 30 30\n+ 80 30 30\n+ 81 30 2\n+ 82 30 30\n+ 83 30 30\n+ 84 30 19\n+ 85 30 0\n+ 86 30 4\n+ 87 30 30\n+ 88 30 0\n+ 89 30 29\n+ 90 30 30\n+ 91 30 25\n+ 92 30 30\n+ 93 30 21\n+ 94 30 0\n+ 95 30 28\n+ 96 30 30\n+ 97 30 0\n+ 98 30 30\n+ 99 30 0\n+100 30 30\n+\n+[pop]=2\n+ 1 30 13\n+ 2 30 3\n+ 3 30 26\n+ 4 30 1\n+ 5 30 30\n+ 6 30 23\n+ 7 30 0\n+ 8 30 30\n+ 9 30 27\n+ 10 30 30\n+ 11 30 10\n+ 12 30 0\n+ 13 30 12\n+ 14 30 30\n+ 15 30 10\n+ 16 30 2\n+ 17 30 30\n+ 18 30 30\n+ 19 30 20\n+ 20 30 30\n+ 21 30 0\n+ 22 30 4\n+ 23 30 30\n+ 24 30 14\n+ 25 30 30\n+ 26 30 0\n+ 27 30 29\n+ 28 30 4\n+ 29 30 1\n+ 30 30 30\n+ 31 30 30\n+ 32 30 1\n+ 33 30 0\n+ 34 30 0\n+ 35 30 29\n+ 36 30 28\n+ 37 30 22\n+ 38 30 30\n+ 39 30 29\n+ 40 30 29\n+ 41 30 15\n+ 42 30 9\n+ 43 30 0\n+ 44 30 29\n+ 45 30 11\n+ 46 30 11\n+ 47 30 14\n+ 48 30 20\n+ 49 30 22\n+ 50 30 30\n+ 51 30 0\n+ 52 30 19\n+ 53 30 30\n+ 54 30 30\n+ 55 30 5\n+ 56 30 27\n+ 57 30 30\n+ 58 30 0\n+ 59 30 1\n+ 60 30 5\n+ 61 30 0\n+ 62 30 30\n+ 63 30 11\n+ 64 30 0\n+ 65 30 30\n+ 66 30 30\n+ 67 30 27\n+ 68 30 30\n+ 69 30 30\n+ 70 30 2\n+ 71 30 0\n+ 72 30 11\n+ 73 30 9\n+ 74 30 0\n+ 75 30 30\n+ 76 30 30\n+ 77 30 15\n+ 78 30 26\n+ 79 30 29\n+ 80 30 30\n+ 81 30 5\n+ 82 30 30\n+ 83 30 30\n+ 84 30 20\n+ 85 30 0\n+ 86 30 14\n+ 87 30 29\n+ 88 30 0\n+ 89 30 26\n+ 90 30 22\n+ 91 30 19\n+ 92 30 30\n+ 93 30 18\n+ 94 30 0\n+ 95 30 12\n+ 96 30 30\n+ 97 30 16\n+ 98 30 30\n+ 99 30 0\n+100 30 30\n+\n+[pop]=3\n+ 1 30 10\n+ 2 30 4\n+ 3 30 17\n+ 4 30 2\n+ 5 30 30\n+ 6 30 17\n+ 7 30 0\n+ 8 30 30\n+ 9 30 18\n+ 10 30 30\n+ 11 30 26\n+ 12 30 0\n+ 13 30 8\n+ 14 30 30\n+ 15 30 10\n+ 16 30 5\n+ 17 30 30\n+ 18 30 30\n+ 19 30 24\n+ 20 30 29\n+ 21 30 0\n+ 22 30 16\n+ 23 30 30\n+ 24 30 26\n+ 25 30 30\n+ 26 30 0\n+ 27 30 13\n+ 28 30 0\n+ 29 30 3\n+ 30 30 29\n+ 31 30 29\n+ 32 30 0\n+ 33 30 0\n+ 34 30 4\n+ 35 30 30\n+ 36 30 30\n+ 37 30 23\n+ 38 30 29\n+ 39 30 19\n+ 40 30 21\n+ 41 30 8\n+ 42 30 8\n+ 43 30 1\n+ 44 30 14\n+ 45 30 27\n+ 46 30 7\n+ 47 30 10\n+ 48 30 22\n+ 49 30 27\n+ 50 30 16\n+ 51 30 6\n+ 52 30 27\n+ 53 30 30\n+ 54 30 30\n+ 55 30 25\n+ 56 30 22\n+ 57 30 30\n+ 58 30 0\n+ 59 30 0\n+ 60 30 0\n+ 61 30 0\n+ 62 30 30\n+ 63 30 27\n+ 64 30 10\n+ 65 30 27\n+ 66 30 30\n+ 67 30 16\n+ 68 30 30\n+ 69 30 27\n+ 70 30 0\n+ 71 30 3\n+ 72 30 15\n+ 73 30 19\n+ 74 30 23\n+ 75 30 30\n+ 76 30 30\n+ 77 30 27\n+ 78 30 10\n+ 79 30 30\n+ 80 30 '..b'17 30 30\n+ 18 30 30\n+ 19 30 29\n+ 20 30 30\n+ 21 30 0\n+ 22 30 26\n+ 23 30 30\n+ 24 30 24\n+ 25 30 27\n+ 26 30 0\n+ 27 30 30\n+ 28 30 0\n+ 29 30 27\n+ 30 30 22\n+ 31 30 30\n+ 32 30 0\n+ 33 30 0\n+ 34 30 0\n+ 35 30 30\n+ 36 30 30\n+ 37 30 26\n+ 38 30 28\n+ 39 30 30\n+ 40 30 30\n+ 41 30 10\n+ 42 30 4\n+ 43 30 0\n+ 44 30 29\n+ 45 30 9\n+ 46 30 5\n+ 47 30 9\n+ 48 30 24\n+ 49 30 21\n+ 50 30 14\n+ 51 30 0\n+ 52 30 11\n+ 53 30 30\n+ 54 30 30\n+ 55 30 22\n+ 56 30 6\n+ 57 30 30\n+ 58 30 0\n+ 59 30 0\n+ 60 30 0\n+ 61 30 10\n+ 62 30 30\n+ 63 30 28\n+ 64 30 0\n+ 65 30 29\n+ 66 30 30\n+ 67 30 25\n+ 68 30 30\n+ 69 30 29\n+ 70 30 0\n+ 71 30 14\n+ 72 30 18\n+ 73 30 18\n+ 74 30 0\n+ 75 30 30\n+ 76 30 30\n+ 77 30 17\n+ 78 30 30\n+ 79 30 30\n+ 80 30 30\n+ 81 30 6\n+ 82 30 30\n+ 83 30 30\n+ 84 30 17\n+ 85 30 0\n+ 86 30 1\n+ 87 30 30\n+ 88 30 0\n+ 89 30 20\n+ 90 30 16\n+ 91 30 20\n+ 92 30 30\n+ 93 30 9\n+ 94 30 2\n+ 95 30 22\n+ 96 30 30\n+ 97 30 1\n+ 98 30 27\n+ 99 30 0\n+100 30 30\n+\n+[pop]=9\n+ 1 30 3\n+ 2 30 28\n+ 3 30 30\n+ 4 30 0\n+ 5 30 30\n+ 6 30 6\n+ 7 30 0\n+ 8 30 30\n+ 9 30 23\n+ 10 30 30\n+ 11 30 1\n+ 12 30 0\n+ 13 30 4\n+ 14 30 30\n+ 15 30 18\n+ 16 30 0\n+ 17 30 29\n+ 18 30 30\n+ 19 30 4\n+ 20 30 16\n+ 21 30 0\n+ 22 30 21\n+ 23 30 30\n+ 24 30 17\n+ 25 30 22\n+ 26 30 0\n+ 27 30 21\n+ 28 30 0\n+ 29 30 28\n+ 30 30 29\n+ 31 30 30\n+ 32 30 0\n+ 33 30 0\n+ 34 30 0\n+ 35 30 28\n+ 36 30 30\n+ 37 30 20\n+ 38 30 27\n+ 39 30 30\n+ 40 30 30\n+ 41 30 0\n+ 42 30 1\n+ 43 30 0\n+ 44 30 26\n+ 45 30 5\n+ 46 30 6\n+ 47 30 12\n+ 48 30 28\n+ 49 30 17\n+ 50 30 6\n+ 51 30 0\n+ 52 30 29\n+ 53 30 30\n+ 54 30 30\n+ 55 30 18\n+ 56 30 28\n+ 57 30 30\n+ 58 30 0\n+ 59 30 0\n+ 60 30 17\n+ 61 30 9\n+ 62 30 29\n+ 63 30 9\n+ 64 30 14\n+ 65 30 30\n+ 66 30 30\n+ 67 30 27\n+ 68 30 30\n+ 69 30 27\n+ 70 30 0\n+ 71 30 12\n+ 72 30 14\n+ 73 30 14\n+ 74 30 3\n+ 75 30 26\n+ 76 30 30\n+ 77 30 15\n+ 78 30 29\n+ 79 30 29\n+ 80 30 17\n+ 81 30 0\n+ 82 30 30\n+ 83 30 30\n+ 84 30 18\n+ 85 30 0\n+ 86 30 8\n+ 87 30 30\n+ 88 30 2\n+ 89 30 29\n+ 90 30 23\n+ 91 30 12\n+ 92 30 30\n+ 93 30 29\n+ 94 30 0\n+ 95 30 23\n+ 96 30 30\n+ 97 30 4\n+ 98 30 30\n+ 99 30 0\n+100 30 30\n+\n+[pop]=10\n+ 1 30 17\n+ 2 30 10\n+ 3 30 29\n+ 4 30 0\n+ 5 30 30\n+ 6 30 28\n+ 7 30 0\n+ 8 30 30\n+ 9 30 25\n+ 10 30 30\n+ 11 30 11\n+ 12 30 0\n+ 13 30 1\n+ 14 30 30\n+ 15 30 0\n+ 16 30 0\n+ 17 30 29\n+ 18 30 30\n+ 19 30 23\n+ 20 30 19\n+ 21 30 0\n+ 22 30 0\n+ 23 30 30\n+ 24 30 29\n+ 25 30 24\n+ 26 30 0\n+ 27 30 12\n+ 28 30 0\n+ 29 30 1\n+ 30 30 23\n+ 31 30 30\n+ 32 30 10\n+ 33 30 0\n+ 34 30 2\n+ 35 30 30\n+ 36 30 30\n+ 37 30 8\n+ 38 30 24\n+ 39 30 29\n+ 40 30 30\n+ 41 30 19\n+ 42 30 17\n+ 43 30 0\n+ 44 30 4\n+ 45 30 27\n+ 46 30 0\n+ 47 30 8\n+ 48 30 9\n+ 49 30 7\n+ 50 30 18\n+ 51 30 0\n+ 52 30 28\n+ 53 30 30\n+ 54 30 30\n+ 55 30 9\n+ 56 30 17\n+ 57 30 30\n+ 58 30 0\n+ 59 30 0\n+ 60 30 5\n+ 61 30 4\n+ 62 30 30\n+ 63 30 22\n+ 64 30 0\n+ 65 30 30\n+ 66 30 30\n+ 67 30 26\n+ 68 30 30\n+ 69 30 24\n+ 70 30 0\n+ 71 30 11\n+ 72 30 9\n+ 73 30 10\n+ 74 30 0\n+ 75 30 30\n+ 76 30 30\n+ 77 30 27\n+ 78 30 19\n+ 79 30 30\n+ 80 30 30\n+ 81 30 5\n+ 82 30 30\n+ 83 30 30\n+ 84 30 10\n+ 85 30 0\n+ 86 30 23\n+ 87 30 30\n+ 88 30 0\n+ 89 30 22\n+ 90 30 9\n+ 91 30 26\n+ 92 30 30\n+ 93 30 28\n+ 94 30 0\n+ 95 30 26\n+ 96 30 30\n+ 97 30 0\n+ 98 30 30\n+ 99 30 0\n+100 30 30\n+\n' |
| b |
| diff -r 000000000000 -r 72f3a333f155 test-data/verif.txt --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/verif.txt Fri Mar 17 15:20:38 2017 -0400 |
| b |
| b'@@ -0,0 +1,1025 @@\n+Summary of parameters and input files.\n+Please check that all is correct while calculation is starting...\n+\n+There are 100 loci.\n+\n+There are 10 populations.\n+\n+Burn in: 50000\n+Thining interval: 10\n+Sample size: 5000\n+Resulting total number of iterations: 100000\n+Nb of pilot runs: 20\n+Length of each pilot run: 5000\n+\n+Allele counts:\n+Pop. 1 locus 1 : 30 26\n+Pop. 1 locus 2 : 30 2\n+Pop. 1 locus 3 : 30 29\n+Pop. 1 locus 4 : 30 0\n+Pop. 1 locus 5 : 30 30\n+Pop. 1 locus 6 : 30 20\n+Pop. 1 locus 7 : 30 0\n+Pop. 1 locus 8 : 30 25\n+Pop. 1 locus 9 : 30 29\n+Pop. 1 locus 10 : 30 30\n+Pop. 1 locus 11 : 30 0\n+Pop. 1 locus 12 : 30 0\n+Pop. 1 locus 13 : 30 30\n+Pop. 1 locus 14 : 30 30\n+Pop. 1 locus 15 : 30 26\n+Pop. 1 locus 16 : 30 0\n+Pop. 1 locus 17 : 30 30\n+Pop. 1 locus 18 : 30 29\n+Pop. 1 locus 19 : 30 27\n+Pop. 1 locus 20 : 30 30\n+Pop. 1 locus 21 : 30 0\n+Pop. 1 locus 22 : 30 0\n+Pop. 1 locus 23 : 30 30\n+Pop. 1 locus 24 : 30 1\n+Pop. 1 locus 25 : 30 28\n+Pop. 1 locus 26 : 30 0\n+Pop. 1 locus 27 : 30 29\n+Pop. 1 locus 28 : 30 0\n+Pop. 1 locus 29 : 30 21\n+Pop. 1 locus 30 : 30 27\n+Pop. 1 locus 31 : 30 29\n+Pop. 1 locus 32 : 30 0\n+Pop. 1 locus 33 : 30 1\n+Pop. 1 locus 34 : 30 2\n+Pop. 1 locus 35 : 30 13\n+Pop. 1 locus 36 : 30 23\n+Pop. 1 locus 37 : 30 6\n+Pop. 1 locus 38 : 30 28\n+Pop. 1 locus 39 : 30 30\n+Pop. 1 locus 40 : 30 30\n+Pop. 1 locus 41 : 30 5\n+Pop. 1 locus 42 : 30 0\n+Pop. 1 locus 43 : 30 0\n+Pop. 1 locus 44 : 30 29\n+Pop. 1 locus 45 : 30 25\n+Pop. 1 locus 46 : 30 1\n+Pop. 1 locus 47 : 30 22\n+Pop. 1 locus 48 : 30 21\n+Pop. 1 locus 49 : 30 24\n+Pop. 1 locus 50 : 30 17\n+Pop. 1 locus 51 : 30 0\n+Pop. 1 locus 52 : 30 13\n+Pop. 1 locus 53 : 30 29\n+Pop. 1 locus 54 : 30 30\n+Pop. 1 locus 55 : 30 13\n+Pop. 1 locus 56 : 30 27\n+Pop. 1 locus 57 : 30 30\n+Pop. 1 locus 58 : 30 0\n+Pop. 1 locus 59 : 30 0\n+Pop. 1 locus 60 : 30 0\n+Pop. 1 locus 61 : 30 1\n+Pop. 1 locus 62 : 30 29\n+Pop. 1 locus 63 : 30 27\n+Pop. 1 locus 64 : 30 18\n+Pop. 1 locus 65 : 30 29\n+Pop. 1 locus 66 : 30 30\n+Pop. 1 locus 67 : 30 29\n+Pop. 1 locus 68 : 30 30\n+Pop. 1 locus 69 : 30 25\n+Pop. 1 locus 70 : 30 0\n+Pop. 1 locus 71 : 30 9\n+Pop. 1 locus 72 : 30 16\n+Pop. 1 locus 73 : 30 25\n+Pop. 1 locus 74 : 30 0\n+Pop. 1 locus 75 : 30 30\n+Pop. 1 locus 76 : 30 30\n+Pop. 1 locus 77 : 30 25\n+Pop. 1 locus 78 : 30 5\n+Pop. 1 locus 79 : 30 30\n+Pop. 1 locus 80 : 30 30\n+Pop. 1 locus 81 : 30 2\n+Pop. 1 locus 82 : 30 30\n+Pop. 1 locus 83 : 30 30\n+Pop. 1 locus 84 : 30 19\n+Pop. 1 locus 85 : 30 0\n+Pop. 1 locus 86 : 30 4\n+Pop. 1 locus 87 : 30 30\n+Pop. 1 locus 88 : 30 0\n+Pop. 1 locus 89 : 30 29\n+Pop. 1 locus 90 : 30 30\n+Pop. 1 locus 91 : 30 25\n+Pop. 1 locus 92 : 30 30\n+Pop. 1 locus 93 : 30 21\n+Pop. 1 locus 94 : 30 0\n+Pop. 1 locus 95 : 30 28\n+Pop. 1 locus 96 : 30 30\n+Pop. 1 locus 97 : 30 0\n+Pop. 1 locus 98 : 30 30\n+Pop. 1 locus 99 : 30 0\n+Pop. 1 locus 100 : 30 30\n+\n+Pop. 2 locus 1 : 30 13\n+Pop. 2 locus 2 : 30 3\n+Pop. 2 locus 3 : 30 26\n+Pop. 2 locus 4 : 30 1\n+Pop. 2 locus 5 : 30 30\n+Pop. 2 locus 6 : 30 23\n+Pop. 2 locus 7 : 30 0\n+Pop. 2 locus 8 : 30 30\n+Pop. 2 locus 9 : 30 27\n+Pop. 2 locus 10 : 30 30\n+Pop. 2 locus 11 : 30 10\n+Pop. 2 locus 12 : 30 0\n+Pop. 2 locus 13 : 30 12\n+Pop. 2 locus 14 : 30 30\n+Pop. 2 locus 15 : 30 10\n+Pop. 2 locus 16 : 30 2\n+Pop. 2 locus 17 : 30 30\n+Pop. 2 locus 18 : 30 30\n+Pop. 2 locus 19 : 30 20\n+Pop. 2 locus 20 : 30 30\n+Pop. 2 locus 21 : 30 0\n+Pop. 2 locus 22 : 30 4\n+Pop. 2 locus 23 : 30 30\n+Pop. 2 locus 24 : 30 14\n+Pop. 2 locus 25 : 30 30\n+Pop. 2 locus 26 : 30 0\n+Pop. 2 locus 27 : 30 29\n+Pop. 2 locus 28 : 30 4\n+Pop. 2 locus 29 : 30 1\n+Pop. 2 locus 30 : 30 30\n+Pop. 2 locus 31 : 30 30\n+Pop. 2 locus 32 : 30 1\n+Pop. 2 locus 33 : 30 0\n+Pop. 2 locus 34 : 30 0\n+Pop. 2 locus 35 : 30 29\n+Pop. 2 locus 36 : 30 28\n+Pop. 2 locus 37 : 30 22\n+Pop. 2 locus 38 : 30 30\n+Pop. 2 locus 39 : 30 29\n+Pop. 2 locus 40 : 30 29\n+Pop. 2 locus 41 : 30 15\n+Pop. 2 locus 42 : 30 9\n+Pop. 2 locus 43'..b'8\n+Pop. 9 locus 49 : 30 17\n+Pop. 9 locus 50 : 30 6\n+Pop. 9 locus 51 : 30 0\n+Pop. 9 locus 52 : 30 29\n+Pop. 9 locus 53 : 30 30\n+Pop. 9 locus 54 : 30 30\n+Pop. 9 locus 55 : 30 18\n+Pop. 9 locus 56 : 30 28\n+Pop. 9 locus 57 : 30 30\n+Pop. 9 locus 58 : 30 0\n+Pop. 9 locus 59 : 30 0\n+Pop. 9 locus 60 : 30 17\n+Pop. 9 locus 61 : 30 9\n+Pop. 9 locus 62 : 30 29\n+Pop. 9 locus 63 : 30 9\n+Pop. 9 locus 64 : 30 14\n+Pop. 9 locus 65 : 30 30\n+Pop. 9 locus 66 : 30 30\n+Pop. 9 locus 67 : 30 27\n+Pop. 9 locus 68 : 30 30\n+Pop. 9 locus 69 : 30 27\n+Pop. 9 locus 70 : 30 0\n+Pop. 9 locus 71 : 30 12\n+Pop. 9 locus 72 : 30 14\n+Pop. 9 locus 73 : 30 14\n+Pop. 9 locus 74 : 30 3\n+Pop. 9 locus 75 : 30 26\n+Pop. 9 locus 76 : 30 30\n+Pop. 9 locus 77 : 30 15\n+Pop. 9 locus 78 : 30 29\n+Pop. 9 locus 79 : 30 29\n+Pop. 9 locus 80 : 30 17\n+Pop. 9 locus 81 : 30 0\n+Pop. 9 locus 82 : 30 30\n+Pop. 9 locus 83 : 30 30\n+Pop. 9 locus 84 : 30 18\n+Pop. 9 locus 85 : 30 0\n+Pop. 9 locus 86 : 30 8\n+Pop. 9 locus 87 : 30 30\n+Pop. 9 locus 88 : 30 2\n+Pop. 9 locus 89 : 30 29\n+Pop. 9 locus 90 : 30 23\n+Pop. 9 locus 91 : 30 12\n+Pop. 9 locus 92 : 30 30\n+Pop. 9 locus 93 : 30 29\n+Pop. 9 locus 94 : 30 0\n+Pop. 9 locus 95 : 30 23\n+Pop. 9 locus 96 : 30 30\n+Pop. 9 locus 97 : 30 4\n+Pop. 9 locus 98 : 30 30\n+Pop. 9 locus 99 : 30 0\n+Pop. 9 locus 100 : 30 30\n+\n+Pop. 10 locus 1 : 30 17\n+Pop. 10 locus 2 : 30 10\n+Pop. 10 locus 3 : 30 29\n+Pop. 10 locus 4 : 30 0\n+Pop. 10 locus 5 : 30 30\n+Pop. 10 locus 6 : 30 28\n+Pop. 10 locus 7 : 30 0\n+Pop. 10 locus 8 : 30 30\n+Pop. 10 locus 9 : 30 25\n+Pop. 10 locus 10 : 30 30\n+Pop. 10 locus 11 : 30 11\n+Pop. 10 locus 12 : 30 0\n+Pop. 10 locus 13 : 30 1\n+Pop. 10 locus 14 : 30 30\n+Pop. 10 locus 15 : 30 0\n+Pop. 10 locus 16 : 30 0\n+Pop. 10 locus 17 : 30 29\n+Pop. 10 locus 18 : 30 30\n+Pop. 10 locus 19 : 30 23\n+Pop. 10 locus 20 : 30 19\n+Pop. 10 locus 21 : 30 0\n+Pop. 10 locus 22 : 30 0\n+Pop. 10 locus 23 : 30 30\n+Pop. 10 locus 24 : 30 29\n+Pop. 10 locus 25 : 30 24\n+Pop. 10 locus 26 : 30 0\n+Pop. 10 locus 27 : 30 12\n+Pop. 10 locus 28 : 30 0\n+Pop. 10 locus 29 : 30 1\n+Pop. 10 locus 30 : 30 23\n+Pop. 10 locus 31 : 30 30\n+Pop. 10 locus 32 : 30 10\n+Pop. 10 locus 33 : 30 0\n+Pop. 10 locus 34 : 30 2\n+Pop. 10 locus 35 : 30 30\n+Pop. 10 locus 36 : 30 30\n+Pop. 10 locus 37 : 30 8\n+Pop. 10 locus 38 : 30 24\n+Pop. 10 locus 39 : 30 29\n+Pop. 10 locus 40 : 30 30\n+Pop. 10 locus 41 : 30 19\n+Pop. 10 locus 42 : 30 17\n+Pop. 10 locus 43 : 30 0\n+Pop. 10 locus 44 : 30 4\n+Pop. 10 locus 45 : 30 27\n+Pop. 10 locus 46 : 30 0\n+Pop. 10 locus 47 : 30 8\n+Pop. 10 locus 48 : 30 9\n+Pop. 10 locus 49 : 30 7\n+Pop. 10 locus 50 : 30 18\n+Pop. 10 locus 51 : 30 0\n+Pop. 10 locus 52 : 30 28\n+Pop. 10 locus 53 : 30 30\n+Pop. 10 locus 54 : 30 30\n+Pop. 10 locus 55 : 30 9\n+Pop. 10 locus 56 : 30 17\n+Pop. 10 locus 57 : 30 30\n+Pop. 10 locus 58 : 30 0\n+Pop. 10 locus 59 : 30 0\n+Pop. 10 locus 60 : 30 5\n+Pop. 10 locus 61 : 30 4\n+Pop. 10 locus 62 : 30 30\n+Pop. 10 locus 63 : 30 22\n+Pop. 10 locus 64 : 30 0\n+Pop. 10 locus 65 : 30 30\n+Pop. 10 locus 66 : 30 30\n+Pop. 10 locus 67 : 30 26\n+Pop. 10 locus 68 : 30 30\n+Pop. 10 locus 69 : 30 24\n+Pop. 10 locus 70 : 30 0\n+Pop. 10 locus 71 : 30 11\n+Pop. 10 locus 72 : 30 9\n+Pop. 10 locus 73 : 30 10\n+Pop. 10 locus 74 : 30 0\n+Pop. 10 locus 75 : 30 30\n+Pop. 10 locus 76 : 30 30\n+Pop. 10 locus 77 : 30 27\n+Pop. 10 locus 78 : 30 19\n+Pop. 10 locus 79 : 30 30\n+Pop. 10 locus 80 : 30 30\n+Pop. 10 locus 81 : 30 5\n+Pop. 10 locus 82 : 30 30\n+Pop. 10 locus 83 : 30 30\n+Pop. 10 locus 84 : 30 10\n+Pop. 10 locus 85 : 30 0\n+Pop. 10 locus 86 : 30 23\n+Pop. 10 locus 87 : 30 30\n+Pop. 10 locus 88 : 30 0\n+Pop. 10 locus 89 : 30 22\n+Pop. 10 locus 90 : 30 9\n+Pop. 10 locus 91 : 30 26\n+Pop. 10 locus 92 : 30 30\n+Pop. 10 locus 93 : 30 28\n+Pop. 10 locus 94 : 30 0\n+Pop. 10 locus 95 : 30 26\n+Pop. 10 locus 96 : 30 30\n+Pop. 10 locus 97 : 30 0\n+Pop. 10 locus 98 : 30 30\n+Pop. 10 locus 99 : 30 0\n+Pop. 10 locus 100 : 30 30\n+\n' |