Модель пропорциональных рисков

Модели пропорциональных рисков представляют собой класс моделей выживания в статистике . Модели выживания связывают время, прошедшее до того, как произойдет какое-либо событие, с одной или несколькими ковариатами , которые могут быть связаны с этим количеством времени. В модели пропорциональных рисков уникальный эффект увеличения ковариаты на единицу мультипликативен по отношению к уровню риска . Например, прием лекарств может вдвое снизить риск возникновения инсульта, а изменение материала, из которого изготовлен компонент, может удвоить риск отказа. Другие типы моделей выживания, такие как модели ускоренного времени отказа, не демонстрируют пропорциональных рисков. Модель ускоренного времени отказа описывает ситуацию, когда биологическая или механическая история жизни события ускоряется (или замедляется).

Предыстория

Модели выживания можно рассматривать как состоящие из двух частей: лежащая в их основе базовая функция риска , часто обозначаемая $\lambda _{0}(t)$, описывающее, как риск события в единицу времени меняется с течением времени на базовых уровнях ковариат; и параметры эффекта, описывающие, как опасность меняется в зависимости от поясняющих ковариат. Типичный медицинский пример может включать такие ковариаты, как назначение лечения, а также характеристики пациента, такие как возраст на момент начала исследования, пол и наличие других заболеваний на момент начала исследования, чтобы уменьшить вариабельность и/или контролировать возникновение искажающих факторов.

Условие пропорциональных рисков утверждает, что ковариаты мультипликативно связаны с опасностью. Например, в простейшем случае стационарных коэффициентов лечение лекарством может, скажем, вдвое снизить риск для субъекта в любой момент времени. $t$, хотя базовый уровень опасности может варьироваться. Однако обратите внимание, что это не удваивает продолжительность жизни субъекта; точное влияние ковариат на время жизни зависит от типа $\lambda _{0}(t)$. Ковариата ; не ограничивается бинарными предикторами в случае непрерывной ковариаты $x$обычно предполагается, что опасность реагирует экспоненциально; каждая единица увеличивается в $x$ приводит к пропорциональному масштабированию опасности.

Модель Кокса

Введение

Сэр Дэвид Кокс заметил, что если предположение о пропорциональных рисках справедливо (или предполагается, что оно справедливо), то можно оценить параметр(ы) эффекта, обозначенный $\beta _{i}$ ниже, без учета полной функции опасности. Такой подход к данным о выживаемости называется применением модели пропорциональных рисков Кокса . иногда сокращается до модели Кокса или модели пропорциональных рисков . Однако Кокс также отметил, что биологическая интерпретация предположения о пропорциональных рисках может быть довольно сложной.

Пусть X _i = ( X _{i 1} , …, X _ip ) будут реализованными значениями p ковариат для субъекта i . Функция риска для модели пропорциональных рисков Кокса имеет вид

$${\begin{aligned}\lambda (t|X_{i})&=\lambda _{0}(t)\exp(\beta _{1}X_{i1}+\cdots +\beta _{p}X_{ip})\\&=\lambda _{0}(t)\exp(X_{i}\cdot \beta )\end{aligned}}$$

Это выражение дает функцию риска в момент времени t для субъекта i с вектором ковариат (независимыми переменными) X _i . Обратите внимание, что между субъектами базовый риск $\lambda _{0}(t)$ идентичен (не зависит от i ). Единственная разница между опасностями субъектов заключается в базовом масштабном коэффициенте. $\exp(X_{i}\cdot \beta )$.

Почему это называется «пропорциональным»

Для начала предположим, что у нас есть только одна ковариата, $x$, и, следовательно, один коэффициент, $\beta _{1}$. Рассмотрим эффект увеличения $x$ на 1:

$${\begin{aligned}\lambda (t|x+1)&=\lambda _{0}(t)\exp(\beta _{1}(x+1))\\&=\lambda _{0}(t)\exp(\beta _{1}x+\beta _{1})\\&={\Bigl (}\lambda _{0}(t)\exp(\beta _{1}x){\Bigr )}\exp(\beta _{1})\\&=\lambda (t|x)\exp(\beta _{1})\end{aligned}}$$

Мы видим, что увеличение ковариаты на 1 масштабирует исходную опасность на константу $\exp(\beta _{1})$. Немного переставив ситуацию, мы видим, что:

