2015-07-21 33 views
1

我正在比較wilcoxon等級總和測試的列數據對,並且我得到了大多數比較的完全相同的pvalue。我想知道從數據來判斷我是否犯了錯誤,或者一切都沒問題。以下是一些比較。Wilcoxon等級總和結果:相同的pvalue一遍又一遍嗎?

這是我用

wtresult<-wilcox.test(datachunk[,i],datachunk[,(i+1)],paired=FALSE) 

在這裏,與上面所使用的數據結果進行比較。

X1  X2  X3      
339.53 354.11 435.56 425.34 434.64 436.08 
X1 X2 X3     
312.1 282.2 281.6 NA NA NA 

Wilcoxon rank sum test 

data: datachunk[, i] and datachunk[, (i + 1)] 
W = 18, p-value = 0.02381 
alternative hypothesis: true location shift is not equal to 0 

X1  X2  X3      
161.21 150.01 183.47 201.51 234.70 321.00 
X1 X2 X3     
501.0 520.1 500.7 NA NA NA 

Wilcoxon rank sum test 

data: datachunk[, i] and datachunk[, (i + 1)] 
W = 0, p-value = 0.02381 
alternative hypothesis: true location shift is not equal to 0 

X1  X2  X3      
247.79 159.64 192.00 262.86 403.33 336.21 
X1 X2 X3     
60.33 66.04 55.23 NA NA NA 

Wilcoxon rank sum test 

data: datachunk[, i] and datachunk[, (i + 1)] 
W = 18, p-value = 0.02381 
alternative hypothesis: true location shift is not equal to 0 

X1 X2 X3     
17.12 15.83 16.88 17.61 18.97 45.92 
X1 X2 X3     
321.8 329.7 334.4 NA NA NA 

回答

2

測試是有點「粗碎塊」的觀測的小號碼,以便如果有一個邊界的情況下(所有的第一個參數的值的比的第二個參數的值或反之亦然更大),你會得到相同的p -values和W統計數據都是0或其他數字(取決於n)。

如需更詳細的答案,我們需要查看您的數據,或者您需要同意查看其他所有我們可以看到的數據。

下面是一個代碼示例,顯示的行爲我談論

i <- 1 
datachunk <- mtcars[1:5,] 
wilcox.test(datachunk[,i],datachunk[,(i+1)],paired=FALSE) 

i <- 2 
wilcox.test(datachunk[,i],datachunk[,(i+1)],paired=FALSE) 
1

由於在每一個測試你的例子,從一個組的所有值比其他組的所有值越大,簽署無論實際數值如何,分配給每個差異的排名都是相同的。所以你的測試統計和p值都是一樣的。

所以在數學上它是有道理的,但我會考慮如果測試的應用對於如此小的樣本量有意義。

相關問題