pneumo <- transform(pneumo, let = log(exposure.time))
fit1 <- vglm(cbind(normal, mild, severe) ~ let, acat,
trace = TRUE, data = pneumo)
pfit1 <- profile(fit1, trace = FALSE)
if (FALSE) {
vplot.profile(pfit1)
vpairs.profile(pfit1)
}
Run the code above in your browser using DataLab