[求助]R语言中的Newton-Raphson算法

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

有哪位大虾指导一下小妹，如何在R中实现Newton-Raphson算法！！！ 我好弱弱呀  万分感谢

[此贴子已经被作者于2007-11-1 8:40:45编辑过]

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Newton PHS RAP NEW APH 语言 算法

范例取自Option Pricing Models and Volatility Using Excel VBA

Page-8 Newton-Raphson Method

f(x) = x^2-7*x+10

f'(x) = 2*x-7

R code:

fun<-function(x){x^2-7*x+10}

dfun<-function(x){2*x-7}

max.it <- 50

tol <- 0.00001

current.x = 1   #start value

for(i in 1:max.it)

{

x=current.x

current.x=x-fun(x)/dfun(x)

dff=current.x-x

if(abs(dff) < tol) break

x=current.x

#print(x)

}

x     #x=2

 

谢谢了，楼上的大虾。由于我的东西远比您提供的公式麻烦，还有嵌套，所以如果解决不了还得向您请教，先谢过了 呵呵

我也想知道  cox 风险比例模型 用newton-raphson  那些软件可以搞定 我很弱智的  再次谢过了

谢谢楼上

谢谢