Выполнить регрессию пропорциональных рисков Кокса на lightbulb набор данных, содержащий смоделированное время жизни лампочек. Первый столбец данных по лампочке содержит срок службы (в часах) двух различных типов лампочек. Второй столбец содержит двоичную переменную, указывающую, является ли колба флуоресцентной или накаливающей; 0 означает, что лампа является флуоресцентной, а 1 означает, что она является накаливающей. Третий столбец содержит информацию о цензуре, где 0 указывает, что лампочка наблюдалась до отказа, а 1 указывает, что наблюдение было подвергнуто цензуре.

Подберите модель пропорциональных рисков Кокса для срока службы лампочек с учетом цензуры. Переменная предиктора - тип колбы.

load lightbulb
coxMdl = fitcox(lightbulb(:,2),lightbulb(:,1), ...
    'Censoring',lightbulb(:,3));

Постройте график базовой функции выживания как функции времени t, что означает вероятность отказа лампочки через некоторое время t. По умолчанию базовая линия рассчитывается для среднего значения предиктора, которое в этом случае равно mean(lightbulb(:,2)) = 0.5.

plotSurvival(coxMdl)

Постройте график выживания флуоресцентных ламп (предиктор = 0) и ламп накаливания (предиктор = 1).

plotSurvival(coxMdl,[0;1])

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

Загрузить coxModel данные. (Эти смоделированные данные создаются в примере объекта модели пропорциональных рисков Кокса.) Модель с именем coxMdl имеет три уровня стратификации (1, 2 и 3) и предиктор X с тремя категориальными значениями (1, 1/20 и 1/100).

load coxModel

Постройте график выживания для X = 1 на трех уровнях стратификации.

c1 = categorical(1);
X = [c1;c1;c1];
stratification = [1;2;3];
plotSurvival(coxMdl,X,stratification)
xlim([1,30])

Загрузить coxModel данные. (Эти смоделированные данные создаются в примере объекта модели пропорциональных рисков Кокса.) Модель с именем coxMdl имеет три уровня стратификации (1, 2 и 3) и предиктор X с тремя категориальными значениями (1, 1/20 и 1/100).

load coxModel

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

h = figure;
axes1 = axes('Parent',h);

Постройте график функции выживания для X значение предиктора categorical(1) и уровень стратификации 3. Этот уровень стратификации представляет постоянную степень опасности. Таким образом, при логарифмировании результирующий график выживания должен быть близок к прямой линии. Постройте график для времени от 1 до 30.

oo = categorical(1);
plotSurvival(coxMdl,axes1,oo,3,'Time',linspace(1,30,300));
axes1.YScale = 'log';