корреляция и р-значение с помощью функции cor.test - PullRequest
0 голосов
/ 06 сентября 2018

Этот запрос дает мне конечную ошибку наблюдения, пожалуйста, помогите мне решить эту проблему В настоящее время я пытаюсь использовать приведенный ниже код, мне нужно найти корреляцию и значение p между переменной (genolocmean) и переменной (genomean), присутствующей в кадре данных genolocmerge.

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

1 Ответ

0 голосов
/ 13 сентября 2018

Вы пробовали использовать функции карты?

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), добавив индекс, чтобы получить только то, что вам нужно.

Без рабочего набора данных я не могу проверить, будет ли этот код работать на самом деле. Если это не помогло, предоставьте образец набора данных, чтобы я мог проверить различные идеи.

...