Commits

Giovanni Marco Dall'Olio  committed 6a7f489

REFACT: keeping the .pos-1/2/3 and .hap-1/2/3 files in the results

  • Participants
  • Parent commits f193bf2

Comments (0)

Files changed (3)

File bin/model.test

-1	6.45998e-09
-16729	4.42110e-05
-16730	6.45998e-09
-16911	5.31905e-05
-16912	6.45998e-09
-51362	9.48934e-05
-51363	6.45998e-09
-60527	1.45550e-04
-60528	6.45998e-09
-85004	8.78432e-05
-85005	6.45998e-09
-87698	1.32578e-05
-87699	6.45998e-09
-88711	3.49265e-05
-88712	6.45998e-09
-88712	1.24941e-05
-88713	6.45998e-09
-89408	3.43308e-05
-89409	6.45998e-09
-89681	2.55479e-07
-89682	6.45998e-09
-90012	1.47019e-04
-90013	6.45998e-09
-104709	1.01626e-04
-104710	6.45998e-09
-105208	1.83174e-07
-105209	6.45998e-09
-107640	3.98121e-05
-107641	6.45998e-09
-145604	6.00750e-05
-145605	6.45998e-09
-167600	6.87432e-05
-167601	6.45998e-09
-175218	2.09023e-05
-175219	6.45998e-09
-175221	6.36380e-05
-175222	6.45998e-09
-178655	2.60928e-06
-178656	6.45998e-09
-183210	1.31393e-03
-183211	6.45998e-09
-183211	7.70500e-05
-183212	6.45998e-09
-204580	3.16075e-04
-204581	6.45998e-09
-216817	2.07928e-05
-216818	6.45998e-09
-228706	6.96908e-07
-228707	6.45998e-09
-243700	6.56951e-09
-243701	6.45998e-09
-255981	3.53365e-05
-255982	6.45998e-09
-256009	1.13852e-05
-256010	6.45998e-09
-287518	1.22591e-03
-287519	6.45998e-09
-331685	2.79672e-06
-331686	6.45998e-09
-375987	1.64384e-05
-375988	6.45998e-09
-384173	8.39366e-05
-384174	6.45998e-09
-391846	5.45144e-05
-391847	6.45998e-09
-394428	6.43799e-05
-394429	6.45998e-09
-406108	2.87408e-06
-406109	6.45998e-09
-406132	5.62938e-05
-406133	6.45998e-09
-433952	6.17438e-07
-433953	6.45998e-09
-434766	7.63706e-06
-434767	6.45998e-09
-487699	3.94078e-06
-487700	6.45998e-09
-487742	2.18766e-05
-487743	6.45998e-09
-493291	2.32957e-05
-493292	6.45998e-09
-493292	1.18569e-05
-493293	6.45998e-09
-496832	9.33022e-06
-496833	6.45998e-09
-504549	1.30071e-04
-504550	6.45998e-09
-517359	1.15190e-05
-517360	6.45998e-09
-520065	4.88849e-04
-520066	6.45998e-09
-523296	1.89838e-07
-523297	6.45998e-09
-523321	7.06767e-07
-523322	6.45998e-09
-531745	2.81392e-07
-531746	6.45998e-09
-571656	7.69714e-09
-571657	6.45998e-09
-583915	1.39634e-05
-583916	6.45998e-09
-586542	2.54590e-06
-586543	6.45998e-09
-587909	6.54054e-06
-587910	6.45998e-09
-590569	3.57423e-05
-590570	6.45998e-09
-592792	6.66699e-05
-592793	6.45998e-09
-603621	2.24611e-05
-603622	6.45998e-09
-603749	3.71147e-06
-603750	6.45998e-09
-605962	2.59705e-04
-605963	6.45998e-09
-606529	4.64871e-06
-606530	6.45998e-09
-606610	7.90841e-06
-606611	6.45998e-09
-629641	1.17098e-04
-629642	6.45998e-09
-629643	2.66746e-04
-629644	6.45998e-09
-656850	7.58523e-06
-656851	6.45998e-09
-675666	1.81962e-07
-675667	6.45998e-09
-681100	5.12769e-05
-681101	6.45998e-09
-686530	1.00930e-06
-686531	6.45998e-09
-686735	1.24703e-04
-686736	6.45998e-09
-686737	2.30948e-05
-686738	6.45998e-09
-686846	2.71776e-03
-686847	6.45998e-09
-687231	2.78948e-07
-687232	6.45998e-09
-689058	1.42564e-07
-689059	6.45998e-09
-689059	1.88465e-07
-689060	6.45998e-09
-689110	1.42073e-07
-689111	6.45998e-09
-694605	3.42260e-05
-694606	6.45998e-09
-742360	1.64878e-06
-742361	6.45998e-09
-745600	2.27287e-04
-745601	6.45998e-09
-750009	2.87847e-06
-750010	6.45998e-09
-750055	5.65489e-05
-750056	6.45998e-09
-750069	2.04456e-04
-750070	6.45998e-09
-778006	1.00709e-03
-778007	6.45998e-09
-790308	3.51371e-04
-790309	6.45998e-09
-790526	2.01022e-04
-790527	6.45998e-09
-833081	1.42186e-04
-833082	6.45998e-09
-850751	8.83702e-08
-850752	6.45998e-09
-861474	5.13581e-05
-861475	6.45998e-09
-876815	3.89720e-04
-876816	6.45998e-09
-878334	1.09240e-04
-878335	6.45998e-09
-883059	2.03735e-04
-883060	6.45998e-09
-884735	4.46816e-05
-884736	6.45998e-09
-934088	2.57428e-06
-934089	6.45998e-09
-934098	4.92019e-05
-934099	6.45998e-09
-943301	1.60342e-05
-943302	6.45998e-09
-943335	3.16324e-07
-943336	6.45998e-09
-946814	1.71728e-04
-946815	6.45998e-09
-949438	9.80147e-07
-949439	6.45998e-09
-949867	1.36662e-04
-949868	6.45998e-09
-949868	1.34890e-04
-949869	6.45998e-09
-950281	6.92339e-09
-950282	6.45998e-09
-950797	3.72067e-05
-950798	6.45998e-09
-952415	1.73353e-04
-952416	6.45998e-09
-958675	2.11599e-05
-958676	6.45998e-09
-967506	6.75444e-05
-967507	6.45998e-09
-972218	1.41983e-06
-972219	6.45998e-09
-972259	7.87384e-07
-972260	6.45998e-09
-972276	3.71420e-05
-972277	6.45998e-09
-972332	1.08166e-06
-972333	6.45998e-09
-976619	9.38260e-05
-976620	6.45998e-09
-976622	1.24012e-05
-976623	6.45998e-09
-977328	2.75143e-07
-977329	6.45998e-09
-977517	3.04172e-04
-977518	6.45998e-09
+1	1.11785e-09
+32352	3.20602e-07
+32353	1.11785e-09
+32605	2.57085e-05
+32606	1.11785e-09
+32702	3.96388e-08
+32703	1.11785e-09
+113716	4.05470e-05
+113717	1.11785e-09
+113985	6.27382e-07
+113986	1.11785e-09
+114003	9.44954e-06
+114004	1.11785e-09
+114502	3.05542e-05
+114503	1.11785e-09
+114565	5.13360e-06
+114566	1.11785e-09
+114567	1.49376e-08
+114568	1.11785e-09
+121680	1.12505e-06
+121681	1.11785e-09
+124266	4.18912e-08
+124267	1.11785e-09
+124267	1.97316e-05
+124268	1.11785e-09
+124274	2.23378e-05
+124275	1.11785e-09
+126292	3.46695e-05
+126293	1.11785e-09
+128722	5.26568e-05
+128723	1.11785e-09
+190676	1.34152e-04
+190677	1.11785e-09
+193139	9.18896e-06
+193140	1.11785e-09
+210106	4.51419e-07
+210107	1.11785e-09
+220170	9.84752e-08
+220171	1.11785e-09
+224722	3.08265e-09
+224723	1.11785e-09
+226186	8.55001e-07
+226187	1.11785e-09
+230160	4.54641e-05
+230161	1.11785e-09
+233599	1.18761e-09
+233600	1.11785e-09
+263867	1.11816e-09
+263868	1.11785e-09
+282603	8.43617e-06
+282604	1.11785e-09
+285301	1.80774e-05
+285302	1.11785e-09
+300016	2.34012e-08
+300017	1.11785e-09
+305280	4.04126e-06
+305281	1.11785e-09
+391588	3.51822e-06
+391589	1.11785e-09
+391989	1.07549e-07
+391990	1.11785e-09
+399924	1.55095e-05
+399925	1.11785e-09
+399968	7.37114e-06
+399969	1.11785e-09
+429257	7.79562e-06
+429258	1.11785e-09
+433256	9.32120e-06
+433257	1.11785e-09
+433261	7.26441e-07
+433262	1.11785e-09
+436385	1.97113e-09
+436386	1.11785e-09
+535277	9.51835e-07
+535278	1.11785e-09
+540292	1.46342e-05
+540293	1.11785e-09
+545817	4.59579e-05
+545818	1.11785e-09
+545860	3.69360e-05
+545861	1.11785e-09
+558007	2.48117e-05
+558008	1.11785e-09
+563470	1.66989e-08
+563471	1.11785e-09
+597097	2.03823e-05
+597098	1.11785e-09
+597121	2.32422e-05
+597122	1.11785e-09
+621309	9.01308e-06
+621310	1.11785e-09
+621311	2.84448e-05
+621312	1.11785e-09
+621313	3.58263e-09
+621314	1.11785e-09
+626225	1.92955e-07
+626226	1.11785e-09
+626948	2.47917e-06
+626949	1.11785e-09
+628011	7.80181e-06
+628012	1.11785e-09
+638271	3.00856e-08
+638272	1.11785e-09
+641279	1.32259e-06
+641280	1.11785e-09
+648580	9.94222e-06
+648581	1.11785e-09
+665368	3.54096e-07
+665369	1.11785e-09
+665410	2.79807e-09
+665411	1.11785e-09
+665680	1.57637e-08
+665681	1.11785e-09
+668109	2.12424e-05
+668110	1.11785e-09
+669673	6.98756e-06
+669674	1.11785e-09
+687439	3.25445e-06
+687440	1.11785e-09
+709375	3.46588e-06
+709376	1.11785e-09
+727712	5.72463e-06
+727713	1.11785e-09
+727914	2.59992e-05
+727915	1.11785e-09
+746244	2.31289e-06
+746245	1.11785e-09
+750930	3.74379e-05
+750931	1.11785e-09
+750949	6.01106e-05
+750950	1.11785e-09
+752446	3.09426e-05
+752447	1.11785e-09
+755179	1.90160e-07
+755180	1.11785e-09
+756667	5.08071e-08
+756668	1.11785e-09
+759247	3.60395e-06
+759248	1.11785e-09
+759311	2.25346e-07
+759312	1.11785e-09
+788165	1.01966e-05
+788166	1.11785e-09
+788166	3.91968e-07
+788167	1.11785e-09
+820741	1.69614e-05
+820742	1.11785e-09
+821012	1.38424e-07
+821013	1.11785e-09
+821047	5.85124e-07
+821048	1.11785e-09
+869119	3.90101e-09
+869120	1.11785e-09
+870074	1.36517e-06
+870075	1.11785e-09
+872005	1.61104e-04
+872006	1.11785e-09
+872175	1.06815e-05
+872176	1.11785e-09
+909867	1.58756e-08
+909868	1.11785e-09
+910068	2.76280e-08
+910069	1.11785e-09
+937382	7.97989e-05
+937383	1.11785e-09
+944646	2.44659e-06
+944647	1.11785e-09
+957765	2.35302e-09
+957766	1.11785e-09
+958531	3.65970e-06
+958532	1.11785e-09
+958952	6.18972e-07
+958953	1.11785e-09
+966864	2.94370e-06
+966865	1.11785e-09
+969087	1.30390e-06
+969088	1.11785e-09
+971990	7.18968e-09
+971991	1.11785e-09
+980990	1.63003e-06
+980991	1.11785e-09
 # per bp per generation
 mutation_rate 1.5e-8
 
