Документация

Введение

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

Анализ выживания состоит из параметрических, полупараметрических, и непараметрических методов. Можно использовать их, чтобы оценить обычно используемые меры в исследованиях выживания, оставшемся в живых и функциях опасности, сравнить их для различных групп и оценить отношение переменных предикторов ко времени выживания. Некоторые статистические вероятностные распределения описывают времена выживания хорошо. Обычно используемые распределения экспоненциальны, Weibull, логарифмически нормальны, Берр и распределения Бирнбаума-Сондерса. Функции Statistics and Machine Learning Toolbox™ ecdf и ksdensity вычислите эмпирическое и оценки плотности ядра cdf, совокупной опасности и функций оставшегося в живых. coxphfit соответствует Cox пропорциональная модель опасностей к данным.

Цензурирование

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

Во время исследования выживания или индивидуум, как наблюдают, приводит к сбою во время T, или наблюдение относительно того индивидуума прекращает во время c. Затем наблюдение является min (T, c), и переменная I _c индикатора показывает, подвергается ли индивидуум цензуре или нет. Вычисления для опасности и функций оставшегося в живых должны быть настроены с учетом цензурирования. Statistics and Machine Learning Toolbox функционирует, такие как ecdf, ksdensity, coxphfit, и mle объясните цензурирование.

Данные

Данные о выживании обычно состоят из времени, пока мероприятие не происходит и информация о цензурировании для каждого индивидуума или компонента. Следующая таблица показывает фиктивное время безработицы индивидуумов в 6-месячном исследовании. Два индивидуума правы подвергнутый цензуре (обозначенный значением цензурирования 1). Один индивидуум был все еще безработным после 24-й недели, когда законченное исследование. Контакт с другим подвергнутым цензуре индивидуумом был потерян в конце 21-й недели.

Данные о выживании могут также включать количество отказов в определенное время (число раз, конкретное время выживания или отказа наблюдалось). Следующая таблица показывает симулированное время до спады светодиодов до 70% его полного уровня светоотдачи, в часах, в ускоренном жизненном тесте.

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

Функция оставшегося в живых

Функция оставшегося в живых является вероятностью выживания как функция времени. Это также вызвано функция выживания. Это дает вероятность, что время выживания индивидуума превышает определенное значение. Поскольку кумулятивная функция распределения, F (t), является вероятностью, что время выживания меньше чем или равно данному моменту времени, функция выживания для непрерывного распределения, S (t), является дополнением кумулятивной функции распределения:

S (t) = 1 – F (t).

Функция оставшегося в живых также связана с функцией опасности. Если данные имеют функцию опасности, h (t), то функция оставшегося в живых

который соответствует

где H (t) является совокупной функцией опасности.

x = 0:0.1:200;
figure()
plot(x,1-cdf('Burr',x,50,3,1))
xlabel('Failure time');
ylabel('Survival probability');

load readmissiontimes

[ReadmissionTime Censored]

ans = 100×2

     5     1
     3     1
    19     0
    17     0
     9     0
    16     0
     4     0
     2     0
     3     0
    15     0
      ⋮

ecdf(ReadmissionTime,'censoring',Censored,'function','survivor')

Функция опасности

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

где Δt является очень маленьким временным интервалом. Показатель риска, поэтому, иногда называется условной интенсивностью отказов. Функция опасности всегда принимает положительное значение. Однако эти значения не соответствуют вероятностям и могут быть больше 1.

Функция опасности связана с функцией плотности вероятности, f (t), кумулятивная функция распределения, F (t) и функция оставшегося в живых, S (t), можно следующим образом:

который также эквивалентен

Так, если вы знаете форму функции выживания, можно также вывести соответствующую функцию опасности.

x = 0:1:200;
Burrhazard = pdf('Burr',x,50,3,1)./(1-cdf('Burr',x,50,3,1));
figure()
plot(x,Burrhazard)
xlabel('Failure time');
ylabel('Hazard rate');

figure
ax1 = subplot(3,1,1);
x1 = 0:0.05:10;
hazard1 = pdf('wbl',x1,3,0.6)./(1-cdf('wbl',x1,3,0.6));
plot(x1,hazard1,'color','b')
set(ax1,'Ylim',[0 0.6]);
legend(ax1,'a=3, b=0.6');

ax2 = subplot(3,1,2);
x2 = 0:0.05:10;
hazard2 = pdf('wbl',x2,9,4)./(1-cdf('wbl',x2,9,4));
plot(x2,hazard2,'color','r')
set(ax2,'Ylim',[0 0.6]);
legend(ax2,'a=9, b=4','location','southeast');

ax3 = subplot(3,1,3);
x3 = 0:0.05:10;
hazard3 = pdf('wbl',x3,2.5,1)./(1-cdf('wbl',x3,2.5,1));
plot(x3,hazard3,'color','g')
set(ax3,'Ylim',[0 0.6]);
legend(ax3,'a=2.5, b=1');