Bạn có thể xem xét một bài kiểm tra hoán vị.
median.test <- function(x,y, NREPS=1e4) {
z <- c(x,y)
i <- rep.int(0:1, c(length(x), length(y)))
v <- diff(tapply(z,i,median))
v.rep <- replicate(NREPS, {
diff(tapply(z,sample(i),median))
})
v.rep <- c(v, v.rep)
pmin(mean(v < v.rep), mean(v>v.rep))*2
}
set.seed(123)
n1 <- 100
n2 <- 200
## the two samples
x <- rnorm(n1, mean=1)
y <- rexp(n2, rate=1)
median.test(x,y)
Cung cấp giá trị p 2 mặt là 0,1112, đây là một minh chứng cho việc kiểm tra trung bình không hiệu quả như thế nào khi chúng ta không kháng cáo với bất kỳ xu hướng phân phối nào.
Nếu chúng ta sử dụng MLE, 95% CI cho trung vị cho bình thường chỉ có thể được lấy từ giá trị trung bình vì giá trị trung bình là trung vị trong phân phối bình thường, do đó, từ 1,00 đến 1,18. CI 95% cho trung vị cho số mũ có thể được đóng khung là , theo phương pháp delta là 0,63 đến 0,80. Do đó, xét nghiệm Wald có ý nghĩa thống kê ở mức 0,05 nhưng xét nghiệm trung bình thì không.log(2)/X¯