$${\frac {\lambda (t|x+1)}{\lambda (t|x)}}=\exp(\beta _{1})$$

Правая часть постоянна во времени (ни один член не имеет $t$ в нем). Эти отношения, $x/y={\text{constant}}$, называется пропорциональной зависимостью .

В более общем плане рассмотрим два предмета, i и j , с ковариатами $X_{i}$ и $X_{j}$ соответственно. Рассмотрим соотношение их опасностей:

$${\begin{aligned}{\frac {\lambda (t|X_{i})}{\lambda (t|X_{j})}}&={\frac {\lambda _{0}(t)\exp(X_{i}\cdot \beta )}{\lambda _{0}(t)\exp(X_{j}\cdot \beta )}}\\&={\frac {{\cancel {\lambda _{0}(t)}}\exp(X_{i}\cdot \beta )}{{\cancel {\lambda _{0}(t)}}\exp(X_{j}\cdot \beta )}}\\&=\exp((X_{i}-X_{j})\cdot \beta )\end{aligned}}$$

Правая часть не зависит от времени, как единственного фактора, зависящего от времени. $\lambda _{0}(t)$, был отменен. Таким образом, соотношение опасностей двух субъектов является постоянным, т.е. опасности пропорциональны.

Отсутствие термина перехвата

Часто в моделях регрессии используется термин-перехват (также называемый постоянным термином или термином смещения). В модели Кокса его нет, потому что базовый риск, $\lambda _{0}(t)$, занимает его место. Давайте посмотрим, что произойдет, если мы все равно добавим термин перехвата, обозначаемый $\beta _{0}$:

$${\begin{aligned}\lambda (t|X_{i})&=\lambda _{0}(t)\exp(\beta _{1}X_{i1}+\cdots +\beta _{p}X_{ip}+\beta _{0})\\&=\lambda _{0}(t)\exp(X_{i}\cdot \beta )\exp(\beta _{0})\\&=\left(\exp(\beta _{0})\lambda _{0}(t)\right)\exp(X_{i}\cdot \beta )\\&=\lambda _{0}^{*}(t)\exp(X_{i}\cdot \beta )\end{aligned}}$$

где мы переопределили $\exp(\beta _{0})\lambda _{0}(t)$ стать новой базовой опасностью, $\lambda _{0}^{*}(t)$. Таким образом, базовая опасность включает в себя все части опасности, которые не зависят от ковариат субъектов, включая любые пересекающиеся термины (которые по определению постоянны для всех субъектов).

Вероятность уникальных моментов времени

Кокса Частичная вероятность , показанная ниже, получается путем использования оценки Бреслоу базовой функции риска, подстановки ее к полной вероятности и последующего наблюдения, что результат является произведением двух факторов. Первым фактором является показанная ниже частичная вероятность, при которой базовый риск «отменился». Второй фактор свободен от коэффициентов регрессии и зависит от данных только посредством шаблона цензурирования . Таким образом, влияние ковариат, оцененных с помощью любой модели пропорциональных рисков, можно представить как коэффициенты рисков .

Вероятность того, что наблюдаемое событие произойдет для субъекта i в момент времени Y _i, можно записать как:

$$L_{i}(\beta )={\frac {\lambda (Y_{i}\mid X_{i})}{\sum _{j:Y_{j}\geq Y_{i}}\lambda (Y_{i}\mid X_{j})}}={\frac {\lambda _{0}(Y_{i})\theta _{i}}{\sum _{j:Y_{j}\geq Y_{i}}\lambda _{0}(Y_{i})\theta _{j}}}={\frac {\theta _{i}}{\sum _{j:Y_{j}\geq Y_{i}}\theta _{j}}},$$

где θ _j = exp( X _j ⋅ β ), а суммирование ведется по множеству субъектов j , где событие не произошло до момента Y _i (включая самого субъекта i ). Очевидно, 0 < L _i (β) ≤ 1. Это частичная вероятность : влияние ковариат можно оценить без необходимости моделирования изменения опасности с течением времени.

Рассматривая испытуемых так, как если бы они были статистически независимы друг от друга, общая вероятность всех реализованных событий — это следующая частичная вероятность, где возникновение события обозначается C _i = 1:

$$L(\beta )=\prod _{i:C_{i}=1}L_{i}(\beta ).$$

Соответствующее логарифмическое частичное правдоподобие равно

