Здравствуйте, гость ( Вход | Регистрация )
10.03.2008 - 07:36
Сообщение
#1
|
|
|
Группа: Пользователи Сообщений: 1219 Регистрация: 13.01.2008 Из: Челябинск Пользователь №: 4704 |
Идеальные данные для анализа выживаемости - когда точно известно сколько человек прожил, например, после операции и когда умер. В этом случае цензурированных наблюдений нет. Другой крайний случай - когда все наблюдения цензурированные и дальнейшая судьба пациентов неизвестна. Например один прожил больше года, другой - больше трех. В этом случае может оказаться, что больше года - это 5, а больше трех - это 4. Поэтому, насколько я понимаю, сравнить выживаемость в двух группах где все наблюдения цензурированные невозможно в принципе. А какова допустимая доля цензурированных наблюдений в выборке? Существуют ли какие-то обоснованные или негласные правила? Полазил в и-нете, заглянул в книжки - пока ответа не нашел, хотя везде рассматриваются примеры где полные данные заметно преобладают над цензурированными. Или может считать цензурированными только точно живых на момент анализа, а всех потерявшихся считать умершими в интервале между двумя осмотрами, как прочитал в одной статье?
|
|
|
![]() |
![]() |
![]() |
1.12.2008 - 02:09
Сообщение
#2
|
|
|
Группа: Пользователи Сообщений: 1325 Регистрация: 27.11.2007 Пользователь №: 4573 |
"Все не так уж просто"
Безусловно, не просто!!! Как и любой прогноз в медицине. Ведь мы его проводим не для того, чтобы он непременно был реализован по нашим моделям, а для предотвращения нежелательного прогноза. Т.е. найти значимые факторы и воздействовать на них. Да, возможно это можно назвать упражнениями, но когда в базе данных более 3,5 тыс наблюдений (101 больной наблюдается более 20 лет), можно создать обучающую выборку из 917 больных и помоделировать. Конечно, помоделировать с целью оценить различные подходы к лечению, с учетом известных или новых факторов риска. Начать действительно нужно с предположения о форме функции выживаемости. В модуле таблицы выживаемости я это могу сделать, причем оценки будут зависеть от выбранного интервала, выбираю 6 мес. Предлагается сравнение с 4 видами распределения Exponential, Linear Hazard, Gompertz, Weibull. Параметры оценки для этих видов очень близки, соответственно р=0,87, р=0,86, р=0,90, р=0,97. Т.е. подходит любое, но мне ближе Weibull, выбираю и получаю таблицы выживания для интервалов по 6 мес. Но в последующем, при создании модели никакие программы меня не спрашивают, что мне ближе. Т.е. вопросов больше, чем ответов. Насчет плохих привычек, я вообще не поняла?., что с чем совмещали раньше. Грустно только, что прогнозы сбываются и в контрольной выборке. |
|
|
![]() |
![]() |
1.12.2008 - 21:57
Сообщение
#3
|
|
|
Группа: Пользователи Сообщений: 1013 Регистрация: 4.10.2006 Пользователь №: 1933 |
"Насчет плохих привычек, я вообще не поняла?., что с чем совмещали раньше. Грустно только, что прогнозы сбываются и в контрольной выборке. Это я имел ввиду, что модель пропорционального риска не предназначена для оценки функции выживаемости. Она получает коэффициенты, повышающие или понижающие риск, а затем используется эмпирическая фнукция выживаемости Каплана-Мейера, которая и меняется в соответствии с коэффициентами модели пропорционального риска. Оценка функции выживаемости - это задача AFTM (той самой модели Вейбулла), однако она должна иметь определенный теоретический вид. Соответственно, писать, что "функцию выживаемости оценивали по модели Кокса" - неправильно, ибо она ее не оценивает. Макисмум, что можно написать - "функцию выживаемости анализировали при помощи коррекции эмпирической функции выживаемости (КМ) коэффициентами модели пропорционального риска". Что же касается численности, то сколько бы ни было вначале, важно на основании какого количества пациентов проводится расчет эмпирической функции в правом углу графика. Если, скажем, было 2000 человек, но 1990 умерли в первый год, а оставшиеся 10 медленно вымирали следующие 5 лет, то считать эмпирическую функцию выживаемости адекватной больше, чем через год после начала исследования нельзя. С AFTM все проще. Поскольку для оценки параметров используется весь массив данных, то фнукция выживаемости справедлива для всего диапазона данных (надо только подбирать распределения, для примера выше это явно не экспонента и не вейбулл, скорее распределение Гомперца - оно позволяет моделировать эффект выздоровления). Вообщем все, что я хотел сказать, это то, что модель Кокса должна использоваться для того, для чего она была разработана - анализ и сравнение факторов риска, а не для того, для чего она не разрабатывалась (изучение функции выживаемости). |
|
|
![]() |
![]() |
nokh Допустимая доля цензурированных наблюдений в анализе выживаемости 10.03.2008 - 07:36
DrgLena Естественно, что анализ выживаемости не проводят,... 12.03.2008 - 23:55
nokh Спасибо за внимание к моей проблеме. С анализом вы... 13.03.2008 - 00:51
плав Цитата(nokh @ 13.03.2008 - 00:51) Сп... 13.03.2008 - 12:58

