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

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

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

Просмотр базовой линии по умолчанию для подогнанной модели.

defaultBaseline = coxMdl.Baseline

defaultBaseline = 0.5000

Вычислить коэффициент опасности колбы накаливания (1) относительно этой базовой линии.

defaultHazard = hazardratio(coxMdl,1)

defaultHazard = 10.6238

Вычислить коэффициент опасности колбы накаливания относительно люминесцентной колбы (0).

relHazard = hazardratio(coxMdl,1,'Baseline',0)

relHazard = 112.8646

Степень опасности колбы накаливания, по оценкам, более чем в 100 раз превышает степень опасности люминесцентной колбы.

Создание модели Кокса из readmissiontimes данные. В этих данных 0 указывает на пациента мужского пола, и 1 указывает на пациентку женского пола.

load readmissiontimes
coxMdl = fitcox([Age,Sex,Weight],ReadmissionTime,'Censoring',Censored);

Рассчитайте относительную опасность 40-летнего мужчины весом 200 фунтов относительно базовой опасности.

hazard = hazardratio(coxMdl,[40 0 200])

hazard = 4.3112

Рассчитайте опасность этого же мужчины относительно 50-летней женщины весом 150 фунтов.

hazard2 = hazardratio(coxMdl,[40 0 200],'Baseline',[50 1 150])

hazard2 = 5.2053

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

load coxModel

Найти коэффициент опасности предикторного значения categorical(1) и уровень стратификации 3 в отношении исходного уровня.

X = categorical(1);
stratification = 3;
hazard = hazardratio(coxMdl,X,stratification)

hazard = 12.7096

Вычислите отношение относительно базовой линии 0.

hazard = hazardratio(coxMdl,X,stratification,'Baseline',0)

hazard = 95.5127

Вычислить отношение categorical(1/100) предиктор относительно базовой линии 0.

X = categorical(1/100);
hazard = hazardratio(coxMdl,X,stratification,'Baseline',0)

hazard = 1