基于生存分析的元件
基于生存分析的元件
基于 EventData 创建一个 ReliabilityDistribution，并计算平均初次出故障时间.
data = RandomVariate[WeibullDistribution[3 / 2, 30], 50];
data2 = RandomVariate[ExponentialDistribution[(1/10)], 50];
ev = Table[If[i > 20, 4, i], {i, data}];
ev2 = Table[If[i > 40, 40, i], {i, data2}];
c = Table[Boole[i == 4], {i, ev}];
c2 = Table[Boole[i == 40], {i, ev2}];
{edata, edata2} = {EventData[ev, c], EventData[ev2, c2]};{Subscript[𝒟, x], Subscript[𝒟, y]} = {EstimatedDistribution[EventData[50], ExponentialDistribution[λ]],
EstimatedDistribution[EventData[50], WeibullDistribution[a, b]]};Subscript[𝒟, sys] = ReliabilityDistribution[x∧y, {{x, Subscript[𝒟, x]}, {y, Subscript[𝒟, y]}}];NExpectation[t, tSubscript[𝒟, sys]]