我在网上找到了替代方法,有人直接写了程序
COR.test=function(X,R){options(digits=4) #求F检验的p值,为矩阵形式,R为简单相关系数矩阵,X为数据集
n=dim(X)[2] #得到p矩阵的阶数
p=diag(0,n) #n阶零矩阵
for(i in 1:n)
{
for(j in 1:n){
f=R[i,j]^2/((1-R[i,j]^2)/(dim(X)[1]-2)) #用F检验对相关系数作显著性检验
p[i,j]=1-pf(f,1,dim(X)[1]-2) #用F检验计算p值
}
}
p
}
COR.test(julei1,cor(julei1))
程序来源:https://wenku.baidu.com/view/5a24ef34f61fb7360a4c65a0.html