ml 最大似然法出问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

命令：
program MyOLS3
args lfn beta sigma
quietly {
replace ‘lfn’ = ln(normalden($ML_y1, ‘beta’, ‘sigma’))
}
end
ml model lf MyOLS3 (y = eduyear1 experi experi2  nonruralspeci marriage health     middle   west) /sigma, tech(nr 5  bfgs 5 dfp 5  bhhh 5)
ml check
ml  search
ml max
                                从ml  check 开始，出现问题：

Test 1:  Calling MyOLS3 to check if it computes log likelihood and
         does not alter coefficient vector...
       opt__eval_cycle():  3301  subscript invalid
             opt__eval():     -  function returned error
    _optimize_evaluate():     -  function returned error
       _mopt__evaluate():     -  function returned error
     mopt__check_test1():     -  function returned error
      _moptimize_check():     -  function returned error
            Mopt_check():     -  function returned error
                 <istmt>:     -  function returned error
r(3301);

                                   

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：最大似然法 最大似然 coefficient Likelihood EFFICIENT

1. use "http://www.ats.ucla.edu/stat/data/hsbdemo"
   
2. 
   
3. 
   
4. /****************************/
   
5. /*    linear regression     */
   
6. /****************************/
   
7. 
   
8. capture program drop lfols
   
9. program lfols
   
10.   version 10.1
    
11.   args lnf xb lnsigma
    
12.   local y "$ML_y1"
    
13.   quietly replace `lnf' = ln(normalden(`y', `xb',exp(`lnsigma')))
    
14. end
    
15. 
    
16. ml model lf lfols (xb: write = read female) (lnsigma:)
    
17. ml maximize
    
18. display exp([lnsigma]_cons)

复制代码

这是一个极大似然的例子，可以参考下
有一本书专门用Stata估计极大似然的，可以找找，里面有例子

ywh19860616 发表于 2015-7-13 17:31
这是一个极大似然的例子，可以参考下
有一本书专门用Stata估计极大似然的，可以找找，里面有例子

非常感谢你的回复，我用了你的例子来ml 我的数据结果成了， 疑问是：
能提供书名吗？呵呵 或着较好的网页 ，我汗颜 居然没搜到，搜到的那个就是我输入这里的还不能成功，再次感谢！！！

好想~~~ 发表于 2015-7-13 19:28
非常感谢你的回复，我用了你的例子来ml 我的数据结果成了， 疑问是：
能提供书名吗？呵呵 或着较好的网页 ...

Maximum Likelihood Estimation with Stata

sadata.dta (265.38 KB)

2015-7-14 11:16:43 上传

switchr 例子的数据文件

ywh19860616 发表于 2015-7-14 10:13
Maximum Likelihood Estimation with Stata

我找到书看了点。。。模模糊糊的。

好想~~~ 发表于 2015-7-14 11:13
sadata.dta
非常感谢你的帮助！！
movestay 是内生转换回归模型，即outcome eq中含有内生的虚拟变量， ...

你还是在想做你之前问的那个模型吗？你是做什么题目呢？

夏目贵志 发表于 2015-7-15 08:05
你还是在想做你之前问的那个模型吗？你是做什么题目呢？

我自己弄晕了 。这几天的新理解是:movestay中的regime项必须是所有外生变量加上工具变量或者是必须所有外生变量，而我的情况是在regime 中会在所有外生中选择一些，也许不能包括所有外生，这样导致该命令我不能使用。switchr 命令允许不是所有的外生变量 ，但是也处了问题，我将发帖求教。

好想~~~ 发表于 2015-7-15 11:59
我自己弄晕了 。这几天的新理解是:movestay中的regime项必须是所有外生变量加上工具变量或者是必须所有外 ...

如果你确定想做你之前问的那个模型，最好的方法就是自己写likelihood function。