作者qwer098 (^^)
看板Statistics
标题[问题]请问哪里有错(R程式)?
时间Thu Jan 22 19:16:43 2009
若Y_t=(y_t,....,y_(t-p+1))'
现在想算Y_t*Y_(t-1)'的和(t=p+1~n)
程式为
a<-array(0,c(98,2,2))
b<-matrix(0,2,2)
y<-c(1:100)
p=2
n=100
for(i in p+1:n)
{
for(j in i:i-p+1)
{
for(k in i:i-p+1)
{
a[i-p,i-j+1,i-k+1]<-y[j-1]*y[k-1]
}
}
}
for(i in p+1:n)
{
b<-a[i,,]+b
}
R跑出来的结果有误....
它会说错误在a[i - p, i - j + 1, i - k + 1] <- y[j - 1] * y[k - 1] :
下标超出边界
错误: 下标超出边界
请大家帮我看一下我该如何修改...谢谢
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 125.224.50.52
1F：→ sleepdragon:没看完你的程式码 01/22 21:38
2F：→ sleepdragon:不过 for()里面要改成( i in (p+1):n) 01/22 21:39