作者andrew43 (討厭有好心推文後刪文者)
看板Statistics
標題Re: [程式] R 的wilcox.test及統計問題
時間Thu Jul 10 13:44:57 2014
※ 引述《innocent817 (Shawn純真)》之銘言:
: 標題: [程式] R 的wilcox.test及統計問題
: 時間: Wed Jul 9 21:59:53 2014
:
:
: [軟體程式類別]:
: R
:
: [程式問題]:
: wilcox.test
:
: [軟體熟悉度]:
: 低(1~3個月)
:
: [問題敘述]:
: 不好意思,小的半路出家,因為論文需求"可能"需要用到wilcoxon signed-rank test因此想請問版上的大大們幾個問題
:
: 1,統計問題
: 我現在有多組成對的排序資料,是根據同一筆資料(10組方案,3個準則)使用不同多準則方法得到的方案排序
: 如下所示:
: 1 2 3 4 5 6 7 8 9 10
: ================================
: A :7 1 4 10 2 6 3 8 5 9
: B1:7 1 4 10 2 6 3 8 5 9
:
: 1 2 3 4 5 6 7 8 9 10
: ================================
: A :7 1 4 10 2 6 3 8 5 9
: B2:7 1 4 10 2 5 3 8 6 9
:
: 1 2 3 4 5 6 7 8 9 10
: ================================
: A :7 1 4 10 2 6 3 8 5 9
: B3:3 1 9.5 9.5 4 7 8 5 2 6
:
: 以A,B3來看,就是第一筆資料在A方法排第7名,在B3方法排第3名
: 第二筆資料在A方法排第1名,在B3方法排第1名
: 第三筆資料在A方法排第4名,在B3方法排第9.5名(3,4共享第9名,(9+10/2=9.5))
:
: 想請問像以上的資料可以使用wilcoxon檢定嗎?
: 因為我在網路上看的範例多為使用不同方法得到的原始資料,而非排序過的數據
: 這樣是否可行?因為從片面了解到wilcoxon可以用來比較兩個方法出來的結果是否相同
:
:
: --
: ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 163.14.45.36
: ※ 文章網址: http://www.ptt.cc/bbs/Statistics/M.1404914395.A.0E5.html
: → andrew43:看不懂你的資料是什麼。 07/09 22:22
: → andrew43:是不是說明清楚資料怎麼產生的? 07/09 22:25
: → innocent817:先生成十組三維資料,i in 1:10 X[i]<-runif(3,0,1) 07/09 23:27
: → innocent817:然後使用一些多準則方法進行排序,上方的1:10代表方案 07/09 23:29
: → innocent817:A1可視為rank(sum(x[i])),B1則是rank(max(x[i])) 07/09 23:31
: → innocent817:簡單來說是要看方案在不同方法後的排序是否有差異 07/09 23:34
: ※ 編輯: innocent817 (163.14.45.36), 07/09/2014 23:40:46
: → innocent817:修改了一下前面文章的表達方式~ 07/09 23:41
: → andrew43:如果10筆資料彼此獨立,那可以採用配對設計。 07/09 23:59
: → andrew43:看起來你在R中沒有採用配對設計,注意一下。 07/09 23:59
: → andrew43:不過我想,如果是20個資料做rank會比10個資料做rank還 07/10 00:01
: → andrew43:要有利。我還是對你這種排序方法進行這個檢驗感到懷疑。 07/10 00:02
: → innocent817:不好意思,可以請問懷疑的點在那嗎?因為我對這個方法 07/10 12:05
: → innocent817:不熟悉,有任何不可行的考量都可以說出來讓我瞭解 07/10 12:05
我舉幾個反駁你的例子看看,
但請先簡單地了解wilcox.test的運算過程才行。
已知有 A B 二種方法,有 1--5 號五種資料。
1.
我假設 B 方法普遍比 A 方法優良,
但因為你排序的方法是方法內自己排序,
所以二個方法下的排序是完全相同的。
於是,wilcox.test 的結果告訴你二方法一致。
2.
我假設 B 方法和 A 方法幾乎等效,
只差在 A 方法的最佳資料是 B 方法的最差資料,其它關係不變。
例於下面資料
資料
方法 1 2 3 4 5
================================
A 1 2 3 4 5
B 5 1 2 3 4
於是,wilcox.test 的結果將傾向 B 方法小於 A 方法。
不知道你看出來你這樣的排序方法拿去做這種檢驗的不當之處了嗎?
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 122.117.37.172
※ 文章網址: http://www.ptt.cc/bbs/Statistics/M.1404971101.A.15A.html
※ 編輯: andrew43 (122.117.37.172), 07/10/2014 13:50:19
→ innocent817:1的部分我可以理解,可是2還是不太清楚,換句話來說 07/10 14:53
→ innocent817:是因為我方法內的排序將我的數值全部侷限在1:10之間 07/10 14:53
→ innocent817:所以其總和其實是一致的嗎?還是這樣的理解也不對? 07/10 14:55
→ andrew43:和總和沒關係,就只是你的資料不適合這個方法。 07/10 15:05
→ innocent817:是因為wilcox是用來判別A方法是否大於、小於或等於B方 07/10 15:16
→ innocent817:法,而不是用來判斷這兩個方法是否有顯著差異(因為一 07/10 15:17
→ innocent817:定有差異,指示是否顯著),而我的資料有某種規律與差 07/10 15:18
→ innocent817:異等限制,所以不適合wilcox的測試? 07/10 15:19
→ andrew43:應該是你說的沒錯。你可以看清楚wilcox.test的運算過程, 07/10 15:20
→ andrew43:應該可以有所體會。 07/10 15:20
推 innocent817:瞭解了,非常謝謝你的解說!太感謝拉~ 07/10 15:25