-#recomb_file /lustre/scratch103/sanger/igr/simulations/model.test_82
+#recomb_file /lustre/scratch103/sanger/igr/simulations/model.test_63
 gene_conversion_rate 4.5e-9
 
 # population info
 
 pop_event change_size "african pop size" 3 17000 12500
 
-pop_event sweep "selective sweep" 2 0 0.0064 0.5 0.4
+pop_event sweep "selective sweep" 1 0 0.008 0.5 0.6
  

File generate_cosi_simulations.py

         print('cp %s %s' %(oldMapFile, newMapFile ))
         print('cp %s %s' %(oldHapFile, newHapFile ))
 #
-        print( 'rm ' + 'out*' )
-        print( 'rm ' + coreName + '*')
+        # Move out.pos-1/2/3 and out.hap-1/2/3 files
+        print('mv %s %s' % ("out.hap-1", outputFolder + 'coalesc_sim' + population + '_neutral_' + str(scenario_id) + '.hap-1')) 
+        print('mv %s %s' % ("out.hap-2", outputFolder + 'coalesc_sim' + population + '_neutral_' + str(scenario_id) + '.hap-2')) 
+        print('mv %s %s' % ("out.hap-3", outputFolder + 'coalesc_sim' + population + '_neutral_' + str(scenario_id) + '.hap-3')) 
+        print('mv %s %s' % ("out.pos-1", outputFolder + 'coalesc_sim' + population + '_neutral_' + str(scenario_id) + '.pos-1')) 
+        print('mv %s %s' % ("out.pos-2", outputFolder + 'coalesc_sim' + population + '_neutral_' + str(scenario_id) + '.pos-2')) 
+        print('mv %s %s' % ("out.pos-3", outputFolder + 'coalesc_sim' + population + '_neutral_' + str(scenario_id) + '.pos-3')) 
+
+#        print( 'rm ' + 'out*' )
+#        print( 'rm ' + coreName + '*')