R Under development (unstable) (2024-03-07 r86063 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.

> ## test use of replicates in svyby, svycontrast
> 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

> 
> data(api)
> dclus1<-svydesign(id=~dnum, weights=~pw, data=apiclus1, fpc=~fpc)
> rclus1<-as.svrepdesign(dclus1)
> 
> meanlogs_without<-svyby(~log(enroll),~stype,svymean, design=rclus1,covmat=TRUE)
> c_without<-svycontrast(meanlogs_without, quote(exp(E-H)))
> vcov(c_without)
           contrast
contrast 0.06862947
> 
> meanlogs_with<-svyby(~log(enroll),~stype,svymean, design=rclus1,covmat=TRUE,return.replicates=TRUE)
> 
> c_with<-svycontrast(meanlogs_with, quote(exp(E-H)))
> 
> v_with<- vcov(rclus1, c_with$replicates)
> 
> r<- attr(meanlogs_with, "replicates")
> vr_with<-vcov(rclus1,exp(r[,1]-r[,2]))
> 
> stopifnot(all.equal(as.numeric(v_with),as.numeric(vr_with)))
> stopifnot(all.equal(as.numeric(v_with),as.numeric(vcov(c_with))))
> 
> proc.time()
   user  system elapsed 
   1.09    0.10    1.18