使用cor.test函数进行相关和p值

Question

该查询给了我有限的观察误差，请帮助我解决这个问题 目前，我正在尝试下面的代码，我需要找到genolocmerge数据帧中存在的变量（genolocmean）和变量（genomean）之间的相关性和p值。

genoloccorr <- genolocmerge %>% group_by(br_field_id) %>% 
           do(tidy(cor.test(.$genolocmean,.$genomean)))

Answer 1

您是否尝试过使用地图功能？

genoloccorr <- genolocmerge %>% group_by(br_field_id) %>% summarise(result=map2(genolocmean,genomean,~cor.test(.x,.y)))

这应该返回一个数据帧，其中的列为br_field_id，列表列的输出为cor.test()。如果只需要p值和相关性，则可以通过附加索引以仅提取所需内容来修改~cor.test(.x,.y)。

没有有效的数据集，我无法测试此代码是否真正起作用。如果没有，请提供样本数据集，以便我可以测试不同的想法。