Fit a Linear Mixed Effects Model
DESCRIPTION:
A method for the generic function lme for objects inheriting from class lmList.
USAGE: 
lme.lmList(fixed, random = fixed, cluster, data = sys.parent(),
           start = list(D = NULL, theta = NULL, alpha = NULL,
           delta = NULL), est.method = c("RML", "ML"), re.block,
           re.structure = c("unstructured", "diagonal", "identity",
           "compsymm", "ar1"), re.paramtr = c("matrixlog", "cholesky",
           "spherical", "logcholesky", "givens"),
           serial.structure = c("identity", "ar1", "ar1.continuous",
           "compsymm", "ar2", "ma1", "ma2"), serial.covariate = NULL,
           serial.covariate.transformation = c("rank.within.cluster",
           "none", "round", "global.rank"), var.function = c("identity",
           "power", "expon", "cte.power"), var.covariate = NULL,
           var.estimate = T, na.action, control)
REQUIRED ARGUMENTS:
fixed:
an object inheriting from class lmList.
OPTIONAL ARGUMENTS:
random:
cluster:
data:
start:
these arguments, though included in the calling sequence to be consistent with the generic function call, should not be supplied in the call to lme.lmList(), as they are obtained from object. Initial values given to these arguments will be overwritten with the values obtained from fixed.
others:
see the documentation on lme.formula
VALUE:
Generic functions such as print, plot and summary have methods to show the results of the fit. See lme.object for the components of the fit. The functions residuals, coefficients and fitted.values can be used to extract some of its components.
SEE ALSO:
lmList , lme.formula , lme.object .
EXAMPLES: 
Ovary.lmList <- lmList(follicles ~ sin(2 * pi * Time) +
                       cos(2 * pi * Time),
                       data = Ovary, cluster = ~ Mare)


ovary.lmList.fit <- lme(Ovary.lmList)


ovary.lmList.fit


# Returns the following:
Call:
  Fixed: follicles ~ sin(2 * pi * Time) + cos(2 * pi * Time)
 Random: follicles ~ sin(2 * pi * Time) + cos(2 * pi * Time)
Cluster:  ~ Mare
   Data: Ovary


Variance/Covariance Components Estimate(s):


  Structure: unstructured
  Parametrization: matrixlog
  Standard Deviation(s) of Random Effect(s)
 (Intercept) sin(2 * pi * Time) cos(2 * pi * Time)
    3.229274           2.092799           1.066919
 Correlation of Random Effects
                   (Intercept) sin(2 * pi * Time)
sin(2 * pi * Time) -0.5697627
cos(2 * pi * Time) -0.8015247   0.1783619


 Cluster Residual Variance: 9.114855


Fixed Effects Estimate(s):
  (Intercept) sin(2 * pi * Time) cos(2 * pi * Time)
    12.18591          -3.296679         -0.8731317


Number of Observations: 308
Number of Clusters: 11