Регрессия Кокса |
Здравствуйте, гость ( Вход | Регистрация )
Регрессия Кокса |
9.01.2012 - 23:18
Сообщение
#1
|
|
Группа: Пользователи Сообщений: 59 Регистрация: 23.12.2011 Пользователь №: 23383 |
Немного запутался. Подскажите, пожалуйста.
Есть две группы больных. Основная и сравнения. Для обеих групп определяю факторы риска смерти по методу регрессии Кокса. Получаю коэффициент регресии и его ст.ошибку, хи-квадрат, р-значение, отношение шансов и 95%. Так вот, в основной группе фактором риска является высеваемая флора (метод высоко эффективен при грамотрицательной флоре) - это закономерно и ожидаемо. Больных с грамотрицательной флорой умерло на самом деле меньше. В группе сравнения летальность при грам-, грам+, грам+/- примерно одинаковая и флора фактором риска не является. Могу ли я сказать, что у больных основной группы достоверно снизилась вероятность умереть при грам- сепсисе? Как выразить это численно? 2. Еще вопрос по регрессии. Можно ли одновременно задавать в качестве предикторных переменных переменные, которые тесно коррелируют между собой. Как я понял из литературы - это крайне нежелательно. А что тогда делать? Просто у меня в анализе выжимаемости регрессией Кокса факторами риска являются как локализация первичного очага, так и высеваемая микрофлора. Разумеется, часто (хоть и не всегда) определенной локализации процесса соответвует определенная конкретная микрофлора. Что более важно и удобно с клинической точки зрения - флора или локализация - сказать сложно. У меня четыре градации признака "локализация" и три градации признака "флора". Что делать? Использовать простой регрессионный анализ отдельно для каждого фактора? Сообщение отредактировал Alex_Z - 10.01.2012 - 00:48 |
|
10.01.2012 - 11:27
Сообщение
#2
|
|
Группа: Пользователи Сообщений: 1325 Регистрация: 27.11.2007 Пользователь №: 4573 |
В какой программе проводите кокс- регрессию?
|
|
10.01.2012 - 11:30
Сообщение
#3
|
|
Группа: Пользователи Сообщений: 59 Регистрация: 23.12.2011 Пользователь №: 23383 |
SPSS v.17
|
|
10.01.2012 - 11:51
Сообщение
#4
|
|
Группа: Пользователи Сообщений: 1325 Регистрация: 27.11.2007 Пользователь №: 4573 |
SPSS v.17 Уже хорошо. А в каком виде анализировалась переменная локализация, ведь фактором риска является определенная локализация относительно другой. Мне пришлось рецензировать работу в которой локализация кодировалась 1,2,3,4,5 и загонялась в Statistica как количественная переменная. С флорой - тоже. Возможно, две группы это условность, а при анализе выживаемости метод лечения, или тот признак , который разбил на две группы может быть предиктором, как и другие факторы. |
|
10.01.2012 - 12:25
Сообщение
#5
|
|
Группа: Пользователи Сообщений: 59 Регистрация: 23.12.2011 Пользователь №: 23383 |
Как я понял, в регрессии Кокса определяется как вырастет риск при увеличении предикторной переменной на единицу значения. Так? Или как увеличится риск в случае дихотомической переменной - да/нет.
А как быть, если необходимо проанализировать влияние номинальных переменных (я так понимаю, закодировав локализацию я получил такую шкалу) с тремя - четырьмя возможными значениями? Сообщение отредактировал Alex_Z - 10.01.2012 - 13:31 |
|
10.01.2012 - 14:04
Сообщение
#6
|
|
Группа: Пользователи Сообщений: 59 Регистрация: 23.12.2011 Пользователь №: 23383 |
О! SPSS можно выделить часть ковариант как категориальные... Это решение?
http://translate.google.ru/translate?hl=ru....htm&anno=2 Вот тут, если я правильно понял, сказано, что SPSS может ввести "фиктивные" переменные, если обозначить их как категориальные. Т.е. я смогу использовать данные "как есть"? "Dummy variables. Whereas SPSS will create dummy variables automatically if a variable is declared categorical suing the Categorical button (see the figure above), this must be done explicitly in Stata" Сообщение отредактировал Alex_Z - 10.01.2012 - 15:00 |
|
10.01.2012 - 15:12
Сообщение
#7
|
|
Группа: Пользователи Сообщений: 59 Регистрация: 23.12.2011 Пользователь №: 23383 |
|
|
10.01.2012 - 20:04
Сообщение
#8
|
|
Группа: Пользователи Сообщений: 1114 Регистрация: 10.04.2007 Пользователь №: 4040 |
Немного запутался. Подскажите, пожалуйста. Есть две группы больных. Основная и сравнения. Для обеих групп определяю факторы риска смерти по методу регрессии Кокса. Очевидно, имеется в виду модель пропорциональных рисков Кокса. Возьмите модуль анализа выживаемости и анализа данных типа времени жизни Survival Analysis программного обеспечения AtteStat. Там сделан подробный обзор и дается возможность расчета. Правда, работает только под Windows (версия для современных операционных систем появится еще не скоро). Ebsignasnan prei wissant Deiws ainat! As gijwans! Sta ast stas arwis!
|
|
10.01.2012 - 20:44
Сообщение
#9
|
|
Группа: Пользователи Сообщений: 1325 Регистрация: 27.11.2007 Пользователь №: 4573 |
Прежде чем пользоваться любой программой, нужно понимать принцип метода и особенности работы с переменными разного типа.
Для того, чтобы понять как логистическая или кокс регрессия работает с категориальными переменными, поймите сначала что такое относительный риск, грибы есть (1), грибов нет (0). Как вы трактуете коэффициент? Что значит 1,521? А теперь посмотрите, как программа предложила вам закодировать имеющиеся локализации, в отчете этого модуля есть такая табличка, а не только та, что вы привели. Коэффициенты будут зависеть от выбора, который вы сделаете при работе с категориальными переменными. А что AtteStat предлагает для номинальных переменных, нужно предварительно создать бинарные для каждой категории локализации? Или все же предикторы могут быть только количественными? |
|
10.01.2012 - 21:05
Сообщение
#10
|
|
Группа: Пользователи Сообщений: 59 Регистрация: 23.12.2011 Пользователь №: 23383 |
DrgLena
Я так понимаю, что при обнаружении грибов, шанс (или риск) того, что наступил событие 1 (смерть) увеличится в 1,521 раз. Или на 152,1% Если я правильно понял, что Exp(B) - это отношение шансов. Верно? Вы про эту таблицу? Игорь, спасибо, постараюсь уже разбираться с SPSS. А то вообще мозги закипят. Сообщение отредактировал Alex_Z - 10.01.2012 - 21:06 |
|
10.01.2012 - 21:35
Сообщение
#11
|
|
Группа: Пользователи Сообщений: 1325 Регистрация: 27.11.2007 Пользователь №: 4573 |
Вот, что получил... А почему SPSS не показывает коэффициент для "Причина" и "Флорп_бак". Теперь вам понятно, почему для последних категорий этих двух переменных коэффициент не нужен? Относительно этих категорий вы риски и считаете, а могли бы и другой выбор сделать, зависит от ваших предварительных знаний о роли каждой категории. А чем все же две группы отличаются, почему вы создаете отдельные модели для каждой? |
|
10.01.2012 - 21:41
Сообщение
#12
|
|
Группа: Пользователи Сообщений: 1325 Регистрация: 27.11.2007 Пользователь №: 4573 |
|
|
10.01.2012 - 22:00
Сообщение
#13
|
|
Группа: Пользователи Сообщений: 59 Регистрация: 23.12.2011 Пользователь №: 23383 |
"Нет, не правильно "в" это не "на". ПО вашему, если риск увеличивается в 1 раз, то это на 100%?"
Да, понял ошибку. А если коэффициент меньше, одного, то снижается в 0,ХХХ раз? "А чем все же две группы отличаются, почему вы создаете отдельные модели для каждой?" Группа один - основная, там применен исследуемый метод лечения (эффективен при грамотрицательном сепсисе). Группа два - сравнения. То есть, если я вас правильно понял, больных обеих групп нужно объеденить и оценить факторы риска для больных обеих групп одновременно? И НЕпроведение процедуры будет фактором риска? Просто я хотел показать, что 1) процедура эффективна - оцениваем клиническую динамику, которая достоверно лучше, чем в группе сравнения; 2) построить кривые Каплан-Меера - показать, что выживаемость увеличилась в основной группе (сноска*); 3) показать, что при применении этого метода лечения у больных вид высеваемых бактерий - фактор риска, для группы сравнения - нет, т.к. у них шансы умереть от любой флоры примерно равны. Мой алгоритм неверен? Т.е. для номинативных факторов, в которых больше двух градаций, выбрать "нулевой" уровень, относительно которого будет идти сравнение? МДа...тогда надо думать как поступить... * я так понимать, это можно сделать тестами лог-ранк и Бреслоу? А не скажете чем они отличаются? Я где-то читал, что один из них более отражает достоверность на начальных этапах, а другой на поздних. Это верно? А где примерная граница применения этих тестов? Не могу никак снова найти этот источник. |
|
10.01.2012 - 23:17
Сообщение
#14
|
|
Группа: Пользователи Сообщений: 1325 Регистрация: 27.11.2007 Пользователь №: 4573 |
Анализ зависит от того что планировалось, как формировались группы. Если они отличались только методом лечения, то его эффективностью и будет знак и оценка коэффициента кокс регрессии, который будет учитывать и флору и локализацию и другие факторы.
К-М кривые могут демонстрировать различия. Анализируя кривые вы можете выбрать срок наблюдения, когда различия будут особенно выражены и доказательства направить на эти сроки. Относительно тестов, Гехана-Вилкоксона, рекомендуют в тех случаях, когда различия кривых более выражены в начальный период наблюдения, а также в случаях, когда не выполняется модель пропорциональных рисков. Когда риски пропорциональны рекомендуется Кокс и Лог-ранг тесты. Критерий Кокса более чувствителен к обнаружению различий на концах распределения. Лог-ранговый рекомендуется, когда число смертей мало. Но эти тесты в Statistica, а что в SPSS ?посмотрите сами, вам до оценки моделей еще далеко, разберитесь с кодировкой предикторов, можно снизить число градаций путем их объединения до ?плохих? и ?хороших? относительно выживания, флора тоже может быть подвержена редуцированию, кроме того могут быть созданы переменные объединяющие и локализацию и флору , например зловредные с плохой локализацией и альтернатива, как вам подсказывают ваши наблюдения. |
|
10.01.2012 - 23:55
Сообщение
#15
|
|
Группа: Пользователи Сообщений: 59 Регистрация: 23.12.2011 Пользователь №: 23383 |
Спасибо большое за разъяснения и советы. Буду много думать...
Ну дело в том, что на самом деле, грамотрицательная флора чаще всего выявляется при локализации инфекции в трансплантате (мочевая система). И эффективность при этом наиболее высока, это факт. Но все равно хочется отдельно анализировать эти два фактора. А если сделать новые факторы? Скажем, по флоре - только грам"-" (да/нет), "смешанная" (да/нет), грам"+" (да/нет). И таким же образом поступить с локализацией? Можно так поступить? А, наверное, я не так выразился. При оценке связи между номинальными признаками используется коэффициент не корреляции, а сопряженности? Достаточно ли будет при описании связи привести сам коэффициент и значимость? Сообщение отредактировал Alex_Z - 12.01.2012 - 10:20 |
|