來自ttest命令的bootstrap估計的語法

Question

我試圖演示小樣本數據的各種測試的特性。我想證明t檢驗的表現，t檢驗以及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) 

每個的在順序程序的錯誤是：需要

無效

1. 表達式列表表達：_b

2. 'E（P）' 評價爲全樣本中缺失

這些結果不與用於bootstrap前綴的文檔是一致的。 e或r類對象的規格是否有問題並且ttest？

編輯：

瞭解，現在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) 
    display r(p) 

不返回p值。

Answer 1

ttest是r級命令並將其重試存儲在r()中。您似乎期望它能將結果保存在e()中，如電子類命令。規範是後者適合模型; 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)