StataのnlからCIの下限と上限を返します

Question

stataでnlコマンドを実行した後、どのようにして95％Confidence Intervalの下限と上限を得ることができますか？

点推定値を得るために、あなただけの操作を行います。

global point_est = _b[/beta] 

と標準誤差は次のように保存することができます。上下のCIの境界を取得するにはどのように

global std_err = _se[/beta] 

？そう、その正常ではない：

95％CI-アッパー= point_est + 1.96 * std_err任意の助け

ありがとう！！

Answer 1

あなたは間違っています：nlは、すべての一般的な推定コマンドとして、対称信頼区間を生成します。 （ブートストラップ信頼区間の一部は非対称であり、logisticやproportionのようなexp、log、logitスケールを自然に許すコマンドでは、これらのスケールでCIを生成し、非対称なものに逆変換することができます）係数と標準誤差を計算し、引用した数式を使って信頼区間を作成します。 help estcomを参照してください。

. sysuse auto, clear 
. nl (price = exp({a} + {b}*mpg)) 

------------------------------------------------------------------------------ 
     price |  Coef. Std. Err.  t P>|t|  [95% Conf. Interval] 
-------------+---------------------------------------------------------------- 
      /a | 9.882335 .2007597 49.22 0.000  9.482128 10.28254 
      /b | -.0569546 .0105899 -5.38 0.000 -.0780651 -.035844 
------------------------------------------------------------------------------ 

. di reldif(_b[/a] - invt(e(df_r),0.975)*_se[/a], 9.482128) 
2.582e-09 

. di reldif(_b[/a] + invt(e(df_r),0.975)*_se[/a], 10.28254) 
1.944e-07 

また、絶対にする必要がない限り、globalマクロを使用しないでください。これはプログラマがa "smell"と呼ぶものです。彼らの唯一の目的は、他の意味のある方法で渡すことができないプログラム間で物を渡すことです。構文解析を伴う構文システムと戻り値のシステムは、ほとんどのニーズを処理します。代わりに、グローバルはおそらく何らかの種類のプレースホルダです。