R Under development (unstable) (2024-03-08 r86070 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2024 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

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.

> ## simulated data with three-stage sample at phase 1, SRS at phase 2
> ## motivated by dietary biomarker substudy in HCHS
> 
> library(survey)
Loading required package: grid
Loading required package: Matrix
Loading required package: survival

Attaching package: 'survey'

The following object is masked from 'package:graphics':

    dotchart

> 
> load("simdata1.RData")
> 
> twophase.full = twophase(id=list(~block+house+ind,~1),
+                          strata=list(~strat,NULL),
+                          probs=list(~P.block+P.house+P.ind,NULL),
+                          subset=~phase2,
+                          data=simdata1,method='full')
> 
> twophase.approx = twophase(id=list(~block+house+ind,~1),
+                          strata=list(~strat,NULL),
+                          probs=list(~P.block+P.house+P.ind,NULL),
+                          subset=~phase2,
+                          data=simdata1,method='approx')
> 
> twophase.rep = twophase(id=list(~block,~1),
+                          strata=list(~strat,NULL),
+                          probs=list(~I(P.block*P.house*P.ind),NULL),
+                          subset=~phase2,
+                          data=simdata1,method='full')
> 
> 
> twophase.repapprox = twophase(id=list(~block,~1),
+                          strata=list(~strat,NULL),
+                          probs=list(~I(P.block*P.house*P.ind),NULL),
+                          subset=~phase2,
+                          data=simdata1,method='approx')
> 
> 
> svymean(~age, twophase.full)
      mean     SE
age 41.116 0.8114
> svymean(~age, twophase.approx)
      mean     SE
age 41.116 0.6693
> svymean(~age, twophase.rep)
      mean    SE
age 41.116 0.784
> svymean(~age, twophase.repapprox)
      mean     SE
age 41.116 0.6592
> 
> proc.time()
   user  system elapsed 
   1.67    0.17    1.82