由买买提看人间百态

boards

本页内容为未名空间相应帖子的节选和存档,一周内的贴子最多显示50字,超过一周显示500字 访问原贴
Statistics版 - 请教matlab non-linear optimization的问题
相关主题
MCMC算法的Posterior Std. 一定是减少的么?在线急等:做maximum likelihood estimation,用optimization怎么都的不出正确的数值
res log pseudo-likelihood function用R的optim Hessian=T option得到负的varaince正常吗
R help.折腾一天,解决不了。求教: R function
请问proc genmod, sas问题请教 C++ 里面 optmization 函数及其用法
any one did EM to find MLE of mixed effects model in R请问如何在R中写recursive function?
问一个R问题再问multinomial logit
为什么有些函数在我的R中无法执行呢? R help!!!
{包子请教} 如何决定用什么样的link function?一般如何证明MCMC算法得到的sample是来自targeted distribution的呢?
相关话题的讨论汇总
话题: iteration话题: function话题: dfunc话题: value话题: maxlik
进入Statistics版参与讨论
1 (共1页)
D******6
发帖数: 6211
1
我需要用一个non-linear optimization function来estimate一个model,近似于logit
model,用maximum likelihood estimation。
如果用maxlik,其实这是一个min函数,就是名字叫max,从这个网站下载的,http://www.spatial-econometrics.com/
结果有点奇怪,不是指我的model参数奇怪,而是maxlik输出的结果如下:
CONVERGENCE CRITERIA MET: Change in Objective Function
iteration function value dfunc
0 2301.06004609 1000.00000000
iteration function value dfunc
1 1966.10277862 0.17036610
iteration function value dfunc
2 1965.50818039 0.00030252
iteration function value dfunc
3 1967.81330103 -0.00117141
1.如果用结果result.f查看,就给最后一个iteration的function value,可是明明是
iteration 2的function value才是最小的啊。不止这一个,我随便运行了很多maxlik
,它总是给最后一个iteration作为optimal,可是,有的时候很明显第一个iteration
的value是optimum啊。我糊涂了。是不是我解读错了输出结果?
2.我明白dfunc value接近0是极值,因为那时一阶导数吧。那么如果存在一个情况是前
面的iteration的function value很小,可是导数还很大,说明在往上走,那么即使后
来走到了一个极值,也该回头比比看看是不是最小吧,算法不会这么差吧,一点不回头?
如有解答,十分感谢!
1 (共1页)
进入Statistics版参与讨论
相关主题
一般如何证明MCMC算法得到的sample是来自targeted distribution的呢?any one did EM to find MLE of mixed effects model in R
SAS问题,关于@和@@的区别问一个R问题
请教Base 70 中的一题 关于@为什么有些函数在我的R中无法执行呢?
提高R速度的一些tips{包子请教} 如何决定用什么样的link function?
MCMC算法的Posterior Std. 一定是减少的么?在线急等:做maximum likelihood estimation,用optimization怎么都的不出正确的数值
res log pseudo-likelihood function用R的optim Hessian=T option得到负的varaince正常吗
R help.折腾一天,解决不了。求教: R function
请问proc genmod, sas问题请教 C++ 里面 optmization 函数及其用法
相关话题的讨论汇总
话题: iteration话题: function话题: dfunc话题: value话题: maxlik