Moving average process in r

Moving Averages in R Sejauh pengetahuan saya, R tidak memiliki fungsi built-in untuk menghitung moving averages. Dengan menggunakan fungsi filter, kita dapat menulis fungsi pendek untuk moving averages: Kita kemudian dapat menggunakan fungsi pada data: mav (data), atau mav (data, 11) jika kita ingin menentukan jumlah titik data yang berbeda. Dari pada default 5 plotting works seperti yang diharapkan: plot (mav (data)). Selain jumlah titik data yang rata-rata, kita juga dapat mengubah argumen sisi fungsi filter: sisi kedua menggunakan kedua sisi, sisi1 hanya menggunakan nilai masa lalu. Share this: Post navigation Comment navigation Comment navigation2.1 Moving Average Models (model MA) Model deret waktu yang dikenal dengan model ARIMA mungkin termasuk istilah autoregresif dan atau istilah rata-rata bergerak. Pada minggu ke 1, kita belajar istilah autoregressive dalam model time series untuk variabel x t adalah nilai lag dari x t. Misalnya, lag 1 autoregressive term adalah x t-1 (dikalikan dengan koefisien). Pelajaran ini mendefinisikan istilah rata-rata bergerak. Istilah rata-rata bergerak dalam model time series adalah kesalahan masa lalu (dikalikan dengan koefisien). Misalkan (wt overset N (0, sigma2w)), yang berarti bahwa w t identik, didistribusikan secara independen, masing-masing dengan distribusi normal memiliki mean 0 dan varian yang sama. Model rata-rata bergerak urutan 1, dilambangkan dengan MA (1) adalah (xt mu wt theta1w) Model rata-rata bergerak urutan 2 yang dinotasikan dengan MA (2) adalah (xt mu wt theta1w theta2w) Model rata-rata pergerakan harga th q th , Dilambangkan dengan MA (q) adalah (xt mu wt theta1w theta2w titik thetaqw) Catatan. Banyak buku teks dan program perangkat lunak menentukan model dengan tanda negatif sebelum persyaratan. Ini tidak mengubah sifat teoritis umum dari model, meskipun ia membalik tanda aljabar dari nilai koefisien perkiraan dan (unsquared) terms dalam formula untuk ACF dan varians. Anda perlu memeriksa perangkat lunak Anda untuk memverifikasi apakah tanda negatif atau positif telah digunakan untuk menuliskan model perkiraan dengan benar. R menggunakan tanda-tanda positif pada model dasarnya, seperti yang kita lakukan di sini. Sifat Teoritis dari Seri Waktu dengan Model MA (1) Perhatikan bahwa satu-satunya nilai nol pada ACF teoritis adalah untuk lag 1. Semua autokorelasi lainnya adalah 0. Jadi sampel ACF dengan autokorelasi signifikan hanya pada lag 1 adalah indikator dari model MA (1) yang mungkin. Bagi siswa yang tertarik, bukti sifat ini adalah lampiran untuk handout ini. Contoh 1 Misalkan model MA (1) adalah x t 10 w t .7 w t-1. Dimana (wt overset N (0,1)). Dengan demikian koefisiennya 1 0,7. ACF teoritis diberikan oleh sebidang ACF berikut. Plot yang baru saja ditunjukkan adalah ACF teoritis untuk MA (1) dengan 1 0,7. Dalam prakteknya, contoh biasanya akan memberikan pola yang jelas. Dengan menggunakan R, kita mensimulasikan n 100 nilai sampel menggunakan model x t 10 w t .7 w t-1 dimana w t iid N (0,1). Untuk simulasi ini, plot time series dari data sampel berikut. Kami tidak tahu banyak dari plot ini. Contoh ACF untuk data simulasi berikut. Kita melihat lonjakan pada lag 1 diikuti oleh nilai-nilai non-signifikan secara umum untuk kelambatan masa lalu 1. Perhatikan bahwa sampel ACF tidak sesuai dengan pola teoritis MA yang mendasarinya, yaitu bahwa semua autokorelasi untuk kelambatan 1 akan menjadi 0 Sampel yang berbeda akan memiliki sampel ACF yang sedikit berbeda yang ditunjukkan di bawah ini, namun kemungkinan memiliki fitur luas yang sama. Sifat Teori dari Seri Waktu dengan Model MA (2) Untuk model MA (2), sifat teoretis adalah sebagai berikut: Perhatikan bahwa satu-satunya nilai tak-nol dalam ACF teoritis adalah untuk lags 1 dan 2. Autokorelasi untuk kelambatan yang lebih tinggi adalah 0 Jadi, sampel ACF dengan autokorelasi signifikan pada kelambatan 1 dan 2, namun autokorelasi yang tidak signifikan untuk kelambatan yang lebih tinggi mengindikasikan model MA (2) yang mungkin. Iid N (0,1). Koefisiennya adalah 0,5 dan 0,3. Karena ini adalah MA (2), ACF teoritis akan memiliki nilai tak-nol hanya pada kelambatan 1 dan 2. Nilai dari kedua autokorelasi tak-nol adalah sebidang ACF teoritis berikut. Seperti yang hampir selalu terjadi, sampel data tidak akan berperilaku sangat sempurna seperti teori. Kami mensimulasikan n 150 nilai sampel untuk model x t 10 w t .5 w t-1, 3 w t-2. Dimana w t iid N (0,1). Plot seri waktu dari data berikut. Seperti halnya plot seri waktu untuk data sampel MA (1), Anda tidak tahu banyak tentangnya. Contoh ACF untuk data simulasi berikut. Pola ini khas untuk situasi di mana model MA (2) mungkin berguna. Ada dua lonjakan yang signifikan secara statistik pada lag 1 dan 2 diikuti oleh nilai non-signifikan untuk kelambatan lainnya. Perhatikan bahwa karena kesalahan sampling, sampel ACF tidak sesuai dengan pola teoritisnya. ACF untuk General MA (q) Model Properti dari model MA (q) secara umum adalah bahwa ada autokorelasi tak-nol untuk q lags pertama dan autokorelasi 0 untuk semua lags gt q. Non-keunikan hubungan antara nilai 1 dan (rho1) pada MA (1) Model. Dalam model MA (1), untuk nilai 1. Timbal balik 1 1 memberikan nilai yang sama untuk Sebagai contoh, gunakan 0,5 untuk 1. Dan kemudian gunakan 1 (0.5) 2 untuk 1. Anda akan mendapatkan (rho1) 0,4 dalam kedua contoh. Untuk memenuhi batasan teoritis yang disebut invertibilitas. Kami membatasi model MA (1) untuk memiliki nilai dengan nilai absolut kurang dari 1. Pada contoh yang diberikan, 1 0,5 akan menjadi nilai parameter yang diijinkan, sedangkan 1 10,5 2 tidak akan. Keterbacaan model MA Model MA dikatakan dapat dibalikkan jika secara aljabar setara dengan model AR tak berhingga yang terkuak. Dengan konvergensi, berarti koefisien AR turun menjadi 0 saat kita bergerak mundur. Invertibilitas adalah pembatasan yang diprogram dalam perangkat lunak time series yang digunakan untuk memperkirakan koefisien model dengan persyaratan MA. Ini bukan sesuatu yang kita periksa dalam analisis data. Informasi tambahan tentang batasan invertibilitas untuk model MA (1) diberikan dalam lampiran. Catatan Teori Lanjutan Untuk model MA (q) dengan ACF tertentu, hanya ada satu model yang dapat dibalik. Kondisi yang diperlukan untuk invertibilitas adalah bahwa koefisien memiliki nilai sedemikian rupa sehingga persamaan 1- 1 y-. - q y q 0 memiliki solusi untuk y yang berada di luar lingkaran unit. Kode R untuk Contoh-Contoh Pada Contoh 1, kami menyusun model teoritis ACF dari model x t 10 w t. 7w t-1. Dan kemudian disimulasikan n 150 nilai dari model ini dan diplotkan time series sampel dan sampel ACF untuk data simulasi. Perintah R yang digunakan untuk merencanakan ACF teoritis adalah: acfma1ARMAacf (mac (0.7), lag. max10) 10 lag dari ACF untuk MA (1) dengan theta1 0.7 lags0: 10 menciptakan variabel yang disebut lags yang berkisar dari 0 sampai 10. plot (Lags, acfma1, xlimc (1,10), ylabr, typeh, ACF utama untuk MA (1) dengan theta1 0.7) abline (h0) menambahkan sumbu horizontal ke plot Perintah pertama menentukan ACF dan menyimpannya dalam objek Bernama acfma1 (pilihan nama kami). Perintah plot (perintah ke-3) cenderung tertinggal dibandingkan nilai ACF untuk lags 1 sampai 10. Parameter ylab memberi label sumbu y dan parameter utama menempatkan sebuah judul pada plot. Untuk melihat nilai numerik ACF cukup gunakan perintah acfma1. Simulasi dan plot dilakukan dengan perintah berikut. Xcarima. sim (n150, list (mac (0.7))) Simulasikan n 150 nilai dari MA (1) xxc10 menambahkan 10 untuk membuat mean 10. Simulasi default berarti 0. plot (x, typeb, mainSimulated MA (1) data) Acf (x, xlimc (1,10), mainACF untuk data sampel simulasi) Pada Contoh 2, kami merencanakan teoritis ACF dari model xt 10 wt .5 w t-1, 3 w t-2. Dan kemudian disimulasikan n 150 nilai dari model ini dan diplotkan time series sampel dan sampel ACF untuk data simulasi. Perintah R yang digunakan adalah acfma2ARMAacf (mac (0,5,0,3), lag. max10) acfma2 lags0: 10 plot (lags, acfma2, xlimc (1,10), ylabr, typeh, ACF utama untuk MA (2) dengan theta1 0.5, Theta20.3) abline (h0) xcarima. sim (n150, list (mac (0.5, 0.3))) xxc10 plot (x, typeb, seri Simulated MA (2)) acf (x, xlimc (1,10) MainACF untuk simulasi MA (2) Data) Lampiran: Bukti Sifat MA (1) Bagi siswa yang berminat, berikut adalah bukti sifat teoritis model MA (1). Varians: teks (teks teks (xt) teks (wt theta1 w) 0 teks (wt) teks (theta1w) sigma2w theta21sigma2w (1theta21) sigma2w) Bila h 1, ungkapan sebelumnya 1 w 2. Untuk h 2, ungkapan sebelumnya 0 Alasannya adalah bahwa, dengan definisi independensi wt. E (w k w j) 0 untuk setiap k j. Selanjutnya, karena meannya 0, E (w j w j) E (w j 2) w 2. Untuk seri waktu, Terapkan hasil ini untuk mendapatkan ACF yang diberikan di atas. Model MA yang dapat dibalik adalah salah satu yang dapat ditulis sebagai model AR berurutan tak terbatas yang menyatu sehingga koefisien AR bertemu ke 0 saat kita bergerak mundur jauh waktu. Nah, mendifinisikan invertibilitas model MA (1). Kita kemudian mengganti hubungan (2) untuk w t-1 dalam persamaan (1) (3) (zt wt theta1 (z-theta1w) dengan theta1z-theta2w) Pada waktu t-2. Persamaan (2) menjadi Kami kemudian mengganti hubungan (4) untuk w t-2 dalam persamaan (3) (zt wt theta1 z-theta21w wta theta1z-theta21w wt theta1z - theta12z theta31w) Jika kita melanjutkan ( Tak terbatas), kita akan mendapatkan model AR tak berhingga (zt wt theta1 z - theta21z theta31z - theta41z dots) Namun perlu dicatat bahwa jika 1 1, koefisien yang mengalikan kelambatan z akan meningkat (tak terbatas) jika kita bergerak kembali waktu. Untuk mencegah hal ini, kita membutuhkan 1 lt1. Ini adalah kondisi untuk model MA (1) yang dapat dibalik. Model MA Order Tak Terhingga Pada minggu ke 3, perhatikan bahwa model AR (1) dapat dikonversi menjadi model MA pesanan tak terbatas: (xt - mu wt phi1w phi21w titik phik1 w titik jumlah phij1w) Penjumlahan istilah white noise masa lalu ini diketahui. Sebagai representasi kausal AR (1). Dengan kata lain, x t adalah tipe MA khusus dengan jumlah tak terhingga dari istilah yang masuk ke masa lalu. Ini disebut MA tak terbatas atau MA (). Urutan terbatas MA adalah AR tak berhingga dan urutan terbatas AR adalah MA tak terbatas. Ingat di Minggu 1, kami mencatat bahwa persyaratan untuk AR stasioner (1) adalah bahwa 1 lt1. Mari menghitung Var (x t) dengan menggunakan representasi kausal. Langkah terakhir ini menggunakan fakta dasar tentang deret geometris yang membutuhkan (phi1lt1) jika rangkaiannya menyimpang. NavigationMoving Average Contoh ini mengajarkan cara menghitung moving average dari deret waktu di Excel. Rata-rata bergerak digunakan untuk memperlancar penyimpangan (puncak dan lembah) agar mudah mengenali tren. 1. Pertama, mari kita lihat rangkaian waktu kita. 2. Pada tab Data, klik Analisis Data. Catatan: cant menemukan tombol Analisis Data Klik disini untuk memuat add-in Analisis ToolPak. 3. Pilih Moving Average dan klik OK. 4. Klik pada kotak Input Range dan pilih range B2: M2. 5. Klik di kotak Interval dan ketik 6. 6. Klik pada kotak Output Range dan pilih sel B3. 8. Plot grafik nilai-nilai ini. Penjelasan: karena kita mengatur interval ke 6, rata-rata bergerak adalah rata-rata dari 5 titik data sebelumnya dan titik data saat ini. Akibatnya, puncak dan lembah dihaluskan. Grafik menunjukkan tren yang semakin meningkat. Excel tidak bisa menghitung moving average untuk 5 poin data pertama karena tidak ada cukup data point sebelumnya. 9. Ulangi langkah 2 sampai 8 untuk interval 2 dan interval 4. Kesimpulan: Semakin besar interval, semakin puncak dan lembah dihaluskan. Semakin kecil interval, semakin dekat rata-rata bergerak ke titik data aktual.