[R-sig-ME] Generalized least squares program
smckinney at bccrc.ca
Wed Jul 16 19:29:16 CEST 2008
Have you tried the predict() function?
You should be able to hand off new data
to your saved fitted model object.
for an example.
From: r-sig-mixed-models-bounces at r-project.org on behalf of Gang Chen
Sent: Wed 7/16/2008 8:21 AM
To: r-sig-mixed-models at r-project.org
Subject: [R-sig-ME] Generalized least squares program
I'm using function gls in nlme to run some extended linear model with
generalized least squares many many times in a loop with the same
design matrix X and same residual covariance structure but different
values for the response variable. A lot of time is wasted in
reformulating the same model structures and overhead matrix
calculations. Such a waste makes the runtime range from one day to
one month. Are there some slots available in gls that would allow me
to avoid the unnecessary waste? If not, is there a counterpart of gls
for extended linear model with generalized least squares available in
lme4 that would give such slots?
R-sig-mixed-models at r-project.org mailing list
More information about the R-sig-mixed-models