L'énoncé du problème est que je veux pouvoir comparer la survie de 2 populations et donner un intervalle de confiance autour de l'estimation ponctuelle. JMP peut tracer un intervalle de confiance simultané ou par paire autour d'une population donnée, mais ne peut pas donner d'intervalle de confiance autour de la différence. Par conséquent, je suggère à tous les membres de mon entreprise d'utiliser R ou Python pour effectuer une analyse de survie, car le but de comparer la survie est de quantifier la différence et l'incertitude de la différence.
Le RMST est un moyen de comparer la survie de 2 populations et d'en dériver des intervalles de confiance.
https://cran.r-project.org/web/packages/survRM2/vignettes/survRM2-vignette3-2.html
Je suggère à JMP d'implémenter RMST.
Merci pour votre scénario. Cela entre les intervalles de confiance autour de chaque population indépendamment. Je demande un intervalle de confiance autour de la différence entre 2 populations, ce qui n'est pas la même chose et est souvent mal utilisé. Pour des exemples, veuillez consulter les Fig 1c et Fig 3c dans cet article : https://doi.org/10.1038/nmeth.2974
La perspective Comparer les groupes de la plate-forme Life Distribution peut vous donner le résultat souhaité pour comparer la capacité de survie entre les groupes.Essayez d'exécuter ce script qui compare deux groupes de distributions de durée de vie et, à l'aide de l'option Estimer la probabilité, compare les taux de survie à des intervalles spécifiques (par exemple, 100, 200, 500 et 1000 cycles de temps).
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" )}
)
)
);
J'ai une implémentation fonctionnelle en python, ici: https://github.com/bayesfactor/lifelines/blob/f447b266549ee666eef37df69469473ba5073076/lifelines/statistics.py#L409
J'espère que cela vous aidera à mettre en œuvre dans JMP. Pouvez-vous s'il vous plaît fournir une mise à jour sur l'emplacement de cette fonctionnalité dans votre file d'attente de développement ? Merci
