大家好。我希望使用 dplyr和/或broom包 从相关矩阵中获取p值,同时测试多个变量。我知道其他方法,但对我来说,dplyr似乎更简单、更直观。此外,dplyr需要将每个变量相关起来以获得特定的p值,这使得该过程更加容易和快速。
我已经查看了其他链接,但它们对这个问题没有用 (示例1,示例2,示例3)。 当我使用这段代码时,会报告相关系数。然而,P值不会报告。
agreg_base_tipo_a %>%
dplyr::select(S2.RT, BIS_total, IDATE, BAI, ASRS_total) %>%
do(as.data.frame(cor(., method="spearman", use="pairwise.complete.obs")))
请查看以下可复制的代码并将其转换为Word文档:
set.seed(1164)
library(tidyverse)
ds <- data.frame(id=(1) ,a=rnorm(10,2,1), b=rnorm(10,3,2), c=rnorm(5,1,05))
ds %>%
select(a,b,c) %>%
do(as.data.frame(cor(., method="spearman", use="pairwise.complete.obs")))
rcorr
很不错的选项。 - akrun