bbmle (version 1.0.25.1)

as.data.frame.profile.mle2: convert profile to data frame

Description

converts a profile of a fitted mle2 object to a data frame

Usage

# S3 method for profile.mle2
as.data.frame(x, row.names=NULL,
optional=FALSE, ...)

Value

a data frame with columns

param

name of parameter being profiled

z

signed square root of the deviance difference from the minimum

parameter values

named par.vals.parname

focal

value of focal parameter: redundant, but included for plotting convenience

Arguments

x

a profile object

row.names

row names (unused)

optional

unused

...

unused

Author

Ben Bolker

Examples

  ## use as.data.frame and lattice to plot profiles
  x <- 0:10
  y <- c(26, 17, 13, 12, 20, 5, 9, 8, 5, 4, 8)
  library(bbmle)
  LL <- function(ymax=15, xhalf=6) {
        -sum(stats::dpois(y, lambda=ymax/(1+x/xhalf), log=TRUE))
    }
  ## uses default parameters of LL
  fit1 <- mle2(LL)
  p1 <- profile(fit1)
  d1 <- as.data.frame(p1)
  library(lattice)
  xyplot(abs(z)~focal|param,data=d1,
    subset=abs(z)<3,
    type="b",
    xlab="",
    ylab=expression(paste(abs(z),
    " (square root of ",Delta," deviance)")),
    scale=list(x=list(relation="free")))