data(qardl_sim)
fit <- qardl(y ~ x1 + x2, data = qardl_sim, tau = c(0.25, 0.50, 0.75), p = 2, q = 2)
wald_results <- qardl_wald(fit)
print(wald_results)
# Pairwise tests
wald_pairwise <- qardl_wald(fit, pairwise = TRUE)
print(wald_pairwise)
Run the code above in your browser using DataLab