<meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />
Qualcuno sa di uno script JMP per eseguire la procedura Hodges-Lehmann sui dati? Alcuni documenti standard europei richiedono che l'analisi di Hodges-Lehmann venga eseguita su tipi specifici di dati. I documenti mostrano come fare Hodges-Lehmann a mano. Posso ottenere l'output desiderato utilizzando il test Wilcoxon a 1 campione in Minitab sebbene richieda molte spiegazioni su quali parti dell'uscita Minitab corrispondono alle informazioni HL desiderate. Questi documenti raccomandano specificatamente StatXact o SAS di fare Hodges-Lehmann. Non ho SAS, ma ho JMP. Speravo di trovare uno script JMP per fare Hodges-Lehmann. Il test che sto eseguendo è un test unilaterale che la mediana di un set di dati è un valore specifico. E devo usare Hodges-Lehmann.
Grazie!
Bruce
Il metodo manuale che ho trovato proviene da un documento standard europeo. Qui è un link sul documento. Non sono sicuro delle informazioni sul copyright o pubblicheremo un'immagine della pagina appropriata o due. Sono stato in grado di trovare un algoritmo generalizzato Qui . Ho anche trovato alcuni passaggi SAS in Questo documento e Questo uno.
Credo di poter citare un paragrafo del documento europeo. Dopo che la mediana è stata trovata, la procedura seguente è descritta per un set di dati di dimensione n = 20.
"Le differenze medie a coppie che non superano la mediana ... sono calcolati Dalla tabella E.5 dei valori critici per le coppie di ranghi firmati di Wilcoxon testate la voce per n = 20 e un livello di significatività 0,025 unilaterale viene trovato il valore critico di 52. Quindi c = 52 + 1 = 53. Le differenze a coppie sono ordinate in ordine decrescente. Il valore 53 ° è 0,11 [NOTA: questo è dai dati di esempio nel documento]. quindi il limite di confidenza superiore unilaterale del 97,5% di Hodges-Lehmann per la differenza di Lg Rs tra RP e PP è 0,11, che è inferiore al margine di inferiorità concordato di 0,6. Pertanto, l'ipotesi di inferiorità della PP viene respinta e si può concludere che la preparazione del test PP non è inferiore a RP ".
Il documento europeo contiene una piccola tabella di valori critici di Wilcoxon, ma non sono stato in grado di trovare una tabella più completa da nessuna parte.
Grazie!
Bruce
bwwhite,
Hai un link per fare l'HL a mano? Inoltre, hai R installato sullo stesso computer di JMP?