在用R做了两次循环时,结果只能输出第一列。程序为f
or(i in 1:10)
for(j in 1:10){
a<-pcor.test(rd[,i],rd[,j],rd[,-c(i,j)],use = c("mat"))
partial[i][j]<-a$estimate
结果如下
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,] 1.00000000 NA NA NA NA NA NA NA NA NA
[2,] 0.57819745 NA NA NA NA NA NA NA NA NA
[3,] 0.32088018 NA NA NA NA NA NA NA NA NA
[4,] -0.33626631 NA NA NA NA NA NA NA NA NA
[5,] -0.34627386 NA NA NA NA NA NA NA NA NA
[6,] 0.02538884 NA NA NA NA NA NA NA NA NA
[7,] -0.16991169 NA NA NA NA NA NA NA NA NA
[8,] 0.14118949 NA NA NA NA NA NA NA NA NA
[9,] 0.27155041 NA NA NA NA NA NA NA NA NA
[10,] -0.06197557 NA NA NA NA NA NA NA NA NA
>