2017-03-22 3 views
0

シーズンは24時間ごとですが、季節のARIMAモデルに合わせたいと思います。 しかし、24時間の季節をRに含めるにはどうすればよいですか?季節のARIMAモデルをRに入れるには

arima(y, order=c(0,0,2), seasonal=c(0,0,5), method = "ML")

しかし、私は正しいだ場合には、ARIMA(0,0,2)(0,0,5)_12モデルですので、私は助けを得ることを望む:これまでのところ、私は次のことを試してみましたその代わりにARIMA(0,0,2)(0,0,5)_24モデルを作成します。

+0

「シーズンは24時間ごとです」とはどういう意味ですか?あなたは毎日の中で繰り返しパターンを探していますか?シリーズの測定頻度はどのくらいですか? – ulfelder

+0

このリンク[https://www.otexts.org/fpp/8/9]があなたに役立つことを願って –

+0

あなたのデータセットの頻度は何ですか? – C8H10N4O2

答えて

1

period=seasonal=list(order=..., period=...)に含める必要があります。観測値が時間単位の場合は、period=24Lを使用します。 1秒あたりの場合、period=24*60*60などを使用してください。

例。

# reproducible example! 
# download file from: 
# https://trends.google.com/trends/explore?date=now%207-d&q=stackoverflow 
df <- read.csv('multiTimeline.csv', skip=3, header=FALSE, stringsAsFactors = FALSE) 
names(df) <- c('Time','Searches') 
df$Time <- as.POSIXlt.character(df$Time, tz='UTC',format = '%Y-%m-%dT%H') 

plot(df, type='l') 

m1 <- arima(x = df$Searches, 
      order = c(0L,0L,2L), 
      seasonal=list(order=c(0L,0L,5L), period=24L) 
) 

> m1 

Call: 
arima(x = df$Searches, order = c(0L, 0L, 2L), seasonal = list(order = c(0L, 
    0L, 5L), period = 24L)) 

Coefficients: 
     ma1  ma2 sma1 sma2  sma3 sma4 sma5 intercept 
     1.0827 0.6160 0.6155 0.1403 -0.1472 0.0104 0.6807 52.1477 
s.e. 0.0631 0.0566 0.2305 0.2005 0.1445 0.2210 0.2176  2.4497 

sigma^2 estimated as 35.69: log likelihood = -575.94, aic = 1169.88 

?arimaを見る

seasonal ARIMAモデル、 の季節一部の仕様プラス(frequency(x)デフォルト)period。これは listとコンポーネントorderperiodである必要がありますが、長さ3の数値ベクトルの 数値ベクトルの指定は、 の順で適切なリストになります。

関連する問題