请问如何获得相关系数矩阵的行列名？

2013-6-10 12:24:06

我也在学习R，试着把输出的矩阵转换成一个数据框，然后用attach(),names()函数调用。
x1 <- 1:10
x2 <- 2:11
x3 <- 10:1
df<- data.frame(x1,x2,x3)
data.cor <- cor(df)
attach(df)
p <- 0.5
for(i in 1:3){
for(j in 1:3)
{
   if(data.cor[i,j] >= p)
   {
   points <- c(names(df),names(df[j]))
   print(points)
   }
}
}

希望能有所帮助。

扫码加我拉你入群

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

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

2013-6-10 14:28:44

Ben910128 发表于 2013-6-10 12:24
我也在学习R，试着把输出的矩阵转换成一个数据框，然后用attach(),names()函数调用。
x1

多谢老兄回答了。但是我试了，不行，直接提示rcor.test的结果不能转化为data frame.

扫码加我拉你入群

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

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

2013-6-11 08:08:12

如果输出结果是矩阵，直接用colnames()就好

扫码加我拉你入群

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

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

2013-6-11 10:10:19

hugebear 发表于 2013-6-11 08:08
如果输出结果是矩阵，直接用colnames()就好

输出结果确实是矩阵，我用class()看了的，返回值是matrix。我是4个变量的相关分析，得到的是一个4*4的矩阵，但是最上面的一行和最左边的一列都是变量名。我现在想首先从矩阵中找到符合条件的相关系数（比如>0.5的），然后获取对应的变量名（就是说满足条件的值是哪两个变量的相关关系），然后以一定的格式输出。

扫码加我拉你入群

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

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

2013-6-11 13:33:09

把数据和你的程序运行结果贴出来胜于千言万语。。。。

扫码加我拉你入群

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

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

点击查看更多内容…

2013-6-11 19:50:54

danbaidong 发表于 2013-6-10 14:28
多谢老兄回答了。但是我试了，不行，直接提示rcor.test的结果不能转化为data frame.

那直接调用最初的数据框呢，是不是可以试试？

扫码加我拉你入群

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

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

2013-6-12 14:16:42

hugebear 发表于 2013-6-11 13:33
把数据和你的程序运行结果贴出来胜于千言万语。。。。

#读入数据#请注意把附件数据的扩展名改为csv，因为论坛不让上传csv格式的文件，故楼主改成了xls
testdata <- read.csv("test.csv")
testdata
#ltm包及相关检验
library(ltm)
#临时数据框
testmat <- cbind(height = testdata$height, shoe = testdata$shoe, hand = testdata$hand)
#相关检验
rcor.test(testmat)
以上代码执行后会得到如下结果：

   height shoe hand
height  *****  0.867  0.574
shoe <0.001  *****  0.516
hand 0.001  0.004  *****

upper diagonal part contains correlation coefficient estimates
lower diagonal part contains corresponding p-values

我想要的输出形式是：
height和shoe显著相关，相关系数是0.867；
height和hand显著相关，相关系数是0.574；
shoe和hand显著相关，相关系数是0.516。
因为现在变量少还好，如果变量很多，即便是相关分析得到的表格装得下，看起来也不方便，如果可以用文字的形式把显著相关的变量及其系数列出来，似乎更有利于信息的传达，尤其是给非专业人士看的时候。故有此问。希望大家可以指点在下一二。

附件列表

test.xls

大小:525 Bytes

马上下载

data for the post

扫码加我拉你入群

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

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

2013-6-12 14:18:17

Ben910128 发表于 2013-6-11 19:50
那直接调用最初的数据框呢，是不是可以试试？

我对这个问题作了补充，欢迎继续讨论，谢谢！

扫码加我拉你入群

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

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

2013-6-12 17:37:40

danbaidong 发表于 2013-6-12 14:18
我对这个问题作了补充，欢迎继续讨论，谢谢！

我明白楼主要做一件什么事情了。我觉得楼主是不是在处理基于相关系数的变量分类问题是产生的这个想法？不过我觉得，当数据比较多时，用字符串的方式输出结果没有以列表的方式输出简洁。不过按照楼主的要求给出代码，最后我是直接做了一个函数，更具有普遍性。应该这样就可以吧。

复制代码

扫码加我拉你入群

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

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

2013-6-12 22:55:08

I am sorry I am at school so I can't use Chinese. Actually, you don't need to convert the matrix into a
dataframe, since the matrix also has "dimnames" attribute.

复制代码

扫码加我拉你入群

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

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

2013-6-12 23:35:51

Ben910128 发表于 2013-6-12 17:37
我明白楼主要做一件什么事情了。我觉得楼主是不是在处理基于相关系数的变量分类问题是产生的这个想法？不 ...

多谢老兄回复！我试着运行了一下，确实可以出来我想要的结果。不过还有两点不明，提出来供商讨：
首先是第09行，为什么要把输入的矩阵转化为数据框呢？
其次是，老兄这里使用的R自带的cor()函数来实现相关系数的计算，但是这个函数只能返回变量的相关系数，并不能对相关性是否具有统计意义进行假设检验（也就是不能返回p-value）。从网上查到ltm包中的rcor.test函数可以同时计算相关系数和p-value，分别在输出的矩阵的上三角和下三角表示，所以我最初的想法就是首先判断p-value是否小于0.1，如果是的话，那就认为是显著相关的，然后再按照我上面提到的格式输出。我在实现的过程中，就发现rcor.test的返回值的类型是“rcor.test”(用class()函数判断)，而且也没有办法转化成矩阵或者数据框等类型，不知道你是否遇到过类似的情况？谢谢

扫码加我拉你入群

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

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

2013-6-12 23:49:12

danbaidong 发表于 2013-6-12 23:35
多谢老兄回复！我试着运行了一下，确实可以出来我想要的结果。不过还有两点不明，提出来供商讨：
首先是 ...

1-关于转换为data.frame的问题。因为我发现转换为数据框之后，names()函数才有用，并且能解决这个问题，所以我就转换了。当然我本人也觉得这不是一个好方法，但是个方法。哈哈~
2-p-value直接用cor.test()函数就可以得到了，我没有用过rcor.test()所以我也不是很清楚你说的那个情况。
我是个新手，接触R才一个星期，所以很多问题我也是很模糊的。刚刚看到楼上的哥们给出他的方法，你要不试试？

扫码加我拉你入群

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

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

2013-6-13 00:51:37

hugebear 发表于 2013-6-12 22:55
I am sorry I am at school so I can't use Chinese. Actually, you don't need to convert the matrix int ...

多谢仁兄！
不过ltm包中的rcor.test函数的返回值类型是rcor.test类型，而非矩阵或者数据框，也就是说根本原因还是在rcor.test函数本身。

扫码加我拉你入群

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

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

2013-6-13 02:16:58

danbaidong 发表于 2013-6-13 00:51
多谢仁兄！
不过ltm包中的rcor.test函数的返回值类型是rcor.test类型，而非矩阵或者数据框，也就是说根本 ...

自己转换一下就好，思想已经都写给你了。要学会灵活变通。

扫码加我拉你入群

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

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

2013-6-13 07:11:53

好吧，对于你这个具体问题，就再给一个程序吧。
注意虽然R的class千变万化，但本质上都是list, 要学会用unclass来查看庙里面究竟住了什么和尚。

复制代码

扫码加我拉你入群

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

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