Здравствуйте, гость ( Вход | Регистрация )
29.03.2012 - 16:00
Сообщение
#1
|
|
|
Группа: Пользователи Сообщений: 49 Регистрация: 7.04.2010 Пользователь №: 15366 |
Добрый день! Помогите разобраться. Имеются данные для анализа выживаемости:
время (месяцы), событие (рецидив или нет), независимые ковариаты 3 переменных - 1. дифференцировка - бинарная (низкодифференированныя и умеренно+высокодифференцированная) 2. уровень РЭА до операции 3. количество метастазов Они остались в результате применения регрессионного анализа Кокса в СПСС. Т.е. я понимаю что можно сделать выводы, что такие то факторы значимо влияют на продолжительность безрецидивной жизни, какой из них более значим. Могу ли я применить это как-то в плане прогнозирования продолжительности безрецидивной жизни? (как например в логистической регрессии). Это пожалуй основной вопрос. Цель исследования не только выявить исходные факторы, влияющие на продолжительность безрецидивного течения, но и иметь возможность прогнозировать это время для любого конкретного больного. Привожу формулу h(t)=h0(t) *exp(b1z1+...+bmzm) В СПСС есть функция сохранить: функция риска - это h(t)? Для 2 больных проживших 1 мес у одного 2,1 а другого 0,6. О чем это говорит? Максимум 4,3 минимум 0,2 - это много или мало. Как мне найти h0(t) - базовая функция риска при нулевых ковариатах (нулевые это равные 0?) и нужно ли это? С экспоненциальной регрессией более понятно, но у меня время жизни распределено не экспоненциально. В Attestat сделала подгонку распределения: получилось логлогистическое. Что мне с ним делать?! |
|
|
![]() |
![]() |
![]() |
2.04.2012 - 10:30
Сообщение
#2
|
|
|
Группа: Пользователи Сообщений: 49 Регистрация: 7.04.2010 Пользователь №: 15366 |
Добрый день! Спасибо большое за указанное направление. Попытки освоить R пока ни к чему не ведут, срубаюсь на этапе загрузки данных, но буду пытаться. Тут на форуме рекомендовали новую книжку. Может она поможет. Сразу приношу извинения за глупые вопросы.
Как я поняла из ветки: для нахождения h(o) строится график функции риска для средних значений ковариат в определенные точки времени и по графику находится h(0)? Но формула то должна быть. В результате построения регрессии Кокса у меня получилось 4 значимых переменных: 1. дифференцировка (0 - умеренная+высокая, 1-низкая) коэффициент - 0,62 среднее 0,61 2. уровень РЭА до операции коэффициент 0,006 среднее 149,78 3. кол-во образований малого диаметра 0,66 среднее 1,58 4. кол-во образований большого диаметра 0,70 среднее 0,67 нахожу p(i)=2,03 h(0)(t=1)=exp(-2,03)=0,13 это для единицы времени т.е. на 1 месяц, а дальше как например для 3, 6 месяцев и так далее? Еще вопрос S(t)=1-h(t), т.е. если выживаемость рассчитана для конкретного больного по нашим данным вычесть из единицы мы получим риск? Но я где-то читала что риск может быть больше 1, я неправильно поняла? Что мы прогнозируем с помощью модели Кокса предполагаемое время жизни определенного больного или вероятность умереть ( в нашем случае иметь МТС) в определенный момент времени, как я понимаю второй вариант. А нельзя ли как-то предсказывать время жизни? Например больной Н. имеет низкодифференцированную опухоль уровень РЭА 20, по 1 большому и маленькому мтс, подставляем в формулу данные и получаем что только через 12 месяцев мы имеем высокую вероятность получить прогрессию. В СПСС сохраняет параметры Функция риска, функция выживаемости. Они в сумме не должны равняться 1? В посте я столкнулась с тем что модели Кокса не применяются для прогнозирования по нравственным соображениям почему? |
|
|
![]() |
![]() |
mamalita регрессия Кокса 29.03.2012 - 16:00
p2004r [вдруг пригодится] вот целая ветка как в R проводи... 29.03.2012 - 17:05
mamalita Извиняюсь, данные загрузились, дальше двигалась по... 2.04.2012 - 11:15
TheThing Цитата(mamalita @ 2.04.2012 - 11:15)... 2.04.2012 - 12:51
mamalita Спасибо! Дело пошло дальше. Тупо повторяю кома... 2.04.2012 - 16:02
p2004r Цитата(mamalita @ 2.04.2012 - 16:02)... 2.04.2012 - 20:31
p2004r У Кокса
plot(survfit(model)) - shows the survivor... 2.04.2012 - 21:32
mamalita У меня новая проблема технического характера: все ... 3.04.2012 - 10:50
p2004r Цитата(mamalita @ 3.04.2012 - 10:50)... 3.04.2012 - 13:05
mamalita Надо просто найти где вычислялся v9 и вычислить. О... 4.04.2012 - 08:50
TheThing Цитата(mamalita @ 4.04.2012 - 08:50)... 4.04.2012 - 09:31
p2004r Цитата(mamalita @ 4.04.2012 - 08:50)... 4.04.2012 - 10:42
mamalita С английским дружо плохо, но посмотрю, просто руки... 12.04.2012 - 14:06
YVR Цитата(mamalita @ 12.04.2012 - 16:06... 12.04.2012 - 14:20
p2004r Цитата(mamalita @ 12.04.2012 - 14:06... 14.04.2012 - 23:09
mamalita Добрый день! Остроумненько. Прилагаю данные. Д... 13.04.2012 - 14:26
YVR Цитата(mamalita @ 13.04.2012 - 16:26... 13.04.2012 - 17:44
mamalita [С переменными: количество метастазов не совсем яс... 13.04.2012 - 18:29
YVR Цитата(mamalita @ 13.04.2012 - 20:29... 13.04.2012 - 19:47
p2004r Цитата(mamalita @ 13.04.2012 - 18:29... 14.04.2012 - 23:21
YVR Попробовал анализировать исходные данные. Выводы: ... 14.04.2012 - 19:37
p2004r Цитата(YVR @ 14.04.2012 - 19:37) Поп... 14.04.2012 - 23:25
p2004r Скопируйте сюда R текст своего анализа. Будет прод... 14.04.2012 - 23:23
mamalita RE: регрессия Кокса 17.04.2012 - 07:15
mamalita Уважаемый YVR! Рецидивы через три года были де... 17.04.2012 - 07:36
mamalita Уважаемый р2004r! у меня проблемы с сохранение... 17.04.2012 - 07:48
p2004r Цитата(mamalita @ 17.04.2012 - 07:48... 17.04.2012 - 11:29
mamalita А где посмотреть сами формулы? Могу ли я это сдела... 19.04.2012 - 10:42
p2004r Цитата(mamalita @ 19.04.2012 - 10:42... 27.04.2012 - 20:00
YVR Цитата(mamalita @ 17.04.2012 - 09:36... 19.04.2012 - 12:47
mamalita А что это за формула и вообще какой метод? И почем... 23.04.2012 - 13:19
YVR Цитата(mamalita @ 23.04.2012 - 15:19... 23.04.2012 - 16:54
mamalita И все таки как же мне найти h0? 23.04.2012 - 13:20
mamalita И все-таки какой метод и что за формула? Вопрос пр... 24.04.2012 - 13:21
YVR Цитата(mamalita @ 24.04.2012 - 15:21... 24.04.2012 - 14:24
p2004r Цитата(YVR @ 24.04.2012 - 14:24) Мат... 27.04.2012 - 20:15
YVR Цитата(p2004r @ 27.04.2012 - 22:15) ... 19.05.2012 - 20:48
DrgLena Цитата(YVR @ 19.05.2012 - 20:48) О ч... 20.05.2012 - 12:14
YVR Цитата(DrgLena @ 20.05.2012 - 14:14)... 20.05.2012 - 13:03
DrgLena Да вы хоть поняли, что у всех больных уже были опе... 20.05.2012 - 13:36
YVR Цитата(DrgLena @ 20.05.2012 - 15:36)... 20.05.2012 - 13:48
DrgLena Цитата(YVR @ 20.05.2012 - 13:48) Вам... 20.05.2012 - 17:43
YVR Цитата(DrgLena @ 20.05.2012 - 19:43)... 20.05.2012 - 18:05
p2004r Цитата(YVR @ 20.05.2012 - 18:05) Эпи... 20.05.2012 - 18:21
YVR Цитата(p2004r @ 20.05.2012 - 20:21) ... 20.05.2012 - 20:16
p2004r Цитата(YVR @ 20.05.2012 - 20:16) Ну ... 20.05.2012 - 21:42![]() ![]() |