comp: Prediction for 'lmer' objects

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

Usage

comp(fit, form = terms(fit), varname, varpattern = "", data = getData(fit), ...)

Arguments

fit
form
varname
varpattern
data
...

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 (fit, form = terms(fit), varname, varpattern = "", data = getData(fit), 
    ...) 
{
    getData <- function(x) UseMethod("getData")
    getData.lme <- function(x) x$data
    getData.lmer <- function(x) x@frame
    disp(form)
    disp(dim(data))
    model.mat <- model.matrix(form, data, ...)
    disp(dim(model.mat))
    ret <- data[rownames(model.mat), , drop = F]
    fe <- fixef(fit)
    effect.names <- grep(varpattern, names(fe), value = T)
    ret$comp <- c(model.mat[, effect.names] %*% fe[effect.names])
    attr(ret, "predictors") <- effect.names
    ret
  }

gmonette/spida documentation built on May 17, 2019, 7:25 a.m.