Programy R

### PRÓBA LPBZ WARSTWOWANA PO WYLOSOWANIU - ESTYMACJA ŚREDNIEJ

###############################################################
# estymacja wartosci sredniej w populacji
# na podstawie próby lpbz
# WARSTWOWANEJ PO WYLOSOWANIU
# dane:
# - pierwsza kolumna to wartości zmiennej badanej z próby
# - druga kolumna to identyfikatory przynalezności do warstw

warpowyl_sr<-function(dane,Nh_ciag) {
n<-nrow(dane)
N<-sum(Nh_ciag)

wh_ciag<-Nh_ciag/N
ys<-dane[,1]
warstwa<-dane[,2]
warstwa<-as.factor(warstwa)

nh_ciag<-c()
ysrsh_ciag<-c()
vgsh_ciag<-c()

for (i in 1:nlevels(warstwa)) {
ysh<-ys[warstwa==levels(warstwa)[i]]
nh_ciag[i]<-length(ysh)
ysrsh_ciag[i]<-mean(ysh)
vgsh_ciag[i]<-var(ysh)
}
est<-sum(wh_ciag*ysrsh_ciag)
eD2<-((N-n)*sum(wh_ciag*vgsh_ciag)/(N*n)) + ((N-n)*sum((1-wh_ciag)*vgsh_ciag)/(N*n^2))
eD<-sqrt(eD2) #wartość estymatora średniego błędu szacunku
egam<-100*eD/est #wartość estymatora średniego względnego błędu szacunku
wynik<-matrix(c(est,eD2,eD,egam),ncol=1)

########### poniższe 3 linie formatujące wynik końcowy mogą być pominięte na kolokwium
#wynik<-formatC(wynik,format="f")
nazwa<-list(c("estymator","estymator wariancji","estymator średniego błędu szacunku","estymator średniego względnego błędu szacunku (w %)"),c("wartość"))
dimnames(wynik)<-nazwa
###########
wynik
}