我尝试创建一个上三角矩阵,对角线中的条目为零,但我没有这样做。
这是我的代码:
w1 <- 0.5
pi <- matrix(w1,5,5)
pi[upper.tri(pi)] <- 0
pi[diag(pi)] <- 0
pi
我希望我的矩阵如下:
> pi
[,1] [,2] [,3] [,4] [,5]
[1,] 0.0 0.0 0.0 0.0 0.0
[2,] 0.5 0.0 0.0 0.0 0.0
[3,] 0.5 0.5 0.0 0.0 0.0
[4,] 0.5 0.5 0.5 0.0 0.0
[5,] 0.5 0.5 0.5 0.5 0.0
然而,我得到了这个:
> pi
[,1] [,2] [,3] [,4] [,5]
[1,] 0.5 0.5 0.5 0.5 0.5
[2,] 0.5 0.5 0.5 0.5 0.5
[3,] 0.5 0.5 0.5 0.5 0.5
[4,] 0.5 0.5 0.5 0.5 0.5
[5,] 0.5 0.5 0.5 0.5 0.5
我的错误在哪里?
然后,如何从我的矩阵(上三角形,零对角线)顺序选择或访问非零元素,即元素接连?
请帮忙吗?