L'affermazione del problema è che voglio essere in grado di confrontare la sopravvivenza di 2 popolazioni e fornire un intervallo di confidenza attorno alla stima puntuale. JMP può tracciare un intervallo di confidenza simultaneo o a coppie attorno a una data popolazione, ma non può fornire un intervallo di confidenza attorno alla differenza. Pertanto, suggerisco a tutti nella mia azienda di utilizzare R o Python per eseguire l'analisi di sopravvivenza, perché lo scopo del confronto della sopravvivenza è quantificare la differenza e l'incertezza nella differenza.
RMST è un modo per confrontare la sopravvivenza tra 2 popolazioni e ricavare intervalli di confidenza.
https://cran.r-project.org/web/packages/survRM2/vignettes/survRM2-vignette3-2.html
Suggerisco l'implementazione di JMP RMST.
Grazie per la tua sceneggiatura. Ciò immette gli intervalli di confidenza attorno a ciascuna popolazione in modo indipendente. Sto chiedendo un intervallo di confidenza attorno alla differenza tra 2 popolazioni, che non è la stessa cosa e comunemente usata in modo improprio. Per esempi, vedere Fig 1c e Fig 3c in questo documento: https://doi.org/10.1038/nmeth.2974
La prospettiva Confronta gruppi della piattaforma Life Distribution potrebbe essere in grado di darti il risultato desiderato per confrontare la sopravvivenza tra i gruppi.Prova a eseguire questo script che confronta due gruppi di distribuzioni di vita e utilizzando l'opzione Stima probabilità, confronta i tassi di sopravvivenza a intervalli specifici (es. 100, 200, 500 e 1000 cicli di tempo).
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" )}
)
)
);
Ho un'implementazione funzionante in Python, qui: https://github.com/bayesfactor/lifelines/blob/f447b266549ee666eef37df69469473ba5073076/lifelines/statistics.py#L409
Spero che questo ti aiuti a implementare in JMP. Puoi fornire un aggiornamento su dove si trova questa funzione nella tua coda di sviluppo? Grazie
