Вы пробовали использовать функции карты?
genoloccorr <- genolocmerge %>% group_by(br_field_id) %>%
summarise(result=map2(genolocmean,genomean,~cor.test(.x,.y)))
Это должно вернуть фрейм данных со столбцом для br_field_id
и список-столбец с выводом cor.test()
. Если вам нужны только р-значение и корреляция, вы можете изменить ~cor.test(.x,.y)
, добавив индекс, чтобы получить только то, что вам нужно.
Без рабочего набора данных я не могу проверить, будет ли этот код работать на самом деле. Если это не помогло, предоставьте образец набора данных, чтобы я мог проверить различные идеи.