parser <- add_argument(parser, "--z_bias", help='how is the classifier biased?', default=1.5)
parser <- add_argument(parser, "--Bxy", help='coefficient of x on y', default=0.1)
parser <- add_argument(parser, "--Bzy", help='coeffficient of z on y', default=-0.1)
+parser <- add_argument(parser, "--B0", help='coeffficient of z on y', default=-0.1)
parser <- add_argument(parser, "--outcome_formula", help='formula for the outcome variable', default="y~x+z")
parser <- add_argument(parser, "--proxy_formula", help='formula for the proxy variable', default="w_pred~y+z")
args <- parse_args(parser)
-B0 <- 0
+B0 <- args$B0
Bxy <- args$Bxy
Bzy <- args$Bzy