# NOT RUN {
library(CIMTx)
set.seed(1)
idata = data_gen(n = 50, ratio =1,scenario = 1)
trt_ind <- as.numeric(idata$trtdat$trt_ind)
all_vars <- idata$trtdat[, -1] #exclude treatment indicator
y <- idata$Yobs
iptw_multiTrt_ate(y=y, trt = trt_ind,SL.library = c("SL.glm"),
trim_alpha = 0.05, method = "IPTW-Logistics-Trim")
causal_multi_treat(y = y,trt = trt_ind,
method = "IPTW-Logistics", estimand = "ATE")
# }
Run the code above in your browser using DataLab