Taylor-Robinson, David
2011-Feb-23 11:57 UTC
[R] can you extract predictions for gaussian process in using lme, can you fit matern class in lme
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