如何替换矩阵的下/上三角形元素？

Question

我的问题：

Amat <- diag(4)

我想用我选择的值替换Amat（即Amat[2,1]，Amat[3,1]，Amat[3,2]等）的所有下三角值（例如） NA）。

显然我不想逐个替换每个元素。

你能用一个命令告诉我最有效的方法吗？

Answer 1

upper.tri的文档中详细记录了这一点。

Amat[upper.tri(Amat)] <- NA
Amat
#      [,1] [,2] [,3] [,4]
# [1,]    1   NA   NA   NA
# [2,]    0    1   NA   NA
# [3,]    0    0    1   NA
# [4,]    0    0    0    1

当然，Amat[lower.tri(Amat)] <- NA会将下三角形转换为NA s。

Answer 2

您要找的是lower.tri和upper.tri吗？

这些功能在R base。