Supondo que as duas distribuições da Laplace tenham a mesma variação,
a) o teste da razão de verossimilhança envolveria uma estatística de teste como:
L=∏ni=112τ^exp(−|xi−μ^|τ^)∏n1i=112τ^1exp(−|xi−μ^1|τ^1)⋅∏ni=n1+112τ^2exp(−|xi−μ^2|τ^2)
Tomando logs, cancelando / simplificando e multiplicando por .−2
−2l=2(nlog(τ^)−n1log(τ^1)−n2log(τ^2)) (onde )l=log(L)
onde , o desvio absoluto médio da mediana na amostra combinada e , o desvio absoluto médio da mediana na amostra .τ^=mτ^i=mii
De acordo com o teorema de Wilks, isso é distribuído assintoticamente como sob o valor nulo; portanto, para um teste de 5%, você rejeitaria se ultrapassasse ,.χ213.84
Experimentos de simulação sugerem que o teste é anticonservador em amostras pequenas (a probabilidade de rejeição é um pouco maior que a nominal), mas em cerca de n = 100, parece ser pelo menos razoável (você fica na ordem de 5,3% - 5,4% taxa de rejeição abaixo do nulo para um teste nominal de 5%, por exemplo; para , parece estar mais próximo de 5,25%).n1,n2>300
b) Também esperamos que seja uma boa estatística de teste (onde representa o mediana da amostra ); se eu não cometer um erro lá, em amostras grandes como a sua, ela será normalmente normalmente distribuída sob o nulo, com média 0 e variação 1, em que pode ser baseado no quadrado da desvio médio absoluto a partir da média da amostra combinada, , apesar de esperar que, na prática, tendem a funcionar melhor baseando-se em uma amostra média ponderada das duas a amostra 's
.μ~1−μ~2v√μ~v=2τ^2(1n1+1n2)τ^2m2m2i†
† (Edit: simulação sugere que a aproximação normal é boa, mas o cálculo da variação não está correto acima; posso ver qual é o problema agora, mas ainda preciso corrigi-lo. A versão de permutação deste teste (consulte o item (c)) ainda deve estar bem).
c) Outra alternativa seria realizar um teste de permutação com base em uma das estatísticas acima. (Uma das respostas aqui fornece um esboço de como implementar o teste de permutação para uma diferença nas medianas.)
d) Você sempre pode fazer um teste de Wilcoxon / Mann-Whitney; será consideravelmente mais eficiente do que tentar usar um teste t no Laplace.
e) Melhor que (d) para os dados de Laplace seria o teste mediano de Mood; embora seja frequentemente recomendado nos livros, ao lidar com os dados da Laplace, eles mostrarão um bom poder. Espero que ele tenha poder semelhante à versão de permutação do teste assintótico da diferença de medianas (um dos testes mencionados em (c)).
A pergunta aqui fornece uma implementação R que usa um teste de Fisher, mas esse código pode ser adaptado para usar um teste do qui-quadrado (o que eu sugeriria mesmo em amostras moderadas); Como alternativa, há um código de exemplo para ele (não como uma função) aqui .
O teste mediano é discutido na Wikipedia aqui , embora não com muita profundidade (a tradução alemã vinculada tem um pouco mais de informação). Alguns livros sobre não paramétricos discutem isso.