Supponendo che entrambe le distribuzioni di Laplace abbiano la stessa varianza,
a) il test del rapporto di verosimiglianza implicherebbe una statistica di test come:
L=∏ni=112τ^exp(−|xi−μ^|τ^)∏n1i=112τ^1exp(−|xi−μ^1|τ^1)⋅∏ni=n1+112τ^2exp(−|xi−μ^2|τ^2)
Prendere i registri, cancellare / semplificare e moltiplicare per .−2
−2l=2(nlog(τ^)−n1log(τ^1)−n2log(τ^2)) (dove )l=log(L)
dove , la deviazione assoluta media dalla mediana nel campione combinato e , la deviazione assoluta media dalla mediana nel campione . τ i=miiτ^=mτ^i=mii
Secondo Wilks' teorema questo è asintoticamente distribuito come sotto l'ipotesi nulla, quindi per un test del 5% che ci si rifiuta se che ha superato . 3,84χ213.84
Gli esperimenti di simulazione suggeriscono che il test è anticonservativo su campioni di piccole dimensioni (la probabilità di rifiuto è leggermente superiore a quella nominale), ma di circa n = 100, sembra essere almeno ragionevole (si arriva all'ordine del 5,3% - 5,4% tasso di rifiuto sotto il valore nullo per un test nominale del 5%, ad esempio; per sembra essere più vicino al 5,25%).n1,n2>300
b) Ci aspetteremmo anche che sarebbe una buona statistica di test (dove rappresenta il mediana campione ); se non ho commesso un errore, in grandi campioni come il tuo sarebbe approssimativamente distribuito normalmente sotto il valore null, con media 0 e varianza 1, dove potrebbe essere basato sul quadrato del deviazione assoluta media dalla media nel campione combinato, , anche se mi aspetto che in pratica tenderebbe a funzionare meglio basandosi su una media ponderata del campione dei due campioni 's
. ~μv=2τ2(1μ~1−μ~2v√μ~ τ 2m2m 2 i †v=2τ^2(1n1+1n2)τ^2m2m2i†
† (Modifica: la simulazione suggerisce che l'approssimazione normale va bene ma il calcolo della varianza non è corretto sopra; posso vedere qual è il problema ora ma devo ancora risolverlo. La versione di permutazione di questo test (vedi elemento (c)) dovrebbe ancora andare bene).
c) Un'altra alternativa sarebbe quella di eseguire un test di permutazione basato su una delle statistiche di cui sopra. (Una delle risposte qui fornisce uno schema di come implementare il test di permutazione per una differenza nelle mediane.)
d) Puoi sempre fare un test Wilcoxon / Mann-Whitney; sarà notevolmente più efficiente rispetto al tentativo di utilizzare un test t sul Laplace.
e) Meglio di (d) per i dati di Laplace sarebbe il test mediano di Mood; mentre spesso raccomandato contro nei libri, quando si tratta di dati Laplace mostrerà un buon potere. Mi aspetto che avrebbe un potere simile alla versione di permutazione del test asintotico della differenza nelle mediane (uno dei test menzionati in (c)).
La domanda qui fornisce un'implementazione R che utilizza un test di Fisher, ma quel codice può essere adattato per utilizzare invece un test chi-quadro (che suggerirei anche in campioni moderati); in alternativa c'è un codice di esempio per questo (non come una funzione) qui .
Il test mediano è discusso qui su Wikipedia , anche se non in modo molto approfondito (la traduzione tedesca collegata ha un po 'più di informazioni). Alcuni libri su non parametrici ne parlano.