# NOT RUN {
n <- 100
outcome <- rnorm(n, 0, 1)
trt <- rbinom(n, 1, 0.5)
trt_pxy <- matrix(rnorm(n, 0, 1), ncol = 1)
out_pxy <- matrix(rnorm(n, 0, 1), ncol = 1)
covariates <- matrix(rnorm(n, 0, 1), ncol = 1)
pcl_object <- pcl(outcome, trt, trt_pxy, out_pxy, covariates)
fit <- pclfit(pcl_object)
# }
Run the code above in your browser using DataLab