$$\ell (\beta )=\sum _{i:C_{i}=1}\left(X_{i}\cdot \beta -\log \sum _{j:Y_{j}\geq Y_{i}}\theta _{j}\right).$$

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

Частичная функция оценки :

$$\ell ^{\prime }(\beta )=\sum _{i:C_{i}=1}\left(X_{i}-{\frac {\sum _{j:Y_{j}\geq Y_{i}}\theta _{j}X_{j}}{\sum _{j:Y_{j}\geq Y_{i}}\theta _{j}}}\right),$$

а матрица Гессе частичного логарифмического правдоподобия равна

$$\ell ^{\prime \prime }(\beta )=-\sum _{i:C_{i}=1}\left({\frac {\sum _{j:Y_{j}\geq Y_{i}}\theta _{j}X_{j}X_{j}^{\prime }}{\sum _{j:Y_{j}\geq Y_{i}}\theta _{j}}}-{\frac {\left[\sum _{j:Y_{j}\geq Y_{i}}\theta _{j}X_{j}\right]\left[\sum _{j:Y_{j}\geq Y_{i}}\theta _{j}X_{j}^{\prime }\right]}{\left[\sum _{j:Y_{j}\geq Y_{i}}\theta _{j}\right]^{2}}}\right).$$

Используя эту оценочную функцию и матрицу Гессе, частичную вероятность можно максимизировать с помощью алгоритма Ньютона-Рафсона . Обратная матрица Гессе, оцененная по оценке β , может использоваться в качестве приблизительной дисперсионно-ковариационной матрицы для оценки и использоваться для получения приблизительных стандартных ошибок для коэффициентов регрессии.

Вероятность, когда существуют равные времена

Было предложено несколько подходов для решения ситуаций, в которых имеются связи во временных данных. Метод Бреслоу описывает подход, при котором описанная выше процедура используется без изменений, даже при наличии связей. Альтернативным подходом, который, как считается, дает лучшие результаты, является метод Эфрона . Пусть t _j обозначает уникальные моменты времени, пусть H _j обозначает набор индексов i таких, что Y _i = t _j ​​и C _i = 1, и пусть m _j = | Ч _Дж |. Подход Эфрона максимизирует следующую частичную вероятность.

$$L(\beta )=\prod _{j}{\frac {\prod _{i\in H_{j}}\theta _{i}}{\prod _{\ell =0}^{m_{j}-1}\left[\sum _{i:Y_{i}\geq t_{j}}\theta _{i}-{\frac {\ell }{m_{j}}}\sum _{i\in H_{j}}\theta _{i}\right]}}.$$

Соответствующее логарифмическое частичное правдоподобие равно

$$\ell (\beta )=\sum _{j}\left(\sum _{i\in H_{j}}X_{i}\cdot \beta -\sum _{\ell =0}^{m_{j}-1}\log \left(\sum _{i:Y_{i}\geq t_{j}}\theta _{i}-{\frac {\ell }{m_{j}}}\sum _{i\in H_{j}}\theta _{i}\right)\right),$$

функция оценки

$$\ell ^{\prime }(\beta )=\sum _{j}\left(\sum _{i\in H_{j}}X_{i}-\sum _{\ell =0}^{m_{j}-1}{\frac {\sum _{i:Y_{i}\geq t_{j}}\theta _{i}X_{i}-{\frac {\ell }{m_{j}}}\sum _{i\in H_{j}}\theta _{i}X_{i}}{\sum _{i:Y_{i}\geq t_{j}}\theta _{i}-{\frac {\ell }{m_{j}}}\sum _{i\in H_{j}}\theta _{i}}}\right),$$

и матрица Гессе есть

$$\ell ^{\prime \prime }(\beta )=-\sum _{j}\sum _{\ell =0}^{m_{j}-1}\left({\frac {\sum _{i:Y_{i}\geq t_{j}}\theta _{i}X_{i}X_{i}^{\prime }-{\frac {\ell }{m_{j}}}\sum _{i\in H_{j}}\theta _{i}X_{i}X_{i}^{\prime }}{\phi _{j,\ell ,m_{j}}}}-{\frac {Z_{j,\ell ,m_{j}}Z_{j,\ell ,m_{j}}^{\prime }}{\phi _{j,\ell ,m_{j}}^{2}}}\right),$$

где

