Saya ingin menentukan apakah ada perbedaan dalam nilai-p rata-rata antara dua kelompok. Untuk melakukan ini, saya melakukan tes peringkat-jumlah Wilcoxon (data tidak terdistribusi normal). Sejauh ini bagus. Akhirnya, saya ingin menghitung ukuran efek yang sesuai. Sayangnya, R tidak menyediakan ini. Itu juga tidak memberikan nilai az dimana ukuran efek dapat dengan mudah dihitung dengan menggunakan: effect size = z / sqrt (N)
berikut ini beberapa contoh kode R:
a=rep(0:1,each=20) #grouping variable
b=c(rnorm(20, .03,.01), rnorm(20, .02, .009)) #vector of p-values
d=cbind(a,b)
test = wilcox.test(b ~ a, data = d) #perform Wilcoxon rank-sum test
test
Apakah ada yang tahu cara mendapatkan ukuran efek?
Jawaban:
Estimator yang sesuai dengan uji Wilcoxon adalah estimator Hodges-Lehmann; dikembalikan dengan
wilcox.test
menggunakanconf.int=TRUE
opsi, di bawah "perbedaan lokasi".Sebagai contoh Anda:
Untuk lebih lanjut tentang Wilcoxon dan asumsi di baliknya, dan apa yang sebenarnya diuji, dan penduga nonparametrik lainnya, dokumen ini (mungkin) membantu: www.stat.umn.edu/geyer/old03/5102/notes/rank.pdf
sumber
wilcox.test(b~a,data=d, conf.int=TRUE)$estimate / sqrt(20)
benar?Dapatkan z untuk formula Anda
dan hitung ukuran efek dengan formula Anda, atur N ke 40
sumber
coin::wilcoxsign_test
fungsi R. Juga, apakah Anda mengacu pada formula OP,