получение p-значений для парной корреляции (dplyr) - PullRequest
Новая
       83

получение p-значений для парной корреляции (dplyr)

1 голос
/ 09 апреля 2020

Я использую код ниже, чтобы получить корреляции между моей зависимой переменной и ответом на вопросник (для разных уровней разных условий). 

BREAK %>%
    group_by(condition, valence) %>%
    summarize(COR=cor(rt, positive_focused_cognitiveER)) %>%
    ungroup()

Это дает мне корреляции и их направления (+/- ). Однако я хотел бы знать, значимы ли эти корреляции. Есть ли способ просто добавить строку в код, который я уже должен получить p-значения? Или другой простой код? (Мне не нужны причудливые вещи, только цифры)

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

Заранее спасибо за любые советы! :)

1 Ответ

1 голос
/ 09 апреля 2020

Вы можете вычислить p-значения с помощью stats::cor.test: 

BREAK %>%
        group_by(condition, valence) %>%
        summarize(COR = stats::cor.test(rt, positive_focused_cognitiveER)$estimate,
                  pval = stats::cor.test(rt, positive_focused_cognitiveER)$p.value
                  ) %>%
        ungroup()
...