# NOT RUN {
res <- HeckmanEM(y, x, w, cc, nu = 4, family = "Normal", error = 1e-05, iter.max = 500,
im = TRUE, criteria = TRUE)
HeckmanEM.envelope(res, ylab="Normalized Quantile Residuals",xlab="Standard normal quantile",
line="quartile", col=c(20,1), pch=19, ylim = c(-5,4))
# }
Run the code above in your browser using DataCamp Workspace