- Apakah ada teknik pemodelan seperti LOESS yang memungkinkan nol, satu, atau lebih diskontinuitas, di mana waktu diskontinuitas tidak diketahui apriori?
- Jika ada teknik, apakah ada implementasi yang ada di R?
r
regression
curve-fitting
change-point
loess
Jeromy Anglim
sumber
sumber
Jawaban:
Kedengarannya seperti Anda ingin melakukan beberapa perubahan changepoint diikuti oleh perataan independen dalam setiap segmen. (Deteksi bisa online atau tidak, tetapi aplikasi Anda tidak mungkin online.) Ada banyak literatur tentang ini; Pencarian internet berbuah.
Saya belum mencari-cari implementasi R (saya telah membuat kode di Mathematica beberapa waktu lalu) tetapi akan sangat menghargai referensi jika Anda menemukannya.
sumber
lakukan dengan regresi garis terputus koencker, lihat halaman 18 dari sketsa ini
http://cran.r-project.org/web/packages/quantreg/vignettes/rq.pdf
Menanggapi komentar terakhir Whuber:
Pengukur ini didefinisikan seperti ini.
, x ( i ) ≥ x ( i - 1 )x∈R ,x(i)≥x(i−1)∀i
,ei:=yi−βix(i)−β0
, z - = maks ( - z , 0 ) ,z+=max(z,0) z−=max(−z,0)
, λ ≥ 0τ∈(0,1) λ≥0
memberikan kuantil yang diinginkan (yaitu dalam contoh, τ = 0,9 ). λ mengarahkan jumlah breakpoint: untuk λ besar estimator ini menyusut menjadi tidak ada break point (sesuai dengan estimator regresi kuantil linier klasikla).τ τ=0.9 λ λ
Quantile Smoothing Splines Roger Koenker, Pin Ng, Stephen Portnoy Biometrika, Vol. 81, No. 4 (Desember, 1994), hlm. 673-680
PS: ada kertas kerja acess terbuka dengan nama yang sama oleh orang yang sama tapi itu bukan hal yang sama.
sumber
Berikut adalah beberapa metode dan paket R terkait untuk mengatasi masalah ini
Estimasi panjang gelombang wavelet dalam regresi memungkinkan untuk ketidaktepatan. Anda dapat menggunakan paket wavethresh di R.
Banyak metode berbasis pohon (tidak jauh dari gagasan wavelet) berguna ketika Anda memiliki diskrititas. Maka paket treethresh, paket pohon!
Dalam berbagai metode " kemungkinan maksimum lokal " ... antara lain: Karya Pozhel dan Spokoiny: Penimbangan bobot adaptif (pengemasan paket) Pekerjaan oleh Catherine Loader: package locfit
Saya kira setiap kernel lebih lancar dengan bandwidth yang bervariasi secara lokal , tetapi saya tidak tahu paket R untuk itu.
catatan: Saya tidak benar-benar mendapatkan apa perbedaan antara LOESS dan regresi ... apakah ini gagasan bahwa dalam alrgoritma LOESS harus "on line"?
sumber
Seharusnya dimungkinkan untuk membuat kode solusi dalam R menggunakan fungsi regresi non-linear nls, b splines (fungsi bs dalam paket spline, misalnya) dan fungsi ifelse.
sumber