2013-12-06 32 views
4

我想用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. 

任何想法來解決這個問題?

回答

1

我不知道這是不是你的問題的根源,但同時在成功和不成功的nl的說法,估計的標準偏差,exp([lnsigma]_b[_cons]),只有在第二個參數劃分只有最後[Tier*]_b[race]詞,然後binormal(這可能是因爲你在第一個參數中將SD約束爲1)。這裏猜測你在第二個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 
+0

謝謝Steve。但是,那不是問題。該公式是'binormal(x1b *,x2b */sigma,rho)'。所以,我的表情是正確的。 – Metrics

+0

是的,但在你的問題中寫的第二個參數的形式是_k_ y2 * + _d_(y2 */sigma):只有部分被sigma分開。 –

+0

真棒史蒂夫。現在,問題通過更正解決。非常感謝。 – Metrics

相關問題