悬赏边界logistic的运算程序或者操作命令

求边界logistic的运算程序或者操作命令，我将送与100论坛币，聊表谢意。最好使用Eviews、或matlab、或 SAS、或stata

模型的表达式为：y =c(1)/(1+exp(-c(2)-c(3)*zx5-c(4)*zx8-c(5)*zx10))
请采用最大似然估计去估计模型的参数，
数据及参考文献见下附件，如果能运算出来，请与我联系,或者通过论坛和我联系也可以。我一定表达我的谢意。

y77 发表于 2011-5-2 14:11 我把我的数据也更新为.dta格式，使用您提供的程序在运行时，出现如下错误，请帮助

*你的因变量y是字符型而非实数型，运行以下可得结果：

use boundlogit,clear
g z=real(y)
cap pr drop blogitml
pr blogitml
args lnl xb w
qui replace `lnl'=ln(`w'/(1+exp(-`xb')))*($ML_y1==1)+ln(1-`w'/(1+exp(-`xb')))*($ML_y1==0)
end
ml mod lf blogitml (xb:z=x1-x3)(w:z=)
ml max

y77 发表于 2011-4-29 22:01 求边界logistic的运算程序或者操作命令，我将送与100论坛币，聊表谢意。最好使用Eviews、或matlab、或 SAS、或stata
模型的表达式为：y =c(1)/(1+exp(-c(2)-c(3)*zx5-c(4)*zx8-c(5)*zx10))
请采用最大似然估计去估计模型的参数

*stata程序：设y是哑元因变量，（非截距）自变量x1-x5。

cap pr drop blogitml
pr blogitml
args lnl xb w
qui replace `lnl'=ln(`w'/(1+exp(-`xb')))*($ML_y1==1)+ln(1-`w'/(1+exp(-`xb')))*($ML_y1==0)
end
ml mod lf blogitml (xb:y=x1-x5)(w:y=)
ml max

非常感谢sungmoo，我把我的数据也更新为.dta格式，使用您提供的程序在运行时，出现如下错误，请帮助。非常感谢！
cap pr drop blogitml
.
. pr blogitml
  1.
. args lnl xb w
  2.
. qui replace `lnl'=ln(`w'/(1+exp(-`xb')))*($ML_y1==1)+ln(1-`w'/(1+exp(-`xb')))*(
> $ML_y1==0)
  3.
. end
.
. ml mod lf blogitml (xb:y=x1-x3)(w:y=)
no observations
r(2000);
.
. ml max
you must issue -ml model- first
r(198);
.
.
.

非常感谢！我已经将您的答案设置为最佳答案，再次谢谢您！

boundlogit及satata实现,再次感谢sungmoo的帮助！

非常感谢上面两位朋友的帮助。。。谢谢！

这个测试数据是不收敛？