tests/testSUMMER/SUMMER.R

## CRAN won't allow this test unless INLA is a listed dependency
## which is impossible, so the test is hidden away here.

library(survey)

data("DemoData2",package="SUMMER")
data("DemoMap2", package="SUMMER")

if(require("INLA",quietly=TRUE)){
INLA::inla.setOption(num.threads="1:1")
 library(survey)
 des0 <- svydesign(ids = ~clustid+id, strata = ~strata,
                  weights = ~weights, data = DemoData2, nest = TRUE)
 Xmat <- aggregate(age~region, data = DemoData2, FUN = mean)

 cts.cov.res <- svysmoothArea(tobacco.use ~ age, 
                          domain = ~region,
                          design = des0,
                          adj.mat = DemoMap2$Amat, 
                          X.domain = Xmat,
                          pc.u = 1,
                          pc.alpha = 0.01,
                          pc.u.phi = 0.5,
                          pc.alpha.phi = 2/3)
print(cts.cov.res)
plot(cts.cov.res)
summary(cts.cov.res)
}
bschneidr/fastsurvey documentation built on March 13, 2024, 11:12 a.m.