]> code.communitydata.science - ml_measurement_error_public.git/blob - simulations/Makefile
update simulation and mle code
[ml_measurement_error_public.git] / simulations / Makefile
1
2 SHELL=bash
3
4 Ns=[1000, 2000, 4000]
5 ms=[200, 400, 800]
6 seeds=[$(shell seq -s, 1 250)]
7 explained_variances=[0.1]
8
9 all:remembr.RDS remember_irr.RDS
10
11 srun=srun -A comdata -p compute-bigmem --time=6:00:00 --mem 4G -c 1
12
13
14 joblists:example_1_jobs example_2_jobs example_3_jobs
15
16 # test_true_z_jobs: test_true_z.R simulation_base.R
17 #       grid_sweep.py --command "Rscript test_true_z.R" --arg_dict '{"N":${Ns},"m":${ms}, "seed":${seeds}, "outfile":["test_true_z.feather"], "y_explained_variancevari":${explained_variances}, "Bzx":${Bzx}}' --outfile test_true_z_jobsb
18
19 # test_true_z.feather: test_true_z_jobs 
20 #       rm -f test_true_z.feather
21 #       sbatch --wait --verbose --array=1-3000 run_simulation.sbatch 0 test_true_z_jobs
22 #       sbatch --wait --verbose --array=3001-6001 run_simulation.sbatch 0 test_true_z_jobs
23
24
25 example_1_jobs: 01_two_covariates.R simulation_base.R
26         grid_sweep.py --command "Rscript 01_two_covariates.R" --arg_dict '{"N":${Ns},"m":${ms}, "seed":${seeds}, "outfile":["example_1.feather"], "y_explained_variance":${explained_variances}, "Bzx":[0.1]}' --outfile example_1_jobs
27
28 example_1.feather: example_1_jobs 
29         rm -f example_1.feather
30         sbatch --wait --verbose --array=1-$(shell cat example_1_jobs | wc -l) run_simulation.sbatch 0 example_1_jobs
31 #       sbatch --wait --verbose --array=3001-6001 run_simulation.sbatch 0 example_1_jobs
32
33 example_2_jobs: 02_indep_differential.R simulation_base.R
34         grid_sweep.py --command "Rscript 02_indep_differential.R" --arg_dict '{"N":${Ns},"m":${ms}, "seed":${seeds}, "outfile":["example_2.feather"],"y_explained_variance":${explained_variances},  "Bzy":[-0.3],"Bxy":[0.3],"Bzx":[0.3], "outcome_formula":["y~x+z"], "proxy_formula":["w_pred~y*z*x"]}' --outfile example_2_jobs
35
36 example_2.feather: example_2_jobs 
37         rm -f example_2.feather
38         sbatch --wait --verbose --array=1-$(shell cat example_2_jobs | wc -l) run_simulation.sbatch 0 example_2_jobs
39 #       sbatch --wait --verbose --array=3001-6001 run_simulation.sbatch 0 example_2_jobs
40
41 # example_2_B_jobs: example_2_B.R
42 #       grid_sweep.py --command "Rscript example_2_B.R" --arg_dict '{"N":${Ns},"m":${ms}, "seed":${seeds}, "outfile":["example_2_B.feather"]}' --outfile example_2_B_jobs
43
44 # example_2_B.feather: example_2_B_jobs
45 #       rm -f example_2_B.feather
46 #       sbatch --wait --verbose --array=1-3000 run_simulation.sbatch 0 example_2_B_jobs
47
48 example_3_jobs: 03_depvar.R simulation_base.R
49         grid_sweep.py --command "Rscript 03_depvar.R" --arg_dict '{"N":${Ns},"m":${ms}, "seed":${seeds}, "outfile":["example_3.feather"], "y_explained_variance":${explained_variances}}' --outfile example_3_jobs
50
51 example_3.feather: example_3_jobs
52         rm -f example_3.feather 
53         sbatch --wait --verbose --array=1-$(shell cat example_3_jobs | wc -l)  run_simulation.sbatch 0 example_3_jobs
54
55 example_4_jobs: 04_depvar_differential.R simulation_base.R
56         grid_sweep.py --command "Rscript 04_depvar_differential.R" --arg_dict '{"N":${Ns},"m":${ms}, "seed":${seeds}, "outfile":["example_4.feather"], "y_explained_variance":${explained_variances}}' --outfile example_4_jobs
57
58 example_4.feather: example_4_jobs
59         rm -f example_4.feather 
60         sbatch --wait --verbose --array=1-$(shell cat example_4_jobs | wc -l)  run_simulation.sbatch 0 example_4_jobs
61
62
63 remembr.RDS:example_1.feather example_2.feather example_3.feather example_4.feather plot_example.R plot_dv_example.R
64         rm -f remembr.RDS
65         ${srun} Rscript plot_example.R --infile example_1.feather --name "plot.df.example.1"
66         ${srun} Rscript plot_example.R --infile example_2.feather --name "plot.df.example.2"
67         ${srun} Rscript plot_dv_example.R --infile example_3.feather --name "plot.df.example.3"
68         ${srun} Rscript plot_dv_example.R --infile example_4.feather --name "plot.df.example.4"
69
70
71 irr_Ns = ${Ns}
72 irr_ms = ${ms}
73 irr_seeds=${seeds}
74 irr_explained_variances=${explained_variances}
75
76 example_5_jobs: 05_irr_indep.R irr_simulation_base.R
77         grid_sweep.py --command "Rscript 05_irr_indep.R" --arg_dict '{"N":${irr_Ns},"m":${irr_ms}, "seed":${irr_seeds}, "outfile":["example_5.feather"], "y_explained_variance":${irr_explained_variances}}' --outfile example_5_jobs
78
79 example_5.feather:example_5_jobs
80         rm -f example_5.feather
81         sbatch --wait --verbose --array=1-$(shell cat example_5_jobs | wc -l)  run_simulation.sbatch 0 example_5_jobs
82
83
84 example_6_jobs: 06_irr_dv.R irr_dv_simulation_base.R
85         grid_sweep.py --command "Rscript 06_irr_dv.R" --arg_dict '{"N":${irr_Ns},"m":${irr_ms}, "seed":${irr_seeds}, "outfile":["example_6.feather"], "y_explained_variance":${irr_explained_variances}}' --outfile example_6_jobs
86
87 example_6.feather:example_6_jobs
88         rm -f example_6.feather
89         sbatch --wait --verbose --array=1-$(shell cat example_6_jobs | wc -l)  run_simulation.sbatch 0 example_6_jobs
90
91 remember_irr.RDS: example_5.feather example_6.feather plot_irr_example.R plot_irr_dv_example.R
92         rm -f remember_irr.RDS
93         ${srun} Rscript plot_irr_example.R --infile example_5.feather --name "plot.df.example.5"
94         ${srun} Rscript plot_irr_dv_example.R --infile example_6.feather --name "plot.df.example.6"
95
96 clean:
97         rm *.feather
98         rm -f remembr.RDS
99         rm -f example_*_jobs
100 #       sbatch --wait --verbose --array=3001-6001 run_simulation.sbatch 0 example_2_B_jobs
101
102 # example_2_B_mecor_jobs:
103 #       grid_sweep.py --command "Rscript example_2_B_mecor.R" --arg_dict '{"N":${Ns},"m":${ms}, "seed":${seeds}, "outfile":["example_2_B_mecor.feather"]}' --outfile example_2_B_mecor_jobs
104
105 # example_2_B_mecor.feather:example_2_B_mecor.R example_2_B_mecor_jobs
106 #       rm -f example_2_B_mecor.feather
107 #       sbatch --wait --verbose --array=1-3000 run_simulation.sbatch 0 example_2_B_mecor_jobs
108 #       sbatch --wait --verbose --array=3001-6001 run_simulation.sbatch 0 example_2_B_mecor_jobs
109
110
111

Community Data Science Collective || Want to submit a patch?