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

## -----------------------------------------------------------------------------
# library(allofus)
# library(tidyverse)
# 
# svy_vars <- aou_survey(questions = c(1585838, 1586135), question_output = c("gender", "birthplace"))

## ----echo = FALSE, eval = TRUE------------------------------------------------
cli::cli_warn(c("No cohort provided.", ">" = "Pulling survey data for entire All of Us cohort."))

## -----------------------------------------------------------------------------
# count(svy_vars, gender, birthplace)

## -----------------------------------------------------------------------------
# cohort <- svy_vars %>%
#   filter(gender == "Woman", birthplace == "Other")
# 
# cohort

## -----------------------------------------------------------------------------
# t2dm <- aou_concept_set(cohort,
#   concepts = c(201826, 4193704),
#   domains = "condition", output = "indicator",
#   concept_set_name = "t2dm"
# )

## -----------------------------------------------------------------------------
# metformin <- aou_concept_set(cohort,
#   concepts = c(40164929, 40164897),
#   domains = "drug", output = "count",
#   start_date = "gender_date", concept_set_name = "metformin"
# )

## -----------------------------------------------------------------------------
# a1c <- aou_concept_set(cohort,
#   concepts = c(3004410, 3005673),
#   domains = "measurement", output = "all", start_date = "gender_date"
# )

## -----------------------------------------------------------------------------
# a1c

## -----------------------------------------------------------------------------
# t2dm_self <- aou_survey(cohort, questions = 43529932, question_output = "t2dm_survey")

## -----------------------------------------------------------------------------
# combined_data <- reduce(list(cohort, t2dm, metformin, a1c, t2dm_self),
#   aou_join,
#   type = "left"
# )

## -----------------------------------------------------------------------------
# final_data <- combined_data %>%
#   group_by(person_id, birthplace, t2dm, metformin, t2dm_survey) %>%
#   summarize(max_a1c = max(value_as_number, na.rm = TRUE), .groups = "drop")

## -----------------------------------------------------------------------------
# final_data