Displaying 1 result from an estimated 1 matches for "mat_ht".
Did you mean:
pat_ht
2006 Jun 01
1
setting the random-effects covariance matrix in lme
...ss the model would have the following form (in hierarchical notation)
Yi|bi,k ~ N(XiB+Zibi, sigmak*Ident)
bi|k ~ N(0, Dk)
K~Bernoulli(p)
I can obtain different sigmas (sigma0 and sigma1 based on the factor 'dx') using the weights option in the call to lme:
lme(fixed = height ~ -1 + bage + mat_ht + pat_ht + dx + time:dx + time2:dx ,
weights=varIdent(form=~1|dx),
random = ~ time + time2 |subject, data = pilot, na.action=na.omit)
but I cannot seem to be able to get two different matrices for the random effects.
I'm particularly interested in obta...