likSurface {Pmisc} | R Documentation |
Evalues multivariable profile likelihoods
likSurface(fit, se = seq(0, 3, len = 10))
fit |
model fit from, for example, glm |
se |
number of standard errors from MLE of each parameter |
List which includes an array of likelihoods
# from glm help file counts <- c(18,17,15,20,10,20,25,13,12) outcome <- gl(3,1,9) treatment <- gl(3,3) glm.D93 <- glm(counts ~ outcome, family = poisson()) x = likSurface(glm.D93) image(x$parameters[,1], x$parameters[,2], apply(x$logLik, c(2,3), min))