D******6
发帖数: 6211 | 1 【 以下文字转载自 Statistics 讨论区 】
发信人: Doha2006 (花猫), 信区: Statistics
标 题: 请教matlab non-linear optimization的问题
发信站: BBS 未名空间站 (Fri Jan 7 23:39:25 2011, 美东)
我需要用一个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很小，可是导数还很大，说明在往上走，那么即使后
来走到了一个极值，也该回头比比看看是不是最小吧，算法不会这么差吧，一点不回头？
如有解答，十分感谢！ |
w********e
发帖数: 962 | |
D******6
发帖数: 6211 | 3 我。。。又一个人虚晃我一枪。。。
我盼老大给我点提示，以前老大好像什么stata sas都用过的。
【在 w********e 的大作中提到】
: 这能闻能捂的。。。哈。
|
f*****s
发帖数: 3868 | 4 看来小二给人的印象就是无所不能的 ....
小二, 到底是谁杀的肯尼迪?
【在 D******6 的大作中提到】
: 我。。。又一个人虚晃我一枪。。。
: 我盼老大给我点提示，以前老大好像什么stata sas都用过的。
|
r****y
发帖数: 26819 | 5 豪门问的问题就是豪门
换我肯定不问这个
【在 f*****s 的大作中提到】
: 看来小二给人的印象就是无所不能的 ....
: 小二, 到底是谁杀的肯尼迪?
|
f*****s
发帖数: 3868 | 6 其实我特想问空姐到底是多少尺码
【在 r****y 的大作中提到】
: 豪门问的问题就是豪门
: 换我肯定不问这个
|
r****y
发帖数: 26819 | 7 口是心非也
【在 f*****s 的大作中提到】
: 其实我特想问空姐到底是多少尺码
|
f*****s
发帖数: 3868 | 8 诛心啊.....
要不你问一个我看看?
【在 r****y 的大作中提到】
: 口是心非也
|
r****y
发帖数: 26819 | 9 豪门空姐你都问了，我只好问小二将来是嫁女给你家呢还是从你家娶儿媳了。
【在 f*****s 的大作中提到】
: 诛心啊.....
: 要不你问一个我看看?
|
f*****s
发帖数: 3868 | 10 这得看他生什么出来啊
【在 r****y 的大作中提到】
: 豪门空姐你都问了，我只好问小二将来是嫁女给你家呢还是从你家娶儿媳了。
|
|
|
r****y
发帖数: 26819 | 11 有话头。你这边呢？:-)
【在 f*****s 的大作中提到】
: 这得看他生什么出来啊
|
f*****s
发帖数: 3868 | 12 他那边还没动静呢, 等他那里出动静了, 我这里可能早就品种齐全了也没准啊:)
【在 r****y 的大作中提到】
: 有话头。你这边呢？:-)
|
r****y
发帖数: 26819 | 13 哈哈！赞，传说中的“总有一款适合你”。。。
【在 f*****s 的大作中提到】
: 他那边还没动静呢, 等他那里出动静了, 我这里可能早就品种齐全了也没准啊:)
|
f*****s
发帖数: 3868 | 14 有女娃不能嫁给小二家 .... 会给娃压力的
除非他真的在墙街挖到了好多钱 ....
【在 r****y 的大作中提到】
: 哈哈！赞，传说中的“总有一款适合你”。。。
|
r****y
发帖数: 26819 | 15 哈哈，咋知小二不能有不肖之子，或者你家女娃不听父母之言。
【在 f*****s 的大作中提到】
: 有女娃不能嫁给小二家 .... 会给娃压力的
: 除非他真的在墙街挖到了好多钱 ....
|
f*****s
发帖数: 3868 | 16 压根就不让见, 哈哈
【在 r****y 的大作中提到】
: 哈哈，咋知小二不能有不肖之子，或者你家女娃不听父母之言。
|
r****y
发帖数: 26819 | 17 狠。豪门深闺，金砖砸开。
【在 f*****s 的大作中提到】
: 压根就不让见, 哈哈
|
f*****s
发帖数: 3868 | 18 恩, 就要看金砖够不够厚了
【在 r****y 的大作中提到】
: 狠。豪门深闺，金砖砸开。
|
D******6
发帖数: 6211 | 19 玩了，我说话得罪人了，高手知道也不告诉我了。。。
【在 f*****s 的大作中提到】
: 看来小二给人的印象就是无所不能的 ....
: 小二, 到底是谁杀的肯尼迪?
|
D******6
发帖数: 6211 | 20 你以前不也是解答过我sas的问题，这次也帮帮忙？要多大的包子？
【在 r****y 的大作中提到】
: 豪门问的问题就是豪门
: 换我肯定不问这个
|
