library(DeclareDesign)
library(ggplot2)
dat <- fabricate(
N = 1000,
A = rnorm(N),
B = rnorm(N),
Z = complete_rs(N),
Y = A*Z + rnorm(N))
# note: remove num.threads = 1 to use more processors
estimates <- causal_forest_handler(data = dat, covariate_names = c("A", "B"), num.threads = 1)
ggplot(data = estimates, aes(A, pred)) + geom_point()
Run the code above in your browser using DataLab