计算关联并写入数据

Question

我试图循环一个数据帧，并计算x，y的每一个观测值的相关性，并存储这是一个新的框架。

sub_dat文件包含日期、硫酸盐、硝酸盐、ID的列。

for (s in sub_dat$ID) {
    corrData <- rbind(cor(sub_dat$sulfate, sub_dat$nitrate))

}

corrData

我得到的只是一个值，而不是矩阵。我已经删除了NA，调试时没有错误。有什么想法吗？

Answer 1

我不太确定您的数据是什么样子，但我假设如下：

             x           y group
1   0.08474419  0.16459807     1
2  -0.52390488 -0.52922622     1
3   0.12789339 -0.37239925     1
4  -0.49442433 -0.51763038     1
5   0.22440165 -0.73524373     1
6  -0.13219889  0.06979758     2
7   2.60532503  1.15088422     2
8   0.49880770  0.73113421     2
9   0.05612694  0.44701820     2
10  0.29456281  0.76804105     2

使用library(dplyr)，您可以做类似于df %>% group_by(group) %>% summarise(corr = cor(x, y))的事情

输出：

      group      corr
1     1          0.2093626
2     2          0.8392375