Displaying 1 result from an estimated 1 matches for "var_est".
Did you mean:
largest
2007 May 25
0
Help with complex lme model fit
...<- lme(fixed = Y1 ~ X1 + X2 + X3 + X4 + X5 + m_sum,
random = pdBlocked(list(pdIdent(~1), pdIdent (~ 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 o...