R version 3.1.0 (2014-04-10) -- "Spring Dance"
Copyright (C) 2014 The R Foundation for Statistical Computing
Platform: x86_64-apple-darwin13.1.0 (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> ## from Takahiro Tsuchiya
> library(survey)
Attaching package: 'survey'
The following object is masked from 'package:graphics':
dotchart
> kigyo<-read.table(tmp<-textConnection(" obs uriage srs.w pps.w
+ 1 1 15 100 20
+ 2 2 143 100 200
+ 3 3 21 100 11
+ 4 4 51 100 25
+ 5 5 337 100 550
+ 6 6 50 100 30
+ 7 7 274 100 250
+ 8 8 145 100 100
+ 9 9 15 100 10
+ 10 10 86 100 55
+ ",open="r"),header=TRUE)
> close(tmp)
> des.srs <- svydesign(ids=~1, weights=~srs.w, data=kigyo)
> (res.srs <- svymean(~uriage, des.srs, deff=TRUE))
mean SE DEff
uriage 113.700 35.626 1.0101
> (SE(res.srs)^2) / ((1-10/1000) * coef(svyvar(~uriage, des.srs)) / 10)
uriage
uriage 1.010101
>
> (tres.srs <- svytotal(~uriage, des.srs, deff=TRUE))
total SE DEff
uriage 113700 35626 1.0101
> (SE(tres.srs)^2) / (1000^2 * (1-10/1000) * coef(svyvar(~uriage, des.srs)) / 10)
uriage
uriage 1.010101
>
>
> des.pps <- svydesign(ids=~1, weights=~pps.w, data=kigyo)
> (res.pps <- svymean(~uriage, des.pps, deff='replace'))
mean SE DEff
uriage 243.914 48.752 1.9741
> (SE(res.pps)^2) / (coef(svyvar(~uriage, des.pps)) / 10)
uriage
uriage 1.974067
> (tres.pps <- svytotal(~uriage, des.pps, deff='replace'))
total SE DEff
uriage 305136 184965 18.157
> (N.hat <- sum(weights(des.pps)))
[1] 1251
> (SE(tres.pps)^2) / (N.hat^2 * coef(svyvar(~uriage, des.pps)) / 10)
uriage
uriage 18.15669
>
> proc.time()
user system elapsed
0.163 0.021 0.189
Generated by dwww version 1.15 on Sat May 18 13:13:02 CEST 2024.