Ben böyle görünüyor sıfırlar bir sürü veri kümesi var:
set.seed(1)
x <- c(rlnorm(100),rep(0,50))
hist(x,probability=TRUE,breaks = 25)
Yoğunluğu için bir çizgi çizmek istiyorum, ancak density()
fonksiyon x'in negatif değerlerini hesaplayan hareketli bir pencere kullanıyor.
lines(density(x), col = 'grey')
Bir density(... from, to)
argüman var, ancak bunlar sadece hesaplamayı kesiyor gibi görünüyor, pencereyi değiştirmiyor, böylece 0'daki yoğunluk, aşağıdaki çizimde görülebileceği gibi verilerle tutarlı olacak şekilde:
lines(density(x, from = 0), col = 'black')
(enterpolasyon değiştiyse, siyah çizginin 0'da gri çizgiden daha yüksek yoğunluğa sahip olmasını beklerdim)
Bu fonksiyona, yoğunluğun sıfırda daha iyi hesaplanmasını sağlayacak alternatifler var mı?