El planteamiento del problema es que quiero poder comparar la supervivencia de 2 poblaciones y dar un intervalo de confianza alrededor de la estimación puntual. JMP puede trazar un intervalo de confianza simultáneo o por pares alrededor de cualquier población determinada, pero no puede dar un intervalo de confianza alrededor de la diferencia. Por lo tanto, sugiero a todos en mi empresa que utilicen R o Python para realizar análisis de supervivencia, porque el objetivo de comparar la supervivencia es cuantificar la diferencia y la incertidumbre en la diferencia.
RMST es una forma de comparar la supervivencia entre dos poblaciones y derivar intervalos de confianza.
https://cran.r-project.org/web/packages/survRM2/vignettes/survRM2-vignette3-2.html
Sugiero que JMP implemente RMST.
Gracias por tu guión. Eso ingresa intervalos de confianza alrededor de cada población de forma independiente. Estoy pidiendo un intervalo de confianza alrededor de la diferencia en 2 poblaciones, que no es lo mismo y comúnmente se usa mal. Para ver ejemplos, consulte la Fig. 1c y la Fig. 3c de este documento: https://doi.org/10.1038/nmeth.2974
La perspectiva Comparar grupos de la plataforma Life Distribution puede brindarle el resultado deseado para comparar la capacidad de supervivencia entre grupos.Intente ejecutar este script que compara dos grupos de distribuciones de vida y, utilizando la opción Estimar probabilidad, compara las tasas de supervivencia en intervalos específicos (por ejemplo, 100, 200, 500 y 1000 ciclos de tiempo).
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" )}
)
)
);
Tengo una implementación de trabajo en python, aquí: https://github.com/bayesfactor/lifelines/blob/f447b266549ee666eef37df69469473ba5073076/lifelines/statistics.py#L409
Espero que esto te ayude a implementar en JMP. ¿Puede proporcionar una actualización sobre dónde se encuentra esta característica en su cola de desarrollo? Gracias