Версия для печати темы

Нажмите сюда для просмотра этой темы в обычном формате

Форум врачей-аспирантов _ Медицинская статистика _ Анализ взаимодействий в регрессионном анализе Кокса

Автор: mix3d 31.10.2014 - 14:42

Здравствуйте,

помогите, пожалуйста, разобраться или наведите на источник, понятный для клинициста.

Определяем независимые предикторы неблагоприятного исхода (регрессионный анализ Кокса). Сначала отобрали возможные кандидаты в однофакторной регрессии с р<0,2. Далее ввели их в многофакторную модель.

Вопрос: требуется учесть взаимодействия, в SPSS такая возможность реализована, но не уверен как ее применить.

Например, если в конечной модели у меня 3 независимых предиктора + пол и возраст - незначимые, но использованы для adjustment, следует ли анализ взаимодействия проводить просто перемножив 3 этих независимых предиктора все вместе (1*2*3), или попарно (1*2, 1*3, 2*3), или еще и возраст, и пол (1*2*3*пол*возраст, 1*2*пол, 1*2*возраст, 1*2*пол*возраст и т.д.).

Стоит ли учитывать все возможные комбинации, или только исходя из физиологической/патофизиологической целесообразности?

Далее, включать ли каждое произведение (двух предикторов) по отдельности в исходную многофакторную модель, или только произведение (т.е. по сути провести сначала однофакторную регрессию для произведений, отобрать возможные предикторы по р<0,02, а затем подставлять их в многофакторную модель).

Стоит ли довериться stepwise или пробовать комбинации вручную?

преимущественно в зарубежных статьях часто наблюдаю "no interactions were found", но что за этим скрыто, и как авторы проводили анализ, никогда не указывается.

спасибо!

PS.: Должно ли распределение предикторов (и остатков, и как их получить в SPSS или Statistica) в регрессии Кокса соответствовать нормальному?

Автор: p2004r 31.10.2014 - 15:22

Определение влияет или нет данная переменная (в том числе в сочетании с другими) на точность прогноза самостоятельный раздел анализа данных -- http://en.wikipedia.org/wiki/Feature_selection .

Например можно использовать library(Boruta) в R. Поиском на форуме легко находиться пример.

Автор: mix3d 31.10.2014 - 22:05

Спасибо за ответ, не уверен, что это то, что нужно. Возможно, если нужно выбрать предикторы из большого количества переменных, да, но не в моем случае. Как я уже писал, всего выявлено 3 независимых предиктора (результаты согласованы по полу и возрасту), вопрос в том, как правильно и с наименьшими потерями времени, выполнить анализ взаимодействий.

Автор: p2004r 1.11.2014 - 19:26

Цитата(mix3d @ 31.10.2014 - 22:05)

1) Вот вместо "выявлено" и примените.

2) А анализ взаимодействия принципиально проводится по всем переменным системы. Критерий "влиятельные для данного показателя и следовательно составляют систему" не пройдет, там имеет место "сеть взаимодействий".

Например library(bnlearn) можно применить для начала.

Автор: mix3d 1.11.2014 - 21:24

Цитата(p2004r @ 1.11.2014 - 19:26)

1) Вот вместо "выявлено" и примените.

Например library(bnlearn) можно применить для начала.

k sozhaleniju, ne vladeju etimi metodami, i net vremeni razbirat'sja

"Simplicity is the best"

Цитата(p2004r @ 1.11.2014 - 19:26)

А анализ взаимодействия принципиально проводится по всем переменным системы.

Vot, sobstvenno, i otvet, "analiz vseh vozmozhnyh kombinacij"

V ljubom sluchae, bol'shoe spasibo za sovet, mozhet byt', v buduschem udastsja vniknut'