[問題]請問哪裡有錯(R程式)?

看板Statistics作者 (^^)時間17年前 (2009/01/22 19:16), 編輯推噓0(002)
留言2則, 1人參與, 最新討論串1/3 (看更多)
若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

01/22 21:38, , 1F
沒看完你的程式碼
01/22 21:38, 1F

01/22 21:39, , 2F
不過 for()裡面要改成( i in (p+1):n)
01/22 21:39, 2F
文章代碼(AID): #19U5MSO7 (Statistics)
文章代碼(AID): #19U5MSO7 (Statistics)