如何计算Stata中单尾检验的p值?
我有以下模型:ln(MPG{I})=\beta{0}+\beta{1}权重{I}+\beta{1}外来{I}+\beta{3}外来{I}*权重{I}+\varepsilon{I,j} 我想使用test命令来测试$\beta_{3}>0.5$in STATA上的系数 我使用了以下代码并获得了此结果: 测试1.外来#c.重量=0.5001 (1)1.国外#c.重量=0.5001如何计算Stata中单尾检验的p值?,stata,Stata,我有以下模型:ln(MPG{I})=\beta{0}+\beta{1}权重{I}+\beta{1}外来{I}+\beta{3}外来{I}*权重{I}+\varepsilon{I,j} 我想使用test命令来测试$\beta_{3}>0.5$in STATA上的系数 我使用了以下代码并获得了此结果: 测试1.外来#c.重量=0.5001 (1)1.国外#c.重量=0.5001 F( 1, 70) = 4.9e+07 Prob > F = 0.0000 所
F( 1, 70) = 4.9e+07
Prob > F = 0.0000
所以我们拒绝空值,因为p值非常小
但问题是,这是一个双尾测试
我的目标是得到这个左尾检验的t检验值,然后存储它。然后用t检验来计算它的p值
在计算p值之后,我决定是否拒绝null。我确信我会拒绝空值,并且p值会非常小。只是需要一些帮助,弄清楚如何用正确的方式编写代码
编辑:我已尝试使用以下命令:
lincom _b[forweight] - 0.5
display invttail(71, 0.5)
最后一个命令输出的值为0。这是左侧t检验的p值吗 本节将对此进行介绍
以下是测试系数是否为0.5的相关代码:
sysuse auto, clear
gen ln_mpg = ln(mpg)
regress ln_mpg i.foreign##c.weight
test _b[1.foreign#c.weight]=.5
local sign_wgt = sign(1.foreign#c.weight)
display "Ho: coef <= 0.5 p-value = " ttail(r(df_r),`sign_wgt'*sqrt(r(F)))
display "Ho: coef >= 0.5 p-value = " 1-ttail(r(df_r),`sign_wgt'*sqrt(r(F)))
sysuse自动,清除
发电机ln\u mpg=ln(mpg)
回归ln#mpg i.国外#c.重量
测试_b[1.外来#c.重量]=0.5
本地标志_wgt=标志(1.外来#c.重量)
显示“Ho:coef=0.5 p-value=”1-ttail(r(df_r),`sign_wgt'*sqrt(r(F)))
您是否已经知道如何执行此测试,并且只想学习Stata命令来执行此测试,或者您的问题可能与不知道要使用什么测试、如何计算或如何解释它有关?你能具体说明你需要什么样的帮助吗?@whuber这是第一部分。我知道这个测试,但我想学习Stata命令来完成它。