作者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