Merge pull request #1397 from cmu-delphi/prepare-contingency-historical-backfill

Add older indicators so contingency matches survey API pipeline

Author: krivard

facebook/delphiFacebook/R/contingency_indicators.R

@@ -73,13 +73,17 @@ get_aggs <- function() {
     ## Mask wearing and distancing
     "pct_wearing_mask_5d", "c_mask_often", compute_binary, jeffreys_binary,
     "pct_wearing_mask_7d", "c_mask_often_7d", compute_binary, jeffreys_binary,
+    "pct_others_masked", "c_others_masked", compute_binary, jeffreys_binary,
     "pct_others_masked_public", "c_others_masked_public", compute_binary, jeffreys_binary,
     "pct_others_distanced_public", "c_others_distanced_public", compute_binary, jeffreys_binary,
 
     "pct_cli", "hh_p_cli", compute_household_binary, jeffreys_count,
     "pct_ili", "hh_p_ili", compute_household_binary, jeffreys_count,
     "pct_hh_cmnty_cli", "hh_community_yes", compute_binary, jeffreys_binary,
+    "pct_nohh_cmnty_cli", "community_yes", compute_binary, jeffreys_binary,
     "pct_direct_contact", "c_direct_contact", compute_binary, jeffreys_binary,
+    "pct_avoid_contact", "c_avoid_contact", compute_binary, jeffreys_binary,
+    "pct_avoid_contact_7d", "c_avoid_contact_7d", compute_binary, jeffreys_binary,
 
     # symptoms
     "pct_symp_fever", "symp_fever", compute_binary, jeffreys_binary,
@@ -303,8 +307,6 @@ get_aggs <- function() {
     "pct_race_treated_fairly_healthcare", "b_race_treated_fairly_healthcare", compute_binary, jeffreys_binary,
     "pct_delayed_care_cost", "b_delayed_care_cost", compute_binary, jeffreys_binary,
 
-    "pct_worried_catch_covid", "mh_worried_catch_covid", compute_binary, jeffreys_binary,
-    
     # topics want to learn about
     "pct_want_info_covid_treatment", "i_want_info_covid_treatment", compute_binary, jeffreys_binary,
     "pct_want_info_vaccine_access", "i_want_info_vaccine_access", compute_binary, jeffreys_binary,
@@ -355,6 +357,8 @@ get_aggs <- function() {
     "pct_anxious_7d", "mh_anxious_7d", compute_binary, jeffreys_binary,
     "pct_depressed_7d", "mh_depressed_7d", compute_binary, jeffreys_binary,
     "pct_felt_isolated_7d", "mh_isolated_7d", compute_binary, jeffreys_binary,
+    # wave 11
+    "pct_worried_catch_covid", "mh_worried_catch_covid", compute_binary, jeffreys_binary,
 
     # travel outside state
     # pre-wave 10
@@ -420,7 +424,11 @@ get_aggs <- function() {
     "pct_vaccine_barrier_type_tried", "v_vaccine_barrier_type_tried", compute_binary, jeffreys_binary,
     "pct_vaccine_barrier_none_tried", "v_vaccine_barrier_none_tried", compute_binary, jeffreys_binary,
     "pct_vaccine_barrier_appointment_location_tried", "v_vaccine_barrier_appointment_location_tried", compute_binary, jeffreys_binary,
-    "pct_vaccine_barrier_other_tried", "v_vaccine_barrier_other_tried", compute_binary, jeffreys_binary
+    "pct_vaccine_barrier_other_tried", "v_vaccine_barrier_other_tried", compute_binary, jeffreys_binary,
+
+    # schooling
+    "pct_inperson_school_fulltime", "s_inperson_school_fulltime", compute_binary, jeffreys_binary,
+    "pct_inperson_school_parttime", "s_inperson_school_parttime", compute_binary, jeffreys_binary
   )
 
   aggs <- create_aggs_product(regions, groups, indicators)

facebook/delphiFacebook/R/contingency_variables.R

@@ -1040,6 +1040,36 @@ code_behaviors <- function(input_data, wave) {
   } else {
     input_data$c_direct_contact <- NA
   }
+
+  # avoid_contact
+  # Percentage of respondents that have reported having intentionally avoided
+  # contact with other people all or most of the time
+  # 1 = all of the time, 2 = most of the time, 3 = some of the time,
+  # 4 = none of the time
+  if ("C7" %in% names(input_data)) {
+    input_data$c_avoid_contact <- case_when(
+      input_data$C7 %in% c(1, 2) ~ 1,
+      input_data$C7 %in% c(3, 4) ~ 0,
+      TRUE ~ NA_real_
+    )
+  } else {
+    input_data$c_avoid_contact <- NA_real_
+  }
+
+  # avoid_contact_7d
+  # Percentage of respondents that have reported having intentionally avoided
+  # contact with other people all or most of the time in the last 7 days
+  # 1 = all of the time, 2 = most of the time, 3 = some of the time,
+  # 4 = a little of the time, 5 = none of the time
+  if ("C7a" %in% names(input_data)) {
+    input_data$c_avoid_contact_7d <- case_when(
+      input_data$C7a %in% c(1, 2) ~ 1,
+      input_data$C7a %in% c(3, 4, 5) ~ 0,
+      TRUE ~ NA_real_
+    )
+  } else {
+    input_data$c_avoid_contact_7d <- NA_real_
+  }
 
   return(input_data)
 }