Question

用R语言计算多水平之间的相关性

如图所示，数据是乱写的，想要在R里编写一个程序，分别计算A1和B1、B2、B3到B12两两之间的相关性，再计算A2和B1、B2、B3到B12两两之间的相关性，然后是A3、A4和B1到B12两两之间的相关性，程序应该怎么编？本人是非计算机专业，对于编程只大概了解点皮毛，望高人指教。

Answer 1

先把整个数据读进来，比如赋值给data。data第一列是A1,A2什么的。

A<-as.matrix(data[grep("A",data[,1]),-1])
B<-as.matrix(data[grep("B",data[,1]),-1])
correlation.matrix<-matrix(nrow=nrow(A),ncol=nrow(B))
for(i in 1:nrow(A))
    for(j in 1:nrow(B))
        correlation.matrix[i,j]<-cor(A[i,],B[j,])

结果就在correlation.matrix里了

追问

遇到这种情况，要怎么办？

追答

不要这样啊，老大。。。。这样吧，你先把数据按照你图里的样子写到Excel里，保存成csv文件，记下完整路径（比如是D:\mydata\data.csv），打开R，输入这些：

data<-read.csv("D:/mydata/data.csv",header=FALSE)
#引号里的是路径，自己改，注意把\变成/

A<-as.matrix(data[grep("A",data[,1]),-1])
B<-as.matrix(data[grep("B",data[,1]),-1])
correlation.matrix<-matrix(nrow=nrow(A),ncol=nrow(B))
for(i in 1:nrow(A))
    for(j in 1:nrow(B))
        correlation.matrix[i,j]<-cor(A[i,],B[j,])

追问

下标出界。。

不好意思啊，再厚着脸皮问一下怎么办

追答

你第一步怎么做的啊？为什么会在输入A<-as.matrix(data[grep("A",data[,1]),-1])之前打印出来那么奇怪的东西？

再有问题你把做了什么一步一步说清楚好吗？

追问

不好意思，麻烦帮看一下

追答

好奇怪。。。不可能出界吧。。。私信里说吧。

追问

追答

看一下Excel里的数据

追问