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

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

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 данные. (Эти симулированные данные сгенерированы в примере Cox Пропорциональный Объект модели Опасностей.) Модель под названием 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 данные. (Эти симулированные данные сгенерированы в примере Cox Пропорциональный Объект модели Опасностей.) Модель под названием 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';