我是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
你有什么线索吗?