Displaying 1 result from an estimated 1 matches for "x7_scalar".
Did you mean:
x6_scalar
2007 May 25
0
Help with complex lme model fit
...Ident (~ X6
-
1),
pdIdent (~ X7 - 1), pdIdent(~ rand_mat -1))), data
=
gpdat)
Which should recover the variance components:
var_label var_est
rand_mat_scalar 0.00021983
X6_scalar 0.62314002
X7_scalar 0.03853604
as recovered by GenStat and used to generate the dataset. Instead I
get:
X6 0.6231819
X7 0.05221481
rand_mat 1.377596e-11
However, If I change or drop either of X5 or X6. I then get much
closer
estimates to what is expected. For example:
mod2 <- lme...