W deki wilcox.test () 'in W istatistiki çıktısı U istatistiğiyle aynı mı?


22

Son zamanlarda Mann-Whitney U testi hakkında okudum. Bu testi R'de gerçekleştirmek için aslında bir Wilcoxon testi uygulamanız gerektiği ortaya çıktı!

Sorum şu: wilcox.testR'nin W istatistiği U istatistiğiyle aynı mı?

Yanıtlar:


25

Wilcoxon, genellikle testin orijinal mucidi olarak kabul edilir *, Mann ve Whitney'in yaklaşımı büyük bir adım attı ve istatistiklerin tablolaştırıldığı durumları genişletti. Tercihim, teste her iki katkıyı tanımak için Wilcoxon-Mann-Whitney olarak atıfta bulunmaktır (Mann-Whitney-Wilcoxon da görülür; Ben de umursamıyorum).

* Bununla birlikte, gerçek tablo biraz daha bulutludur, diğer bazı yazarlar da bu kez veya daha önce aynı veya benzer istatistiklerle veya bazı durumlarda testle yakından bağlantılı katkılarda bulunur. En azından kredinin bir kısmı başka bir yere gitmeli.

Wilcoxon testi ve Mann-Whitney U testi, aynı durumlarda her zaman aynı davaları reddetmeleri bakımından eşdeğerdir (ve yardım ettiklerini belirtir); Test istatistiklerinin çoğu sadece bir kayma (ve bazı durumlarda, sadece bir işaret değişikliği) ile değişecektir.

Wilcoxon testi literatürde birden fazla şekilde tanımlanmıştır (ve bu belirsizlik test istatistiğinin orijinal tablolarına dayanmaktadır, bir dakikadan daha fazladır), bu nedenle bir Wilcoxon testinin tartışılmasında dikkat edilmesi gerekir.

En yaygın iki tanım biçimi bu iki gönderide tartışılmaktadır:

Wilcoxon rütbe toplamı testi R

Wilcoxon rütbe toplamı testi için test istatistiklerini hesaplamanın farklı yolları

Özellikle R'de ne olduğunu ele almak için:

wilcox.testR'de kullanılan istatistik yardımda ( ?wilcox.test) tanımlanmakta ve Mann-Whitney U istatistiği ile ilişkinin sorusu burada açıklanmaktadır:

Wilcoxon sıra toplamının ve Mann-Whitney testlerinin tanımları ile ilgili literatür oybirliği değildir.

En yaygın iki tanım, çıkarılan asgari değere sahip ilk örneğin saflarının toplamına tekabül eder: R çıkarır ve S-PLUS, a için m (m + 1) / 2 ile daha büyük bir değer vererek, ilk beden büyüklüğü m. (Wilcoxon'un orijinal makalesinde rütbelerin ayarlanmamış toplamı kullanılmış, ancak sonraki tablolar en aza indirilmiş gibi görünüyor.)

R değeri, tüm çift sayısı olarak hesaplanabilir (x[i], y[j])olan y[j]daha büyük değildir x[i], Mann-Whitney testi en yaygın tanımı.

Bu son cümle, sorunuzun bu yönüne tamamen cevap veriyor - R'nin ortaya koyduğu W sürümü * aynı zamanda U'nun değeridir.

* Örnek 1'deki sıraların toplamı, alabileceği en düşük değeri eksi (yani eksi ).n1(n1+1)2


7

Hem Wilcoxon rütbe toplamı testi hem de Mann-Whitney testi, bağımsız t-testinin parametrik olmayan eşdeğerleridir . Bazı durumlarda R'nin verdiği W sürümü aynı zamanda U'nun da değeridir. Ancak her durumda değil.

Kullandığınızda: wilcox.test(df$var1 ~ df$var2, paired=FALSE)verilen W, U ile aynıdır. Dolayısıyla, Mann-Whitney U istatistiği olarak rapor edebilirsiniz.

Ancak kullandığınızda: wilcox.test(df$var1 ~ df$var2, paired=TRUE)Aslında bir Wilcoxon imzalı rütbe testi gerçekleştiriyorsunuz. Wilcoxon işaretli rütbe testi, bağımlı t-testinin eşdeğeridir .

Kaynak: "Field R kullanarak istatistik bulma" Andy Field (2013)


1
Eğer argümanı kullanırsan paired=TRUE, Wilcoxon-Mann-Whitney değil, imzalı rütbenin belirtilmesi iyi bir fikirdir .
Glen_b -Reinstate Monica

1

Ancak, kodun: wilcox.test(df$var1 ~ df$var2, paired=FALSE)('~' kullanarak) olduğuna dikkat edin.

aşağıdakilerden farklı bir W istatistiği üretecektir: wilcox.test(df$var1, df$var2, paired=FALSE)(',' kullanarak)


Siteye Hoşgeldiniz, @Tilen. Bu, OP'nin sorusuna bir cevap, OP'den veya cevap verenlerden birinden açıklama isteyen bir yorum mu yoksa kendi yeni bir sorunuzu mu amaçladı? Lütfen asıl soruya cevap vermek için sadece "Cevabınız" alanını kullanın. İtibarınız> 50 olduğunda her yerde yorum yapabileceksiniz. Yeni bir sorunuz varsa ASK QUESTION , sayfanın üstündeki griye tıklayın ve orada sorun, o zaman size doğru şekilde yardımcı olabiliriz. Burada yeniyken , yeni kullanıcılar için bilgi sahibi olan turumuza katılmak isteyebilirsiniz .
gung - Reinstate Monica

Çok teşekkürler. Aslında, sanırım hem cevap hem de soru birazdı. Tavsiyeler için teşekkürler. Bunun bir parçası olmayı dört gözle bekliyorum.
Tilen,

'Cevaplar' sadece sorunun doğrudan cevapları içindir. Bunu daha kesin bir cevap haline getirmek için düzenleyebilir misiniz?
gung - Reinstate Monica

Bunun için üzgünüm. Şimdi düzenlemiştim. Şimdi daha iyi görünüyor mu?
Tilen

1
Özre gerek yok. Sanırım cevabınız buraya koyduğunuz şeyle örtük ama bence insanlar satırları okumak zorunda kalacak. Açık ve her şeyi bir yayda bağlayabilir misin? Ancak bunun doğru olduğundan emin değilim. Formül sürümünü kullanmak için bir değerler vektörüne ve bir gruplandırma göstergesine ihtiyacınız vardır. Örneğin, wilcox.test(values~ind, with(df, stack(var1=var1, var2=var2)), paired=FALSE). Bunu yaptığımda, Wher iki şekilde de aynı olur .
gung - Reinstate Monica
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.