data(sim_bounded, package = "unitquantreg")
sim_bounded_curr <- sim_bounded[sim_bounded$family == "uweibull", ]
models <- c("uweibull", "kum", "ulogistic")
lt_fits <- lapply(models, function(fam) {
unitquantreg(formula = y1 ~ x, tau = 0.5, data = sim_bounded_curr,
family = fam)
})
ans <- pairwise.vuong.test(lt = lt_fits)
ans
Run the code above in your browser using DataLab