Re: [R] standard errors for predict.nls?

2008-11-12 Thread Ben Bolker
Bolker > Sent: Tue 11/11/2008 3:09 PM > To: [EMAIL PROTECTED] > Cc: r-help@r-project.org > Subject: Re: [R] standard errors for predict.nls? > > > Here's how far I've gotten. It's a start, but > I can't so far find a way to extract the actual > t

Re: [R] standard errors for predict.nls?

2008-11-12 Thread R Heberto Ghezzo, Dr
nd t(qq) gives us the values obviously quantile can be used with any set of p's [EMAIL PROTECTED] -Original Message- From: [EMAIL PROTECTED] on behalf of Ben Bolker Sent: Tue 11/11/2008 3:09 PM To: [EMAIL PROTECTED] Cc: r-help@r-project.org Subject: Re: [R] standard errors for predict.nls

Re: [R] standard errors for predict.nls?

2008-11-11 Thread Ben Bolker
Here's how far I've gotten. It's a start, but I can't so far find a way to extract the actual thing we want -- the bootstrap confidence intervals on the predictions. I apologize for not taking the time to go read up on the boot library etc. etc. (I *have* RTFM, but I haven't backtracked to th

Re: [R] standard errors for predict.nls?

2008-11-09 Thread Rolf Turner
On 7/11/2008, at 11:33 PM, Christoph Scherber wrote: Dear all, I would like to get standard errors (or confidence intervals) for *predicted* values from an nls fit. I have tried to implement code from p.225 in MASS (bootstrapping a nls fit), but this gives only the confidence intervals o

Re: [R] standard errors for predict.nls?

2008-11-07 Thread Christoph Scherber
Dear all, I would like to get standard errors (or confidence intervals) for *predicted* values from an nls fit. I have tried to implement code from p.225 in MASS (bootstrapping a nls fit), but this gives only the confidence intervals of the parameter estimates, but not an overall confidence in

Re: [R] standard errors for predict.nls?

2008-11-06 Thread Christoph Scherber
Dear Prof Ripley, Am I correct if I use the following code to get c.i.´s for predicted values of the nls fit: puro1<-nls(rate~a*conc/(b+conc), data=Puromycin[1:12,], start=list(a=200, b=1)) #set up nls model # assume only one predicted value is obtained using predict(puro1,list(conc=0.02)):

Re: [R] standard errors for predict.nls?

2008-11-04 Thread Prof Brian Ripley
On Mon, 3 Nov 2008, Ben Bolker wrote: Prof Brian Ripley wrote: Christoph Scherber agr.uni-goettingen.de> writes: Dear all, Is there a way to retrieve standard errors from nls models? The help page tells me that arguments such as se.fit are ignored... Many thanks and best wishes Christoph

Re: [R] standard errors for predict.nls?

2008-11-03 Thread Ben Bolker
Prof Brian Ripley wrote: >> Christoph Scherber agr.uni-goettingen.de> >> writes: >> >>> >>> Dear all, >>> >>> Is there a way to retrieve standard errors from nls models? >>> The help page tells me that arguments >>> such as se.fit are ignored... >>> >>> Many thanks and best wishes >>> Christoph >

Re: [R] standard errors for predict.nls?

2008-11-03 Thread Prof Brian Ripley
On Mon, 3 Nov 2008, Ben Bolker wrote: Christoph Scherber agr.uni-goettingen.de> writes: Dear all, Is there a way to retrieve standard errors from nls models? The help page tells me that arguments such as se.fit are ignored... Many thanks and best wishes Christoph I have written some rea

Re: [R] standard errors for predict.nls?

2008-11-03 Thread Christian Ritz
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

Re: [R] standard errors for predict.nls?

2008-11-03 Thread Ben Bolker
Christoph Scherber agr.uni-goettingen.de> writes: > > Dear all, > > Is there a way to retrieve standard errors from nls models? > The help page tells me that arguments > such as se.fit are ignored... > > Many thanks and best wishes > Christoph I have written some reasonably generic delta-

[R] standard errors for predict.nls?

2008-11-03 Thread Christoph Scherber
Dear all, Is there a way to retrieve standard errors from nls models? The help page tells me that arguments such as se.fit are ignored... Many thanks and best wishes Christoph -- Dr. rer.nat. Christoph Scherber University of Goettingen DNPW, Agroecology Waldweg 26 D-37073 Goettingen Germa