r温度をARIMA时系列分析

篇首语：本文由编程笔记#小编为大家整理，主要介绍了r 温度をARIMA时系列分析相关的知识，希望对你有一定的参考价值。

library(forecast)
# データ読み込み
d = read.csv("temp.csv")
# unixtimeをPOSIXctに変換
tsPOSIX = as.POSIXct(d[,1], origin="1970-01-01")
d = cbind(tsPOSIX, d)
colnames(d)[1] = "ts"
head(d)
tail(d)
# auto.arimaを使って時系列分析
d.arima = auto.arima(d[, 2], stepwise = T, trace=T)
# 100期先まで予測する
d.arima.forecast = forecast(d.arima, level=c(50,95), h = 100)
# 予測値をプロット
plot(d.arima.forecast)