# NOT RUN {
mte_fit <- mte(selection = d ~ x + z, outcome = y ~ x,
method = "localIV", data = toydata)
x <- seq(0.05, 0.95, 0.05)
mte_tilde_p <- eval_mte_tilde(mte_fit, p = x, u = 0.5)$mte_tilde
mte_tilde_u <- eval_mte_tilde(mte_fit, p = 0.5, u = x)$mte_tilde
mprte_tilde_p <- eval_mte_tilde(mte_fit, p = x, u = x)$mte_tilde
out <- cbind(mte_tilde_p, mte_tilde_u, mprte_tilde_p)
matplot(x = x, y = out, type = "l")
# }
Run the code above in your browser using DataLab