## 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)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.