问题陈述是我希望能够比较两个群体的生存率并给出点估计的置信区间。JMP 可以绘制任何给定群体的同步或成对置信区间,但不能给出差异的置信区间。因此,我建议公司里的每个人都使用 R 或 Python 进行生存分析,因为比较生存率的目的是量化差异以及差异中的不确定性。
RMST 是比较两个群体的生存率并得出置信区间的一种方法。
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" )}
)
)
);