$$\phi _{j,\ell ,m_{j}}=\sum _{i:Y_{i}\geq t_{j}}\theta _{i}-{\frac {\ell }{m_{j}}}\sum _{i\in H_{j}}\theta _{i}$$

$$Z_{j,\ell ,m_{j}}=\sum _{i:Y_{i}\geq t_{j}}\theta _{i}X_{i}-{\frac {\ell }{m_{j}}}\sum _{i\in H_{j}}\theta _{i}X_{i}.$$

Обратите внимание, что когда H _j пусто (все наблюдения за время t _j подвергаются цензуре), слагаемые в этих выражениях рассматриваются как ноль.

Примеры

Ниже приведены некоторые примеры применения модели Кокса на практике.

Одна бинарная ковариата

Предположим, что интересующая нас конечная точка — это выживаемость пациентов в течение 5-летнего периода наблюдения после операции. Пациенты могут умереть в течение 5 лет, и мы записываем, когда они умерли, или пациенты могут прожить более 5 лет, и мы записываем только то, что они прожили более 5 лет. Операция была проведена в одной из двух больниц, A или B , и мы хотели бы знать, связано ли расположение больницы с 5-летней выживаемостью. В частности, мы хотели бы знать относительное увеличение (или уменьшение) риска от операции, проведенной в больнице А по сравнению с больницей Б. Предоставляются некоторые (поддельные) данные, где каждая строка представляет пациента: T — как долго пациент находился в больнице. наблюдался до смерти или в течение 5 лет (измеряется в месяцах), а C обозначает, если пациент умер в течение 5-летнего периода. Мы закодировали больницу как двоичную переменную, обозначенную X если из больницы A , 0, если из больницы B. : 1 ,

больница	Х	Т	С
Б	0	60	ЛОЖЬ
Б	0	32	Истинный
Б	0	60	ЛОЖЬ
Б	0	60	ЛОЖЬ
Б	0	60	ЛОЖЬ
А	1	4	Истинный
А	1	18	Истинный
А	1	60	ЛОЖЬ
А	1	9	Истинный
А	1	31	Истинный
А	1	53	Истинный
А	1	17	Истинный

Наша одноковариатная пропорциональная модель Кокса выглядит следующим образом: $\beta _{1}$ представляю эффект больницы, и я индексирую каждого пациента:

$$\overbrace {\lambda (t|X_{i})} ^{\text{hazard for i}}=\underbrace {\lambda _{0}(t)} _{{\text{baseline}} \atop {\text{hazard}}}\cdot \overbrace {\exp(\beta _{1}X_{i})} ^{\text{scaling factor for i}}$$

Используя статистическое программное обеспечение, мы можем оценить $\beta _{1}$ будет 2.12. Коэффициент риска представляет собой экспоненту этой величины, $\exp(\beta _{1})=\exp(2.12)$. Чтобы понять почему, рассмотрим соотношение опасностей, а именно:

$${\frac {\lambda (t|X=1)}{\lambda (t|X=0)}}={\frac {{\cancel {\lambda _{0}(t)}}\exp(\beta _{1}\cdot 1)}{{\cancel {\lambda _{0}(t)}}\exp(\beta _{1}\cdot 0)}}=\exp(\beta _{1})$$

Таким образом, отношение рисков больницы А к больнице Б составляет $\exp(2.12)=8.32$. Отложив на мгновение статистическую значимость, мы можем заявить, что у пациентов в больнице А риск смерти в любой короткий период времени в 8,3 раза выше, чем в больнице Б.

При интерпретации следует сделать несколько важных замечаний:

