Kebingungan tentang kriging

9

Saya sedang membaca artikel wikipedia ini terkait dengan kriging. Saya tidak mengerti bagian ketika mengatakan itu

Kriging menghitung terbaik linear berisi , dari Z ( x 0 ) sehingga kriging varians dari diminimalkan dengan kondisi unbiasedness. Saya tidak mendapatkan derivasi dan juga bagaimana variansnya diminimalkan. Ada saran?Z^(x0)Z(x0)

Khususnya, saya tidak mendapatkan bagian di mana berlaku diminimalkan tunduk pada kondisi tidak memihak.

Saya pikir seharusnya begitu

E [Z '(x0) -Z (x0)] bukan E [Z' (x) -Z (x)] bukan. 'Ini setara dengan topi di artikel wiki. Juga saya tidak mengerti bagaimana kesalahan kriging diturunkan

pengguna31820
sumber
Di mana Anda terpaku pada derivasi?
whuber
Bagian di mana ia menghitung kesalahan kriging dan memaksakan kondisi ketidakberpihakan. Baik untuk mengatakan bahwa kondisi yang tidak memihak berarti ekspektasi penduga dan yang benar adalah sama. Saya telah mengedit posting untuk memasukkan rinciannya.
user31820
Saya pikir Anda benar bahwa ekspresi Wikipedia harus membaca . E[Z(x0)Z(x0)]
Whuber

Jawaban:

13

Misalkan adalah vektor yang diasumsikan memiliki distribusi multivariat dari mean yang tidak diketahui ( μ , μ , ... , μ ) dan matriks varians-kovarians yang diketahui Σ . Kami mengamati ( z 1 , z 2 , ... , z n ) dari distribusi ini dan ingin memprediksi z 0 dari informasi ini menggunakan prediktor linier yang tidak bias:(Z0,Z1,,Zn)(μ,μ,,μ)Σ(z1,z2,,zn) z0

  • Linear berarti prediksi harus berbentuk untuk koefisien λ i yang akan ditentukan. Koefisien-koefisien ini dapat bergantung paling banyak pada apa yang diketahui sebelumnya: yaitu, entri Σ .z0^=λ1z1+λ2z2++λnznλiΣ

Prediktor ini juga dapat dianggap sebagai variabel acak .Z0^=λ1Z1+λ2Z2++λnZn

  • Tidak sesuai artinya harapan sama dengan rata-rata (tidak diketahui) μ .Z0^μ

Menuliskan sesuatu memberikan beberapa informasi tentang koefisien:

μ=E[Z0^]=E[λ1Z1+λ2Z2++λnZn]=λ1E[Z1]+λ2E[Z2]++λnE[Zn]=λ1μ++λnμ=(λ1++λn)μ.

Baris kedua adalah karena linearitas harapan dan sisanya adalah aljabar sederhana. Karena prosedur ini seharusnya berfungsi terlepas dari nilai , jelas koefisien harus dijumlahkan menjadi satu. Menulis koefisien dalam notasi vektor λ = ( λ i ) , ini dapat ditulis dengan rapi 1 λ = 1 .μλ=(λi)1λ=1

Di antara set semua prediktor linier yang tidak bias seperti itu, kami mencari satu yang menyimpang sesedikit mungkin dari nilai sebenarnya , diukur dalam ruang mean square. Ini, sekali lagi, adalah perhitungan. Itu bergantung pada bilinearitas dan simetri kovarians, yang penerapannya bertanggung jawab atas penjumlahan pada baris kedua:

E[(Z0^Z0)2]=E[(λ1Z1+λ2Z2++λnZnZ0)2]=i=1nj=1nλiλjvar[Zi,Zj]2i=1nλivar[Zi,Z0]+var[Z0,Z0]=i=1nj=1nλiλjΣi,j2i=1nλiΣ0,i+Σ0,0.

1λ=1

Zx0,,xnZ(Z(x0),,Z(xn))Zi=Z(xi)n+1xin+1

σOK2Σ(Z0,,Zn)z0

  1. Rata-rata prediksi kami benar.

  2. z0σOKz0

Jauh lebih banyak yang perlu dikatakan sebelum ini dapat diterapkan pada situasi praktis seperti memperkirakan permukaan dari data tepat waktu: kita memerlukan asumsi tambahan tentang bagaimana karakteristik statistik dari proses spasial bervariasi dari satu lokasi ke lokasi lain dan dari satu realisasi ke yang lain (walaupun , dalam praktiknya, biasanya hanya satu realisasi yang akan tersedia). Tetapi penjelasan ini harus cukup untuk mengikuti bagaimana pencarian untuk Predictor Linear Tidak Berfungsi "Terbaik" ("BLUP") mengarah langsung ke sistem persamaan linear.


ΣΣΣdan memprediksi kumpulan nilai di lokasi yang tidak diketahui. Mereka membutuhkan asumsi yang sedikit lebih kuat (normalitas multivarian) untuk mencapai prestasi ini.

whuber
sumber
Ada sebuah situs web di luar sana di mana mereka berteriak menentang kriging dan sepertinya dia memiliki beberapa poin yang valid. Saya pikir paragraf terakhir Anda di sini sangat mencerahkan.
Wayne
@Wayne Ya, Anda bisa tahu apa yang saya bereaksi. Tetapi meskipun kriging telah digunakan sebagai "minyak ular" oleh para konsultan, ia telah banyak digunakan untuk itu, termasuk teori "perubahan dukungan" untuk membandingkan data yang diperoleh dari (katakanlah) sampel kecil suatu media dengan data yang diperoleh dari jauh lebih besar. bagian dari medium itu. Kriging akhirnya berada di bawah pemodelan spatio-temporal paling canggih saat ini. Ini juga merupakan cara yang berguna untuk mengevaluasi proposal alternatif: misalnya, banyak interpolator spasial linier (atau dapat linierisasi), jadi adil untuk membandingkan varians estimasi mereka dengan yang dari kriging.
whuber
1

Kriging hanyalah estimasi kuadrat terkecil untuk data spasial. Karena itu ia menyediakan penaksir linier yang tidak memihak yang meminimalkan jumlah kesalahan kuadrat. Karena tidak bias, MSE = varians estimator dan minimum.

Michael R. Chernick
sumber
Saya tidak mendapatkan bagian menghitung kesalahan kriging. Juga saya bingung dengan varian dan varian kriging. Apa perbedaannya dan apa signifikansi mereka
user31820
@whuber. Terima kasih atas penjelasannya tetapi saya tidak mendapatkan derivasi persamaan ketika Anda menghitung MSE dari nilai yang diprediksi oleh estimasi bias dan estimator sejati. Baris kedua untuk lebih spesifik dalam persamaan itu
user31820
@whuber Juga saya tidak mendapatkan bagian wiki ketika menghitung varians kriging yang mirip dengan yang ada di jawaban Anda. Mereka memiliki hasil yang sama tetapi persyaratan awal berbeda. Bagaimana bisa?
user31820