1 2## lonely PSUs by design 3library(survey) 4data(api) 5## not certainty PSUs by fpc 6ds<-svydesign(id = ~1, weights = ~pw, strata = ~dnum, data = apiclus1) 7summary(ds) 8 9options(survey.lonely.psu="fail") 10try(svymean(~api00,ds)) 11try(svymean(~api00, as.svrepdesign(ds))) 12options(survey.lonely.psu="remove") 13svymean(~api00,ds) 14svymean(~api00, as.svrepdesign(ds)) 15options(survey.lonely.psu="certainty") 16svymean(~api00,ds) 17svymean(~api00, as.svrepdesign(ds)) 18options(survey.lonely.psu="adjust") 19svymean(~api00,ds) 20svymean(~api00, as.svrepdesign(ds)) 21options(survey.lonely.psu="average") 22svymean(~api00,ds) 23svymean(~api00, as.svrepdesign(ds)) 24 25## fpc specified 26fpc<-ifelse(apiclus1$dnum==413, 1,1000) 27ds<-svydesign(id = ~1, weights = ~pw, strata = ~dnum, data = apiclus1,fpc=fpc) 28summary(ds) 29 30options(survey.lonely.psu="fail") 31try(svymean(~api00,ds)) 32svymean(~api00, as.svrepdesign(ds)) 33options(survey.lonely.psu="remove") 34svymean(~api00,ds) 35svymean(~api00, as.svrepdesign(ds)) 36options(survey.lonely.psu="certainty") 37svymean(~api00,ds) 38svymean(~api00, as.svrepdesign(ds)) 39options(survey.lonely.psu="adjust") 40svymean(~api00,ds) 41svymean(~api00, as.svrepdesign(ds)) 42options(survey.lonely.psu="average") 43svymean(~api00,ds) 44svymean(~api00, as.svrepdesign(ds)) 45 46rs<-as.svrepdesign(ds) 47svytotal(~api00,rs) 48SE(svytotal(~api00,subset(rs, dnum==413)))==0 49 50## lonely PSUs after subsetting 51ds<-svydesign(id = ~1, weights = ~pw, strata = ~dnum, data = subset(apiclus1,dnum !=413)) 52ds1<-ds[-31,] 53summary(ds1) 54 55options(survey.lonely.psu="fail") 56svymean(~api00,ds1) 57options(survey.lonely.psu="remove") 58svymean(~api00,ds1) 59options(survey.lonely.psu="certainty") 60svymean(~api00,ds1) 61options(survey.lonely.psu="adjust") 62svymean(~api00,ds1) 63options(survey.lonely.psu="average") 64svymean(~api00,ds1) 65 66## with adjustment 67options(survey.adjust.domain.lonely=TRUE) 68ds<-svydesign(id = ~1, weights = ~pw, strata = ~dnum, data = subset(apiclus1,dnum !=413)) 69ds1<-ds[-31,] 70summary(ds1) 71 72options(survey.lonely.psu="fail") 73try(svymean(~api00,ds1)) 74options(survey.lonely.psu="remove") 75svymean(~api00,ds1) 76options(survey.lonely.psu="certainty") 77svymean(~api00,ds1) 78options(survey.lonely.psu="adjust") 79svymean(~api00,ds1) 80options(survey.lonely.psu="average") 81svymean(~api00,ds1) 82