Синтаксис оценок начальной загрузки из команды ttest

Я пытаюсь продемонстрировать характеристики различных тестов для небольших выборок данных. Я хотел бы продемонстрировать производительность 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)

Ошибки для каждой из процедур по порядку:

  1. #P4#
  2. #P5#
  3. #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-значение.


person Todd D    schedule 18.06.2015    source источник
comment
return list после bootstrap не содержит r(p),\ покажет вам, что r(p) не возвращается bootstrap, поэтому ваша последняя команда выдаст ошибку.   -  person Steve Samuels    schedule 24.06.2015


Ответы (1)


ttest — это команда r-класса, результаты которой сохраняются в 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)
person Roberto Ferrer    schedule 18.06.2015
comment
Это решение не работает, как указано выше в отредактированном запросе. - person Todd D; 19.06.2015
comment
bootstrap имеет собственный файл справки, и вы можете легко получить к нему доступ с помощью help bootstrap. Там вы можете обнаружить, что результат начальной загрузки находится в e(b). Запустите matrix list e(b), чтобы отобразить его. Что касается вашего желания создать переменную, если вы хотите получить промежуточные результаты, отметьте опцию saving() для bootstrap. Если вы хотите сохранить окончательный результат, теперь вы знаете, как получить к нему доступ. См. help postfile о сохранении результатов в наборе данных Stata. Я рекомендую лучше использовать руководства help files и Stata. (Вы можете сделать любой ответ неуместным, достаточно изменив свой вопрос.) - person Roberto Ferrer; 20.06.2015