[R] nls question

Bill Simpson wsi at gcal.ac.uk
Fri Apr 27 11:25:36 CEST 2001


I have a question about passing arguments to the function f that nlm
minimizes.

I have no problems if I do this:
x<-seq(0,1,.1)
y<-1.1*x + (1-1.1) + rnorm(length(x),0,.1)
fn<-function(p)
	{
	yhat<-p*x+(1-p)
	sum((y-yhat)^2)
	}
out<-nlm(fn,p=1.5,hessian=TRUE)

But I would like to define
fn<-function(x,y,p)
        {
        yhat<-p*x+(1-p)
        sum((y-yhat)^2)
        }
so that I can pass the x and y variables along, doing something like
out<-nlm(fn,x=this, y=that,p=1.5,hessian=TRUE)
This version doesn't work of course. I have looked at the examples and
can't figure out how to do what I want. Is it possible?

Thanks very much for any help.

(BTW I know nls() can be used for nonlinear least squares fitting. It
choked on this problem.)

Bill Simpson

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list