用optim函数拟合的结果

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

$par
[1]  0.170097755 -0.005719282  0.005023195 -0.055810510  0.208030681 -0.005023231  0.055810358 -0.208030743
$value
[1] -897069.4
$counts
function gradient
     121       14
$convergence
[1] 0
$message
NULL



请问下面三个
$counts
function gradient
     121       14
$convergence
[1] 0
$message
NULL
分别是什么意思啊？急用啊

扫码加我 拉你入群

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

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

分享1 收藏0 回帖

关键词：Optim OPT Tim Convergence function 函数

convergence = 0 说明收敛了

counts代表调用了fn和gr函数各多少次

message是optim函数返回的一些信息

qoiqpwqr 发表于 2013-5-26 23:09
counts代表调用了fn和gr函数各多少次

message是optim函数返回的一些信息

接版大的话小的来补充下吧，
对于counts
fn某种程度上意味着算了几次收敛，因为一般情况下每次迭代都会用到fn
梯度函数的话就不一定了，看你的算法了，一般用的比较多的拟牛顿类和共轭梯度类都会要用到，但肯定会比fn少很多，当然像单纯型这种就会出现NULL(根本不需要)，就是function的次数会比较那啥。。。
convergence等于恭喜lz，其他的一些情况我也来给lz罗列下吧
当然他一般都是整数
输出1表示达到了你所设置的最大收敛次数，lz可通过control参数设置，默认没记错的话500次或者50次吧忘了
输出10只有当你选择单纯型算法才会碰到（不是线性规划里的单纯形法，收敛速度略慢。。。慎用。。。）
输出51说明使用边界约束BFGS算法时可能出现了问题（警告，出现错误会输出52）
希望对lz有用

qoiqpwqr 发表于 2013-5-26 23:09
counts代表调用了fn和gr函数各多少次

message是optim函数返回的一些信息

不好意思，我再请问一下，这个我是用的optim函数，如下：
res<-optim(params,LL,method="BFGS",hessian=TRUE,control=list(fnscale=-1))
请问后面那个control里面设置的是什么意思？，还有我这个optim是求的最大值还是最小值？

求证1加1 发表于 2013-5-26 23:39
接版大的话小的来补充下吧，
对于counts
fn某种程度上意味着算了几次收敛，因为一般情况下每次迭代都会 ...

不好意思，我再请问一下，这个我是用的optim函数，如下：
res<-optim(params,LL,method="BFGS",hessian=TRUE,control=list(fnscale=-1))
请问后面那个control里面设置的是什么意思？，还有我这个optim是求的最大值还是最小值？

pjq521 发表于 2013-5-28 14:20
不好意思，我再请问一下，这个我是用的optim函数，如下：
res

https://bbs.pinggu.org/forum.php? ... ;page=1#pid19138430

求证1加1 发表于 2013-5-28 14:23
https://bbs.pinggu.org/forum.php?mod=viewthread&tid=2437141&page=1#pid19138430

又来问题了，我这个运行出来的结果为什么会使这样
> f=log(1+exp(-(params[1]*x1+params[2]*x2+params[3]*x3+params[4]*x4+params[5]*x5-params[6]*y1-params[7]*y2-params[8]*y3+params[9])))
> f
numeric(0)

pjq521 发表于 2013-5-28 14:52
又来问题了，我这个运行出来的结果为什么会使这样
> f=log(1+exp(-(params[1]*x1+params[2]*x2+params[3 ...

再怎么样也不可能等于0啊，这是为什么啊

你得先定义函数，不能直接这么定义的。。。这不是matlab也不是maple。。。。