
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 
