cell.default: Confidence ellipse

In gmonette/spida15: Collection of miscellaneous functions for mixed models etc. prepared for SPIDA 2009+ (development version)

Usage

cell.default(model, which.coef, levels = 0.95, Scheffe = FALSE, dfn = 2, center.pch = 19, center.cex = 1.5, segments = 51, xlab, ylab, las = par("las"), col = palette()[2], lwd = 2, lty = 1, add = FALSE, ...)

Arguments

model
which.coef
levels
Scheffe
dfn
center.pch
center.cex
segments
xlab
ylab
las
col
lwd
lty
add
...

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (model, which.coef, levels = 0.95, Scheffe = FALSE, 
    dfn = 2, center.pch = 19, center.cex = 1.5, segments = 51, 
    xlab, ylab, las = par("las"), col = palette()[2], lwd = 2, 
    lty = 1, add = FALSE, ...) 
{
    which.coef <- if (length(coefficients(model)) == 2) 
        c(1, 2)
    else {
        if (missing(which.coef)) {
            if (any(names(coefficients(model)) == "(Intercept)")) 
                c(2, 3)
            else c(1, 2)
        }
        else which.coef
    }
    coef <- coefficients(model)[which.coef]
    xlab <- if (missing(xlab)) 
        paste(names(coef)[1], "coefficient")
    ylab <- if (missing(ylab)) 
        paste(names(coef)[2], "coefficient")
    if (missing(dfn)) {
        if (Scheffe) 
            dfn <- sum(df.terms(model))
        else 2
    }
    dfd <- df.residual(model)
    shape <- vcov(model)[which.coef, which.coef]
    ret <- numeric(0)
    ret <- ell(coef, shape, sqrt(dfn * qf(levels, dfn, dfd)))
    colnames(ret) <- c(xlab, ylab)
    ret
  }

gmonette/spida15 documentation built on May 17, 2019, 7:26 a.m.