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

Question

这个查询给了我有限的观察误差，请帮我解决这个问题。目前我正在尝试下面的代码，我需要找到genolocmerge数据框中存在的变量（genolocmean）和变量（genean）之间的相关性和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)以仅提取您需要的值。

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