This function is a wrapper for lme.formula that allows
Brownian motion, fractional Brownian motion or integrated Ornstein-Uhlenbeck
components to be included in linear mixed models, with related parameter
estimates and confidence intervals returned in their natural parameterisation.
An optional corStruct object describing the within-group
covariance structure. In addition to those available in nlme,
covBM can be used to incorporate a Brownian motion component, covFracBM
can be used to incorporate a fractional Brownian motion component and covIOU
can be used to incorporate an integrated Ornstein-Uhlenbeck process in relation to
a continuous variable.