我试图nlcom
在有限因变量模型中使用Stata(Delta方法)估计边际效应的标准误差。在计算边际效应时,我们通常计算概率,条件和无条件部分。我也能够nlcom
在Stata中使用的有条件和无条件零件都获得标准误差。但是,在使用nlcom
概率时,出现了错误“超出了最大迭代次数” 。我确实在这里查看了解决方案,并且确实解决了一些二进制变量(3)。但是,对于4个变量,我仍然遇到相同的错误。我深入研究问题,发现错误与伪变量为0的部分有关:
\ begin {equation} \ Delta P(y> 0)= P(y> 0 | x = 1)-P(y> 0 | x = 0)\ end {equation}
这是竞速变量(虚拟变量)的示例代码,nlcom
即使考虑了前面提到的网站中建议的改正,虚拟变量为1(无校正)但不适用于0的零件也适用。增加迭代次数并乘以更大的幅度也无济于事:
_nl_1: binormal(1.733373263151318-(.0438395852394023*[Tier1]_b[race])+(1*[Tier1]_b[race]),7.575546508411257-(.0438395852394023*[Tier2]_b[race])+
> (1*[Tier2]_b[race])/exp([lnsigma]_b[_cons]),tanh([tau]_b[_cons]))
------------------------------------------------------------------------------
| Coef. Std. Err. z P>|z| [95% Conf. Interval]
-------------+----------------------------------------------------------------
_nl_1 | .9223902 .0129743 71.09 0.000 .8969611 .9478193
------------------------------------------------------------------------------
_nl_1: binormal(1.733373263151318-(.0438395852394023*[Tier1]_b[race])+(0*[Tier1]_b[race]),7.575546508411257-(.0438395852394023*[Tier2]_b[race])+
> (0*[Tier2]_b[race])/exp([lnsigma]_b[_cons]),tanh([tau]_b[_cons]))*(1e+2)
Maximum number of iterations exceeded.
有解决这个问题的主意吗?
我不知道这是否是您问题的根源,但是在成功和失败的nl
陈述中,估计的标准偏差exp([lnsigma]_b[_cons])
只会除最后[Tier*]_b[race]
一项,然后仅除(的第二个参数),binormal
这可能是因为您限制了SD1
位于第一个参数中)。这是您对第二个nl
语句中的意思的猜测,但这仅仅是一个猜测。
#delim ;
_nl_1:
binormal
((1.733373263151318-(.0438395852394023*[Tier1]_b[race])+(0*[Tier1]_b[race]))
/exp([lnsigma]_b[_cons]),
(7.575546508411257-(.0438395852394023*[Tier2]_b[race])+ (0*[Tier2]_b[race]))
/exp([lnsigma]_b[_cons]),
tanh([tau]_b[_cons])) ;
#delim cr
本文收集自互联网,转载请注明来源。
如有侵权,请联系[email protected] 删除。
我来说两句