1

我使用以下for循环生成了一个矩阵。
现在我正在尝试使用while循环生成相同的矩阵,但不知道该怎么做。
有人能帮忙吗?非常感谢。

i<- 1
j<- 1
m1<- matrix(0, nrow=5, ncol=5)
for(i in 1:5) {
for(j in i:5) {
m1[i,j]<- (j-i)*2+1
}
}
m1
4

2 回答 2

2
i <- 1
j <- 1
m1 <- matrix(0, nrow = 5, ncol = 5)
while (i <= 5) {
    while (j <= 5) {
        m1[i, j] <- (j - i) * 2 + 1
        j <- j + 1
    }
    i <- i + 1
    j <- i
}
m1 
于 2012-10-06T01:55:35.520 回答
2

那这个呢:

> m       <- matrix(0, nrow=5, ncol=5)
> mu      <- upper.tri(m)
> m[mu]   <- apply(which(mu, arr.ind=TRUE), 1, function(x) (x[2] - x[1]) * 2 + 1)
> diag(m) <- 1
> m
     [,1] [,2] [,3] [,4] [,5]
[1,]    1    3    5    7    9
[2,]    0    1    3    5    7
[3,]    0    0    1    3    5
[4,]    0    0    0    1    3
[5,]    0    0    0    0    1

编辑:抱歉,刚刚意识到您不想对任务进行矢量化,而是重写代码以使用while而不是for. 虽然我不知道为什么:)

于 2012-10-06T10:06:18.027 回答