R语言 在相关矩阵/配对图文本中添加p值和/或95% CI

tvz2xvvm  于 2023-01-18  发布在  其他
关注(0)|答案(1)|浏览(175)

是否可以通过在此代码中添加一些内容,在相关性矩阵图的文本中,在提及点估计值的位置旁边添加精确的p值和/或95% CI?

我也欢迎使用其他软件包的解决方案。

ggpairs(baseline2, title="correlogram ")

df

baseline2 <- structure(list(delta_bp = c(-21.5, 7, -17.5, -10.5, -21, -7.5, 
4.5, 3, -9, 9, -22.5, -9.5), delta_bp_05 = c(-21.5, 0, 3, -13.5, 
-13, -4, -16.5, -8, 5, -5, -12, 0.5), delta_bp_10 = c(-26.5, 
1, -6, -10.5, -9, -3, -20.5, -10, 1, -6, -22, -0.5)), row.names = c(NA, 
-12L), class = c("tbl_df", "tbl", "data.frame"), na.action = structure(c(`17` = 17L, 
`27` = 27L, `32` = 32L), class = "omit"))
zbq4xfa0

zbq4xfa01#

您可以通过upper"statistic"选项使用ggally_statistic函数,并提供自己的text_fn。您为text_fn给予的函数接受数据列,并应返回文本字符串以显示。这里是添加到相关性的p值和置信区间。这可以根据您的喜好进行自定义。

tf = function(x,y) {
  ct <- cor.test(x, y)
  sprintf("Corr: %0.2f, p: %0.3f\n95%% CI: (%0.2f to %0.2f)", 
          ct$estimate, ct$p.value, ct$conf.int[1], ct$conf.int[2])
}

ggpairs(baseline2, title="correlogram ", 
        upper=list(continuous=wrap("statistic",text_fn=tf,title=NULL, sep=NULL)))

相关问题