Saya memiliki dataset 20-tahun dari jumlah tahunan kelimpahan spesies untuk satu set poligon (~ 200 poligon kontinu berbentuk tidak teratur). Saya telah menggunakan analisis regresi untuk menyimpulkan tren (perubahan dalam hitungan per tahun) untuk setiap poligon, serta agregasi data poligon berdasarkan batas manajemen.
Saya yakin bahwa ada autokorelasi spasial dalam data, yang pasti akan berdampak pada analisis regresi untuk data agregat. Pertanyaan saya adalah - bagaimana cara menjalankan tes SAC untuk data deret waktu? Apakah saya perlu melihat SAC residu dari regresi saya untuk setiap tahun (global Moran's I)? Atau bisakah saya menjalankan satu tes dengan semua tahun?
Setelah saya uji ya ada SAC, adakah yang mudah untuk mengatasinya? Latar belakang statistik saya minimal dan semua yang saya baca tentang pemodelan spatio-temporal terdengar sangat kompleks. Saya tahu bahwa R memiliki fungsi autocovariate berbobot jarak - apakah ini mudah digunakan?
Saya benar-benar bingung bagaimana cara menilai / menambahkan SAC untuk masalah ini dan akan sangat menghargai setiap saran, tautan, atau referensi. Terima kasih sebelumnya!
Jawaban:
Menurut makalah ini , OLS konsisten dengan adanya autokorelasi spasial, tetapi kesalahan standar tidak benar dan perlu disesuaikan. Solomon Hsiang menyediakan stata dan kode matlab untuk melakukannya. Sayangnya saya tidak terbiasa dengan kode R untuk ini.
Tentu saja ada pendekatan lain untuk masalah semacam ini dalam statistik spasial yang secara eksplisit memodelkan proses spasial. Yang ini hanya mengembang kesalahan standar.
Sayangnya para ahli ekonometrik teoretis tampaknya senang mengaburkan. Kertas yang ditautkan sangat sulit dibaca. Pada dasarnya yang dikatakannya adalah menjalankan regresi apa pun yang Anda inginkan, dan kemudian perbaiki kesalahan standar nanti yaitu: menggunakan kode dari Hsiang. Ruang tidak muncul sampai Anda mencoba memperkirakan varians dari estimator Anda. Secara intuitif, jika semua perbedaannya berdekatan, Anda kurang yakin bahwa perkiraan Anda bukan hanya peninggalan dari beberapa kejutan spasial yang tidak teramati.
Perhatikan bahwa Anda perlu menentukan bandwidth kernel di mana Anda pikir proses spasial mungkin beroperasi.
Jawaban ini pada dasarnya adalah salinan / tempel pengulangan dari jawaban serupa yang saya buat di sini
sumber
Jika masalah ini adalah kesalahan autokorelasi, , maka OLS konsisten tetapi tidak efisien, seperti kata ACD. Ini seperti korelasi serial dalam ekonometrik deret waktu.y=Xβ+u,u=ρWu+ϵ
Tetapi jika ada autodependensi spasial (juga disebut autokorelasi, membingungkan), , maka OLS tidak konsisten. Itu sama dengan variabel bias yang hilang. Jika Anda memiliki kedua masalah tersebut, Anda perlu menggunakan Model Durbin Spasial, .y=ρWy+Xβ+ϵ y=ρWy+Xβ+WXλ+ϵ
The spdep paket untuk R mengandung banyak fungsi yang menghitung bobot spasial matriks, memperkirakan regresi spasial, dan melakukan hal-hal lain. Saya punya banyak pengalaman dengan
lagsarlm
fungsi - fungsi, tetapi lihat dalam dokumentasi paket bahwa adasacsarlm
fungsi yang tampaknya lebih dari apa yang Anda cari.Sejauh aspek temporal dari masalah Anda, asumsi yang Anda buat tentang ketergantungan akan sangat menentukan spesifikasi model Anda. Apakah area Anda berinteraksi satu sama lain secara langsung? Contohnya adalah pasar perdagangan atau perumahan; ekspor dari satu negara sangat bergantung pada impor di negara lain, dan harga jual rumah yang baru dibeli merupakan kontribusi yang sangat penting bagi harga jual rumah di sekitarnya. Dalam hal ini, masuk akal untuk menentukan bobot Anda matriks untuk mengakomodasi ketergantungan ini. Izinkan rumah yang dibeli di waktu menjadi "tetangga" dengan rumah di waktu , tetapi tidak dengan rumah di waktuW t t−1 t+1
sumber