Saya membuat semivariogram di R menggunakan package gstat, variogram()
function. Saya ingin melihat apakah ada autokorelasi spasial dalam residu model saya (kelimpahan spesies sebagai fungsi habitat, melintasi situs yang berjarak beberapa km hingga 900 km terpisah, menggunakan glmm).
Unit saya berada dalam km, dan interpretasi saya adalah bahwa jangkauannya hanya lebih dari 100 km sampai autokorelasi spasial tidak lagi menjadi "masalah". Saya bertanya-tanya apakah seseorang dapat menjelaskan mengapa nugget itu tampak begitu tinggi? Apakah ini berarti bahwa bahkan di lokasi yang sama, masih ada tingkat perbedaan yang relatif tinggi? Atau, apakah variogram bergelombang ini berarti saya harus menyesuaikan jumlah jeda dan jarak jeda saya sampai saya mendapatkan bentuk yang lebih khas?
Untuk menyelidiki sedikit lebih jauh, saya juga menggunakan fungsi variog()
dalam paket geoR, dan digunakan breaks=seq(0,100,10)
, untuk mencoba dan melihat jarak yang lebih dekat (menggunakan titik dan residu model yang sama). Ini menunjukkan bahwa poin terdekat lebih berbeda, yang juga tidak masuk akal. Mungkin ini menunjukkan tidak ada autokorelasi spasial, dan model saya sudah memperhitungkan hal ini.
Saya menemukan sumber yang sangat bagus ini, "Geostats tanpa air mata" , dan di halaman 51 ada beberapa saran bagus tentang pemasangan variograms. Dengan saran ini, yang pertama saya tampaknya memiliki kisaran yang benar. Jadi ini kembali ke pertanyaan pertama - bagaimana saya menafsirkan ini?
sumber
Jawaban:
Ya, efek nugget tinggi (semivarian tinggi pada asalnya) memberi tahu ada lemahnya (atau tidak ada) ketergantungan spasial (autokorelasi) antara data sampel pada jarak kecil. Bisa jadi itu adalah kasus struktur data yang memiliki rentang yang lebih pendek daripada interval pengambilan sampel, tetapi gambar kedua tampaknya mengindikasikan bahwa itu bukan kasusnya.
sumber