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

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

01/22 23:13, , 1F
好像有誤?..y=1~100只是方便說明..實際上是沒規律的數列
01/22 23:13, 1F

01/22 23:45, , 2F
那就把 (i-j)*(i-k) 改成 y[i-j]*y[i-k]
01/22 23:45, 2F
文章代碼(AID): #19U8bfdE (Statistics)
文章代碼(AID): #19U8bfdE (Statistics)