library(sampling)
data(swissmunicipalities)
srs <- srswor(200,nrow(swissmunicipalities))
dat_srs <- swissmunicipalities[srs == 1,]
nrow(dat_srs)
pik <- inclusionprobabilities(swissmunicipalities$Pop020,200)
sum(pik)
pps <- UPtille(pik)
dat_pps <- swissmunicipalities[pps == 1,]
nrow(dat_pps)
numpop <- table(swissmunicipalities[,2])
numpop <- round((numpop/sum(numpop))*200)
stratif <- strata(swissmunicipalities,stratanames=c("REG"),size=numpop,method="srswor")
data_str <- getdata(swissmunicipalities, stratif)
pik <- rep(200/nrow(swissmunicipalities),nrow(swissmunicipalities))
x <- as.matrix(swissmunicipalities[,c(13,14,15,16)])
bal <- samplecube(x,pik,order=1,comment=TRUE)
dat_bal <- swissmunicipalities[bal == 1,]
nrow(dat_bal)

library(survey)
dstratif <- svydesign(id=~1,strata=~REG, probs=~stratif[,3], data= data_str, fpc=~stratif[,3])
totest <- svytotal(~~H00P01+H00P02+H00P03+H00P04, dstratif, deff=TRUE)
totest
colSums(swissmunicipalities[,c(18,19,20,21)])
cv(totest)
coef(totest)
totpop <- c(`(Intercept)`=nrow(swissmunicipalities), Pop020=sum(swissmunicipalities[,13]), Pop2040=sum(swissmunicipalities[,14]), Pop4065=sum(swissmunicipalities[,15]))
dstratifg <- calibrate(dstratif,~Pop020+Pop2040+Pop4065, totpop)
totestg <- svytotal(~H00P01+H00P02+H00P03+H00P04, dstratifg, deff=TRUE)
totestg
cv(totestg)
coef(totestg)
diff <- weights(dstratifg)/weights(dstratif)
summary(diff)
dstratifg <- calibrate(dstratif,~Pop020+Pop2040+Pop4065, calfun="logit",bounds=c(0.5,2),maxit=2,epsilon=1e-10,totpop)
totestg <- svytotal(~H00P01+H00P02+H00P03+H00P04, dstratifg, deff=TRUE)
totestg
cv(totestg)
coef(totestg)
diff <- weights(dstratifg)/weights(dstratif)
summary(diff)