tests/poisson.R
In bschneidr/fastsurvey: Analysis of Complex Survey Samples
## check poisson sampling
library(survey)
data(api)
set.seed(2021-7-15)
apipop$prob<-apipop$api00/1000
insample<-rbinom(nrow(apipop),1,apipop$prob)
apipois<-apipop[insample,]
des<-svydesign(id=~1, prob=~prob, pps=poisson_sampling(apipois$prob), data=apipois)
stopifnot(isTRUE(all.equal(
as.vector(SE(svytotal(~api00,design=des))),
as.vector(sqrt(sum( (apipois$api00*weights(des))^2*(1-apipois$prob))))
)))
bschneidr/fastsurvey documentation built on March 13, 2024, 11:12 a.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
## check poisson sampling
library(survey)
data(api)
set.seed(2021-7-15)
apipop$prob<-apipop$api00/1000
insample<-rbinom(nrow(apipop),1,apipop$prob)
apipois<-apipop[insample,]
des<-svydesign(id=~1, prob=~prob, pps=poisson_sampling(apipois$prob), data=apipois)
stopifnot(isTRUE(all.equal(
as.vector(SE(svytotal(~api00,design=des))),
as.vector(sqrt(sum( (apipois$api00*weights(des))^2*(1-apipois$prob))))
)))
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.