全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 Stata专版
2249 1
2011-10-17
用ml command对正态分布中的mu 和 sigma平方 做最大似然估计
用的是ml model if method
定义了 log density function
语句如下:
program drop _all                                                                                                                
program mynormal                                                                                                                     
args lnf mu sigma                                                                                                       
quietly replace `lnf' = log((1/`sigma')*normalden(($ML_y1 - `mu')/`sigma'))
end


结果在运行do file时显示说
1/ invalid name

请教各位大牛该如何修正? 最神奇的是 同样的语句在老师给的范例中是运行正常的

小弟在此拜谢!
二维码

扫码加我 拉你入群

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

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

全部回复
2011-10-18 06:38:58
*设待检验的变量是x。运行以下:

prog drop _all                                                                                                                 
prog mynormal                                                                                                                     
args lnf mu sigma                                                                                                        
qui replace `lnf' = log((1/`sigma')*normalden(($ML_y1-`mu')/`sigma'))
end

ml mod lf mynormal (x=)(x=)
ml max


*另,end之前一步可简化为:
qui replace `lnf'=ln(normalden($ML_y1,`mu',`sigma'))
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群