likSurface: Likelihood surfaces

View source: R/profL.R

likSurfaceR Documentation

Likelihood surfaces

Description

Evalues multivariable profile likelihoods

Usage

likSurface(fit, se = seq(0, 3, len = 10))

Arguments

fit

model fit from, for example, glm

se

number of standard errors from MLE of each parameter

Value

List which includes an array of likelihoods

Examples

    # 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))


Pmisc documentation built on Feb. 14, 2024, 3 a.m.