使用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) 。

沒有有效的數據集，我將無法測試此代碼是否真正起作用。 如果沒有，請提供樣本數據集，以便我可以測試不同的想法。