从Stata表手工计算p值

Question

我想问一个关于如何在没有t-stat表的情况下计算p值的问题，只需查看该表，就像下面的链接http://faculty.arts.ubc.ca/dwhistler/326UBC/stataHILL.pdf中pdf的第一页一样。例如，如果我不知道值0.062，我怎么能通过查看表格中的其他信息来知道它是0.062？

Answer 1

您需要使用ttail()函数，它返回学生的反向累积t分布，也就是概率T> t：

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

第一个参数38是自由度(样本大小减去参数数量)，第二个参数1.92是感兴趣系数的绝对值除以其标准误差，即t-stat。2的因素来自于Stata正在进行双尾测试的事实。您还可以通过使用存储的DoF

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

您还可以使用Adrian Mander的integrate“手动”进行t密度的积分

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

这会得到0.93761229，但您需要的是Pr(T>|t|)，即1-0.93761229=0.06238771。

Answer 2

如果你看了很多统计学教科书，你会发现一个叫做Z表的表，它会给你Z超过你的检验统计量的概率。该表实际上是正态曲线的累积分布函数。

当人们带着四功能计算器去上学时，统计测试中的一个或多个问题将包括这个Z表的副本，亲爱的学生将不得不插入数字列来找到p值。在您的示例中，您将看到.06和.07之间的测试统计数据，这些手指会指出它更接近.06，并进行线性插值以得出.062。

今天，p值是Stata或SAS将为您计算的值。

这是另一个可能感兴趣的SO问题：How do I calculate a p-value if I have the t-statistic and d.f. (in Perl)?

这是一个关于如何“手动”确定p值的基本页面：http://www.dummies.com/how-to/content/how-to-determine-a-pvalue-when-testing-a-null-hypo.html

下面是如何使用Excel确定p值：http://ms-office.wonderhowto.com/how-to/find-p-value-with-excel-346366/

===EDIT===

我的Stata文本(“使用Stata的微观计量经济学”，修订后的Ed，Cameron & Trivedi)在第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