3

我正在尝试使用 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

我一定做错了什么,但就是看不到。我将不胜感激这方面的任何帮助。

阿涅斯卡

4

0 回答 0