#' @title plot.comp2regr.ecdf
#' @description TODO
#' @export
plot.comp2regr.ecdf <- function(x, ...){
graphics::par(mfrow = c(2, 2))
graphics::plot(c(x$x1, x$x2), c(x$y1, x$y2), type = "n", xlab = "covariate",
ylab = "response", main = "resgression functions")
graphics::points(x$x1, x$y1, col = "red")
graphics::points(x$x2, x$y2, col = "blue")
graphics::lines(sort(x$x1), x$m1x1.hat[order(x$x1)], col = "red")
graphics::lines(sort(x$x1), x$m0x1.hat[order(x$x1)], col = "red", lty = 2)
graphics::lines(sort(x$x2), x$m2x2.hat[order(x$x2)], col = "blue")
graphics::lines(sort(x$x2), x$m0x2.hat[order(x$x2)], col = "blue", lty = 2)
graphics::plot(c(x$x1, x$x2), c(x$y1, x$y2), type = "n", xlab = "covariate",
ylab = "sigma", main = "conditional variance functions")
graphics::lines(sort(x$x1), x$sigma1x1.hat[order(x$x1)], col = "red")
graphics::lines(sort(x$x2), x$sigma2x2.hat[order(x$x2)], col = "blue")
graphics::plot(ecdf(x$eps1.hat), col = "red", main = "F_eps1")
graphics::plot(ecdf(x$eps01.hat), col = "black", add = T)
graphics::plot(ecdf(x$eps2.hat), col = "blue", main = "F_eps2")
graphics::plot(ecdf(x$eps02.hat), col = "black", add = T)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.