d中对于cor的分组

Question

我希望获得大量观察的按组（年）成对相关性。

当不使用for循环时，我得到了我想要的结果，即：

ddply(mydata, .(year), summarise, corr=cor(x, y, use="pairwise.complete.obs"))  

我想要的结果：

1   1 0.8366892
2   2 0.8929666
3   3 0.8364396
4   4 0.6201038
5   5 0.8914541

但是当我使用for循环遍历我的数据集的列时 喜欢：

for (i in 1:length(x))
ddply(mydata, .(year), summarise, corr=cor(x[[i]], y[[i]], use="pairwise.complete.obs"))

我明白了：

 grp     corr
1   1 0.835378
2   2 0.835378
3   3 0.835378
4   4 0.835378
5   5 0.835378

不同年份的平均相关性

我不理解ddply的工作方式吗？