[程式] 用SAS模擬ARMA (附上code 10.06.00:16)

看板Statistics作者 (賴打)時間15年前 (2010/10/05 21:28), 編輯推噓3(303)
留言6則, 3人參與, 最新討論串1/1
[軟體程式類別]: SAS [程式問題]: 不知道該如何模擬 [軟體熟悉度]: 難以說明...... 兩年前的實習課曾經學過一年的SAS 但是老師(+助教)專長在個體計量 SAS有關時間序列部份的code完全沒教 [問題敘述]: 簡而言之就是要模擬一個ARMA process 例如: y(t)=1.5y(t-1)+0.5y(t-2)+ε(t)+0.3ε(t-1) output 100筆數值+圖 有google過看看網路上有沒有類似的說明 但是沒有令我滿意的結果 所以上來問問看 煩請賜教 謝謝~ [code]: y1=0.9; y11=0.9; u1 = 0; do time = -50 to 100; u = rannor( 32565 ); y =1.5*y1+0.5*y11+u+0.3*u1 ; if time > 0 then output; y11=y1; y1=y; u1 = u; end; run; p.s 1.起始值是個很重要的問題嗎? 2.有沒有推薦什麼SAS關於時間序列的網站 -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 203.73.175.130 ※ 編輯: saar 來自: 203.73.175.130 (10/05 22:50) ※ 編輯: saar 來自: 203.73.175.130 (10/06 00:16) ※ 編輯: saar 來自: 203.73.175.130 (10/06 00:51)

10/07 11:37, , 1F
模擬而言起始值不重要 符合WHITE NOISE即可 跑完再踢除就好
10/07 11:37, 1F

10/07 11:40, , 2F
這樣跑是跑不出ARMA(2,2)的 你要先作兩個WHITE NOISE當起始
10/07 11:40, 2F

10/07 11:40, , 3F
按照定義 去一筆一筆填入WHITE NOISE
10/07 11:40, 3F

10/07 22:49, , 4F
謝謝 我暫時先用R解決了 網路上有範例 SAS再研究看看......
10/07 22:49, 4F

10/09 02:55, , 5F
這語法應該沒錯 結果不滿意不表示語法錯誤
10/09 02:55, 5F

10/09 13:55, , 6F
可以跑樣本大一點 再直接用手算確認是否模擬正確 R王道阿
10/09 13:55, 6F
文章代碼(AID): #1CgoVoBh (Statistics)