Выполните регрессию пропорциональных рисков Кокса на 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';