0

私は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) }
}

ただし、最初の行のみを返します。

4

1 に答える 1

0

これを試して。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)
于 2012-10-07T14:52:41.750 に答える