問題陳述是我希望能夠比較兩個群體的生存情況並給出點估計的置信區間。 JMP 可以圍繞任何給定總體繪製同步或成對置信區間,但無法給出差異周圍的置信區間。所以我建議公司裡的大家都用R或python來做生存分析,因為比較生存的點就是量化差異,以及差異中的不確定性。
RMST 是比較 2 個群體的生存率並得出置信區間的一種方法。
https://cran.r-project.org/web/packages/survRM2/vignettes/survRM2-vignette3-2.html
我建議 JMP 實作 RMST。
謝謝你的劇本。 這獨立地輸入了每個群體的置信區間。 我要求圍繞 2 個人群的差異設置一個置信區間,這不是一回事並且經常被誤用。 例子請看本文的Fig 1c和Fig 3c: https://doi.org/10.1038/nmeth.2974
Life Distribution 平台的 Compare Groups 透視圖可能能夠為您提供所需的結果來比較各組的生存能力。嘗試運行此腳本來比較兩組生命分佈並使用“估計概率”選項,比較特定時間間隔(例如 100、200、500 和 1000 個時間週期)的存活率。
Open( "$SAMPLE_DATA/Reliability/Blenders.jmp" );
Life Distribution(
Perspective( Compare Groups ),
Y( :Time Cycles ),
Grouping( :Group ),
Censor( :Censor ),
Censor Code( 1 ),
Confidence Interval Method( Wald ),
Select Distribution( Distribution, Weibull ),
Select Scale( Weibull ),
Tabbed Report( 0 ),
Estimate Probability( 1 ),
Estimate Probability( Compute( [1000, 500, 200, 100] ) ),
SendToReport(
Dispatch(
{"Life Distribution - Compare Groups", "Statistics"},
"Group Homogeneity Tests",
OutlineBox,
{Close( 0 )}
),
Dispatch(
{"Life Distribution - Compare Groups", "Individual Group"},
"",
TabListBox,
{Set Selected( 2 )}
),
Dispatch(
{"Life Distribution - Compare Groups", "Individual Group"},
"",
List Box( 2 ),
{Set Summary Behavior( "Collapse" )}
),
Dispatch(
{"Life Distribution - Compare Groups", "Individual Group"},
"",
List Box( 76 ),
{Set Summary Behavior( "Collapse" )}
)
)
);