code for nonlinear least square using newton raphson method

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

NLS.pdf (102.99 KB, 需要: 10 个论坛币)

2012-12-4 06:50:23 上传

极有用
需要: 10 个论坛币  [购买]

找到这个东东以后我骂了句粗话，早找到就好了，写这个code写得我晕了，function不同的要要稍微改改。

必须跟大家分享了。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏4 回帖

关键词：Nonlinear nonlinea Nonline Square Method method least

精彩部分，免费分享，欲知详情请下载全文。

#Exponential model
fun=function(th,x){
y = th[2]+((th[1]-th[2])*exp(-th[3]*x))
return(y)
}
fun(th,x) #try it out -- should give you fitted values
#Calculate the residual sum of squares
resid=function(th,x,y){
yhat=fun(th,x)
rss=t(y-yhat)%*%(y-yhat)
return(rss)
}
resid(th,x,y) #try it out -- should give a scalar for output
#Calulate the Jacobian matrix
firder=function(th,x){
j1=exp(-th[3]*x)
j2=1-j1
j3=-x*(th[1]-th[2])*exp(-th[3]*x)
j=cbind(j1,j2,j3)
return(j)
}
firder(th,x)
#Calculate the Second Derivatives
secder=function(th,x){
h11=rep(0,times=length(x))
h21=h11
h22=h11
h31=-x*exp(-th[3]*x)
h32=-h31
h33=(x**2)*((th[1]-th[2])*exp(-th[3]*x))
hmat=cbind(h11,h21,h22,h31,h32,h33)
return(hmat)
}
secder(th,x)
#Gradient evaluated at "th"
grad=function(th,x,y){
fd=firder(th,x)
hm=secder(th,x)
yhat=fun(th,x)
r=y-yhat

10个币？

具体什么用处？