パッケージ 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
私は何か間違ったことをしているに違いないのですが、それを見ることができません。これについて何か助けていただければ幸いです。
アニエスカ