Фиттинг непараметрических распределений по данным выживаемости 

Любое поддерживаемое непараметрическое распределение может быть построено непосредственно по цензурированным или усечённым данным. Зададим дважды (слева и справа) некоторые цензурированные данные о времени наступления изучаемых событий.

In[1]:=

X \[ScriptCapitalD] = EventData[{10, 11, 12, 13, 14, 15, 16, 17, 18}, {{4, 0, 0}, {12, 0, 0}, {19, 2, 0}, {24, 15, 1}, {20, 24, 2}, {13, 18, 3}, {3, 14, 2}, {1, 6, 3}, {0, 4, 1}}];

Построим по этим данным непараметрические распределения и найдём их средние значения.

In[2]:=

X edist = EmpiricalDistribution[\[ScriptCapitalD]]; hdist = HistogramDistribution[\[ScriptCapitalD], {.8}]; sdist = SmoothKernelDistribution[\[ScriptCapitalD], .3]; kdist = KernelMixtureDistribution[\[ScriptCapitalD]]; Mean /@ {edist, hdist, sdist, kdist}

Out[3]=

Оценки совокупной интенсивности отказов по дважды цензурированным данным, хранимым в объекте EventData.