如何在R中计算Spearman相关系数的置信区间？

你可以通过以下步骤，使用R的DescTools包中的SpearmanRho函数获取Spearman相关系数及其置信区间：

步骤1：安装并加载DescTools包 #

如果还没安装DescTools包，先执行安装命令：

install.packages("DescTools")

安装完成后加载包：

library(DescTools)

步骤2：准备对应的数据 #

将你Python代码中的x和y数据转换为R的向量：

x <- c(1.76405235, 0.40015721, 0.97873798, 2.2408932, 1.86755799, -0.97727788)
y <- c(2.71414076, 0.2488, 0.87551913, 2.6514917, 2.01160156, 0.47699563)

步骤3：调用SpearmanRho函数计算结果 #

直接调用SpearmanRho函数，默认就会返回置信区间（默认置信水平95%，可通过conf.level参数调整）：

result <- SpearmanRho(x, y)
print(result)

结果解释 #

执行后会输出类似如下的结果：

Spearman's rank correlation rho

data: x and y
S = 8, rho = 0.8571, p-value = 0.09524
alternative hypothesis: true rho is not equal to 0
95 percent confidence interval:
0.1464467 0.9828254
sample estimates:
rho
0.8571429

其中：

• rho是Spearman相关系数
• p-value是检验的p值
• 95 percent confidence interval就是你需要的相关系数置信区间

你也可以单独提取置信区间：

result$conf.int

内容的提问来源于stack exchange，提问作者Simd