evaluating variance functions in nlme
Yes, this probably should go to R-help, and before you do that, I suggest you PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html. It can increase the likelihood you will get a useful reply relatively quickly. I tried the your two calls to "gnls", which of course produced errors for me as I don't have your data. The help page for "gnls" included an example that looked quite similar to me: > fm1 <- gnls(weight ~ SSlogis(Time, Asym, xmid, scal), Soybean, + weights = varPower()) > summary(fm1)> fm1 <- gnls(weight ~ SSlogis(Time, Asym, xmid, scal), Soybean, + weights = varPower()) > summary(fm1) Those results plus "str(fm1)" looked like they might help answer your question. However, I don't understand enough of what you are asking to say. If an answer might still be worth pursuing to you, I suggest you read the posting guide and submit a question following that model to r-help. spencer graves
Nicholas Lewin-Koh wrote:
Hi,
I guess this is a final plea, and maybe this should go to R-help but
here goes.
I am writing a set of functions for calibration and prediction, and to
calculate standard
errors and intervals I need the variance function to be evaluated at new
prediction points.
So for instance
fit<-gnls(Y~SSlogis(foo,Asym,xmid,scal),weights=varPower())
fit2<-gnls(Y~SSlogis(foo,Asym,xmid,scal),weights=varPower(form=~foo))
Now using fit or fit2 I would like to get the variance function
evaluated at new points.
I have played with getCovariateFormula, and looked at Initialize.gnls,
summary etc.
but it is not clear to me how to evaluate the form component, especially
in the case of fit
above where form=~fitted(.), in any safe way.
I can grep for "fitted" in the formula eg.
grep("fitted",deparse(getCovariateFormula(fit$modelStruct$varFunc)))
and try to calculate predicted values from the model for the new points
but how to substitute back in the new terms?
I don't need this problem solved on a platter, I just need to unedrstand
an approach,
because my stabs are failing.
Thanks
Nicholas
______________________________________________ R-devel at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-devel
Spencer Graves, PhD Senior Development Engineer PDF Solutions, Inc. 333 West San Carlos Street Suite 700 San Jose, CA 95110, USA spencer.graves at pdf.com www.pdf.com <http://www.pdf.com> Tel: 408-938-4420 Fax: 408-280-7915