作者wheado (principal component QQ)
看板R_Language
标题[问题] Conjugate Gradient Method in optim
时间Mon Nov 14 15:23:44 2016
[问题类型]:
程式谘询(我想用R 做某件事情,但是我不知道要怎麽用R 写出来)
[软体熟悉度]:
入门(写过其他程式,只是对语法不熟悉)
[问题叙述]:
我想自己练习写一个Nonlinear Conjugate Gradient Method,可是我参考文献
http://www.cs.cmu.edu/%7Equake-papers/painless-conjugate-gradient.pdf
中的 "52页 B4-table" 所写的程式码却跟R内建optim找到的不一样,
很明显不是最小值,试了很久还是不懂哪里出问题,我把我的程式贴在下面。
我想写出来的程式是想跟optim一样只需要
Input 函数 起始点 迭代次数 误差
就可以 output 一个解使得函数是local min
我使用numDeriv套件找到该点的 梯度 跟 hessian
[程式范例]:
http://ideone.com/xCSQ8s
[环境叙述]:
R version 3.3.2 (2016-10-31)
[关键字]:
nonlinear conjugate gradient for minimization
感谢大家的指导跟建议
这真的好困难QQ
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 1.200.224.102
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1479108227.A.C2E.html