1. Более высокий риск смерти в 8,3 раза не означает, что в больнице Б умрет в 8,3 раза больше пациентов: анализ выживаемости изучает, насколько быстро происходят события, а не просто происходят ли они.
2. Более конкретно, «риск смерти» является мерой уровня. Скорость имеет единицы измерения, например, метры в секунду. Однако относительная скорость не имеет значения: велосипед может двигаться в два раза быстрее, чем другой велосипед (эталонный велосипед), без указания каких-либо единиц измерения. Аналогично, риск смерти (скорость смертности) в больнице А в 8,3 раза выше (быстрее), чем риск смерти в больнице Б (контрольная группа).
3. обратная величина, $1/8.32={\frac {1}{\exp(2.12)}}=\exp(-2.12)=0.12$ коэффициент риска больницы B по отношению к больнице A. —
4. Мы не сделали никаких выводов о вероятности выживания между больницами. Это потому, что нам понадобится оценка базового уровня опасности, $\lambda _{0}(t)$, а также наш $\beta _{1}$ оценивать. Однако стандартная оценка модели пропорционального риска Кокса не позволяет напрямую оценить базовый уровень риска.
5. Поскольку мы проигнорировали единственный изменяющийся во времени компонент модели — базовый уровень опасности, наша оценка не зависит от масштаба времени. Например, если бы мы измеряли время в годах, а не в месяцах, мы получили бы ту же оценку.
6. Заманчиво сказать, что именно больница стала причиной разницы в рисках между двумя группами, но поскольку наше исследование не является причинно-следственным (то есть мы не знаем, как были получены данные), мы придерживаемся такой терминологии, как «связанный».

Одна непрерывная ковариата

Чтобы продемонстрировать менее традиционный вариант использования анализа выживания, следующим примером будет экономический вопрос: какова связь между соотношением цены и прибыли (P/E) компании в первую годовщину IPO и ее будущим выживанием? Точнее, если мы рассматриваем «событие рождения» компании как ее первую годовщину IPO, а любое банкротство, продажу, переход в частную собственность и т. д. как событие «смерти» компании, мы хотели бы знать влияние компаний ' Соотношение цена/прибыль на момент их «рождения» (первая годовщина IPO) зависит от их выживания.

Предоставляется (поддельный) набор данных с данными о выживании 12 компаний: T представляет собой количество дней между первой годовщиной IPO и смертью (или датой окончания 01 января 2022 г., если компания не умерла). C означает, умерла ли компания до 1 января 2022 г. или нет. P/E представляет собой соотношение цены и прибыли компании на момент ее первой годовщины IPO.

Ко.	Дата IPO через 1 год	Дата смерти*	С	Т	Цена/прибыль
0	2000-11-05	2011-01-22	Истинный	3730	9.7
1	2000-12-01	2003-03-30	Истинный	849	12.0
2	2011-01-05	2012-03-30	Истинный	450	3.0
3	2010-05-29	2011-02-22	Истинный	269	5.3
4	2005-06-23	2022-01-01	ЛОЖЬ	6036	10.8
5	2000-06-10	2002-07-24	Истинный	774	6.3
6	2011-07-11	2014-05-01	Истинный	1025	11.6
7	2007-09-27	2022-01-01	ЛОЖЬ	5210	10.3
8	2006-07-30	2010-06-03	Истинный	1404	8.0
9	2000-07-13	2001-07-19	Истинный	371	4.0
10	2013-06-10	2018-10-10	Истинный	1948	5.9
11	2011-07-16	2014-08-15	Истинный	1126	8.3

В отличие от предыдущего примера, где использовалась двоичная переменная, в этом наборе данных есть непрерывная переменная P/E; однако модель выглядит похоже:

$$\lambda (t|P_{i})=\lambda _{0}(t)\cdot \exp(\beta _{1}P_{i})$$

где $P_{i}$ представляет собой коэффициент P/E компании. Пропуск этого набора данных через модель Кокса дает оценку значения неизвестного. $\beta _{1}$, что составляет -0,34. Таким образом, оценка всей опасности равна:

$$\lambda (t|P_{i})=\lambda _{0}(t)\cdot \exp(-0.34P_{i})$$

Поскольку базовая опасность, $\lambda _{0}(t)$, не была оценена, всю опасность невозможно рассчитать. Однако рассмотрим соотношение рисков компаний i и j :

$${\begin{aligned}{\frac {\lambda (t|P_{i})}{\lambda (t|P_{j})}}&={\frac {{\cancel {\lambda _{0}(t)}}\cdot \exp(-0.34P_{i})}{{\cancel {\lambda _{0}(t)}}\cdot \exp(-0.34P_{j})}}\\&=\exp(-0.34(P_{i}-P_{j}))\end{aligned}}$$

Все условия справа известны, поэтому можно рассчитать соотношение опасностей между компаниями. Поскольку справа нет члена, зависящего от времени (все члены постоянны), опасности пропорциональны друг другу. Например, отношение рисков компании 5 к компании 2 равно $\exp(-0.34(6.3-3.0))=0.33$. Это означает, что в интервале исследования риск «смерти» компании 5 составляет 0,33 ≈ 1/3 от риска смерти компании 2.

