r - 如何在 R 中制作 APA 就绪的相关表

Question

我想将 R 中的互变量相关表转换为根据 APA 样式格式化的表。该表应具有以下特点

• 只有矩阵的下对角线或上对角线
• 通常只显示带有 2 位小数且没有前导零的相关系数，例如 0.25、-.05
• 用星号表示重要性

Answer 1

我曾经使用Hmisc创建相关表，p-values但没有包这是可行的，特别是如果你想要 Kendall 的 tau 相关性。

#First indicate which columns do you want in your correlation table
cor_cols = c('gender','age','education','expense',
              'learningmotivation', 'reading')


#Slice the data and calculate correlation table:
mt = as.matrix(data[,cor_cols])
library(Hmisc)
myrcorr <- rcorr(mt, type="spearman")

#format numbers
cor_table <-as.matrix( round(myrcorr$r,2))

#remove rownames and colnames
colnames(cor_table)<- rownames(cor_table)<- NULL

#remove leading zeros and add stars
n <- length(cor_cols)
for (c in 1:n) {
  for (r in 1:n) {
      pval <- myrcorr$P[r,c]
      stars <- ifelse(pval < .001, "***", ifelse(pval < .01, "** ", ifelse(pval < .05, "*  ", "   ")))
      coeff = sub("0*\\.",".",cor_table[r,c]) #remove leading zeros
      cor_table[r,c] = paste(c(coeff,stars),collapse='')
  }
}

cor_table[upper.tri(cor_table)] <- '' # erase the upper triangle
diag(cor_table) <- '-' # replace the diagonals by dash(-)