在R中循环，带有用于多边形功能的cumsum

Question

我正在研究R：

的循环

dypol和wnc是1 x 3矩阵，x是1 x 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)