我有一个数据帧如下:
x <- data.frame(Name=c("a", "b","c", "d", "e"),A=(1:5), B=(2:6), C=(7:11), D=c(1,1,1,1,1))
我想得到一个数据帧,包括a vs b的所有pearson系数,a vs c,a vs d,a vs e,b vs a,b vs c,b vs d,b vs e .... etc(不包括自我比较,即aa,bb等) .
换句话说,我想更改以下数据帧
A B C D
a 1 2 7 1
b 2 3 8 1
c 3 4 9 1
d 4 5 10 1
e 5 6 11 1
对于这样的结果数据帧
a b 0.989143436
a c 0.963282002
a d 0.93048421
a e 0.89585579
b c 0.9922798
b d 0.974216034
b e 0.951427811
c d 0.994675789
c e 0.982264673
d e 0.996357429
请善意地指导这样做的有效方法 .
EDIT
感谢Michaelv2的帮助 .
根据建议的代码,我发现结果如下:
X1 X2 value
1 A A 1
2 B A 1
3 C A 1
4 D A NA
5 A B 1
6 B B 1
7 C B 1
8 D B NA
9 A C 1
10 B C 1
11 C C 1
12 D C NA
13 A D NA
14 B D NA
15 C D NA
16 D D 1
错误消息是“警告消息:在cor(x [2:5],方法=”pearson“):标准偏差为零”
在我看来,我可能误用了代码,请你进一步指示如何进一步解决这个问题?谢谢 .
1 回答
您可以使用以下内容:
结果: