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

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

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





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



Идеальные данные для анализа выживаемости - когда точно известно сколько человек прожил, например, после операции и когда умер. В этом случае цензурированных наблюдений нет. Другой крайний случай - когда все наблюдения цензурированные и дальнейшая судьба пациентов неизвестна. Например один прожил больше года, другой - больше трех. В этом случае может оказаться, что больше года - это 5, а больше трех - это 4. Поэтому, насколько я понимаю, сравнить выживаемость в двух группах где все наблюдения цензурированные невозможно в принципе. А какова допустимая доля цензурированных наблюдений в выборке? Существуют ли какие-то обоснованные или негласные правила? Полазил в и-нете, заглянул в книжки - пока ответа не нашел, хотя везде рассматриваются примеры где полные данные заметно преобладают над цензурированными. Или может считать цензурированными только точно живых на момент анализа, а всех потерявшихся считать умершими в интервале между двумя осмотрами, как прочитал в одной статье?
Вернуться в начало страницы
 
+Ответить с цитированием данного сообщения
 
 
Открыть тему
Ответов
DrgLena
сообщение 3.12.2008 - 16:05
Сообщение #2





Группа: Пользователи
Сообщений: 1325
Регистрация: 27.11.2007
Пользователь №: 4573



Плав, спасибо большое за внимание и терпение, с которым вам удается находить необходимые слова и деликатный тон общения с людьми разного уровня подготовки. Вы оказываете реальную помощь людям, которые самостоятельно пытаются разобраться в большом количестве информации в столь интересной области знаний.
Именно с формулой я и пытаюсь разобраться, вы пишете?
В модели пропорционального риска риск (hazard) каждого человека определяется как часть общей опасности:
h_i(t)/h_j(t)=exp(b1*x_i-x_j)+...).

Я же вижу в формулах из разных источников, которые я приводила, что риск для каждого человека, не часть общей опасности, а отношение его риска и базовому, т.е. к среднему в выборке по которой модель построена. И если умножить обе части на h_j(t), то это и будет его риск.
Эту ветку по кокс регрессии не я начала и интерес к ней скромный, мало кто использует этот метод анализа, но я благодаря вашим постам нашла в Statistica весьма полезный модуль, Process Analysis, с прекрасными графическими возможностями, получила и параметр формы и параметр масштаба, а также нашла лучший параметр положения (location), поскольку распределение Вейбулла ограничено слева.
Спасибо!!!
Вернуться в начало страницы
 
+Ответить с цитированием данного сообщения
 
плав
сообщение 3.12.2008 - 22:21
Сообщение #3





Группа: Пользователи
Сообщений: 1013
Регистрация: 4.10.2006
Пользователь №: 1933



Цитата(DrgLena @ 3.12.2008 - 16:05) *
Плав, спасибо большое за внимание и терпение, с которым вам удается находить необходимые слова и деликатный тон общения с людьми разного уровня подготовки. Вы оказываете реальную помощь людям, которые самостоятельно пытаются разобраться в большом количестве информации в столь интересной области знаний.
Именно с формулой я и пытаюсь разобраться, вы пишете?
В модели пропорционального риска риск (hazard) каждого человека определяется как часть общей опасности:
h_i(t)/h_j(t)=exp(b1*x_i-x_j)+...).

Я же вижу в формулах из разных источников, которые я приводила, что риск для каждого человека, не часть общей опасности, а отношение его риска и базовому, т.е. к среднему в выборке по которой модель построена. И если умножить обе части на h_j(t), то это и будет его риск.
Эту ветку по кокс регрессии не я начала и интерес к ней скромный, мало кто использует этот метод анализа, но я благодаря вашим постам нашла в Statistica весьма полезный модуль, Process Analysis, с прекрасными графическими возможностями, получила и параметр формы и параметр масштаба, а также нашла лучший параметр положения (location), поскольку распределение Вейбулла ограничено слева.
Спасибо!!!

Да, формула такая, как я привел. То, что обычно приводится - (h_i(t)=h(t)*...) это просто модель пропорционального риска, а не модель Кокса (как, например, модель Вейбулла). Но тогда, поскольку h(t) в правой части, надо ее описать. А для эмпирической кривой, каждое значение - отдельная точка. Соответственно для оценки h(t) будем использовать уже все точки и степеней свободы (точек, наблюдений) для оценки влияния факторов риска не будет. Альтернативой является описать функицю риска каким-то параметрическим уравнением (например, распределением Вейбулла) - у него всего два (три) параметра. Тогда останется и на оценку факторов риска. А и то и другое вместе - никак. Кокс же предложил анализ, который позволял вообще обойтись без оценки h(t), путем того, что вместо опасности анализировалось отношение опасностей, а уж h_i(t) можно было найти проанализировав влияние факторов риска и построив эмпирическую кривую (со всеми вытекающими проблемами в правой части).
А вообще-то жалко, что методом (у нас) пользуются мало - он мощный, красивый и более информативный, чем, например, логистическая регрессия (но надо иметь время до события, а при принятом у нас "ретроспективном" сборе материла...)

Вернуться в начало страницы
 
+Ответить с цитированием данного сообщения
 

Сообщений в этой теме
- 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


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