上三角矩阵,对角线为零

时间:2017-08-01 10:27:18

标签: r

我尝试创建一个上三角矩阵,对角线中的条目为零,但我没有这样做。

这是我的代码:

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

我的错误在哪里?

然后,如何从我的矩阵(上三角形,零对角线)顺序选择或访问非零元素,即元素接连?

请帮忙吗?

0 个答案:

没有答案