Wolfram
Mathematica
8의 신기능: 비모수 분포, 파생 분포, 포뮬라 분포
◄
이전
|
다음
►
핵심 알고리즘
비모수와 모수 신외성 모델 비교
SurvivalDistribution
과
CensoredDistribution
을 각각 이용하여 로프 제조에 대한 신뢰성 함수의 NPMLE와 MLE를 추정해 봅니다.
In[1]:=
X
data = {100, 43, 100, 100, 63, 26, 100, 100, 100, 100, 100, 100, 100, 69, 58, 99, 95};
In[2]:=
X
pars = FindDistributionParameters[data, CensoredDistribution[{0, 100}, WeibullDistribution[a, b]]]; sf[x_] = SurvivalFunction[ CensoredDistribution[{-Infinity, 100}, WeibullDistribution[a, b]], x] /. pars;
In[3]:=
X
\[ScriptCapitalD] = SurvivalDistribution[Censoring[data, Boole[# >= 100] & /@ data]];
In[4]:=
X
Plot[{SurvivalFunction[\[ScriptCapitalD], x], sf[x]}, {x, 0, 120}, PlotPoints -> 100, Exclusions -> None, Filling -> {1 -> {2}}, PlotStyle -> Thick, Frame -> {{True, False}, {True, False}}, FrameLabel -> (Style[#, FontFamily -> "Verdana", Bold, 12] &) /@ {"Force", "Reliability"}, ImageSize -> 500]
Out[4]=