Я пытаюсь продемонстрировать характеристики различных тестов для небольших выборок данных. Я хотел бы продемонстрировать производительность t-теста, t-теста с начальной оценкой и теста ранговой суммы. Меня интересует получение p-значения для каждого теста на нескольких наборах данных с использованием simulate
. Однако я не могу получить оценки t-теста, используя префикс bootstrap
и команду ttest
.
Данные генерируются:
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)
Ошибки для каждой из процедур по порядку:
#P4#
#P5#
#P6#
Эти результаты не согласуются с документацией для префикса 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-значение.
return list
послеbootstrap
не содержитr(p)
,\ покажет вам, чтоr(p)
не возвращаетсяbootstrap
, поэтому ваша последняя команда выдаст ошибку. - person Steve Samuels   schedule 24.06.2015