我想用Stata(Delta方法)中的nlcom
來估計邊際效應的標準誤差,用於受限因變量模型。我們通常在計算邊際效應時計算概率,有條件和無條件的部分。我還能夠在Stata中使用nlcom
獲得有條件和無條件部分的標準錯誤。然而,當使用nlcom
作爲概率時,我得到了錯誤「超過最大迭代次數」。我看過解決方案here,它確實解決了一些二進制變量(3)。但是,我仍然得到4個變量的相同錯誤。我研究這個問題,並發現,在僞變量爲0錯誤與部分相關聯:錯誤:在Stata中使用nlcom的「超過最大迭代次數」
\開始{方程} \德爾塔P(Y > 0)= P(Y > 0 | X = 1) -P(Y > 0 | x = 0時) \ {端方程}
下面是種族變量(虛擬變量)該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.
任何想法來解決這個問題?
謝謝Steve。但是,那不是問題。該公式是'binormal(x1b *,x2b */sigma,rho)'。所以,我的表情是正確的。 – Metrics
是的,但在你的問題中寫的第二個參數的形式是_k_ y2 * + _d_(y2 */sigma):只有部分被sigma分開。 –
真棒史蒂夫。現在,問題通過更正解決。非常感謝。 – Metrics