从Stata中的nl返回CI下限和上限
在stata中运行nl命令后,如何获得95%置信区间的下限和上限 要获得点估计值,只需执行以下操作:从Stata中的nl返回CI下限和上限,stata,Stata,在stata中运行nl命令后,如何获得95%置信区间的下限和上限 要获得点估计值,只需执行以下操作: global point_est = _b[/beta] 标准误差可以这样保存: global std_err = _se[/beta] 那么,如何获得CI的上下界呢?这不正常,所以: 95%置信区间上限!=点测试+1.96*标准错误 谢谢你的帮助 您错了:nl,因为每个通用估计命令都会产生对称的置信区间。(一些引导置信区间是不对称的,一些自然允许exp、log或logit标度的命令,如lo
global point_est = _b[/beta]
标准误差可以这样保存:
global std_err = _se[/beta]
那么,如何获得CI的上下界呢?这不正常,所以:
95%置信区间上限!=点测试+1.96*标准错误
谢谢你的帮助 您错了:
nl
,因为每个通用估计命令都会产生对称的置信区间。(一些引导置信区间是不对称的,一些自然允许exp、log或logit标度的命令,如logistic
或proportion
,可以在这些标度上生成CI,然后将其反变换为不对称的)。它只需要估计系数和标准误差,并使用您引用的公式生成置信区间。看
此外,不要使用global
宏,除非您必须这样做。这就是程序员所说的。它们唯一的目的是在程序之间传递其他任何有意义的方式都无法传递的内容。带解析的语法系统,加上返回值系统,可以满足大多数需求。相反,globals可能是某种占位符
. 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