我想計算
sum_{i}^{n} x_{i}^2 * lambda^(n-i) / sum_{i}^{n} lambda^(n-i)
在 R 中這樣做我設法做到以下幾點。但是出現了零點。我做錯了什么？
n = 10
x = seq(1,n,1);x
lambda = 0.99
mat = 0
for (i in 2:n) {
mat = (lambda^(n-i)*x[i-1]^2) / (lambda^(n-i))
}
mat
結果必須是每個條目的方差x_{i}
uj5u.com熱心網友回復：
很好，我會試著看看這是否是你需要的@The Red：
n <- 10
x <- seq(1,n,1);x
lambda <- 0.99
mat <- 0
vect_mat <- rep(0, n-1); vect_mat
for (i in 2:n) {
mat <- mat (lambda^(n-i)*x[i-1]^2) / (lambda^(n-i))
vect_mat[i-1] <- mat
}
mat
vect_mat
運行后，結果如下：
> mat
[1] 285
> vect_mat
[1] 1 5 14 30 55 91 140 204 285
如果這是您需要的，請接受答案。
轉載請註明出處，本文鏈接：https://www.uj5u.com/ruanti/427742.html