# NOT RUN {
## generate data
set.seed(1)
n <- 200 # number of observations
p <- 5 # number of covariates
D <- rbinom(n, 1, 0.5) # random treatment assignment
Y <- runif(n) # outcome variable
propensity_scores <- rep(0.5, n) # propensity scores
proxy_BCA <- runif(n) # proxy BCA estimates
proxy_CATE <- runif(n) # proxy CATE estimates
membership <- quantile_group(proxy_CATE) # group membership
## perform BLP
BLP <- BLP(Y, D, propensity_scores, proxy_BCA, proxy_CATE)
## perform GATES
GATES <- GATES(Y, D, propensity_scores, proxy_BCA, proxy_CATE, membership)
## get estimates of the lambda parameters
lambda_parameters(BLP, GATES, proxy_CATE, membership)
# }
Run the code above in your browser using DataLab