r - 使用 R 的非零零假设的相关显着性

Question

我正在测试两个变量之间的相关性:

set.seed(123)
x <- rnorm(20)
y <- x + x * 1:20
cor.test(x, y, method = c("spearman"))

这使:

Spearman's rank correlation rho

data:  x and y 
S = 54, p-value = 6.442e-06
alternative hypothesis: true rho is not equal to 0 
sample estimates:
   rho 
0.9594 

p 值正在检验相关性为零的原假设。是否有一个 R 函数可以让我测试不同的零假设——比如相关性小于或等于 0.3？

Answer 1

您可以使用 bootstrap 计算 rho 的置信区间:

1) Make 函数来提取 cor.test 的估计值(记住放置索引以便引导可以对数据进行采样):

rho <- function(x, y, indices){
  rho <- cor.test(x[indices], y[indices],  method = c("spearman"))
  return(rho$estimate)
}

2) 使用 boot包来引导您的估计:

library(boot)    
boot.rho <- boot(x ,y=y, rho, R=1000)

3)取置信区间:

boot.ci(boot.rho)