# Kod 4.9. Przedział ufności dla oczekiwanego dochodu netto z wykorzystaniem nieparametrycznej estymacji funkcji gęstości
pu<- function (dane, ufn) {
n=length(dane)
zmienna=c()
for (i in 1:n) zmienna=c(zmienna,mean(sample(dane,n,replace=TRUE)))
kde = kdensity(zmienna, kernel = "gaussian", normalized=FALSE)
lewy=min(zmienna)
prawy=max(zmienna)
while (integrate(function(x) kde(x), lower=-Inf, upper=lewy)$value<((1-ufn)/2))
lewy = lewy+0.01
while (integrate(function(x) kde(x), lower=prawy,upper=Inf)$value<((1-ufn)/2))
prawy=prawy-0.01
lewy
prawy
pu
list(lewy=lewy, prawy=prawy, ufn=ufn)
}
n=25
d_s=sample(dochód,n)
pu(d_s,.90)
pu(d_s,.95)[1:2]
t.test(d_s)$conf.int