从1-ibs矩阵创建数字序列(组)

时间:2018-11-21 23:03:49

标签: python r windows matrix

我是R的新手,我正尝试从矩阵生成这种序列(请参见下文),然后估计错误率。只是让您知道,每个数字都有三份,这就是为什么我有X1,X1.1和X1.2

              X1        X1.1  X1.2
X1   0.000000000          NA   NA
X1.1 0.013452613 0.000000000   NA
X1.2 0.008110273 0.007781853    0  


             X73      X73.1   X73.2
X73   0.00000000         NA    NA
X73.1 0.05336268 0.00000000    NA
X73.2 0.04450677 0.05150902     0

以此类推...

我应用了此代码(NB:261个样本)

mat <- 1-ibsmat
mat1 <- mat[-99:-100,-99:-100]

x <-  seq(1,261,3)
err.r <-  rep(NA, length(x))
for (i in 1:(length(x)-1)){
  k <- x[i]
  k1=k+2
  ibx <- mat1[k:k1,k:k1]
  print(ibx)
  err.r[i] <- mean(ibx[lower.tri(ibx)])
}
errorrate <- mean(na.omit(err.r))

一切正常,直到数字开始与其他数字重叠(请参见下面的示例)。

           X91.1     X91.2  X12
X91.1 0.00000000        NA  NA
X91.2 0.01587456 0.0000000  NA
X12   0.10722157 0.1139528   0

你有什么线索吗?

0 个答案:

没有答案