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

## ----message=FALSE, warning=FALSE---------------------------------------------
library(rankingQ)
library(estimatr)
data(identity)

## -----------------------------------------------------------------------------
out_weights <- imprr_weights(
  identity,
  main_q = c("party", "religion", "gender", "race"),
  anc_correct = "anc_correct_identity"
)

out_weights$est_p_random

## -----------------------------------------------------------------------------
lm_robust(
  party ~ 1,
  data = out_weights$results,
  weights = out_weights$results$weights
) |>
  tidy()

## -----------------------------------------------------------------------------
items_df <- data.frame(
  variable = c("party", "religion", "gender", "race"),
  item = c("Party", "Religion", "Gender", "Race")
)

avg_rank(out_weights$results, items = items_df, weight = "weights", raw = FALSE)

## -----------------------------------------------------------------------------
out_boot <- imprr_weights_boot(
  identity,
  main_q = c("party", "religion", "gender", "race"),
  anc_correct = "anc_correct_identity",
  n_bootstrap = 10,
  seed = 123
)

out_boot$est_p_random
subset(out_boot$results, qoi == "average rank")

## -----------------------------------------------------------------------------
subset(out_boot$results, qoi == "pairwise ranking")