Secara umum, gali ke dalam buku teks analisis deret waktu tingkat lanjut (buku pengantar biasanya akan mengarahkan Anda untuk hanya mempercayai perangkat lunak Anda), seperti Analisis Rangkaian Waktu oleh Box, Jenkins & Reinsel. Anda juga dapat menemukan detail tentang prosedur Box-Jenkins dengan googling. Perhatikan bahwa ada pendekatan lain selain Box-Jenkins, misalnya yang berbasis AIC.
Di R, Anda pertama-tama mengonversi data Anda menjadi objek ts
(seri waktu) dan memberi tahu R bahwa frekuensinya adalah 12 (data bulanan):
require(forecast)
sales <- ts(c(99, 58, 52, 83, 94, 73, 97, 83, 86, 63, 77, 70, 87, 84, 60, 105, 87, 93, 110, 71, 158, 52, 33, 68, 82, 88, 84),frequency=12)
Anda dapat memplot fungsi autokorelasi (sebagian):
acf(sales)
pacf(sales)
Ini tidak menyarankan perilaku AR atau MA.
Kemudian Anda cocok dengan model dan memeriksanya:
model <- auto.arima(sales)
model
Lihat ?auto.arima
bantuan. Seperti yang kita lihat, auto.arima
pilih model sederhana (0,0,0), karena tidak melihat tren atau musiman atau AR atau MA dalam data Anda. Akhirnya, Anda dapat memperkirakan dan merencanakan rangkaian waktu dan perkiraan:
plot(forecast(model))
Lihatlah ?forecast.Arima
(perhatikan ibukota A!).
Buku teks online gratis ini merupakan pengantar yang bagus untuk analisis deret waktu dan perkiraan menggunakan R. Sangat direkomendasikan.
Dua hal. Rangkaian waktu Anda adalah bulanan, Anda memerlukan setidaknya 4 tahun data untuk estimasi ARIMA yang masuk akal, sebagaimana tercermin 27 poin tidak memberikan struktur autokorelasi. Ini juga dapat berarti bahwa penjualan Anda dipengaruhi oleh beberapa faktor eksternal, bukannya berkorelasi dengan nilainya sendiri. Coba cari tahu faktor apa yang memengaruhi penjualan Anda dan apakah faktor itu sedang diukur. Kemudian Anda dapat menjalankan regresi atau VAR (Vector Autoregression) untuk mendapatkan perkiraan.
Jika Anda benar-benar tidak memiliki hal lain selain nilai-nilai ini, cara terbaik Anda adalah menggunakan metode smoothing eksponensial untuk mendapatkan perkiraan naif. Perataan eksponensial tersedia dalam R.
Kedua tidak melihat penjualan produk secara terpisah, penjualan dua produk mungkin berkorelasi misalnya peningkatan penjualan kopi dapat mencerminkan penurunan penjualan teh. gunakan informasi produk lainnya untuk meningkatkan perkiraan Anda.
Ini biasanya terjadi dengan data penjualan di ritel atau rantai pasokan. Mereka tidak menunjukkan banyak struktur autokorelasi dalam seri ini. Sementara di sisi lain metode seperti ARIMA atau GARCH biasanya bekerja dengan data pasar saham atau indeks ekonomi di mana Anda biasanya memiliki autokorelasi.
sumber
Ini benar-benar sebuah komentar tetapi melebihi yang diizinkan jadi saya mempostingnya sebagai jawaban semu karena menyarankan cara yang benar untuk menganalisis data deret waktu. .
Fakta yang terkenal tetapi sering diabaikan di sini dan di tempat lain adalah bahwa ACF / PACF teoretis yang digunakan untuk merumuskan model tentatif ARIMA tidak menggunakan Pulsa / Pergeseran Tingkat / Pulsa Musiman / Pulsa Musiman / Tren Waktu Lokal. Selain itu itu menempatkan parameter konstan dan varians kesalahan konstan sepanjang waktu. Dalam hal ini, observasi ke-21 (nilai = 158) dengan mudah ditandai sebagai pencilan / Denyut nadi dan penyesuaian yang disarankan dari -80 menghasilkan nilai yang dimodifikasi dari 78. ACF / PACF yang dihasilkan dari seri yang dimodifikasi menunjukkan sedikit atau tidak ada bukti struktur stokastik (ARIMA). Dalam hal ini operasi berhasil tetapi pasien meninggal. Sampel ACF didasarkan pada kovarians / varians dan varians yang terlalu tinggi / kembung menghasilkan bias ke bawah pada ACF. Prof Keith Ord pernah menyebut ini sebagai "efek Alice in Wonderland"
sumber
Seperti yang telah ditunjukkan oleh Stephan Kolassa, tidak ada banyak struktur dalam data Anda. Fungsi autokorelasi tidak menyarankan struktur ARMA (lihat
acf(sales)
,pacf(sales)
) danforecast::auto.arima
tidak memilih perintah AR atau MA.Namun demikian, perhatikan bahwa nol normalitas dalam residu ditolak pada tingkat signifikansi 5%.
Selain Catatan:
JarqueBera.test
didasarkan pada fungsi yangjarque.bera.test
tersedia dalam pakettseries
.Termasuk outlier aditif pada pengamatan 21 yang terdeteksi dengan
tsoutliers
membuat normalitas dalam residu. Dengan demikian, estimasi intersep dan ramalan tidak terpengaruh oleh pengamatan terluar.sumber