[R] specifying random effects covariance structure in nlme

Chris Stride c@b@@tride @ending from @heffield@@c@uk
Sat Jun 16 13:00:23 CEST 2018


I'm trying to fit a mixed effects exponential decay model, in which I 
have random effects for the initial value (init), the asymptote (asymp), 
and the rate (rate).

The catch is that I'd also like to estimate the correlation between init 
and asymp, but not between init and rate, or asymp and rate.

Now using  random = pdDiag(init + asymp + rate ~ 1) has none of the 
random effects correlated

And using  random = pdSymm(init + asymp + rate ~ 1) has all three of the 
random effects correlated

How do I specify just the correlation I want?



More information about the R-help mailing list