[R] can you extract predictions for gaussian process in using lme, can you fit matern class in lme

[Taylor-Robinson, David]

Dear R users group

I have a couple of questions:

1. I'm fitting a linear mixed effects model using lme (package nlme) with gaussian errors (using exponential correlation structure at the moment): Is it possible to extract predictions for the gaussian process?

exp1 <- corExp(value=c(7,(40/300)),form=~age|id,nugget=T)
exp1 <- Initialize(exp1,d)
mv <- lme(fev1 ~ age, data=d, random= ~ 1|id, method="ML",correlation=exp1)

2. Can you fit an "off the shelf" matern class in lme, and if not has anyone coded one up using the constructor function in lme?

Any advice appreciated

Best wishes

David

Dr David Taylor-Robinson
MRC Population Health Scientist
Department of Health Inequalities and the Social Determinants of Health
University of Liverpool
Whelan Building
Liverpool
L69 3GB

0151 794 4314
dctr at liv.ac.uk