ВалНест Речь идет о цензурированных сулчаях, т.е. известно... 1.06.2008 - 09:16
ВалНест Цитата(nokh @ 13.03.2008 - 01:51) Сп... 1.06.2008 - 09:10
плав Цитата(ВалНест @ 1.06.2008 - 10:10) ... 3.06.2008 - 22:57
DrgLena Вы затронули очень важную тему в анализе данных. Ч... 13.03.2008 - 12:53
DrgLena Обведенная колонка у меня съехала, я хотела отмети... 13.03.2008 - 12:57
плав To DrgLena
Если известно, что человек был жив чере... 13.03.2008 - 13:24
ВалНест Уважаемый ПЛАВ!
У меня проблемы при сравнении ... 1.06.2008 - 08:48
DrgLena Как именно игнориорвать мне и не ясно. Условие вкл... 13.03.2008 - 14:12
DrgLena Вот пример реальный, из той же базы данных, сравни... 13.03.2008 - 14:37
плав Цитата(DrgLena @ 13.03.2008 - 14:37)... 13.03.2008 - 18:38
DrgLena График убежал, повторяю jpg 13.03.2008 - 14:40
nokh Спасибо! На свои вопросы я ответы получил. Вых... 13.03.2008 - 22:55
Igoroshka Позвольте и свои 5 копеек добавить .
Усеченные (ц... 4.08.2008 - 17:58
плав Ну и каким образом это "улучшает" выжива... 5.08.2008 - 13:24
Igoroshka Цитата(плав @ 5.08.2008 - 13:24) Ну ... 5.08.2008 - 14:11
плав Цитата(Igoroshka @ 5.08.2008 - 15:11... 9.08.2008 - 20:14
Igoroshka Цитата(плав @ 9.08.2008 - 20:14) Глу... 12.08.2008 - 12:20
плав Цитата(Igoroshka @ 12.08.2008 - 13:2... 20.08.2008 - 11:10
Игорь Цитата(плав @ 9.08.2008 - 20:14) Глу... 10.08.2008 - 18:58
Игорь По информационным каналам прошло сообщение. Как ра... 1.10.2008 - 15:26
Igoroshka Цитата(Игорь @ 1.10.2008 - 15:26) По... 4.10.2008 - 09:42
nokh Еще информация к размышлению по части анализа выжи... 9.10.2008 - 21:44
Игорь Цитата(nokh @ 9.10.2008 - 22:44) Еще... 10.10.2008 - 13:19
nokh Не стал открывать новую ветку, т.к. вопрос опять п... 14.11.2008 - 08:59
плав Цитата(nokh @ 14.11.2008 - 08:59) Не... 14.11.2008 - 16:43
nokh Большое спасибо, пока все получилось красиво. Но я... 17.11.2008 - 18:28
плав Цитата(nokh @ 17.11.2008 - 18:28) Бо... 21.11.2008 - 23:13
nokh Спасибо. Получается, что если переменная бинарная,... 22.11.2008 - 18:05
DrgLena Да, если предиктор бинарный, то экспоненциальный к... 23.11.2008 - 00:56
nokh Спасибо за разъяснения. Выходит что если предиктор... 27.11.2008 - 19:21
плав Цитата(nokh @ 27.11.2008 - 19:21) 1)... 27.11.2008 - 22:17
DrgLena «Как вы это проделываете?»
Да также ка... 28.11.2008 - 03:30
плав Цитата(DrgLena @ 28.11.2008 - 03:30)... 28.11.2008 - 10:29
DrgLena Мы очевидно для разных целей используем кокс регре... 28.11.2008 - 14:27
плав Цитата(DrgLena @ 28.11.2008 - 14:27)... 28.11.2008 - 17:54
DrgLena Да, мне наконец удалось, посчитать руками, помогла... 29.11.2008 - 04:11
плав Цитата(DrgLena @ 29.11.2008 - 04:11)... 1.12.2008 - 00:53
nokh Спасибо за ответы и интересную дискуссию. Использо... 30.11.2008 - 23:23
плав Цитата(nokh @ 30.11.2008 - 23:23) Сп... 1.12.2008 - 00:56
Игорь Цитата(DrgLena @ 1.12.2008 - 02:09) ... 1.12.2008 - 18:37
DrgLena Я не стала бы возражать уважаемому модератору, но ... 2.12.2008 - 02:37
плав Цитата(DrgLena @ 2.12.2008 - 02:37) ... 2.12.2008 - 22:18
DrgLena Ответ уже нашла
http://www.weibull.com/ 2.12.2008 - 21:37
DrgLena Плав, спасибо большое за внимание и терпение, с к... 3.12.2008 - 16:05
плав Цитата(DrgLena @ 3.12.2008 - 16:05) ... 3.12.2008 - 22:21
DrgLena Да, не используется не только кокс регрессионное м... 3.12.2008 - 23:58![]() ![]() |