→ qwer098:好像有誤?..y=1~100只是方便說明..實際上是沒規律的數列 01/22 23:13
→ KnightX:那就把 (i-j)*(i-k) 改成 y[i-j]*y[i-k] 01/22 23:45
※ 引述《qwer098 (^^)》之銘言：
: 若Y_t=(y_t,....,y_(t-p+1))'
: 現在想算Y_t*Y_(t-1)'的和(t=p+1~n)
: R跑出來的結果有誤....
: 它會說錯誤在a[i - p, i - j + 1, i - k + 1] <- y[j - 1] * y[k - 1] :
: 下標超出邊界
: 錯誤: 下標超出邊界
: 請大家幫我看一下我該如何修改...謝謝
指標有點亂, 應該可以更簡化,
而且既然 y 是 1~100 那也可以用指標來取代,
以下是不負責任的改寫, 請自行驗證：
A <- matrix(sapply(3:100, function(i) {
sapply(1:2, function(j) {
sapply(1:2, function(k) {
(i-j)*(i-k)
})
})
}), c(98, 2, 2))
B <- apply(A, 2:3, sum)
# 註：我程式裡的 A 是你原本的 a, B 是 b
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 218.169.112.41