Форум врачей-аспирантов

Здравствуйте, гость ( Вход | Регистрация )

> Допустимая доля цензурированных наблюдений в анализе выживаемости, а что если 100%?
nokh
сообщение 10.03.2008 - 07:36
Сообщение #1





Группа: Пользователи
Сообщений: 1219
Регистрация: 13.01.2008
Из: Челябинск
Пользователь №: 4704



Идеальные данные для анализа выживаемости - когда точно известно сколько человек прожил, например, после операции и когда умер. В этом случае цензурированных наблюдений нет. Другой крайний случай - когда все наблюдения цензурированные и дальнейшая судьба пациентов неизвестна. Например один прожил больше года, другой - больше трех. В этом случае может оказаться, что больше года - это 5, а больше трех - это 4. Поэтому, насколько я понимаю, сравнить выживаемость в двух группах где все наблюдения цензурированные невозможно в принципе. А какова допустимая доля цензурированных наблюдений в выборке? Существуют ли какие-то обоснованные или негласные правила? Полазил в и-нете, заглянул в книжки - пока ответа не нашел, хотя везде рассматриваются примеры где полные данные заметно преобладают над цензурированными. Или может считать цензурированными только точно живых на момент анализа, а всех потерявшихся считать умершими в интервале между двумя осмотрами, как прочитал в одной статье?
Вернуться в начало страницы
 
+Ответить с цитированием данного сообщения
 
 
Открыть тему
Ответов
DrgLena
сообщение 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



Цитата(DrgLena @ 1.12.2008 - 02:09) *
"Насчет плохих привычек, я вообще не поняла?., что с чем совмещали раньше. Грустно только, что прогнозы сбываются и в контрольной выборке.

Это я имел ввиду, что модель пропорционального риска не предназначена для оценки функции выживаемости. Она получает коэффициенты, повышающие или понижающие риск, а затем используется эмпирическая фнукция выживаемости Каплана-Мейера, которая и меняется в соответствии с коэффициентами модели пропорционального риска. Оценка функции выживаемости - это задача 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
|- - плав   Цитата(DrgLena @ 1.12.2008 - 02:09) ...   1.12.2008 - 21:57
- - Игорь   Цитата(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


Добавить ответ в эту темуОткрыть тему