[R] NLM

Uwe Ligges ligges at statistik.tu-dortmund.de
Sun Sep 27 19:11:27 CEST 2009



Andrew Wang wrote:
> I am trying to understand NLM package, so I generated this data set consisting y and x using 
> 
> y= a + b*x +c*x^2 + N(0,10), with a=3.5,b=4.5,c=5.5
> 
> 
> Given y and x, I am trying to use NLM to have estimates of parameters a, b and c that minimize the least square error  
> 
> my code looks like
>  
> f<- function(y,x,a,b,c) {sum((y-(a+b*x+c*x^2))^2)}
> nlm(f, y, x, a=3, b=4,c=5)
> 
> 
> But it comes up with rediculous result. My understanding of this NLM must be wrong.



See ?nlm for the arguments you have to pass to nlm() and how the 
function needs to be written. There are also some examples on that help 
page.

Uwe Ligges



> Please help!
> 
> Andy
> 
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.




More information about the R-help mailing list