data(api)
apipop$prob<-with(apipop, 200*api00/sum(api00))
insample<-as.logical(rbinom(nrow(apipop),1,apipop$prob))
apipois<-apipop[insample,]
despois<-svydesign(id=~1, prob=~prob, pps=poisson_sampling(apipois$prob), data=apipois)
svytotal(~api00, despois)
## SE formula
sqrt(sum( (apipois$api00*weights(despois))^2*(1-apipois$prob)))
Run the code above in your browser using DataLab