[公告] 時序作業 5/20交
5/20交
#7.4 #7.7
#8.4 #8.5
期中報告5/9交
可以用的指令:
#arima.sim :生成time series data
set.seed(5, kind = NULL) #固定隨機變數的種子
z=arima.sim(n=1000,list(ar=0.3),sd = sqrt(0.1796),n.start=1)
#驗證
set.seed(5, kind = NULL) #固定隨機變數的種子,以防模擬誤差
za=rep(0,1001)
for(i in 2:1001){
u=rnorm(1,mean=0,sd=sqrt(0.1796))
za[i]=0.3*za[i-1]+u
}
#可以比較一下 arima.sim丟2個,所以從第3個以後就一樣了。
# ie za[3] = z[1] , za[4] = z[2] , za[5]=z[3]...
#arima:求參數:
mod=arima(z,order=c(1,0,0))
names(mod) #可以看有那些物件
#[1] "coef" "sigma2" "var.coef" "mask" "loglik" "aic"
#[7] "arma" "residuals" "call" "series" "code" "n.cond"
#[13] "model"
#用錢就可以取出來
mod$coef
# ar1 intercept
#0.27688839 0.01020825
#關於物件的說明請參考 arima的help檔 value的部分
#不過,它沒寫耶! 大家自行研究吧!
acf
acfz=acf(z,lag.max=100) #這樣就會生出100個,不夠的話請自行選擇
#同樣的也可以看看acfz裡面有什麼
pacf
#自行研究…
對了!第六章的作業大家分數都加2分喔!
期中報告要好好做囉~~
--
已經沒有利用價值的助教要去準備考試了!
最好有那個好心人來教一下數統^^
--
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.114.36.142
推
04/30 19:18, , 1F
04/30 19:18, 1F
推
04/30 19:38, , 2F
04/30 19:38, 2F
推
05/01 00:22, , 3F
05/01 00:22, 3F
→
05/01 00:22, , 4F
05/01 00:22, 4F
推
05/01 01:11, , 5F
05/01 01:11, 5F
→
05/01 13:55, , 6F
05/01 13:55, 6F
※ 編輯: tsta 來自: 140.114.36.142 (05/01 14:04)
推
05/01 21:30, , 7F
05/01 21:30, 7F