我正在尝试使用来自package psych（psych_1.6.9）的corr.test，但它似乎在使用method =“spearman”时从cor.test给出不同的p值。相关系数相同，但p值不同。我已经汇总了一些示例代码和输出，如下所示。
library(psych)
set.seed(5)
df<-data.frame("S1"=runif(5, min=1, max=3), "S2"=runif(5, min=1, max=3), "S3"=runif(5, min=1, max=3), "S4"=runif(5, min=1, max=3), "S5"=runif(5, min=1, max=3))
rownames(df)<-c("A","B","C","D","E")
dfm<-as.matrix(df)
x1<-head(dfm,n=3)
x2<-tail(dfm,n=2)
a<-corr.test(t(x1),t(x2),method="spearman",ci=FALSE, adjust="none")
a$p
D E
A 1.0000000 0.8728886
B 0.2847570 0.1881204
C 0.8728886 0.6238377
cor.test(dfm["B",],dfm["D",], method="spearman", alternative="two.sided")
Spearman's rank correlation rho
data: dfm["B", ] and dfm["D", ]
S = 32, p-value = 0.35
alternative hypothesis: true rho is not equal to 0
sample estimates:
rho
-0.6
a$r
D E
A 0.0 0.1
B -0.6 -0.7
C 0.1 -0.3
我必须做错事但却看不到。我很感激任何帮助。
Agnieszka