ordinal: slice.clm – R documentation

Become an expert in R — Interactive courses, Cheat Sheets, certificates and more!

slice.clm

Slice the likelihood of a clm

Description

Slice likelihood and plot the slice. This is usefull for illustrating the likelihood surface around the MLE (maximum likelihood estimate) and provides graphics to substantiate (non-)convergence of a model fit. Also, the closeness of a quadratic approximation to the log-likelihood function can be inspected for relevant parameters. A slice is considerably less computationally demanding than a profile.

Usage

slice(object, ...)

## S3 method for class 'clm'
slice(object, parm = seq_along(par), lambda = 3,
     grid = 100, quad.approx = TRUE, ...)

## S3 method for class 'slice.clm'
plot(x, parm = seq_along(x),
    type = c("quadratic", "linear"), plot.mle = TRUE,
    ask = prod(par("mfcol")) < length(parm) && dev.interactive(), ...)

Arguments

`object`	for the `clm` method an object of class `"clm"`, i.e., the result of a call to `clm`.
`x`	a `slice.clm` object, i.e., the result of `slice(clm.object)`.
`parm`	for `slice.clm` a numeric or character vector indexing parameters, for `plot.slice.clm` only a numeric vector is accepted. By default all parameters are selected.
`lambda`	the number of curvature units on each side of the MLE the slice should cover.
`grid`	the number of values at which to compute the log-likelihood for each parameter.
`quad.approx`	compute and include the quadratic approximation to the log-likelihood function?
`type`	`"quadratic"` plots the log-likelihood function which is approximately quadratic, and `"linear"` plots the signed square root of the log-likelihood function which is approximately linear.
`plot.mle`	include a vertical line at the MLE (maximum likelihood estimate) when `type = "quadratic"`? Ignored for `type = "linear"`.
`ask`	logical; if `TRUE`, the user is asked before each plot, see `par(ask=.)`.
`...`	further arguments to `plot.default` for the plot method. Not used in the slice method.

Value

The slice method returns a list of data.frames with one data.frame for each parameter slice. Each data.frame contains in the first column the values of the parameter and in the second column the values of the (positive) log-likelihood "logLik". A third column is present if quad.approx = TRUE and contains the corresponding quadratic approximation to the log-likelihood. The original model fit is included as the attribute "original.fit".

The plot method produces a plot of the likelihood slice for each parameter.

Author(s)

Rune Haubo B Christensen

Examples

## fit model:
fm1 <- clm(rating ~ contact + temp, data = wine)
## slice the likelihood:
sl1 <- slice(fm1)

## three different ways to plot the slices:
par(mfrow = c(2,3))
plot(sl1)
plot(sl1, type = "quadratic", plot.mle = FALSE)
plot(sl1, type = "linear")

## Verify convergence to the optimum:
sl2 <- slice(fm1, lambda = 1e-5, quad.approx = FALSE)
plot(sl2)

ordinal

Regression Models for Ordinal Data

v2019.12-10

GPL (>= 2)

Authors

Rune Haubo Bojesen Christensen [aut, cre]

Initial release

2019-12-10