В этом примере показано, как оценить и построить график кумулятивной опасности и функций выживших для различных групп.
Загрузите образцы данных.
load('readmissiontimes.mat')
Данные имеют время реадмиссии пациентов с информацией об их поле, возрасте, весе, статусе курения и цензуре. Это смоделированные данные.
Создайте матрицу времени реадмиссии и цензуры для каждого пола.
female = [ReadmissionTime(Sex==1),Censored(Sex==1)]; male = [ReadmissionTime(Sex==0),Censored(Sex==0)];
Постройте график оценки Каплана-Мейера кумулятивной функции распределения для пациентов женского и мужского пола.
figure() ecdf(gca,female(:,1),'Censoring',female(:,2)); hold on [f,x] = ecdf(male(:,1),'Censoring',male(:,2)); stairs(x,f,'--r') hold off legend('female','male','Location','SouthEast')
Сравните функции выживших для пациентов женского и мужского пола.
figure() ax1 = gca; ecdf(ax1,female(:,1),'Censoring',female(:,2),'function','survivor'); hold on [f,x] = ecdf(male(:,1),'Censoring',male(:,2),'function','survivor'); stairs(x,f,'--r') legend('female','male')
Эта цифра показывает, что время реадмиссии для пациентов мужского пола короче, чем для пациентов женского пола.
Подберите распределения Вейбулла для времени реадмиссии пациентов женского и мужского пола.
pd = fitdist(female(:,1),'wbl','Censoring',female(:,2))
pd =
WeibullDistribution
Weibull distribution
A = 12.5593 [10.749, 14.6745]
B = 1.99834 [1.56489, 2.55185]
pd2 = fitdist(male(:,1),'wbl','Censoring',male(:,2))
pd2 =
WeibullDistribution
Weibull distribution
A = 4.63991 [3.91039, 5.50551]
B = 1.94422 [1.48496, 2.54552]
pd2 = fitdist(male(:,1),'wbl','Censoring',male(:,2))
pd2 =
WeibullDistribution
Weibull distribution
A = 4.63991 [3.91039, 5.50551]
B = 1.94422 [1.48496, 2.54552]
Постройте график функций выжившего Вейбулла для пациентов женского и мужского пола на оценочных функциях выжившего.
plot(0:1:25,1-cdf('wbl',0:1:25,12.5593,1.99834),'-.') plot(0:1:25,1-cdf('wbl',0:1:25,4.63991,1.94422),':r') hold off legend('Festimated','Mestimated','FWeibull','MWeibull')
Распределение Вейбулла обеспечивает хорошее соответствие данных.
Оценка кумулятивной функции опасности для полов и подбор кумулятивных функций опасности Вейбулла.
figure() [f,x] = ecdf(female(:,1),'Censoring',female(:,2),... 'function','cumhazard'); plot(x,f) hold on plot(x,cumsum(pdf(pd,x)./(1-cdf(pd,x))),'-.') [f,x] = ecdf(male(:,1),'Censoring',male(:,2),... 'function','cumhazard'); plot(x,f,'--r') plot(x,cumsum(pdf(pd2,x)./(1-cdf(pd2,x))),':r') legend('Festimated','FWeibull','Mestimated','MWeibull',... 'Location','North')
