Как сгруппировать стандартные ошибки с помощью Stata stcrprep

Question

Я использую Stata для моделирования кумулятивных функций инцидентов при наличии конкурирующих рисков. Каждое заклинание заканчивается появлением события («сбой»), но предметы имеют возможность ввести новое заклинание впоследствии. Следовательно, заклинания не являются независимыми, а вложены в предметы.

Стандартная команда stata stcrreg может обрабатывать эту структуру путем моделирования стандартных ошибок, которые сгруппированы на уровне субъекта. Однако мой набор данных огромен (более 3 миллионов наблюдений), и время вычислений огромно. Другие пользователи предложили использовать написанную пользователем программу stcrprep, которая также обладает дополнительными функциями.

Кто-нибудь знает, как сгруппировать стандартные ошибки с помощью stcrprep?

Answer 1

Я нашел решение из информации, доступной в файле справки stcrprep и объясненной в презентации автором команды (http://repec.org/usug2013/lambert.uk13.pdf).

Сложные стандартные ошибки могут быть легко достигнуты следующим образом.

1. Подготовьте данные для анализа выживаемости с помощью stcrprep.
2. Создайте нужный формат, используя stset, указав весовые коэффициенты как весы. Указание их как iweights сделало бы невозможным оценку надежных и кластерных стандартных ошибок.
3. Проведите регрессию конкурирующих рисков с использованием stcox и укажите желаемый вариант ,vce(). Обратите внимание, что благодаря такому способу подготовки данных это даст точно такие же результаты, как stcrreg.

Для заинтересованного читателя, примера оценки риска рака при наличии конкурирующего сердечно-сосудистых заболеваний и другие причины смерти, с цензурой наблюдений и пациентов кластерным в терапии:

use "http://www.stata-journal.com/software/sj4-2/st0059/prostatecancer", clear 
stset time, failure(status=1,2,3) id(id) //specify all possible events 
stcrprep, events(status) keep(treatment age therapy) trans(1) //trans creates failcode (event of interest) 
gen event = failcode == status 
stset tstop [pw=weight_c], failure(event==1) enter(tstart) 
stcox treatment age, vce(cluster therapy)