lme4 (version 1.1-7)

refit: Refit a model with a new response, by maximum likelihood criterion

Description

Refit a model with a different response vector

Usage

refit(object, newresp, ...)
  ## S3 method for class 'merMod':
refit(object, newresp = NULL, rename.response=FALSE, ...)

Arguments

object

a fitted model, usually of class lmerMod, to be refit with a new response

newresp

an (optional) numeric vector providing the new response, of the same length as the original response (see Details for information on NA handling). May also be a data frame with a single numeric column, e.g. as produced by simula

rename.response

when refitting the model, should the name of the response variable in the formula and model frame be replaced with the name of newresp?

...

optional additional parameters. None are used at present.

Value

• an object like x, but fit by maximum likelihood

Details

Refit a model, possibly after modifying the response vector. This could be done using an update method but this approach should be faster because it bypasses the creation of the model representation and goes directly to the optimization step.

Setting rename.response to TRUE may be necessary if one wants to do further operations (such as update) on the fitted model. However, the refitted model will still be slightly different from the equivalent model fitted via update; in particular, the terms component is not updated to reflect the new response variable, if it has a different name from the original.

If newresp has an na.action attribute, then it is assumed that NA values have already been removed from the numeric vector; this allows the results of simulate(object) to be used even if the original response vector contained NA values. Otherwise, the length of newresp must be the same as the original length of the

Examples

## ex. 1: using refit() to fit each column in a matrix of responses
set.seed(101)
Y <- matrix(rnorm(1000),ncol=10)
res <- list()
## combine first column of responses with predictor variables
d <- data.frame(y=Y[,1],x=rnorm(100),f=rep(1:10,10))
## (use check.conv.grad="ignore" to disable convergence checks because we
##  are using a fake example)
## fit first response
fit1 <- lmer(y~x+(1|f),data=d,
             control=lmerControl(check.conv.grad="ignore",
             check.conv.hess="ignore"))
## combine fit to first response with fits to remaining responses
res <- c(fit1,lapply(as.data.frame(Y[,-1]),refit,object=fit1))
## ex. 2: refitting simulated data using data that contain NA values
sleepstudyNA <- sleepstudy
sleepstudyNA$Reaction[1:3] <- NA
fm0 <- lmer(Reaction~Days+(1|Subject),sleepstudyNA)
## the special case of refitting with a single simulation works ...
ss0 <- refit(fm0,simulate(fm0))
## ... but if simulating multiple responses (for efficiency),
## need to use na.action=na.exclude in order to have proper length of data
fm1 <- lmer(Reaction~Days+(1|Subject),sleepstudyNA,na.action=na.exclude)
ss <- simulate(fm1,5)
res2 <- refit(fm1,ss[[5]])