用optim函数拟合的结果

14672

收藏 2013-05-26

$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
分别是什么意思啊？急用啊

扫码加我拉你入群

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

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

全部回复

ntsean

2013-5-26 22:25:15

convergence = 0 说明收敛了

扫码加我拉你入群

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

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

qoiqpwqr

2013-5-26 23:09:22

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

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

扫码加我拉你入群

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

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

求证1加1

2013-5-26 23:39:50

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有用

扫码加我拉你入群

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

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

pjq521

2013-5-28 14:18:26

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是求的最大值还是最小值？

扫码加我拉你入群

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

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

pjq521

2013-5-28 14:20:44

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

扫码加我拉你入群

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

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

点击查看更多内容…

求证1加1

2013-5-28 14:23:33

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

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

扫码加我拉你入群

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

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

pjq521

2013-5-28 14:52:24

求证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:53

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

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

扫码加我拉你入群

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

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

求证1加1

2013-5-28 14:53:49

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

扫码加我拉你入群

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

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

pjq521

2013-5-28 15:00:02

求证1加1 发表于 2013-5-28 14:53
你得先定义函数，不能直接这么定义的。。。这不是matlab也不是maple。。。。

我把我程序发给你看看吧，

复制代码

这是我的程序，但是optim就是不给我迭代，为什么

扫码加我拉你入群

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

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

pjq521

2013-5-28 15:00:38

pjq521 发表于 2013-5-28 15:00
我把我程序发给你看看吧，这是我的程序，但是optim就是不给我迭代，为什么

> cat("rohbar=",(L0-(LL-length(b)))/L0)
rohbar= 0.9764778> print(res)
$par
[1] 0 0 0 0 0 0 0 0 0

$value
[1] 0

$counts
function gradient
1 1

扫码加我拉你入群

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

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

求证1加1

2013-5-28 15:14:54

你再看看有没有输错的，还有就是我不明白为什么要在函数输入里面加d，感觉没这必要，后面没有必要return，直接sum(f)就行了，再者就是你再看看你的数据

扫码加我拉你入群

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

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

pjq521

2013-5-28 16:17:24

求证1加1 发表于 2013-5-28 15:14
你再看看有没有输错的，还有就是我不明白为什么要在函数输入里面加d，感觉没这必要，后面没有必要return，直 ...

误于optim(params, LL, method = "BFGS", hessian = TRUE, control = list(fnscale = -1)) :
non-finite finite-difference value [2]
确实是我数据的格式不对，但是他又出来了个错误，你能告诉我这是什么意思么

扫码加我拉你入群

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

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

qinglingxuanfen

2014-11-14 16:28:04

求助！伽马参数估计出现这情况怎么处理？
> optim(c(1,2),logL,dat=dat)
Error in optim(c(1, 2), logL, dat = dat) :
function cannot be evaluated at initial parameters
In addition: Warning message:
In lapply(X = x, FUN = .Generic, ...) : NaNs produced

扫码加我拉你入群

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

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

qinglingxuanfen

2014-11-14 16:31:18

pjq521 发表于 2013-5-28 16:17
误于optim(params, LL, method = "BFGS", hessian = TRUE, control = list(fnscale = -1)) :
non-fin ...

> optim(c(1,2),logL,dat=dat)
Error in optim(c(1, 2), logL, dat = dat) :
function cannot be evaluated at initial parameters
In addition: Warning message:
In lapply(X = x, FUN = .Generic, ...) : NaNs produced
我这个怎么处理啊求助大神！

扫码加我拉你入群

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

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

dw928762571

2017-10-30 09:46:08

感觉这个对非凸函数没啥用

扫码加我拉你入群

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

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

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群