Stata 从ttest命令进行引导估计的语法

我试图证明小样本数据的各种测试的特征。我想演示一下t-test、带bootstrap估计的t-test和ranksum测试的性能。我对使用

simulate

获得多组数据的每次测试的p值感兴趣。但是，我无法使用

bootstrap

前缀和

ttest

命令获得t检验估计值

数据由以下内容生成：

clear
set obs 60
gen level = abs(rnormal(0,1))
gen group = "A"
replace group = "B" if [_n] >30

bootstrap, reps(100): ttest level, by(group)
bootstrap _b, reps(100): ttest level, by(group)
bootstrap boot_p = e(p), reps(100): ttest level, by(group)

每个程序的错误顺序如下：

需要表达式列表

无效的表达式：\u b

“e（p）”在完整样本中评估为缺失

这些结果与

引导

前缀的文档不一致。e或r类对象的规范和测试是否有问题

编辑：

既然理解了r-class是正确的标量组，我仍然没有根据解决方案中提供的代码生成变量“p”。此外：

clear
set more off

set obs 60

gen level = abs(rnormal(0,1))
gen group = "A"
replace group = "B" if [_n] >30

bootstrap p=r(p), reps(100): ttest level, by(group)
    display r(p)

不返回p值

ttest

是一个r类命令，将其REULT存储在

r（）

中。您似乎期望它将结果保存在

e（）

中，就像e-class命令一样。标准是后者符合模型<代码>测试不在此类别中

双面p值存储在

r（p）

中，如

帮助测试中所示：

clear
set more off

set obs 60

gen level = abs(rnormal(0,1))
gen group = "A"
replace group = "B" if [_n] >30

bootstrap p=r(p), reps(100): ttest level, by(group)

此解决方案的工作方式与上面编辑的查询中所述不同。
bootstrap
有自己的帮助文件，您可以通过
help bootstrap
轻松访问它。在那里，您可以发现引导结果在
e（b）
中。运行
矩阵列表e（b）
以显示它。关于创建变量的愿望，如果需要中间结果，请检查
引导的
选项
saving（）
。如果要保存最终结果，您现在知道如何访问它。有关将结果保存到Stata数据集的信息，请参见帮助postfile。我的建议是更好地使用
帮助文件
和Stata手册。（您可以使任何答案不相关，只需修改您的问题即可。）
返回列表
在
引导后
不包含
r（p）
，\将显示
r（p）
不是由
引导返回的，因此您的最后一个命令将产生错误。