Мейєра (порівняно з методом таблиць життя) полягає в тому, що оцінки не залежать від розбиття часу спостереження на інтервали, тобто від угруповання. Метод розмножувальних оцінок і метод таблиць часів життя призводять, по суті, до однакових результатів, якщо тимчасові інтервали містять, максимум, по одному спостереженню. p align="justify"> Вибір найбільш відповідного розбиття - непросте завдання. Оцінки значень величин, отриманих за методом Каплана-Мейєра не залежать від розбиття часу спостереження на інтервали, а залежать тільки від часу життя кожного окремого пацієнта. Тому досліднику простіше проводити аналіз, та й результати нерідко виявляються якісніше результатів В«звичайногоВ» аналізу. br/>
2. Програма-функція
k2 (lambda, n, m)% вхідні параметри
a = exprnd ((1/lambda), 1, n);% вибірка = sort (a);
a1 = a; = randperm (n);% тимчасова змінна для цензурірованіяm> 0 (t (1: m)) = [];% цензурування
end
for j = 1: (n-m)
r (j) = sum (a> = a1 (j));% число відпрацьованих елементів
end (1) = 1; (1) = 1; j = 1: (n-m)
s (j +1) = (s (j)) * (1 - (1/r (j)));% ф-ія дожиття Каплана-Мейєра
sigma = s * sqrt (sum (t1));% оцінка точності за ф-ле Грінвуда
hold on (a2, sigma) (1) = []; = (t1)./(log (s). ^ 2); = log (-log (s)) -1.96. * sqrt (v); = log (-log (s)) +1.96. * sqrt (v);
c1 = exp (-exp (b2));% нижня межа довірчого інтервалу = exp (-exp (b1));% верхня межа довірчого інтервалу
stairs (a1, c1, 'r-')% зображення нижньої межі довірчого інтервалу
stairs (a1, c2, 'g-')% зображення верхньої межі довірчого інтервалу
xlabel ('t')% підпис осі x
ylabel ('S (t)')% підпис осі y
Обчислення для значень параметрів
? 0.50.50.5n1050150m0515
Графіки точної функції дожиття для значень параметрів ?, n, m та оцінки Каплана-Мейєра з 95% довірчими інтервалами у фіксованих точках, розрахованих за допомогою формули Грінвуда.
. Графіки: S (зображений світло-синім кольором), нижній і верхній інтервали для S (t) (показані жовтим і зеленими кольорами відповідно) і оцінка точності за формулою Грінвуда (показана синім кольором), при значеннях ? = 0.5, n = 10, m = 0 ( рис. 1).
В
Рис. 1
. Графіки: S (зображений світло-синім кольором), нижній і верхні...
