This function generate a Q-Q plot for residuals.
Examples
# Instantiate
x <- rand_uniform()
e <- rand_normal()
test <- vi_model(prm = list(x = x, e = e),
prm_type = list(x = "r", e = "r"),
formula = y ~ 1 + x + 10 * x^2 + e,
null_formula = y ~ x,
alt_formula = y ~ x + I(x^2))
# Residual plot of the null model
test$plot_qq(test$gen(1000, fit_model = TRUE))