用 NA 替换对称矩阵的元素对

Question

我有一个正定对称矩阵。粘贴使用以下代码生成的矩阵：

set.seed(123)
m <- genPositiveDefMat(
    dim = 3, 
    covMethod = "unifcorrmat",
    rangeVar = c(0,1) )

x <- as.matrix(m$Sigma)
  diag(x) <- 1

x

#Output
          [,1]       [,2]       [,3]
[1,]  1.0000000 -0.2432303 -0.4110525
[2,] -0.2432303  1.0000000 -0.1046602
[3,] -0.4110525 -0.1046602  1.0000000

现在，我想通过迭代运行矩阵，并且在每次迭代中我想用 NA 替换对称对。例如，

迭代 1： x[1,2] = x[2,1]

迭代2： x[1,3] = x[3,1]

等等……

我的想法是使用 for 循环检查

原型：

for( r in 1:nrow(x)
for( c in 1:ncol(x)

if x[r,c]=x[c,r]<-NA

else

x[r,c]

我的代码的问题在于第 1 行和第 1 列，值相等，因此设置为 0（这是错误的）。此外，当它不是 NA 时，它就会退出循环。

在这里感谢任何帮助。

谢谢

Answer 1

如果您需要迭代完成替换，您可以使用upper.tri(x)/lower.tri 表示的值的索引来逐对进行替换。这将允许您在每次替换之前/之后将结果传递给函数，例如：

idx <- which(lower.tri(mat), arr.ind=TRUE)
sel <- cbind(
  replace(mat, , seq_along(mat))[ idx ],
  replace(mat, , seq_along(mat))[ idx[,2:1] ]
)
#     [,1] [,2]
#[1,]    2    4 ##each row represents the lower/upper pair
#[2,]    3    7
#[3,]    6    8

for( i in seq_len(nrow(sel)) ) {
    mat[ sel[i,] ] <- NA
    print(mean(mat, na.rm=TRUE))
}
#[1] 0.2812249
#[1] 0.5581359
#[1] 1