## ----setup, include=FALSE-----------------------------------------------------
knitr::opts_chunk$set(collapse = TRUE, comment = "#>")

## -----------------------------------------------------------------------------
library(interSAE)
x <- inter_example()
design <- inter_design(x$survey, "weight", "area", "strata", "psu")
margins <- census_margins(x$margins)
fit <- fit_inter_sae(
  unemployed ~ sex + age + residence + area,
  design, margins,
  family = "binomial",
  domain = c("area", "sex")
)
head(inter_estimates(fit))

## -----------------------------------------------------------------------------
max(margin_diagnostics(fit)$relative_error)
check_identifiability(fit)

## ----eval=requireNamespace("lpSolve", quietly=TRUE)---------------------------
# inter_sensitivity(
#   fit,
#   gamma = c(1, 1.25, 1.5, 2),
#   domain = list(area = "A1", sex = "Female")
# )

## ----eval=FALSE---------------------------------------------------------------
# boot <- inter_bootstrap(fit, R = 199, type = "multiplier", seed = 42)
# boot$summary