likSurface | 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))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.