# \donttest{
data(GasolineYield)
# Fit nested models
fit_restricted <- gkwreg(yield ~ temp, data = GasolineYield, family = "kw")
fit_full <- gkwreg(yield ~ batch + temp, data = GasolineYield, family = "kw")
# Likelihood ratio test
lrtest(fit_restricted, fit_full)
# }
Run the code above in your browser using DataLab