corM=function(X){
n=ncol(X)
M=matrix(0,nrow=n,ncol=n)
for(i in 1:n){
for(j in i:n){
M[i,j]=cor.test(X[,i],X[,j])$p.value
}
}
return(M)
}
我不知道我理解的对不对,你是想求X这个数据矩阵各列之间的相关系数的显著概率的话,应该可以用上面这个函数,它会输出一个上三角矩阵,矩阵的i行j列的数值代表X这个矩阵第i列和第j列的相关系数的显著概率。我没用过cor.test这个函数,随便编的,你看看对不对吧。