将R中的矩阵转换成上/下三角矩阵除以次对角线

Question

我知道使用upper.tri（）和lower.tri（）我可以得到上下三角矩阵除以主对角线。

但是，得到三角矩阵除以次对角线的最快方法是什么？

Answer 1

只需修改lower.tri的代码，例如，

lower.anti.tri<-function(m) col(m)+row(m) > dim(m)[1]+1 

m<-matrix(1:16,4) 
lower.anti.tri(m) 

 
     [,1] [,2] [,3] [,4] 
[1,] FALSE FALSE FALSE FALSE 
[2,] FALSE FALSE FALSE TRUE 
[3,] FALSE FALSE TRUE TRUE 
[4,] FALSE TRUE TRUE TRUE 

m[lower.anti.tri(m)]<-NA 
m 

 
    [,1] [,2] [,3] [,4] 
[1,] 1 5 9 13 
[2,] 2 6 10 NA 
[3,] 3 7 NA NA 
[4,] 4 NA NA NA 

Answer 2

您可以使用apply(mat, 1, rev)水平翻转矩阵状翻转lower.tri()：

m <- matrix(1:16,4) 
m[apply(lower.tri(m), 1, rev)] <- 0 

m 
#  [,1] [,2] [,3] [,4] 
#[1,] 1 5 9 13 
#[2,] 2 6 10 0 
#[3,] 3 7 0 0 
#[4,] 4 0 0 0