#stata
#stata
Вопрос:
Я пытаюсь в Stata создать симуляцию, использующую загрузку моего населения bsample .
Набор stan3 данных, и моя цель-получить процентили времени выживания (10 20 30 40 50). Я знаю, что могут быть и другие способы сделать то же самое, но я хотел бы понять, почему, даже если мой код, кажется, работает, после прохождения 20-го процентиля я получаю все те же результаты.
use https://www.stata-press.com/data/r17/stan3.dta, clear set trace on tempname p tempvar s i postfile `p' int(rep) float(s i) using btsur, replace forval i=10(10)50 { // i are percentiles forval nsim = 1/100 { quietly bsample quietly sts gen `s'= s // estimates the chances of survival in the sample obtained by bootstraping quietly summarize _t if `s'lt;=(100-`i')/100, meanonly post `p' (`nsim') (`r(min)') (`i') // i save r(min) because I want the lowest time possible where survival is 100-`i' drop `s' } } postclose `p' use btsur,clear browse in 80/250 Комментарии:
1. Вы каждый раз проходите цикл после первой загрузки образца из предыдущего загруженного образца, а не из исходных данных. Это озадачивающий выбор, и я подозреваю, что это ошибка в вашем коде. Кроме того, я не могу легко следовать вашей процедуре расчета процентилей. Я бы рекомендовал прямое использование
pctile. Кроме того, естественная процедура, по-видимому, заключалась бы в том, чтобы взять образец начальной загрузки, а затем вычислить из него нужные процентили, а не брать другую выборку начальной загрузки для каждого процентиля.2. Привет, Ник, и спасибо за ваш ответ, идея заключалась в том, чтобы оценить доверительный интервал вокруг каждого процентиля (10 20 30 40 50) начальной загрузки времени выживания. Вот в чем причина этой симуляции. Выше этого ваша точка зрения последовательна. Я изложу пошаговую идею, лежащую в основе кода:-мы берем функцию выживания-мы указываем на 10 20 30 40 50 процентилей времени выживания — мы загружаемся с помощью bsample для создания CI со 100 повторениями — мы выбираем r(min), потому что хотим как можно меньше времени, когда выживание равно 100 — `я»
3. Вы можете запросить результаты сразу для нескольких статистических данных. Глядя на справку
bootstrap, я вижу, что первый пример показывает именно это. Иначе говоря, почему бы не использоватьbootstrapздесь?4. Я лично знаком с Тимом Моррисом, но я не читал эту чушь. Я сомневаюсь, что Тим подпишется на то, что вы, похоже, делаете, особенно на загрузку загруженного образца. Извините, но я думаю, что мне нужно отключиться здесь, так как ваш вопрос продолжает меняться и поднимает слишком много общих и конкретных вопросов.
5. Я не раздражен, и «извините» означает извинение. Теперь я сомневаюсь, что могу сказать что-то еще, что поможет вам, не в последнюю очередь потому, что я не работаю с данными о выживании и не являюсь медицинским статистиком. Это ограничение каждого сайта, форума или списка рассылки, что они лучше всего работают с простыми, конкретными вопросами. Как только вопрос, кажется, превращается в сложную дискуссию, вам нужен эксперт на другом конце и/или какой-то другой вид обучения, обучения или обратной связи.