作者innocent817 (Shawn纯真)
看板Statistics
标题[程式] 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可以用来比较两个方法出来的结果是否相同
2.R的wilcox.test
使用wilcox.test把上述的资料丢进来後出现下列结果：
Wilcoxon signed rank test with continuity correction
data: as.numeric(researching_10[7, 11:20])
V = 27, p-value = 0.6344
alternative hypothesis: true location is not equal to 0
警告讯息：
1: In wilcox.test.default(as.numeric(researching_10[7, 11:20])) :
cannot compute exact p-value with ties
2: In wilcox.test.default(as.numeric(researching_10[7, 11:20])) :
cannot compute exact p-value with zeroes
想请问这样是两个数据不相等的意思？
下面的警告讯息是因为我的样本数太少吗？
还是像A1&B1的情况，相同的数据太多，所以差为0的笔数太多？
由於是半路出家，所以对一些基本的统计问题不太了解
有看过版上一些相关的文章，可是感觉跟我遇到的问题又不太一样
因此想请问版上大家对这个方法的了解，谢谢
[程式范例]:
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 163.14.45.36
※ 文章网址: http://webptt.com/cn.aspx?n=bbs/Statistics/M.1404914395.A.0E5.html
1F：→ andrew43:看不懂你的资料是什麽。 07/09 22:22
2F：→ andrew43:是不是说明清楚资料怎麽产生的? 07/09 22:25
3F：→ innocent817:先生成十组三维资料，i in 1:10 X[i]<-runif(3,0,1) 07/09 23:27
4F：→ innocent817:然後使用一些多准则方法进行排序，上方的1:10代表方案 07/09 23:29
5F：→ innocent817:A1可视为rank(sum(x[i]))，B1则是rank(max(x[i])) 07/09 23:31
6F：→ innocent817:简单来说是要看方案在不同方法後的排序是否有差异 07/09 23:34
※ 编辑: innocent817 (163.14.45.36), 07/09/2014 23:40:46
7F：→ innocent817:修改了一下前面文章的表达方式~ 07/09 23:41
8F：→ andrew43:如果10笔资料彼此独立，那可以采用配对设计。 07/09 23:59
9F：→ andrew43:看起来你在R中没有采用配对设计，注意一下。 07/09 23:59
10F：→ andrew43:不过我想，如果是20个资料做rank会比10个资料做rank还 07/10 00:01
11F：→ andrew43:要有利。我还是对你这种排序方法进行这个检验感到怀疑。 07/10 00:02
12F：→ innocent817:不好意思，可以请问怀疑的点在那吗？因为我对这个方法 07/10 12:05
13F：→ innocent817:不熟悉，有任何不可行的考量都可以说出来让我了解 07/10 12:05
14F：→ innocent817:另外，成对比较是指加上paired=TRUE吧！这个我会加上 07/10 12:07
15F：→ innocent817:去的;样本其实有分10组跟100组，不过还是谢谢大大的 07/10 12:07
16F：→ innocent817:提醒 07/10 12:07