Dear Zelig developers!
I have already complained about the issue of using colnames(model.matrix()) to get the
coefficient names in qi.mixed.R
http://lists.hmdc.harvard.edu/lists/zelig/2008_01/msg00020.html
Recent version of Zelig (3.4.0) again uses
fTermsNames <- colnames(model.matrix(parsefml$fixed, data = D))
which IMHO would be better as
fTermsNames <- dimnames(betas)[[2]]
I noticed the same use for random effects, but I can not suggest the solution, because
I am not sure how "special" random effects are handled, i.e., varying intercept
and varying
slopes model etc. Using the model.matrix() just for getting the names of location
parameters
in the model surely is not the most efficient way, since model.matrix can be a really fat
memory hog for large datasets.
Lep pozdrav / With regards,
Gregor Gorjanc
----------------------------------------------------------------------
University of Ljubljana PhD student
Biotechnical Faculty www:
http://gregor.gorjanc.googlepages.com
Animal Science Department blog:
http://ggorjan.blogspot.com
Groblje 3 mail: gregor.gorjanc <at> bfro.uni-lj.si
SI-1230 Domzale fax: +386 (0)1 72 17 888
Slovenia, Europe tel: +386 (0)1 72 17 861
----------------------------------------------------------------------
-
Zelig Mailing List, served by Harvard-MIT Data Center
Send messages: zelig(a)lists.gking.harvard.edu
[un]subscribe Options:
http://lists.gking.harvard.edu/?info=zelig
Zelig program information:
http://gking.harvard.edu/zelig/