Message-ID: <ef0374eb-d395-4a8d-994e-25ca78c649e9@sheffield.ac.uk>
Date: 2018-06-16T11:01:49Z
From: Dr C B Stride
Subject: random effects covariance matrix specification in nlme
In-Reply-To: <3602e304-c27f-0b79-35fd-65268807a931@highstat.com>
Hi
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?
cheers
Chris