stata - Stata テーブルから手動で p 値を計算する

Question

次のリンクhttp://faculty.arts.ubc. ca/dwhistler/326UBC/stataHILL.pdf . たとえば、値 0.062 がわからない場合、テーブルの他の情報を見て、それが 0.062 であることをどのように知ることができますか?

Answer 1

ttail()逆累積スチューデントの t 分布、別名確率 T > t を返す関数を使用する必要があります。

display ttail(38,abs(_b[_cons]/_se[_cons]))*2

最初の引数 38 は自由度 (サンプル サイズからパラメーターの数を引いたもの) で、2 番目の引数 1.92 は対象係数の絶対値を標準誤差で割った値、つまり t-stat です。係数 2 は、Stata が両側検定を行っているという事実に由来します。保存された DoF を次のように使用することもできます

display ttail(e(df_r),abs(_b[_cons]/_se[_cons]))*2

Adrian Mander の を使用して、「手」で t 密度の積分を行うこともできますintegrate。

ssc install integrate
integrate, f(tden(38,x)) l(-1.92) u(1.92)

これにより 0.93761229 が得られますが、1-0.93761229=0.06238771 である Pr(T>|t|) が必要です。

Answer 2

多くの統計の教科書を見ると、Z テーブルと呼ばれるテーブルが見つかります。これは、Z が検定統計量を超える確率を示します。表は、実際には正規曲線の累積分布関数です。

人々が 4 関数電卓を持って学校に通うと、統計テストの 1 つまたは複数の問題にこの Z テーブルのコピーが含まれ、親愛なる生徒は数値の列を補間して p 値を見つける必要がありました。あなたの例では、.06 と .07 の間のテスト統計が表示され、それらの指はそれが .06 に近いことをタップし、線形補間を行って .062 を見つけます。

現在、p 値は、Stata または SAS が計算するものです。

興味深いかもしれない別の SO の質問があります: t 統計量と df (Perl で) がある場合、p 値を計算するにはどうすればよいですか?

「手で」p値を決定する方法に関する基本的なページは次のとおりです。 -hypo.html

Excel を使用して p 値を決定する方法は次のとおりです: http://ms-office.wonderhowto.com/how-to/find-p-value-with-excel-346366/

===編集===

私の Stata テキスト ("Microeconometrics using Stata", Revised Ed, Cameron & Trivedi) は、p. 402。

* p-values for t(30), F(1,30), Z, and chi(1) at y=2
. scalar y=2
. scalar p_t30 = 2 * ttail(30,y)
. scalar p_f1and30 = Ftail(1,30,y^2)
. scalar p_z = 2 * (1 - normal(y))
. scalar p_chi1 = chi2tail(1,y^2)
. display "p-values" " t(30)=" %7.4f p_t30
p-values t(30) = 0.0546