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