Saya punya beberapa data yang jelas terpotong di sebelah kiri. Saya ingin mencocokkannya dengan estimasi kepadatan yang akan menanganinya dalam beberapa cara alih-alih mencoba memuluskannya.
Metode apa yang diketahui (seperti biasa, dalam R) yang dapat mengatasi ini?
Kode sampel:
set.seed(1341)
x <- c(runif(30, 0, 0.01), rnorm(100,3))
hist(x, br = 10, freq = F)
lines(density(x), col = 3, lwd = 3)
Terima kasih :)
r
pdf
histogram
kernel-smoothing
Tal Galili
sumber
sumber
Jawaban:
Paket logspline untuk R memiliki fungsi oldlogspline yang akan memperkirakan kepadatan menggunakan campuran data yang diamati dan disensor.
sumber
Fungsi kerapatan juga memiliki
from
parameter untuk menunjukkan sisi paling kiri "dari grid di mana kerapatan diperkirakan". Melanjutkan dari contoh di atas:Namun, seperti yang Anda lihat, distribusi ini persis sama tanpa
from
parameter seperti di atas. Itu baru dimulai dari 0, itu saja.sumber