# Kod 3.11. Ocena błędu standardowego ilorazu osób napływających i odpływających z powiatów województwa śląskiego w roku 2012, konstrukcja
przedziału ufności i dla danych rzeczywistych
dane=read.table('http://stat.ue.katowice.pl/powiaty.txt')
dane1=dane[dane$Rok==2012,]
dane1
attach(dane1)
sum(Napływ)/sum(Odpływ)
###
x=Napływ/Odpływ
mean(x)
n=length(x)
boot_s=c()
for (i in 1:n) {
num_s=sample(1:n,replace=TRUE)
boot_s=c(boot_s,sum(Napływ[num_s])/sum(Odpływ[num_s]))
}
mean(boot_s)
sd(boot_s)
###
library(boot)
ratio <- function(d, w) sum(d$Napływ*w) /sum(d$Odpływ *w)
boot(dane1, ratio, R = 999, stype = "w")
b1=boot(dane1, ratio, R = 999, stype = "w")
boot.ci(b1, conf = c(0.90, 0.95), type = c("perc", "norm"))