При интерпретации следует сделать несколько важных замечаний:

1. Отношение рисков – это количество $\exp(\beta _{1})$, что $\exp(-0.34)=0.71$ в приведенном выше примере. Судя по последнему расчету, приведенному выше, это можно интерпретировать как соотношение опасностей между двумя «субъектами», переменные которых отличаются на одну единицу: если $P_{i}=P_{j}+1$, затем $\exp(\beta _{1}(P_{i}-P_{j})=\exp(\beta _{1}(1))$. Выбор «отличаться на одну единицу» удобен, поскольку он точно передает значение $\beta _{1}$.
2. Базовую опасность можно представить, если масштабный коэффициент равен 1, т.е. $P=0$.
   $$\lambda (t|P_{i}=0)=\lambda _{0}(t)\cdot \exp(-0.34\cdot 0)=\lambda _{0}(t)$$
   Можем ли мы интерпретировать базовый риск как риск «базовой» компании, P/E которой равен 0? Такая интерпретация базовой опасности как «опасности для исходного субъекта» несовершенна, поскольку возможно, что ковариата, равная 0, невозможна. В этом приложении коэффициент P/E, равный 0, не имеет смысла (это означает, что цена акций компании равна 0, т. е. они «мертвы»). Более подходящей интерпретацией была бы «опасность, когда все переменные равны нулю».

1. Соблазнительно захотеть понять и интерпретировать такое значение, как $\exp(\beta _{1}P_{i})$ представлять опасность компании. Однако подумайте, что это на самом деле представляет собой: $\exp(\beta _{1}P_{i})=\exp(\beta _{1}(P_{i}-0))={\frac {\exp(\beta _{1}P_{i})}{\exp(\beta _{1}0)}}={\frac {\lambda (t|P_{i})}{\lambda (t|0)}}$. Здесь неявно существует соотношение рисков, сравнивающее риск компании i с воображаемой базовой компанией с 0 P/E. Однако, как объяснялось выше, P/E, равный 0, в этом приложении невозможен, поэтому $\exp(\beta _{1}P_{i})$ в данном примере бессмысленно. Однако соотношения между вероятными опасностями имеют важное значение.

во времени предикторы

Расширение переменных, зависящих от времени, слоев, зависящих от времени, и нескольких событий на одного субъекта может быть включено в формулировку процесса подсчета Андерсена и Гилла. Одним из примеров использования моделей рисков с изменяющимися во времени регрессорами является оценка влияния страхования по безработице на периоды безработицы.

В дополнение к разрешению изменяющихся во времени ковариатов (т.е. предикторов), модель Кокса также может быть обобщена на изменяющиеся во времени коэффициенты. То есть пропорциональный эффект лечения может меняться со временем; например, лекарство может быть очень эффективным, если его ввести в течение одного месяца после заболевания , и со временем стать менее эффективным. Затем можно проверить гипотезу об отсутствии изменения со временем (стационарности) коэффициента. Подробности и программное обеспечение ( пакет R ) доступны у Мартинуссена и Шайке (2006).

В этом контексте можно также упомянуть, что теоретически возможно определить влияние ковариат, используя аддитивные опасности, т.е. указание

$$\lambda (t|X_{i})=\lambda _{0}(t)+\beta _{1}X_{i1}+\cdots +\beta _{p}X_{ip}=\lambda _{0}(t)+X_{i}\cdot \beta .$$

Если такие модели аддитивных опасностей используются в ситуациях, когда целью является максимизация (логарифмического) правдоподобия, необходимо позаботиться о том, чтобы ограничить $\lambda (t\mid X_{i})$ до неотрицательных значений. Возможно, из-за этой сложности такие модели встречаются редко. Если вместо этого целью является метод наименьших квадратов, ограничение неотрицательности строго не требуется.

Определение базовой функции риска

Модель Кокса может быть специализированной, если существует причина предполагать, что базовая опасность имеет определенную форму. В этом случае базовая опасность $\lambda _{0}(t)$ заменяется заданной функцией. Например, если предположить, что функция риска является Вейбулла функцией риска , получается модель пропорциональных рисков Вейбулла .

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

Общий термин «параметрические модели пропорциональных рисков» можно использовать для описания моделей пропорциональных рисков, в которых указана функция риска. Модель пропорциональных рисков Кокса, иногда называют полупараметрической моделью напротив, .

Некоторые авторы используют термин « модель пропорциональных рисков Кокса» даже при указании базовой функции риска. признать долг всего поля перед Дэвидом Коксом.

Термин «модель регрессии Кокса» (без учета пропорциональных рисков ) иногда используется для описания расширения модели Кокса за счет включения факторов, зависящих от времени. Однако такое использование потенциально неоднозначно, поскольку модель пропорциональных рисков Кокса сама по себе может быть описана как регрессионная модель.

с

Существует связь между моделями пропорциональных рисков и моделями регрессии Пуассона , которая иногда используется для соответствия приблизительным моделям пропорциональных рисков в программном обеспечении для регрессии Пуассона. Обычной причиной этого является то, что расчет происходит намного быстрее. Это было более важно во времена более медленных компьютеров, но все еще может быть полезно для особенно больших наборов данных или сложных проблем. Лэрд и Оливье (1981) дайте математические подробности. Они отмечают: «Мы не предполагаем, что [модель Пуассона] верна, а просто используем ее как инструмент для определения вероятности». МакКаллах и Нелдер В книге по обобщенным линейным моделям есть глава, посвященная преобразованию моделей пропорциональных рисков в обобщенные линейные модели .

В многомерной настройке

В многомерном случае, когда количество ковариат p велико по сравнению с размером выборки n, метод LASSO является одной из классических стратегий выбора модели. Тибширани (1997) предложил процедуру Лассо для параметра регрессии пропорционального риска. Оценка Лассо параметра регрессии β определяется как минимизатор противоположности частичного логарифмического правдоподобия Кокса при L ¹ -ограничение типа нормы .

$$\ell (\beta )=\sum _{j}\left(\sum _{i\in H_{j}}X_{i}\cdot \beta -\sum _{\ell =0}^{m_{j}-1}\log \left(\sum _{i:Y_{i}\geq t_{j}}\theta _{i}-{\frac {\ell }{m_{j}}}\sum _{i\in H_{j}}\theta _{i}\right)\right)+\lambda \|\beta \|_{1},$$

В последнее время наблюдается теоретический прогресс по этой теме.

Программные реализации

• Математика : CoxModelFit функция.
• Р : coxph() функция, расположенная в пакете выживания .
• САС : phreg процедура
• Был : stcox команда
• Питон : CoxPHFitter находится в библиотеке «Линии жизни» . phreg в библиотеке статистических моделей.
• SPSS : доступен в разделе «Регрессия Кокса» .
• МАТЛАБ : fitcox или coxphfit функция
• Джулия : Доступно в библиотеке Survival.jl .
• JMP : доступен на платформе Fit Proportional Hazards .
• Призма : доступна в анализах выживания и анализах с несколькими переменными.

См. также

• Модель ускоренного времени отказа
• Правило одного из десяти
• Распределение Вейбулла
• Гипертабастическое распределение

Ссылки

• Багдонавичюс, В.; Левулин, Р.; Никулин, М. (2010). «Критерии соответствия модели Кокса на основе усеченных слева и цензурированных справа данных». Журнал математических наук . 167 (4): 436–443. дои : 10.1007/s10958-010-9929-6 . S2CID   121788950 .
• Кокс, доктор медицинских наук; Оукс, Д. (1984). Анализ данных о выживании . Нью-Йорк: Чепмен и Холл. ISBN  978-0412244902 .
• Коллетт, Д. (2003). Моделирование данных о выживаемости в медицинских исследованиях (2-е изд.). Бока-Ратон: CRC. ISBN  978-1584883258 .
• Гурьеру, Кристиан (2000). «Модели продолжительности» . Эконометрика качественных зависимых переменных . Нью-Йорк: Издательство Кембриджского университета. стр. 284–362. ISBN  978-0-521-58985-7 .
• Певица Джудит Д.; Уиллетт, Джон Б. (2003). «Подбор моделей регрессии Кокса» . Прикладной продольный анализ данных: моделирование изменений и возникновения событий . Нью-Йорк: Издательство Оксфордского университета. стр. 503–542. ISBN  978-0-19-515296-8 .
• Терно, ТМ; Грамбш, премьер-министр (2000). Моделирование данных о выживании: расширение модели Кокса . Нью-Йорк: Спрингер. ISBN  978-0387987842 .