如何重複循環？

Question

這是我的第一個循環：

e13 = rnorm(122) 

y = 0*e13 

y[1]=2 
y[2]=5 

for(k in 3:length(e13)) { 
    y[k]= 0.6 * y[k-1] + e13[k] 
} 

現在我想重複10次，以找到另外10個值相同的Y [k]的一樣，10個不同的y[1]，10個不同的y[2]等。

然後我不需要使用這些數字作爲其他功能。

我嘗試過，但沒有奏效

for (i in 1:20) { 
    for(k in 3:length(e13)) { 
    y[k,i] = 0.6 * y[k-1,i] + e13[k,i] 
    } 
} 

它說：「維打錯了」

現在我的代碼是這樣的：

e13 = rnorm(122) 
y <- array(0, c(length(e13), 20)) 
y[1, ] <- 2 
y[2, ] <- 5 
i<-c(1:20) 
y[1]=2 
y[2]=5 

for (i in 1:20) { 
    for(k in 3:length(e13)) { 
    y[k,i] = 0.6 * y[k-1,i] + e13[k,i] 
    } 
} 

但現在的問題是e13的維度（正態分佈），即時試圖設置正確的維度

Answer 1

他re是一個簡單的解決方案，雖然cbind會使您的數據減慢。如果需要更多速度，則不需要指定2維，則可以將數據矩陣指定爲長度爲length(e13)*(how many times you want to repeat)的向量，在此情況下爲length(e13)*20。然後在第二個循環中，每次循環的重複都會沿着數據矩陣的長度進一步放置。有時候這樣做更簡單，如果你想使用matrix，你最後把它變成矩陣。

datamatrix=NULL 
for(n in 1:20) { 
e13 = rnorm(122) 

y = 0*e13 

y[1]=2 
y[2]=5 

for(k in 3:length(e13)) { 
    y[k]= 0.6 * y[k-1] + e13[k] 
} 
datamatrix=cbind(datamatrix,y)} 

Answer 2

你應該採取一些介紹[R文本（例如，this，this，this，和其他人）好好看看。對你的問題的評論引導你自己解決這個問題。你幾乎在那裏，但這裏是我如何做你想做的事情。

e13 <- matrix(rnorm(122*20), ncol=20) 
y <- array(0, dim(e13)) 
y[1, ] <- 2 
y[2, ] <- 5 

for (i in 1:ncol(y)) { 
    for(k in 3:nrow(y)) { 
    y[k, i] <- 0.6 * y[k-1, i] + e13[k, i] 
    } 
}