Commits

Giovanni Marco Dall'Olio  committed b7de7af

ADD: script used by Irene to launch simulations

  • Participants
  • Parent commits 3075fc0

Comments (0)

Files changed (3)

File bin/model.test

-1	2.25135e-09
-63552	3.71702e-06
-63553	2.25135e-09
-79619	3.80844e-07
-79620	2.25135e-09
-79636	2.49054e-04
-79637	2.25135e-09
-80202	3.82834e-09
-80203	2.25135e-09
-81058	3.40463e-04
-81059	2.25135e-09
-96030	4.09795e-04
-96031	2.25135e-09
-106930	1.64606e-06
-106931	2.25135e-09
-107176	8.03586e-05
-107177	2.25135e-09
-107226	2.25135e-09
-107227	2.25135e-09
-116150	5.57769e-06
-116151	2.25135e-09
-116169	4.87374e-04
-116170	2.25135e-09
-191629	1.89948e-06
-191630	2.25135e-09
-198054	5.71155e-06
-198055	2.25135e-09
-221644	1.71899e-04
-221645	2.25135e-09
-221645	2.14628e-04
-221646	2.25135e-09
-221683	1.34468e-03
-221684	2.25135e-09
-222438	1.02649e-03
-222439	2.25135e-09
-222454	5.13092e-06
-222455	2.25135e-09
-262621	5.95597e-09
-262622	2.25135e-09
-262975	9.58612e-07
-262976	2.25135e-09
-268721	8.42742e-06
-268722	2.25135e-09
-268880	2.39729e-05
-268881	2.25135e-09
-276713	1.60709e-05
-276714	2.25135e-09
-281365	4.29015e-05
-281366	2.25135e-09
-281436	9.90965e-05
-281437	2.25135e-09
-314606	4.58283e-05
-314607	2.25135e-09
-314770	2.46754e-04
-314771	2.25135e-09
-333085	1.52458e-05
-333086	2.25135e-09
-335754	8.10121e-04
-335755	2.25135e-09
-370430	2.68284e-05
-370431	2.25135e-09
-370539	7.72021e-06
-370540	2.25135e-09
-370575	8.75805e-05
-370576	2.25135e-09
-429490	2.26501e-04
-429491	2.25135e-09
-492121	7.39705e-05
-492122	2.25135e-09
-500241	3.16566e-06
-500242	2.25135e-09
-504893	6.03137e-06
-504894	2.25135e-09
-513785	8.96092e-05
-513786	2.25135e-09
-513790	2.58438e-04
-513791	2.25135e-09
-519079	6.04212e-04
-519080	2.25135e-09
-543325	1.55985e-04
-543326	2.25135e-09
-553667	1.04328e-05
-553668	2.25135e-09
-570110	2.08581e-05
-570111	2.25135e-09
-570455	5.34964e-05
-570456	2.25135e-09
-627049	6.65792e-04
-627050	2.25135e-09
-627095	5.45923e-07
-627096	2.25135e-09
-627121	5.50419e-05
-627122	2.25135e-09
-629468	1.65001e-04
-629469	2.25135e-09
-631266	1.11951e-05
-631267	2.25135e-09
-634839	5.23912e-08
-634840	2.25135e-09
-638315	1.20085e-04
-638316	2.25135e-09
-641031	8.93889e-04
-641032	2.25135e-09
-641150	1.33759e-05
-641151	2.25135e-09
-648264	2.92965e-04
-648265	2.25135e-09
-648267	4.63829e-04
-648268	2.25135e-09
-649565	1.24430e-04
-649566	2.25135e-09
-655297	2.13711e-03
-655298	2.25135e-09
-658193	3.02478e-05
-658194	2.25135e-09
-677908	2.67801e-09
-677909	2.25135e-09
-678702	4.49808e-06
-678703	2.25135e-09
-703853	3.09363e-05
-703854	2.25135e-09
-703896	8.70477e-04
-703897	2.25135e-09
-705604	6.43649e-04
-705605	2.25135e-09
-713590	5.68652e-06
-713591	2.25135e-09
-718151	6.40389e-05
-718152	2.25135e-09
-722474	8.25119e-04
-722475	2.25135e-09
-722539	5.47510e-04
-722540	2.25135e-09
-722609	3.63312e-05
-722610	2.25135e-09
-722648	2.61380e-09
-722649	2.25135e-09
-724707	8.20241e-04
-724708	2.25135e-09
-760575	3.84631e-09
-760576	2.25135e-09
-762699	4.65657e-04
-762700	2.25135e-09
-764877	7.76432e-07
-764878	2.25135e-09
-850244	2.78472e-05
-850245	2.25135e-09
-868143	1.12572e-05
-868144	2.25135e-09
-872534	2.12519e-04
-872535	2.25135e-09
-872635	1.94607e-04
-872636	2.25135e-09
-872637	3.35004e-05
-872638	2.25135e-09
-920803	2.42055e-04
-920804	2.25135e-09
-923814	2.35457e-09
-923815	2.25135e-09
-923819	4.57767e-05
-923820	2.25135e-09
-933303	1.53301e-06
-933304	2.25135e-09
-936211	5.53995e-06
-936212	2.25135e-09
-975261	3.69658e-04
-975262	2.25135e-09
+1	3.77941e-10
+29761	1.58486e-05
+29762	3.77941e-10
+42204	5.89969e-07
+42205	3.77941e-10
+42205	7.50514e-06
+42206	3.77941e-10
+44372	4.19899e-06
+44373	3.77941e-10
+49041	6.43195e-06
+49042	3.77941e-10
+57499	2.16034e-06
+57500	3.77941e-10
+65270	1.00436e-06
+65271	3.77941e-10
+106696	5.68301e-08
+106697	3.77941e-10
+108542	1.12141e-05
+108543	3.77941e-10
+108744	4.46872e-06
+108745	3.77941e-10
+150903	4.81381e-06
+150904	3.77941e-10
+156000	1.65204e-07
+156001	3.77941e-10
+156282	1.78179e-08
+156283	3.77941e-10
+162523	1.02660e-07
+162524	3.77941e-10
+172547	5.62925e-07
+172548	3.77941e-10
+184488	9.46549e-09
+184489	3.77941e-10
+184664	1.14035e-05
+184665	3.77941e-10
+211348	1.95265e-09
+211349	3.77941e-10
+224532	4.35420e-07
+224533	3.77941e-10
+250191	6.06932e-07
+250192	3.77941e-10
+250406	5.37313e-10
+250407	3.77941e-10
+254971	2.53684e-05
+254972	3.77941e-10
+255457	1.03111e-09
+255458	3.77941e-10
+257371	1.60828e-06
+257372	3.77941e-10
+259553	9.11885e-07
+259554	3.77941e-10
+259662	7.58679e-06
+259663	3.77941e-10
+270770	5.61837e-09
+270771	3.77941e-10
+271227	1.90773e-05
+271228	3.77941e-10
+271244	1.63797e-06
+271245	3.77941e-10
+271259	9.97257e-06
+271260	3.77941e-10
+289114	3.88320e-07
+289115	3.77941e-10
+289291	1.11963e-05
+289292	3.77941e-10
+289954	7.48615e-07
+289955	3.77941e-10
+290700	8.51531e-07
+290701	3.77941e-10
+298197	5.30204e-06
+298198	3.77941e-10
+319368	1.57674e-06
+319369	3.77941e-10
+323007	8.60956e-08
+323008	3.77941e-10
+323593	4.53830e-08
+323594	3.77941e-10
+327344	1.70941e-07
+327345	3.77941e-10
+333365	1.22971e-09
+333366	3.77941e-10
+333594	3.77941e-10
+333595	3.77941e-10
+333601	6.23437e-07
+333602	3.77941e-10
+335415	2.77678e-07
+335416	3.77941e-10
+339258	3.81915e-06
+339259	3.77941e-10
+341545	2.29996e-06
+341546	3.77941e-10
+341753	2.79587e-06
+341754	3.77941e-10
+349015	9.54847e-09
+349016	3.77941e-10
+349357	1.03575e-06
+349358	3.77941e-10
+362881	8.42847e-07
+362882	3.77941e-10
+404986	1.18154e-06
+404987	3.77941e-10
+405705	7.52449e-07
+405706	3.77941e-10
+405725	5.28603e-10
+405726	3.77941e-10
+412331	3.81084e-10
+412332	3.77941e-10
+429809	6.59553e-09
+429810	3.77941e-10
+430013	1.49192e-05
+430014	3.77941e-10
+437937	4.15220e-10
+437938	3.77941e-10
+447623	4.24340e-07
+447624	3.77941e-10
+447624	2.45066e-06
+447625	3.77941e-10
+462072	1.48698e-08
+462073	3.77941e-10
+470406	4.89214e-08
+470407	3.77941e-10
+470592	3.79427e-10
+470593	3.77941e-10
+471515	9.90847e-06
+471516	3.77941e-10
+478452	4.80144e-06
+478453	3.77941e-10
+499306	5.43030e-09
+499307	3.77941e-10
+519000	4.83018e-07
+519001	3.77941e-10
+521412	2.93323e-06
+521413	3.77941e-10
+522111	5.34041e-06
+522112	3.77941e-10
+523454	1.70283e-06
+523455	3.77941e-10
+529710	9.71163e-07
+529711	3.77941e-10
+529927	1.53416e-06
+529928	3.77941e-10
+532925	1.33904e-06
+532926	3.77941e-10
+581817	1.38178e-07
+581818	3.77941e-10
+581818	7.92382e-06
+581819	3.77941e-10
+590651	1.83720e-06
+590652	3.77941e-10
+606399	7.35884e-08
+606400	3.77941e-10
+610930	1.09446e-05
+610931	3.77941e-10
+617675	2.02582e-09
+617676	3.77941e-10
+621423	3.11311e-07
+621424	3.77941e-10
+624902	1.83868e-06
+624903	3.77941e-10
+627146	7.43833e-07
+627147	3.77941e-10
+630733	2.98900e-08
+630734	3.77941e-10
+638854	6.74400e-08
+638855	3.77941e-10
+639083	4.32968e-06
+639084	3.77941e-10
+658997	5.61351e-08
+658998	3.77941e-10
+659034	2.53464e-07
+659035	3.77941e-10
+659475	1.24170e-07
+659476	3.77941e-10
+659557	1.60565e-06
+659558	3.77941e-10
+673783	3.82368e-08
+673784	3.77941e-10
+674271	5.06022e-07
+674272	3.77941e-10
+690529	6.44553e-07
+690530	3.77941e-10
+692230	9.28690e-07
+692231	3.77941e-10
+701066	4.93853e-10
+701067	3.77941e-10
+701068	8.58739e-07
+701069	3.77941e-10
+701300	2.54325e-06
+701301	3.77941e-10
+701341	1.66425e-05
+701342	3.77941e-10
+712871	2.93964e-06
+712872	3.77941e-10
+714069	2.58602e-08
+714070	3.77941e-10
+714565	4.93259e-07
+714566	3.77941e-10
+722511	9.32584e-06
+722512	3.77941e-10
+733431	1.67456e-05
+733432	3.77941e-10
+736286	5.35142e-06
+736287	3.77941e-10
+801528	5.36447e-08
+801529	3.77941e-10
+801745	2.62432e-05
+801746	3.77941e-10
+804676	1.26204e-06
+804677	3.77941e-10
+805236	1.31477e-06
+805237	3.77941e-10
+808299	6.38127e-07
+808300	3.77941e-10
+842281	1.81167e-07
+842282	3.77941e-10
+842282	1.52172e-08
+842283	3.77941e-10
+842367	9.03718e-07
+842368	3.77941e-10
+842901	3.60543e-06
+842902	3.77941e-10
+843056	7.26718e-07
+843057	3.77941e-10
+874146	3.50732e-07
+874147	3.77941e-10
+887299	4.21014e-07
+887300	3.77941e-10
+947934	3.24056e-08
+947935	3.77941e-10
+948273	4.42803e-10
+948274	3.77941e-10
+955789	1.21364e-08
+955790	3.77941e-10
+955826	4.32215e-06
+955827	3.77941e-10
+957323	1.28761e-06
+957324	3.77941e-10
 # per bp per generation
 mutation_rate 1.5e-8
 
-#recomb_file /lustre/scratch103/sanger/igr/simulations/model.test_77
+#recomb_file /lustre/scratch103/sanger/igr/simulations/model.test_75
 gene_conversion_rate 4.5e-9
 
 # population info
 
 pop_event change_size "african pop size" 3 17000 12500
 
-pop_event sweep "selective sweep" 1 0 0.0064 0.5 0.4
+pop_event sweep "selective sweep" 3 0 0.0172 0.5 1
  

File irene_script.pl

+for ($i=1; $i<=$iterate; $i++){
+       my $resp = `~/software/cosi_1.2/recosim/nfs/users/nfs_i/igr/scripts/cosi/recParams $length`;
+       print $resp;
+       my @recom = ("cp /lustre/scratch103/sanger/igr/simulations/model.test /lustre/scratch103/sanger/igr/test/bestfit_neutral_params_rec_$i");
+       system (@recom);
+       my $output = ("/lustre/scratch103/sanger/igr/test/bestfit_neutral_params_$i");
+       $resp = "~/software/cosi_1.2/coalescent -p $model -o $output";
+       system($resp);
+}