[R] similar limma's contrasts.fit() for lme (mixed effect model) object

[shirley zhang]

Dear R help,

In limma package, contrasts.fit() function is very useful. I am
wondering whether there is a similar function for lme object, which
means given a mixed linear model fit, compute estimated coefficients
and standard errors for a given set of contrasts.

Thanks,
Shirley