Question

我试图展示小数据样本的各种测试的特征。我想展示t-test，t-test和bootstrap估计以及rankum测试的性能。我有兴趣使用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)

按顺序执行的每个程序的错误是：

需要表达式列表
表达式无效：_b
＆＃39; E（P）＆＃39;评估为完整样本中缺失

这些结果与bootstrap前缀的文档不一致。规范e或r类对象和ttest是否存在问题？

编辑：

现在了解r-class是正确的标量组，我仍然不会生成变量＆＃39; p＆＃39;给出解决方案中提供的代码。此外：

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值。

Answer 1

ttest是 r-class 命令，并将其结果存储在r()中。您似乎期望它将结果保存在e()中，就像 e-class 命令一样。规范是后者适合模型; ttest不在此类别中。

双面 p值存储在r(p)中，如help ttest所示：

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)

来自ttest命令的bootstrap估计语法

1 个答案: