如何将Pearson相关矩阵打印到r中的Word

Question

我对R很陌生，我想在word文档中绘制一个奇特的相关矩阵。

这里是我的数据的一个小例子：

provaSO1 <- structure(list(TotalDebt = c(0.637, 0.517, 0.581, 0.785, 0.687,0.703, 0.474), 
Long_ok = c(0.157, 0.121, 0.051, 0.29, 0.153,0.301, 0.102), 
Short_ok = c(0.48, 0.396, 0.531, 0.495, 0.535,0.402, 0.372), 
Size = c(6.184, 6.184, 6.663, 7.302, 6.714, 6.949,7.627), 
ln_Age = c(3.638, 3.664, 3.689, 3.714, 3.738, 3.761, 3.784), 
liquidity = c(0.99, 0.988, 0.995, 0.965, 0.949, 0.949,0.53),
Asset_Tangibility = c(0.005, 0.006, 0.003, 0.023, 0.033, 0.03, 0.457), 
profitability = c(-0.058, -0.202, -0.106, 0.032, 0.03, 0.013, 0.042)), 
row.names = c(NA, 7L), class = "data.frame")

首先我想计算相关矩阵：

c_matrix = cor(provaSO1, method = c("pearson"))

我使用corstars function来增加显著性水平，并且只获得下三角：

corstarsl <- function(x){
require(Hmisc)
x <- as.matrix(x)
R <- rcorr(x)$r
p <- rcorr(x)$P
mystars <- ifelse(p < .001, "***", ifelse(p < .01, "** ", ifelse(p < .05, "* ", " ")))
R <- format(round(cbind(rep(-1.11, ncol(x)), R), 3))[,-1]
Rnew <- matrix(paste(R, mystars, sep=""), ncol=ncol(x))
diag(Rnew) <- paste(diag(R), " ", sep="")
rownames(Rnew) <- colnames(x)
colnames(Rnew) <- paste(colnames(x), "", sep="")
Rnew <- as.matrix(Rnew)
Rnew[upper.tri(Rnew, diag = TRUE)] <- ""
Rnew <- as.data.frame(Rnew)

Rnew <- cbind(Rnew[1:length(Rnew)-1])
return(Rnew)
}

最后，我绘制了我的表格并导出到excel:

matrix_correlation <- corstarsl(provaSO1)
library(xtable)
corr_print <- xtable(matrix_correlation)
library(rio) 
library(openxlsx)
export(corr_print,"correlation_ok.xlsx")

我想知道是否有更优雅的方法来打印Pearson相关矩阵并将其放入Word文档(例如，回归表的Stargazer(。

提前谢谢。。

Answer 1

您可以使用R的modelsummary包将关联表直接保存到Word文档中。(免责声明：我是作者。(使用datasummary_correlation函数及其output参数，我们得到了：

library(modelsummary)
datasummary_correlation(provaSO1,
output = "correlation_table.docx")