Now, doing python compute_stats.py [treefile] will do replicates, with parameters drawn from Yi-Ming's posteriors in recap/data/posterior.txt. The number of each type of replicate is hardcoded, and they are output to files like:
-rw-r--r-- 1 peter peter 26756 Feb 15 22:50 param_grid/sim_runs/run000003/sim_2377345170879_171455_286145_1.pairstats.csv
-rw-r--r-- 1 peter peter 2926 Feb 15 22:50 param_grid/sim_runs/run000003/sim_2377345170879_171455_286145_1.stats.csv
-rw-r--r-- 1 peter peter 26754 Feb 15 22:50 param_grid/sim_runs/run000003/sim_2377345170879_171455_286145_0.pairstats.csv
-rw-r--r-- 1 peter peter 2926 Feb 15 22:49 param_grid/sim_runs/run000003/sim_2377345170879_171455_286145_0.stats.csv
-rw-r--r-- 1 peter peter 26654 Feb 15 22:49 param_grid/sim_runs/run000003/sim_2377345170879_171455_956804_1.pairstats.csv
-rw-r--r-- 1 peter peter 2902 Feb 15 22:49 param_grid/sim_runs/run000003/sim_2377345170879_171455_956804_1.stats.csv
-rw-r--r-- 1 peter peter 26649 Feb 15 22:49 param_grid/sim_runs/run000003/sim_2377345170879_171455_956804_0.pairstats.csv
-rw-r--r-- 1 peter peter 2906 Feb 15 22:49 param_grid/sim_runs/run000003/sim_2377345170879_171455_956804_0.stats.csv
-rw-r--r-- 1 peter peter 26735 Feb 15 22:49 param_grid/sim_runs/run000003/sim_2377345170879_650513_57938_1.pairstats.csv
-rw-r--r-- 1 peter peter 2929 Feb 15 22:49 param_grid/sim_runs/run000003/sim_2377345170879_650513_57938_1.stats.csv
-rw-r--r-- 1 peter peter 26760 Feb 15 22:49 param_grid/sim_runs/run000003/sim_2377345170879_650513_57938_0.pairstats.csv
-rw-r--r-- 1 peter peter 2925 Feb 15 22:49 param_grid/sim_runs/run000003/sim_2377345170879_650513_57938_0.stats.csv
-rw-r--r-- 1 peter peter 26637 Feb 15 22:48 param_grid/sim_runs/run000003/sim_2377345170879_650513_512413_1.pairstats.csv
-rw-r--r-- 1 peter peter 2910 Feb 15 22:48 param_grid/sim_runs/run000003/sim_2377345170879_650513_512413_1.stats.csv
-rw-r--r-- 1 peter peter 26634 Feb 15 22:48 param_grid/sim_runs/run000003/sim_2377345170879_650513_512413_0.pairstats.csv
-rw-r--r-- 1 peter peter 2901 Feb 15 22:48 param_grid/sim_runs/run000003/sim_2377345170879_650513_512413_0.stats.csv
... but, it seems to work!
Next, we need to figure out if there's a better way to write out the information.
Now, doing
python compute_stats.py [treefile]will do replicates, with parameters drawn from Yi-Ming's posteriors inrecap/data/posterior.txt. The number of each type of replicate is hardcoded, and they are output to files like:... but, it seems to work!
Next, we need to figure out if there's a better way to write out the information.