作者KnightX (........)
看板Statistics
标题Re: [问题]请问哪里有错(R程式)?
时间Thu Jan 22 22:57:43 2009
※ 引述《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
1F：→ qwer098:好像有误?..y=1~100只是方便说明..实际上是没规律的数列 01/22 23:13
2F：→ KnightX:那就把 (i-j)*(i-k) 改成 y[i-j]*y[i-k] 01/22 23:45