set.seed(123)
n <- 500
x <- rchisq(n, df = 4) - 4
e <- rchisq(n, df = 3) - 3
y <- 0.5 * x + e
d <- data.frame(x, y)
result <- dda.resdist(y ~ x, pred = "x", data = d,
B = 50, conf.level = 0.90, prob.trans = TRUE)
print(result)
Run the code above in your browser using DataLab