全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
11561 3
2013-06-13
用R软件对一个矩阵进行friedman检验,矩阵是一个12*5的矩阵,12个处理,5个区组,显示不出结果

> d<-read.delim("F:data/book6.txt")
> d
      X2011     X2010     X2009      X2008    X2007
1  25.33432  15.60497  21.97791  18.191992 66.50104
2  33.85897 236.29722  10.16129  18.320611 55.95238
3  43.67464 258.29805  15.93794  25.709220 45.36082
4  34.96316 172.00949  76.84564 -12.780488 70.26578
5  40.95074 197.58707  20.03167  28.223350 41.93084
6  30.58526  17.53873  24.01747  39.068826 75.80071
7  53.23715 243.37123  18.27731  27.272727 57.80591
8  17.48572 446.22836  63.41463   8.753316 81.25000
9  87.75134 334.76947  17.59425  25.733634 30.29412
10 27.61323 902.49183  39.14141  36.082474 43.34975
11 19.93419 333.28439  50.70028  39.453125 89.62963
12 30.93760 222.15606 -34.10014  57.569296 58.44595
> friedman.test(d)
Error in friedman.test.default(d) :
  argument "groups" is missing, with no default
>

但是要是把数据矩阵转置以后就能出结果,只是自由度,结果都不是期望的结果,求高所点拨。。。
r初学。。。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2013-6-15 07:32:23
你的数据格式不对,你读入的是data frame格式,friedman.test要求的数据是向量或矩阵。
把你的数据转换成矩阵就OK了
e <- as.matrix(d)
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2015-4-28 21:28:19
那么请问老师,如果这组数据做出来P<0.05,我如何再进行下一步的两两检验呢?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2015-4-28 23:03:17
你会对12个样本的5个区组做friedman测试,说明每个样本都进行了重复测试,而是是非参数检验,所以做pairwise test的时候只需要取两行数据,做pairwise.wilcox.test,而且必须加上paired=TRUE,因为是重复测量
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群