data("tobin", package = "survival")
# selection double hurdle model
model110i <- mhurdle(durable ~ age | quant | 0, tobin, dist = "n")
# double-hurdle p-tobit
model011i <- mhurdle(durable ~ 0 | quant | age, tobin, dist = "n")
# Vuong test for strictly non-nested models
vuongtest(model011i, model110i)
Run the code above in your browser using DataLab