library(survey)
data(api)
srs_design <- svydesign(id=~1, fpc=~fpc, data=apisrs)
srs_design
svytotal(~enroll, srs_design)
svymean(~enroll, srs_design)
nofpc <- svydesign(id=~1, weights=~pw, data=apisrs)
svytotal(~enroll, nofpc)
strat_design <- svydesign(id=~1, strata=~stype, fpc=~fpc, data=apistrat)
svytotal(~enroll, strat_design)
svytotal(~stype, strat_design)
svytotal(~api00+api99+stype, strat_design)
svyquantile(~api00,design=strat_design,quantiles=c(0.25,0.5,0.75),ci=TRUE)
table(apisrs$sch.wide,apisrs$awards)
tab <- svytotal(~interaction(sch.wide,awards), strat_design)
tab
ftab<-ftable(tab,rownames=list(Obiettivo=c("No","Yes"),Premio=c("No","Yes")))
ftab
emerg_high<- subset(strat_design,emer>20)
emerg_low <- subset(strat_design,emer==0)
svymean(~enroll,emerg_high)
svymean(~enroll,emerg_low)
svyby(~api00+api99,~stype,strat_design,svytotal)
svyby(~api00+api99,~awards,strat_design,svytotal)
clus2_design <- svydesign(id=~dnum+snum, fpc=~fpc1+fpc2, data=apiclus2)
svytotal(~api00+api99, clus2_design)

data(election)
election$votes<-with(election, Bush+Kerry+Nader)
election$p<-40*election$votes/sum(election$votes)
library(sampling)
insample<-UPtille(election$p)
ppsample <- election[insample==1,] 
ppsample$wt<-1/ppsample$p
pps_design<-svydesign(id=~1, weight=~wt, data=ppsample)
svytotal(~Bush+Kerry+Nader,pps_design,deff=TRUE)