logLik1.lme: Calculates contribution of one subject to the log-likelihood...
In nlmeU: Datasets and utility functions enhancing functionality of nlme package.

Description Usage Arguments Details Value Author(s) Examples

This is method for logLik1() generic function.

1 2	## S3 method for class 'lme' logLik1(modfit, dt1, dtInit)

`modfit`	an `lme` object representing model fitted using maximum likelihood.
`dt1`	a data frame with data for one subject, for whom the log-likelihood function is to be evaluated
`dtInit`	an optional auxiliary data frame.

Calculates profile likelihood (with beta profiled out) for *one* subject. Data with *one* level of grouping only. correlation component in modelStruct not implemented.

numeric scalar value representing contribution of a given subject to the overall log-likelihood returned by logLik() function applied to lme object defined by modfit argument.

Andrzej Galecki and Tomasz Burzykowski

require(nlme)
lm3.form <- visual ~ visual0 + time + treat.f
(fm16.5ml <-                                # M16.5
   lme(lm3.form,
   random = list(subject = pdDiag(~time)),
   weights = varPower(form = ~time),
   data = armd, method = "ML"))
 df1 <- subset(armd, subject == "1")       # Panel R20.7
 logLik1(fm16.5ml, df1)