我有一个数据集描述了属于不同（1）种，（2）性别和（3）生殖阶段的个体的几个行为变量。 我的目标是创建一系列循环，这些循环允许我计算每个行为变量（例如＆＃34;隐藏＆＃34;下面）与另一个感兴趣变量之间的Spearman相关系数和p值（＆＃34;嵌套＆ ＃34;以下）每个品种每个阶段的性别。 我最终想要的是一个如下所示的数据框：
Species Sex Stage Spearman_corr_coef p-value behaviour
PO m virgin 0.78 0.01 hiding
PO f virgin 0.54 0.09 hiding
PO m parent 0.03 0.71 hiding
PO f parent 0.14 0.59 hiding
BW m virgin 0.71 0.02 hiding
BW f virgin 0.79 0.04 hiding
... ... ... ... ... ...
以下是如何创建类似于我的假数据集：
ID<-as.factor(rep(seq(from=1, to=40), each=2))
hiding <- rnorm(80)
attacking <- rnorm(80)
nesting <- rnorm(80)
sex <- as.factor(rep(c("m","f"),each=2, times=20))
species <- as.factor(rep(c("A", "B"), each=40))
stage<-as.factor(rep(c("virgin", "parent"), times=80))
df <- data.frame(ID, sex, species, stage, nesting, hiding, attacking)
我想关联＆＃34;嵌套＆＃34;的值。 （感兴趣的变量）到＆＃34;隐藏＆＃34;和＃34;攻击＆＃34;，每个阶段每个性别使用Spearman相关性测试。然后，我想提取相关系数和p值，并将它们全部集成到一个数据帧中（如上所示）。 提前谢谢！