핵심 알고리즘

좌, 우, 구간 중도 절단 데이터 분석

우변 중도 절단 데이터와 두 그룹에서 10일때의 존속 생존 중간값을 구하기 위해 SurvivalDistribution을 이용하여 두 개의 카플란 마이어를 추정 비교해 봅니다.

In[1]:=

In[2]:=

In[3]:=

sm = SurvivalFunction[\[ScriptD]1, #1] &; snm =

SurvivalFunction[\[ScriptD]2, #1] &;

cPlot[t_, cens_, sf_, col_] :=

Block[{loc = Extract[t, Position[cens, 1]]},

Show[Table[

Graphics[{col, Line[{{i, sf[i] - .02}, {i, sf[i] + .02}}]}], {i,

loc}]]];

tbl = Grid[

Join[{Text[Style[#, Bold, FontFamily -> "Verdana"]] & /@ {

"\!\(\*SubscriptBox[\(t\), \(i\)]\)",

"\!\(\*SubscriptBox[\(P\), \

\(T\)]\)[t>\!\(\*SubscriptBox[\(t\), \(i\)]\)]",

"\!\(\*SubscriptBox[\(P\), \

\(C\)]\)[t>\!\(\*SubscriptBox[\(t\), \(i\)]\)]"}},

Table[{Text[Style[ToString@i, FontFamily -> "Verdana"]],

PaddedForm[

Probability[t > i, t \[Distributed] \[ScriptD]1], {4, 3}],

PaddedForm[

Probability[t > i, t \[Distributed] \[ScriptD]2], {4, 3}]}, {i,

Range[0, 45, 5]}]], Alignment -> Right,

Dividers -> {{False, True, False}, {False, True}},

ItemStyle -> {Automatic,

Automatic, {{{2, -1}, {2, 2}} ->

Directive[Blend[{Black, Red}],

FontFamily -> "Verdana"], {{2, -1}, {3, 3}} ->

Directive[Blend[{Black, Blue}], FontFamily -> "Verdana"]}}];

Grid[{{tbl, ,

Show[Plot[{sm[x], snm[x]}, {x, 0, 60}, PlotRange -> {0, 1},

Exclusions -> None,

PlotStyle -> {Blend[{Black, Red}], Blend[{Black, Blue}]},

Filling -> {1 -> {{2}, Automatic}}, ImageSize -> 325],

cPlot[m, cm, sm, Blend[{Black, Red}]],

cPlot[nm, cnm, snm, Blend[{Black, Blue}]]]}}]

Out[3]=

In[4]:=

Out[4]=

In[5]:=

Out[5]=