我通过rcorr计算了编码编码的相关性。现在我想通过此代码过滤和提取结果,但是有一个错误。为什么代码有问题?
代码:
library(Hmisc)
coef_BS<- matrix(data = res5$r, nrow =20 ,ncol=20)
pval_BS<- matrix(data = res5$P, nrow =20 ,ncol=20)
for(i in 1:20)
{
for (j in 1:20)
{
if (i>=j)
{
rcor_out<-rcorr(c.num[,i], c.nc.num[,j], type="spearman")
coef_BS[i,j] <- rcor_out$r[1, 2]
pval_BS[i,j] <- rcor_out$P[1, 2]
}
}
print(i)
}
错误是:
Error in c.nc.num[, j] : subscript out of bounds