# Kod 4.8. Przedział ufności dla wartości oczekiwanej i współczynnika zmienności
library(boot)
n=25
d_s=sample(dochód,25)
# Przedział ufności dla wartości oczekiwanej
t.test(d_s,conf.level=0.95)$conf.int[1:2]
d_s.mean <- mean(d_s)
d_s.n <- nrow(d_s)
d_s.v <- d_s.mean^2/d_s.n
norm.ci(t0 = d_s.mean, var.t0 = d_s.v)
### Dla współczynnika zmienności
t.test(d_s,conf.level=0.95)$conf.int[1:2]
wsp.zm <- function(x, i) sd(x[i])/mean(x[i])
d_s.boot <- boot(d_s, wsp.zm, R = 999, sim = "ordinary")
d_s.boot
boot.ci(d_s.boot, conf=0.95,type = c("norm", "basic", "perc"))