Tobit回帰を実行するときに、Stataにp値の3桁以上を表示させる簡単な方法を知っている人はいますか?
通常、Stataはp値がまたはであると報告しますが、たとえば、または.001
。.065
をもっと多くの桁で表示したいと思います。.0011123
.065320
明確にするために、回帰テーブルでのデータの生成方法を(必然的に)変更したくありません。
興味のあるp値に対してStataがより多くの桁を表示できるようにしたいだけです。
Tobit回帰を実行するときに、Stataにp値の3桁以上を表示させる簡単な方法を知っている人はいますか?
通常、Stataはp値がまたはであると報告しますが、たとえば、または.001
。.065
をもっと多くの桁で表示したいと思います。.0011123
.065320
明確にするために、回帰テーブルでのデータの生成方法を(必然的に)変更したくありません。
興味のあるp値に対してStataがより多くの桁を表示できるようにしたいだけです。
tobitコマンドを次のようにフォローアップします
est tab, p(%12.10g)
(例えば)。これは、かなり古いバージョンのStataでも機能するはずです。少し簡単ではありませんが、独自の出力プロシージャを作成することです。
Stata 11.1ではset pformat
、係数テーブルのp値の出力形式を指定するコマンドが導入されました。( STATAについてはわかりませんが、1980年代に廃止されたと思います)。
ちなみに、ここではなく、Statalistでこのような完全にStata固有の質問をする方がよいでしょう。
多くの場合、内部名でp値を知っていれば、最高の精度を得ることができます。私は通常、真剣に使用するほぼすべてのコマンドを入力するreturn list
か、ereturn list
その後に入力してから、p値を含むスカラーのように見える結果を取得します。e(p)
r(p)
e(p_chi2)
トービット回帰の後、コマンドを使用してtest
、帰無仮説x1=0からp値を取得できます。
sysuse auto
tobit weight trunk length headroom, ll(1500)
test trunk
結果はr(p)
return list
tobit
コマンドヘルプファイルの最初の例を使用すると、次のようになります。
. sysuse auto, clear
. generate wgt = weight / 1000
. tobit mpg wgt, ll(17)
Tobit regression Number of obs = 74
LR chi2(1) = 72.85
Prob > chi2 = 0.0000
Log likelihood = -164.25438 Pseudo R2 = 0.1815
------------------------------------------------------------------------------
mpg | Coef. Std. Err. t P>|t| [95% Conf. Interval]
-------------+----------------------------------------------------------------
wgt | -6.87305 .7002559 -9.82 0.000 -8.268658 -5.477442
_cons | 41.49856 2.05838 20.16 0.000 37.39621 45.6009
-------------+----------------------------------------------------------------
/sigma | 3.845701 .3663309 3.115605 4.575797
------------------------------------------------------------------------------
18 left-censored observations at mpg <= 17
56 uncensored observations
0 right-censored observations
次の結果から返される結果から任意のp値を簡単に取得できますr()
。
. matrix list r(table)
r(table)[9,3]
model: model: sigma:
wgt _cons _cons
b -6.8730504 41.498557 3.8457011
se .70025591 2.0583803 .36633085
t -9.8150552 20.160782 .b
pvalue 5.610e-15 1.471e-31 .b
ll -8.2686584 37.396211 3.1156048
ul -5.4774424 45.600903 4.5757975
df 73 73 73
crit 1.9929971 1.9929971 1.9929971
eform 0 0 0
そして、それに応じてフォーマットします。
. matrix results = r(table)
. display %18.17f results[4,1]
0.00000000000000561
詳細については、Stataのコマンドプロンプトhelp format
から入力してください。