迴路中的R與cumsum多邊形功能

Question

我的工作循環中R：

dypol和wnc是1×3矩陣和x是1通過-100矩陣。 我想循環返回一個3×100的矩陣（每列累計）。我有這個：

For (i in 1:100) { 
i=dypol*t(x^2)-dypol+wnc 
{yi = cumsum(i) } 
} 

但它只返回第一行。

Answer 1

試試這個。我已轉換的循環成sapply呼叫（這有效地遍歷x），然後對所得到的矩陣的列施加cumsum：

x  <- 1:100 
wnc <- c(0.123, 0.263, 0.223) 
dypol <- c(.05, .30, .02) 

Z <- t(sapply(x, function(x)dypol * x^2 - dypol + wnc)) 
apply(Z, 2, cumsum)