Programy R

### LOSOWANIE WARSTWOWE - WARUNKOWA MINIMALIZACJA WARIANCJI

###############################################################
# dane Z PRÓBKI WSTĘPNEJ
# - pierwsza kolumna to wartości zmiennej badanej z próby
# - druga kolumna to identyfikatory przynalezności do warstw

war_min_wariancji<-function(dane,Nh_ciag,kh_ciag,K) {
n<-nrow(dane)
N<-sum(Nh_ciag)

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

vgsh_ciag<-c()

for (i in 1:nlevels(warstwa)) {
ysh<-ys[warstwa==levels(warstwa)[i]]
vgsh_ciag[i]<-var(ysh)
}

c0<-K/(sum(wh_ciag*sqrt(vgsh_ciag*kh_ciag)))
nh_ciag<-c0*wh_ciag*sqrt(vgsh_ciag/kh_ciag)
#sprawdzenie warunku 3.17
n0<-K/min(kh_ciag)
spr<-min(Nh_ciag)>n0

list("warunek spełniony"=spr,"optymalne liczebności prób"=nh_ciag)

}