如何在r中使用模型摘要包时向相关性表添加p值?

pn9klfpd  于 2023-04-18  发布在  其他
关注(0)|答案(1)|浏览(126)

我能够创建一个包含13个变量的相关性表,但根据我的任务,我需要在此相关性表中包含p值。我使用了模型摘要包中的data summary_correlation,并能够创建一个发布就绪表。我在网上找不到任何有助于添加p值的东西。
我希望我能通过发布我的代码得到更多的帮助,但数据是研究敏感的。
我希望有一种方法可以将p值列添加到表的末尾。

ni65a41a

ni65a41a1#

没有内置的方法来添加p值,但datasummary_correlation()允许用户提交任意相关函数。这使得它相对容易破解。例如:

library(modelsummary)
library(correlation)

dat = mtcars[1:5, 1:5]

fun = function(x) {
    out = correlation(x)
    stars = c("*" = .2, "**" = .15, "***" = .05)
    p = modelsummary:::make_stars(out$p, stars)
    out$r = sprintf("%.2f%s", out$r, p)
    out = as.matrix(out)
    return(out)
}

datasummary_correlation(dat, method = fun)
每加仑圆柱体分散马力德拉特
每加仑1个-0.98***-.86- 点95**点51口径
圆柱体-0.98***1个点八九口径点九一口径-.59
分散-.86点八九1个0.90-.87
马力- 点95**点九一口径0.901个-.60
德拉特点51口径-.59-.87-.601个

相关问题