[R] standard errors for predict.nls?

Christian Ritz ritz at life.ku.dk
Mon Nov 3 16:15:34 CET 2008


Dear Christoph,

using the package 'alr3' it's not difficult!

Have a look at the following example:

## Fitting a Michaelis-Menten model
Puromycin.m1<-nls(rate~a*conc/(b+conc), data=Puromycin[1:12,], start=list(a=200, b=1))

library(alr3)

## Predictions (with standard errors) at concentrations 0.02, 0.06, and 011
delta.method(Puromycin.m1, "a*0.02/(b+0.02)")
delta.method(Puromycin.m1, "a*0.06/(b+0.06)")
delta.method(Puromycin.m1, "a*0.11/(b+0.11)")


Hope you can make it work for your own application!

Christian



More information about the R-help mailing list