请教各位大虾，假如我的数据是这样的：
data <- data.frame(y = rnorm(30), x = rnorm(30), grp = sample(c("east", "center", "west"), 30, replace = T))

在grp的每个水平下，分析x和y的相关性，应该怎么写程序呢？

cor.test(x,y,data=data[grp=="east",])
cor.test(x,y,data=data[grp=="center",])
cor.test(x,y,data=data[grp=="west",])

这样好像不行，因为cor.test命令好像不能指定data的名字

如果先生成subset，然后再做cor.test倒是可以，不过这样就要逐一的去生成subset，当分组变量有很多个，比如grp1,grp2,grp3.。。的时候，这种方法就不太可靠了

请问有好的方法吗？

比如f<-function(){
?????
}
进一步的，如果计算得出了结果，怎么分别保存对应的相关系数和p值