1 Overview

This report reproduces and extends the analyses in portugal.R using the uploaded Qualtrics CSV file.

The CSV appears to include a first row containing question text, followed by respondent data. The report therefore stores the first row as a question-label lookup and then removes it from the analytic data.

Where possible, model names are matched to the corresponding question labels in the CSV. For example, model Q32_1_bin is labelled with the text from Q32_1.

`%||%` <- function(x, y) {
  if (is.null(x) || length(x) == 0 || is.na(x) || x == "") y else x
}

clean_question_text <- function(x) {
  x %>%
    as.character() %>%
    str_replace_all("\\s+", " ") %>%
    str_trim()
}

strip_derived_suffix <- function(x) {
  x %>%
    str_remove("_bin$") %>%
    str_remove("_ord_num$") %>%
    str_remove("_ord$")
}

p_stars <- function(p) {
  case_when(
    is.na(p) ~ "",
    p < .001 ~ "***",
    p < .01 ~ "**",
    p < .05 ~ "*",
    p < .1 ~ "†",
    TRUE ~ ""
  )
}

safe_aov_tidy <- function(data, outcome, predictor) {
  f <- as.formula(paste(outcome, "~", predictor))
  out <- tryCatch(
    aov(f, data = data) %>% broom::tidy(),
    error = function(e) tibble(term = predictor, df = NA_real_, statistic = NA_real_, p.value = NA_real_)
  )

  out %>%
    mutate(outcome = outcome, predictor = predictor) %>%
    dplyr::select(outcome, predictor, term, df, statistic, p.value)
}

safe_glm <- function(formula, data) {
  tryCatch(
    glm(formula, data = data, family = binomial(link = "logit"), na.action = na.omit),
    error = function(e) structure(list(error = conditionMessage(e)), class = "model_error")
  )
}

safe_polr <- function(formula, data) {
  tryCatch(
    MASS::polr(formula, data = data, method = "logistic", Hess = TRUE, na.action = na.omit),
    error = function(e) structure(list(error = conditionMessage(e)), class = "model_error")
  )
}

is_model_error <- function(x) inherits(x, "model_error")

safe_tidy_model <- function(model, outcome, exponentiate = TRUE) {
  empty_result <- function(error_message = NA_character_) {
    tibble(
      outcome = outcome,
      term = NA_character_,
      estimate = NA_real_,
      std.error = NA_real_,
      statistic = NA_real_,
      p.value = NA_real_,
      conf.low = NA_real_,
      conf.high = NA_real_,
      model_error = error_message
    )
  }

  if (is_model_error(model)) {
    return(empty_result(model$error))
  }

  tidy_result <- tryCatch(
    broom::tidy(model, conf.int = TRUE, exponentiate = exponentiate),
    error = function(e) empty_result(conditionMessage(e))
  )

  # `broom::tidy()` does not always return p-values for MASS::polr models.
  # When a z/t statistic is available, approximate a two-sided p-value.
  if (!"p.value" %in% names(tidy_result)) {
    tidy_result$p.value <- if ("statistic" %in% names(tidy_result)) {
      2 * pnorm(abs(tidy_result$statistic), lower.tail = FALSE)
    } else {
      NA_real_
    }
  }

  required_numeric_cols <- c("estimate", "std.error", "statistic", "p.value", "conf.low", "conf.high")
  for (col in required_numeric_cols) {
    if (!col %in% names(tidy_result)) {
      tidy_result[[col]] <- NA_real_
    }
  }

  if (!"term" %in% names(tidy_result)) {
    tidy_result$term <- NA_character_
  }

  if (!"model_error" %in% names(tidy_result)) {
    tidy_result$model_error <- NA_character_
  }

  tidy_result %>%
    mutate(outcome = outcome) %>%
    dplyr::select(outcome, term, estimate, std.error, statistic, p.value, conf.low, conf.high, model_error)
}

2 Load data

raw_data <- readr::read_csv(
  params$data_file,
  col_types = cols(.default = col_character()),
  show_col_types = FALSE
)

question_labels <- raw_data %>%
  slice(1) %>%
  pivot_longer(everything(), names_to = "variable", values_to = "question") %>%
  mutate(question = clean_question_text(question))

question_lookup <- setNames(question_labels$question, question_labels$variable)

qlabel <- function(variable) {
  base_variable <- strip_derived_suffix(variable)
  question_lookup[[base_variable]] %||% variable
}

dataframe <- raw_data %>%
  slice(-1) %>%
  mutate(across(everything(), ~ na_if(.x, "")))

n_respondents <- nrow(dataframe)
n_variables <- ncol(dataframe)

tibble(
  respondents = n_respondents,
  variables = n_variables
) %>%
  knitr::kable(caption = "Analytic dataset dimensions after removing the Qualtrics question-text row.")
Analytic dataset dimensions after removing the Qualtrics question-text row.
respondents variables
1221 191

3 Variable preparation

The original script expects some multi-select fields to be exported as separate columns. In this CSV, Q34 is stored as one comma-separated multi-select field. Because one of the response options itself contains a comma, this report avoids splitting on commas and instead searches for each known response option using fixed-string matching.

q34_options <- c(
  "I cannot leave work early to avoid travelling",
  "I cannot work from home to avoid travelling",
  "I cannot avoid doing outdoor tasks or physical exertion during extreme heat because of work responsibilities",
  "I cannot avoid working in an exposed area during extreme heat or thunderstorms",
  "Difficulty leaving my home in the event of a flood or wildfire warning because I depend on others for transport",
  "I cannot avoid travelling because of caring responsibilities for one or more vulnerable people",
  "Difficulty finding a cool place to stay during extreme heat",
  "Difficulty finding sheltered place to stay during heavy rain",
  "I would not feel comfortable cancelling trips or travel to visit friends or family if avoiding unnecessary travel due to bad weather was recommended",
  "It is difficult for me or other household members to adopt the suggested measures because of mobility problems",
  "I cannot or do not like to open windows at night to let cool air in",
  "I do not have the physical strength to secure and protect all objects and equipment on my property, and have no one to help me do so",
  "Other response provided"
)

q34_lookup <- setNames(q34_options, paste0("Q34_", seq_along(q34_options)))
question_lookup <- c(question_lookup, q34_lookup)

for (i in seq_along(q34_options)) {
  new_var <- paste0("Q34_", i)

  if (i < 13) {
    dataframe[[new_var]] <- if_else(
      !is.na(dataframe$Q34) & str_detect(dataframe$Q34, fixed(q34_options[[i]])),
      q34_options[[i]],
      NA_character_
    )
  } else {
    dataframe[[new_var]] <- if_else(
      "Q34_13_TEXT" %in% names(dataframe) & !is.na(dataframe$Q34_13_TEXT),
      dataframe$Q34_13_TEXT,
      NA_character_
    )
  }
}

q26_cols <- paste0("Q26_", 1:13)
q31_cols <- paste0("Q31_", 1:7)
q32_cols <- paste0("Q32_", 1:4)
q34_cols <- paste0("Q34_", 1:13)

present_cols <- function(cols) cols[cols %in% names(dataframe)]

outcome_vars <- present_cols(c(q31_cols, q32_cols, q34_cols, q26_cols))

# Create one combined carer variable from any Q63-style columns in the export.
# This does not attempt to preserve Q63 labels; the final variable is simply called `carer`.
q63_cols <- grep("^Q63($|[._])", names(dataframe), value = TRUE)

if (length(q63_cols) > 0) {
  q63_response_matrix <- dataframe %>%
    transmute(
      across(
        all_of(q63_cols),
        ~ !is.na(.x) &
          str_trim(.x) != "" &
          !str_detect(.x, regex("none|no caring|not applicable|prefer not", ignore_case = TRUE))
      )
    )

  dataframe$carer <- if_else(
    rowSums(as.data.frame(q63_response_matrix), na.rm = TRUE) > 0,
    1L,
    0L
  )
} else {
  dataframe$carer <- NA_integer_
}

recoded_data <- dataframe %>%
  mutate(
    gender = case_when(
      Q3 == "Male" ~ "Male",
      Q3 == "Female" ~ "Female",
      TRUE ~ NA_character_
    ),
    gender = factor(gender, levels = c("Male", "Female")),

    educ = case_when(
      Q65 %in% c("None", "Basic education, 1st, 2nd, or 3rd cycle") ~ "Primary",
      Q65 %in% c(
        "Secondary education",
        "Post-secondary non-higher education (Technological Specialisation Courses)",
        "Post-secondary non-tertiary education (Technological Specialization Courses)"
      ) ~ "Secondary",
      Q65 %in% c(
        "Bachelor's degree or equivalent", "Bachelor’s degree or equivalent",
        "Master's degree", "Master’s degree", "Doctorate"
      ) ~ "Tertiary",
      Q65 %in% c("Other", "Prefer not to answer") ~ NA_character_,
      TRUE ~ NA_character_
    ),
    educ = factor(educ, levels = c("Primary", "Secondary", "Tertiary")),

    income = case_when(
      Q64 %in% c("870", "992", "1070", "1136", "1187") ~ "Low",
      Q64 %in% c("1787", "2078", "2432", "3233") ~ "Middle",
      Q64 %in% c("5547", "20221", "Greater than 20221") ~ "High",
      Q64 == "Prefer not to answer" ~ "Prefer not to answer",
      TRUE ~ NA_character_
    ),
    income = factor(income, levels = c("Low", "Middle", "High", "Prefer not to answer")),

    urbanicity = case_when(
      Q61 == "Urban" ~ "Urban",
      Q61 == "Suburban" ~ "Suburban",
      Q61 == "Semi-rural" ~ "Semi-rural",
      Q61 == "Rural" ~ "Rural",
      Q61 == "Other/prefer not to answer" ~ NA_character_,
      TRUE ~ NA_character_
    ),
    urbanicity = factor(urbanicity, levels = c("Urban", "Suburban", "Semi-rural", "Rural")),

    carer = factor(carer, levels = c(0, 1), labels = c("Non-carer", "Carer")),

    age = readr::parse_number(as.character(Q2)),
    age = if_else(age %in% c(5, 28399), NA_real_, age),

    Exp1Type = factor(Exp1Type),
    Exp2Type = factor(Exp2Type)
  ) %>%
  mutate(
    across(
      all_of(present_cols(q26_cols)),
      ~ case_when(
        .x %in% c("Much less frequent") ~ 1,
        .x %in% c("Slightly less frequent", "Somewhat less frequent") ~ 2,
        .x %in% c("Neither more nor less frequent") ~ 3,
        .x %in% c("Slightly more frequent", "Somewhat more frequent") ~ 4,
        .x %in% c("Much more frequent") ~ 5,
        TRUE ~ NA_real_
      ),
      .names = "{.col}_ord"
    )
  ) %>%
  mutate(
    across(
      ends_with("_ord"),
      ~ ordered(
        .x,
        levels = c(1, 2, 3, 4, 5),
        labels = c(
          "Much less frequent",
          "Somewhat less frequent",
          "Neither more nor less frequent",
          "Somewhat more frequent",
          "Much more frequent"
        )
      )
    ),
    across(
      ends_with("_ord"),
      ~ as.numeric(.x),
      .names = "{.col}_num"
    ),
    across(
      all_of(present_cols(q31_cols)),
      ~ case_when(
        .x == "Yes" ~ 1,
        .x %in% c("No", "I am not sure", "Not sure", "Not applicable") ~ 0,
        TRUE ~ NA_real_
      ),
      .names = "{.col}_bin"
    ),
    across(
      all_of(present_cols(q32_cols)),
      ~ case_when(
        .x == "Yes" ~ 1,
        .x %in% c("No", "Not applicable") ~ 0,
        TRUE ~ NA_real_
      ),
      .names = "{.col}_bin"
    ),
    across(
      all_of(present_cols(q34_cols)),
      ~ case_when(
        is.na(.x) ~ 0,
        .x == "" ~ 0,
        TRUE ~ 1
      ),
      .names = "{.col}_bin"
    )
  )

dataframe <- recoded_data

analysis_vars <- c("age", "gender", "educ", "income", "carer", "urbanicity", "Exp1Type", "Exp2Type")

dataframe %>%
  summarise(across(all_of(analysis_vars), ~ sum(is.na(.x)))) %>%
  pivot_longer(everything(), names_to = "variable", values_to = "missing_n") %>%
  knitr::kable(caption = "Missingness in analysis variables.")
Missingness in analysis variables.
variable missing_n
age 10
gender 3
educ 10
income 30
carer 0
urbanicity 7
Exp1Type 0
Exp2Type 0

4 Sample descriptives

descriptive_table <- bind_rows(
  dataframe %>% count(variable = "gender", value = as.character(gender), name = "n"),
  dataframe %>% count(variable = "educ", value = as.character(educ), name = "n"),
  dataframe %>% count(variable = "income", value = as.character(income), name = "n"),
  dataframe %>% count(variable = "carer", value = as.character(carer), name = "n"),
  dataframe %>% count(variable = "urbanicity", value = as.character(urbanicity), name = "n")
) %>%
  group_by(variable) %>%
  mutate(percent = round(100 * n / sum(n), 1)) %>%
  ungroup()

descriptive_table %>%
  arrange(variable, desc(n)) %>%
  knitr::kable(caption = "Categorical sample descriptives.")
Categorical sample descriptives.
variable value n percent
carer Carer 686 56.2
carer Non-carer 535 43.8
educ Secondary 608 49.8
educ Tertiary 505 41.4
educ Primary 98 8.0
educ NA 10 0.8
gender Female 613 50.2
gender Male 605 49.5
gender NA 3 0.2
income Middle 700 57.3
income Low 358 29.3
income Prefer not to answer 86 7.0
income High 47 3.8
income NA 30 2.5
urbanicity Urban 606 49.6
urbanicity Suburban 346 28.3
urbanicity Rural 133 10.9
urbanicity Semi-rural 129 10.6
urbanicity NA 7 0.6 
age_summary <- dataframe %>%
  summarise(
    n = sum(!is.na(age)),
    mean = mean(age, na.rm = TRUE),
    sd = sd(age, na.rm = TRUE),
    median = median(age, na.rm = TRUE),
    min = min(age, na.rm = TRUE),
    max = max(age, na.rm = TRUE)
  )

age_summary %>%
  mutate(across(where(is.numeric), ~ round(.x, 2))) %>%
  knitr::kable(caption = "Age summary.")
Age summary.
n mean sd median min max
1211 45.34 15.28 45 1 88 
descriptive_table %>%
  filter(!is.na(value)) %>%
  ggplot(aes(x = fct_reorder(value, percent), y = percent)) +
  geom_col() +
  coord_flip() +
  facet_wrap(~ variable, scales = "free_y") +
  scale_y_continuous(labels = label_percent(scale = 1)) +
  labs(x = NULL, y = "Percent", title = "Sample characteristics") +
  theme_minimal()

5 Univariate outcome results

outcome_choice_percents <- dataframe %>%
  dplyr::select(dplyr::all_of(outcome_vars)) %>%
  pivot_longer(cols = everything(), names_to = "outcome", values_to = "response") %>%
  mutate(
    response = if_else(is.na(response) | response == "", "Missing / blank", as.character(response)),
    question = map_chr(outcome, qlabel)
  ) %>%
  count(outcome, question, response, name = "n") %>%
  group_by(outcome) %>%
  mutate(percent = round(100 * n / sum(n), 1)) %>%
  ungroup() %>%
  arrange(outcome, desc(percent))

outcome_choice_percents %>%
  dplyr::select(outcome, question, response, n, percent) %>%
  knitr::kable(caption = "Response distributions for outcome variables.")
Response distributions for outcome variables.
outcome question response n percent
Q26_1 To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme heat Somewhat more frequent 609 49.9
Q26_1 To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme heat Much more frequent 392 32.1
Q26_1 To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme heat Neither more nor less frequent 138 11.3
Q26_1 To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme heat Somewhat less frequent 64 5.2
Q26_1 To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme heat Much less frequent 18 1.5
Q26_10 To what extent do you feel the following have become more common in Portugal over your lifetime? - Floods Somewhat more frequent 530 43.4
Q26_10 To what extent do you feel the following have become more common in Portugal over your lifetime? - Floods Much more frequent 428 35.1
Q26_10 To what extent do you feel the following have become more common in Portugal over your lifetime? - Floods Neither more nor less frequent 177 14.5
Q26_10 To what extent do you feel the following have become more common in Portugal over your lifetime? - Floods Somewhat less frequent 52 4.3
Q26_10 To what extent do you feel the following have become more common in Portugal over your lifetime? - Floods Much less frequent 34 2.8
Q26_11 To what extent do you feel the following have become more common in Portugal over your lifetime? - Poor air quality Somewhat more frequent 488 40.0
Q26_11 To what extent do you feel the following have become more common in Portugal over your lifetime? - Poor air quality Much more frequent 337 27.6
Q26_11 To what extent do you feel the following have become more common in Portugal over your lifetime? - Poor air quality Neither more nor less frequent 292 23.9
Q26_11 To what extent do you feel the following have become more common in Portugal over your lifetime? - Poor air quality Somewhat less frequent 68 5.6
Q26_11 To what extent do you feel the following have become more common in Portugal over your lifetime? - Poor air quality Much less frequent 36 2.9
Q26_12 To what extent do you feel the following have become more common in Portugal over your lifetime? - Landslides Somewhat more frequent 523 42.8
Q26_12 To what extent do you feel the following have become more common in Portugal over your lifetime? - Landslides Neither more nor less frequent 296 24.2
Q26_12 To what extent do you feel the following have become more common in Portugal over your lifetime? - Landslides Much more frequent 289 23.7
Q26_12 To what extent do you feel the following have become more common in Portugal over your lifetime? - Landslides Somewhat less frequent 71 5.8
Q26_12 To what extent do you feel the following have become more common in Portugal over your lifetime? - Landslides Much less frequent 42 3.4
Q26_13 To what extent do you feel the following have become more common in Portugal over your lifetime? - Sea swell Somewhat more frequent 530 43.4
Q26_13 To what extent do you feel the following have become more common in Portugal over your lifetime? - Sea swell Much more frequent 341 27.9
Q26_13 To what extent do you feel the following have become more common in Portugal over your lifetime? - Sea swell Neither more nor less frequent 284 23.3
Q26_13 To what extent do you feel the following have become more common in Portugal over your lifetime? - Sea swell Somewhat less frequent 42 3.4
Q26_13 To what extent do you feel the following have become more common in Portugal over your lifetime? - Sea swell Much less frequent 24 2.0
Q26_2 To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme cold Somewhat more frequent 497 40.7
Q26_2 To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme cold Neither more nor less frequent 311 25.5
Q26_2 To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme cold Much more frequent 209 17.1
Q26_2 To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme cold Somewhat less frequent 157 12.9
Q26_2 To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme cold Much less frequent 47 3.8
Q26_3 To what extent do you feel the following have become more common in Portugal over your lifetime? - Heavy rain Somewhat more frequent 565 46.3
Q26_3 To what extent do you feel the following have become more common in Portugal over your lifetime? - Heavy rain Much more frequent 397 32.5
Q26_3 To what extent do you feel the following have become more common in Portugal over your lifetime? - Heavy rain Neither more nor less frequent 183 15.0
Q26_3 To what extent do you feel the following have become more common in Portugal over your lifetime? - Heavy rain Somewhat less frequent 59 4.8
Q26_3 To what extent do you feel the following have become more common in Portugal over your lifetime? - Heavy rain Much less frequent 17 1.4
Q26_4 To what extent do you feel the following have become more common in Portugal over your lifetime? - Strong wind Somewhat more frequent 569 46.6
Q26_4 To what extent do you feel the following have become more common in Portugal over your lifetime? - Strong wind Much more frequent 378 31.0
Q26_4 To what extent do you feel the following have become more common in Portugal over your lifetime? - Strong wind Neither more nor less frequent 206 16.9
Q26_4 To what extent do you feel the following have become more common in Portugal over your lifetime? - Strong wind Somewhat less frequent 51 4.2
Q26_4 To what extent do you feel the following have become more common in Portugal over your lifetime? - Strong wind Much less frequent 17 1.4
Q26_5 To what extent do you feel the following have become more common in Portugal over your lifetime? - Snow Neither more nor less frequent 535 43.8
Q26_5 To what extent do you feel the following have become more common in Portugal over your lifetime? - Snow Somewhat more frequent 228 18.7
Q26_5 To what extent do you feel the following have become more common in Portugal over your lifetime? - Snow Much less frequent 203 16.6
Q26_5 To what extent do you feel the following have become more common in Portugal over your lifetime? - Snow Somewhat less frequent 189 15.5
Q26_5 To what extent do you feel the following have become more common in Portugal over your lifetime? - Snow Much more frequent 66 5.4
Q26_6 To what extent do you feel the following have become more common in Portugal over your lifetime? - Ice Neither more nor less frequent 542 44.4
Q26_6 To what extent do you feel the following have become more common in Portugal over your lifetime? - Ice Somewhat more frequent 226 18.5
Q26_6 To what extent do you feel the following have become more common in Portugal over your lifetime? - Ice Somewhat less frequent 204 16.7
Q26_6 To what extent do you feel the following have become more common in Portugal over your lifetime? - Ice Much less frequent 173 14.2
Q26_6 To what extent do you feel the following have become more common in Portugal over your lifetime? - Ice Much more frequent 76 6.2
Q26_7 To what extent do you feel the following have become more common in Portugal over your lifetime? - Winter storms Somewhat more frequent 542 44.4
Q26_7 To what extent do you feel the following have become more common in Portugal over your lifetime? - Winter storms Much more frequent 315 25.8
Q26_7 To what extent do you feel the following have become more common in Portugal over your lifetime? - Winter storms Neither more nor less frequent 230 18.8
Q26_7 To what extent do you feel the following have become more common in Portugal over your lifetime? - Winter storms Somewhat less frequent 78 6.4
Q26_7 To what extent do you feel the following have become more common in Portugal over your lifetime? - Winter storms Much less frequent 56 4.6
Q26_8 To what extent do you feel the following have become more common in Portugal over your lifetime? - Summer storms Neither more nor less frequent 464 38.0
Q26_8 To what extent do you feel the following have become more common in Portugal over your lifetime? - Summer storms Somewhat more frequent 416 34.1
Q26_8 To what extent do you feel the following have become more common in Portugal over your lifetime? - Summer storms Much more frequent 142 11.6
Q26_8 To what extent do you feel the following have become more common in Portugal over your lifetime? - Summer storms Somewhat less frequent 118 9.7
Q26_8 To what extent do you feel the following have become more common in Portugal over your lifetime? - Summer storms Much less frequent 81 6.6
Q26_9 To what extent do you feel the following have become more common in Portugal over your lifetime? - Wildfires Much more frequent 588 48.2
Q26_9 To what extent do you feel the following have become more common in Portugal over your lifetime? - Wildfires Somewhat more frequent 453 37.1
Q26_9 To what extent do you feel the following have become more common in Portugal over your lifetime? - Wildfires Neither more nor less frequent 134 11.0
Q26_9 To what extent do you feel the following have become more common in Portugal over your lifetime? - Wildfires Somewhat less frequent 26 2.1
Q26_9 To what extent do you feel the following have become more common in Portugal over your lifetime? - Wildfires Much less frequent 20 1.6
Q31_1 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Bottled water for three days Yes 815 66.7
Q31_1 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Bottled water for three days No 343 28.1
Q31_1 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Bottled water for three days I am not sure 60 4.9
Q31_1 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Bottled water for three days Missing / blank 3 0.2
Q31_2 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Food that will keep (non-perishable) for three days Yes 956 78.3
Q31_2 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Food that will keep (non-perishable) for three days No 189 15.5
Q31_2 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Food that will keep (non-perishable) for three days I am not sure 75 6.1
Q31_2 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Food that will keep (non-perishable) for three days Missing / blank 1 0.1
Q31_3 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Flashlight Yes 1003 82.1
Q31_3 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Flashlight No 179 14.7
Q31_3 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Flashlight I am not sure 33 2.7
Q31_3 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Flashlight Missing / blank 6 0.5
Q31_4 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Charger or portable battery, for example to charge a mobile phone during a power outage Yes 784 64.2
Q31_4 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Charger or portable battery, for example to charge a mobile phone during a power outage No 395 32.4
Q31_4 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Charger or portable battery, for example to charge a mobile phone during a power outage I am not sure 38 3.1
Q31_4 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Charger or portable battery, for example to charge a mobile phone during a power outage Missing / blank 4 0.3
Q31_5 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Batteries Yes 989 81.0
Q31_5 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Batteries No 181 14.8
Q31_5 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Batteries I am not sure 48 3.9
Q31_5 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Batteries Missing / blank 3 0.2
Q31_6 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Battery-powered or hand-crank radio Yes 602 49.3
Q31_6 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Battery-powered or hand-crank radio No 572 46.8
Q31_6 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Battery-powered or hand-crank radio I am not sure 44 3.6
Q31_6 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Battery-powered or hand-crank radio Missing / blank 3 0.2
Q31_7 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - First aid kit Yes 809 66.3
Q31_7 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - First aid kit No 351 28.7
Q31_7 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - First aid kit I am not sure 52 4.3
Q31_7 In an emergency causing interruption to utilities and other important services, do you have the following items at home? - First aid kit Missing / blank 9 0.7
Q32_1 In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Essential medicines for household members who need them, lasting at least three days Yes 968 79.3
Q32_1 In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Essential medicines for household members who need them, lasting at least three days No 154 12.6
Q32_1 In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Essential medicines for household members who need them, lasting at least three days Not applicable 98 8.0
Q32_1 In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Essential medicines for household members who need them, lasting at least three days Missing / blank 1 0.1
Q32_2 In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Baby food for baby/babies in your household lasting at least three days Not applicable 668 54.7
Q32_2 In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Baby food for baby/babies in your household lasting at least three days Yes 310 25.4
Q32_2 In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Baby food for baby/babies in your household lasting at least three days No 240 19.7
Q32_2 In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Baby food for baby/babies in your household lasting at least three days Missing / blank 3 0.2
Q32_3 In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Pet food lasting at least three days Yes 758 62.1
Q32_3 In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Pet food lasting at least three days Not applicable 299 24.5
Q32_3 In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Pet food lasting at least three days No 161 13.2
Q32_3 In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Pet food lasting at least three days Missing / blank 3 0.2
Q32_4 In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Spare glasses or contact lenses Yes 522 42.8
Q32_4 In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Spare glasses or contact lenses No 394 32.3
Q32_4 In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Spare glasses or contact lenses Not applicable 296 24.2
Q32_4 In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Spare glasses or contact lenses Missing / blank 9 0.7
Q34_1 I cannot leave work early to avoid travelling Missing / blank 904 74.0
Q34_1 I cannot leave work early to avoid travelling I cannot leave work early to avoid travelling 317 26.0
Q34_10 It is difficult for me or other household members to adopt the suggested measures because of mobility problems Missing / blank 1095 89.7
Q34_10 It is difficult for me or other household members to adopt the suggested measures because of mobility problems It is difficult for me or other household members to adopt the suggested measures because of mobility problems 126 10.3
Q34_11 I cannot or do not like to open windows at night to let cool air in Missing / blank 1067 87.4
Q34_11 I cannot or do not like to open windows at night to let cool air in I cannot or do not like to open windows at night to let cool air in 154 12.6
Q34_12 I do not have the physical strength to secure and protect all objects and equipment on my property, and have no one to help me do so Missing / blank 997 81.7
Q34_12 I do not have the physical strength to secure and protect all objects and equipment on my property, and have no one to help me do so I do not have the physical strength to secure and protect all objects and equipment on my property, and have no one to help me do so 224 18.3
Q34_13 Other response provided Missing / blank 1181 96.7
Q34_13 Other response provided None 5 0.4
Q34_13 Other response provided nenhum 3 0.2
Q34_13 Other response provided . 1 0.1
Q34_13 Other response provided Having difficulty taking the pets 1 0.1
Q34_13 Other response provided Having financial difficulties and living in a car 1 0.1
Q34_13 Other response provided I can leave anyway; there is no reason holding me back 1 0.1
Q34_13 Other response provided I have no problem implementing the previous solutions 1 0.1
Q34_13 Other response provided I would not be able to protect all my belongings, e.g. PC, camera… 1 0.1
Q34_13 Other response provided I would not leave without taking my dogs and cat with me. 1 0.1
Q34_13 Other response provided Lack of transport conditions 1 0.1
Q34_13 Other response provided No gosto de sair, portanto quanto pior o clima, melhor o meu dia. Eu amo estar em casa. 1 0.1
Q34_13 Other response provided No trabalho 1 0.1
Q34_13 Other response provided No vejo nenhuma impossibilidade 1 0.1
Q34_13 Other response provided None caso se aplica. 1 0.1
Q34_13 Other response provided None destes 1 0.1
Q34_13 Other response provided None destes motivos se aplica 1 0.1
Q34_13 Other response provided None dos acima referidos 1 0.1
Q34_13 Other response provided None dos anteriores 1 0.1
Q34_13 Other response provided None motivo 1 0.1
Q34_13 Other response provided None se aplica a mim 1 0.1
Q34_13 Other response provided Nonea das anteriores 1 0.1
Q34_13 Other response provided Not applicable 1 0.1
Q34_13 Other response provided Not being able to gather the pets quickly and bring them with me 1 0.1
Q34_13 Other response provided Nothing a pontar 1 0.1
Q34_13 Other response provided Nothing mencionado se aplica 1 0.1
Q34_13 Other response provided On that occasion I will certainly be at home 1 0.1
Q34_13 Other response provided Sem motivos de impedimento 1 0.1
Q34_13 Other response provided So faria alguma coisa quando os membros da familis estivessem todos juntos 1 0.1
Q34_13 Other response provided Unforeseen situations that can happen during events of this nature 1 0.1
Q34_13 Other response provided depending on the time available, it would be possible to carry out some protective work 1 0.1
Q34_13 Other response provided n/a 1 0.1
Q34_13 Other response provided nenhuma 1 0.1
Q34_13 Other response provided not applicable 1 0.1
Q34_13 Other response provided there is no one to provide help 1 0.1
Q34_2 I cannot work from home to avoid travelling Missing / blank 813 66.6
Q34_2 I cannot work from home to avoid travelling I cannot work from home to avoid travelling 408 33.4
Q34_3 I cannot avoid doing outdoor tasks or physical exertion during extreme heat because of work responsibilities Missing / blank 1032 84.5
Q34_3 I cannot avoid doing outdoor tasks or physical exertion during extreme heat because of work responsibilities I cannot avoid doing outdoor tasks or physical exertion during extreme heat because of work responsibilities 189 15.5
Q34_4 I cannot avoid working in an exposed area during extreme heat or thunderstorms Missing / blank 1033 84.6
Q34_4 I cannot avoid working in an exposed area during extreme heat or thunderstorms I cannot avoid working in an exposed area during extreme heat or thunderstorms 188 15.4
Q34_5 Difficulty leaving my home in the event of a flood or wildfire warning because I depend on others for transport Missing / blank 1054 86.3
Q34_5 Difficulty leaving my home in the event of a flood or wildfire warning because I depend on others for transport Difficulty leaving my home in the event of a flood or wildfire warning because I depend on others for transport 167 13.7
Q34_6 I cannot avoid travelling because of caring responsibilities for one or more vulnerable people Missing / blank 1081 88.5
Q34_6 I cannot avoid travelling because of caring responsibilities for one or more vulnerable people I cannot avoid travelling because of caring responsibilities for one or more vulnerable people 140 11.5
Q34_7 Difficulty finding a cool place to stay during extreme heat Missing / blank 993 81.3
Q34_7 Difficulty finding a cool place to stay during extreme heat Difficulty finding a cool place to stay during extreme heat 228 18.7
Q34_8 Difficulty finding sheltered place to stay during heavy rain Missing / blank 1030 84.4
Q34_8 Difficulty finding sheltered place to stay during heavy rain Difficulty finding sheltered place to stay during heavy rain 191 15.6
Q34_9 I would not feel comfortable cancelling trips or travel to visit friends or family if avoiding unnecessary travel due to bad weather was recommended Missing / blank 1145 93.8
Q34_9 I would not feel comfortable cancelling trips or travel to visit friends or family if avoiding unnecessary travel due to bad weather was recommended I would not feel comfortable cancelling trips or travel to visit friends or family if avoiding unnecessary travel due to bad weather was recommended 76 6.2 
outcome_choice_percents %>%
  group_by(outcome) %>%
  slice_max(order_by = percent, n = 1, with_ties = FALSE) %>%
  ungroup() %>%
  mutate(outcome_label = paste0(outcome, ": ", str_trunc(question, 70))) %>%
  ggplot(aes(x = fct_reorder(outcome_label, percent), y = percent)) +
  geom_col() +
  coord_flip() +
  scale_y_continuous(labels = label_percent(scale = 1)) +
  labs(
    x = NULL,
    y = "Percent",
    title = "Most common response for each outcome"
  ) +
  theme_minimal()

6 ANOVAs by experimental condition

These models test each outcome against Exp1Type and Exp2Type separately.

anova_outcomes <- c(
  paste0(q31_cols, "_bin"),
  paste0(q32_cols, "_bin"),
  paste0(q34_cols, "_bin"),
  paste0(q26_cols, "_ord_num")
) %>%
  present_cols()

anova_exp1_results <- map_dfr(anova_outcomes, ~ safe_aov_tidy(dataframe, .x, "Exp1Type"))
anova_exp2_results <- map_dfr(anova_outcomes, ~ safe_aov_tidy(dataframe, .x, "Exp2Type"))

anova_results <- bind_rows(anova_exp1_results, anova_exp2_results) %>%
  filter(term %in% c("Exp1Type", "Exp2Type")) %>%
  mutate(
    question = map_chr(outcome, qlabel),
    p.value = round(p.value, 4),
    statistic = round(statistic, 3),
    sig = p_stars(p.value)
  ) %>%
  dplyr::select(outcome, question, predictor, df, statistic, p.value, sig)

anova_results %>%
  arrange(p.value) %>%
  knitr::kable(caption = "One-way ANOVA results, sorted by p-value.")
One-way ANOVA results, sorted by p-value.
outcome question predictor df statistic p.value sig
Q32_3_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Pet food lasting at least three days Exp2Type 5 2.768 0.0171 *
Q34_10_bin It is difficult for me or other household members to adopt the suggested measures because of mobility problems Exp1Type 8 1.937 0.0512 †
Q34_9_bin I would not feel comfortable cancelling trips or travel to visit friends or family if avoiding unnecessary travel due to bad weather was recommended Exp1Type 8 1.866 0.0616 †
Q26_8_ord_num To what extent do you feel the following have become more common in Portugal over your lifetime? - Summer storms Exp2Type 5 2.044 0.0700 †
Q26_12_ord_num To what extent do you feel the following have become more common in Portugal over your lifetime? - Landslides Exp2Type 5 1.936 0.0856 †
Q26_11_ord_num To what extent do you feel the following have become more common in Portugal over your lifetime? - Poor air quality Exp2Type 5 1.833 0.1035
Q31_4_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Charger or portable battery, for example to charge a mobile phone during a power outage Exp2Type 5 1.832 0.1037
Q26_5_ord_num To what extent do you feel the following have become more common in Portugal over your lifetime? - Snow Exp2Type 5 1.702 0.1313
Q34_7_bin Difficulty finding a cool place to stay during extreme heat Exp1Type 8 1.534 0.1409
Q26_4_ord_num To what extent do you feel the following have become more common in Portugal over your lifetime? - Strong wind Exp1Type 8 1.525 0.1437
Q31_4_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Charger or portable battery, for example to charge a mobile phone during a power outage Exp1Type 8 1.489 0.1567
Q26_9_ord_num To what extent do you feel the following have become more common in Portugal over your lifetime? - Wildfires Exp1Type 8 1.385 0.1987
Q34_7_bin Difficulty finding a cool place to stay during extreme heat Exp2Type 5 1.440 0.2070
Q34_8_bin Difficulty finding sheltered place to stay during heavy rain Exp2Type 5 1.412 0.2169
Q26_13_ord_num To what extent do you feel the following have become more common in Portugal over your lifetime? - Sea swell Exp1Type 8 1.330 0.2241
Q31_3_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Flashlight Exp2Type 5 1.382 0.2282
Q34_5_bin Difficulty leaving my home in the event of a flood or wildfire warning because I depend on others for transport Exp1Type 8 1.309 0.2347
Q26_6_ord_num To what extent do you feel the following have become more common in Portugal over your lifetime? - Ice Exp2Type 5 1.309 0.2576
Q26_7_ord_num To what extent do you feel the following have become more common in Portugal over your lifetime? - Winter storms Exp1Type 8 1.245 0.2692
Q31_1_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Bottled water for three days Exp2Type 5 1.234 0.2908
Q32_3_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Pet food lasting at least three days Exp1Type 8 1.189 0.3018
Q31_7_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? - First aid kit Exp1Type 8 1.180 0.3074
Q26_11_ord_num To what extent do you feel the following have become more common in Portugal over your lifetime? - Poor air quality Exp1Type 8 1.163 0.3183
Q34_10_bin It is difficult for me or other household members to adopt the suggested measures because of mobility problems Exp2Type 5 1.167 0.3233
Q26_10_ord_num To what extent do you feel the following have become more common in Portugal over your lifetime? - Floods Exp1Type 8 1.121 0.3462
Q26_10_ord_num To what extent do you feel the following have become more common in Portugal over your lifetime? - Floods Exp2Type 5 1.111 0.3526
Q34_11_bin I cannot or do not like to open windows at night to let cool air in Exp2Type 5 1.103 0.3572
Q26_9_ord_num To what extent do you feel the following have become more common in Portugal over your lifetime? - Wildfires Exp2Type 5 1.100 0.3586
Q32_2_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Baby food for baby/babies in your household lasting at least three days Exp1Type 8 1.093 0.3649
Q34_5_bin Difficulty leaving my home in the event of a flood or wildfire warning because I depend on others for transport Exp2Type 5 0.962 0.4397
Q32_2_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Baby food for baby/babies in your household lasting at least three days Exp2Type 5 0.942 0.4526
Q34_1_bin I cannot leave work early to avoid travelling Exp2Type 5 0.925 0.4639
Q32_4_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Spare glasses or contact lenses Exp2Type 5 0.917 0.4687
Q26_7_ord_num To what extent do you feel the following have become more common in Portugal over your lifetime? - Winter storms Exp2Type 5 0.914 0.4709
Q31_7_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? - First aid kit Exp2Type 5 0.826 0.5311
Q34_1_bin I cannot leave work early to avoid travelling Exp1Type 8 0.876 0.5359
Q26_4_ord_num To what extent do you feel the following have become more common in Portugal over your lifetime? - Strong wind Exp2Type 5 0.818 0.5371
Q26_3_ord_num To what extent do you feel the following have become more common in Portugal over your lifetime? - Heavy rain Exp2Type 5 0.814 0.5398
Q31_1_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Bottled water for three days Exp1Type 8 0.866 0.5444
Q34_4_bin I cannot avoid working in an exposed area during extreme heat or thunderstorms Exp2Type 5 0.791 0.5559
Q34_2_bin I cannot work from home to avoid travelling Exp1Type 8 0.810 0.5936
Q32_1_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Essential medicines for household members who need them, lasting at least three days Exp1Type 8 0.798 0.6040
Q26_12_ord_num To what extent do you feel the following have become more common in Portugal over your lifetime? - Landslides Exp1Type 8 0.794 0.6084
Q26_1_ord_num To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme heat Exp2Type 5 0.719 0.6093
Q26_8_ord_num To what extent do you feel the following have become more common in Portugal over your lifetime? - Summer storms Exp1Type 8 0.785 0.6158
Q34_12_bin I do not have the physical strength to secure and protect all objects and equipment on my property, and have no one to help me do so Exp1Type 8 0.773 0.6270
Q34_12_bin I do not have the physical strength to secure and protect all objects and equipment on my property, and have no one to help me do so Exp2Type 5 0.678 0.6405
Q31_6_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Battery-powered or hand-crank radio Exp1Type 8 0.755 0.6431
Q34_8_bin Difficulty finding sheltered place to stay during heavy rain Exp1Type 8 0.748 0.6486
Q26_2_ord_num To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme cold Exp2Type 5 0.635 0.6732
Q32_1_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Essential medicines for household members who need them, lasting at least three days Exp2Type 5 0.611 0.6912
Q31_3_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Flashlight Exp1Type 8 0.685 0.7055
Q34_4_bin I cannot avoid working in an exposed area during extreme heat or thunderstorms Exp1Type 8 0.680 0.7092
Q31_2_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Food that will keep (non-perishable) for three days Exp2Type 5 0.582 0.7135
Q34_13_bin Other response provided Exp2Type 5 0.541 0.7453
Q26_13_ord_num To what extent do you feel the following have become more common in Portugal over your lifetime? - Sea swell Exp2Type 5 0.512 0.7675
Q32_4_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Spare glasses or contact lenses Exp1Type 8 0.589 0.7880
Q34_6_bin I cannot avoid travelling because of caring responsibilities for one or more vulnerable people Exp2Type 5 0.464 0.8030
Q31_5_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Batteries Exp2Type 5 0.456 0.8089
Q26_3_ord_num To what extent do you feel the following have become more common in Portugal over your lifetime? - Heavy rain Exp1Type 8 0.555 0.8153
Q26_1_ord_num To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme heat Exp1Type 8 0.551 0.8184
Q34_3_bin I cannot avoid doing outdoor tasks or physical exertion during extreme heat because of work responsibilities Exp1Type 8 0.525 0.8383
Q31_6_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Battery-powered or hand-crank radio Exp2Type 5 0.408 0.8435
Q31_5_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Batteries Exp1Type 8 0.492 0.8624
Q34_2_bin I cannot work from home to avoid travelling Exp2Type 5 0.350 0.8821
Q26_5_ord_num To what extent do you feel the following have become more common in Portugal over your lifetime? - Snow Exp1Type 8 0.438 0.8986
Q31_2_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? - Food that will keep (non-perishable) for three days Exp1Type 8 0.431 0.9027
Q34_9_bin I would not feel comfortable cancelling trips or travel to visit friends or family if avoiding unnecessary travel due to bad weather was recommended Exp2Type 5 0.277 0.9256
Q34_11_bin I cannot or do not like to open windows at night to let cool air in Exp1Type 8 0.360 0.9412
Q34_3_bin I cannot avoid doing outdoor tasks or physical exertion during extreme heat because of work responsibilities Exp2Type 5 0.233 0.9481
Q26_6_ord_num To what extent do you feel the following have become more common in Portugal over your lifetime? - Ice Exp1Type 8 0.344 0.9487
Q34_13_bin Other response provided Exp1Type 8 0.259 0.9786
Q34_6_bin I cannot avoid travelling because of caring responsibilities for one or more vulnerable people Exp1Type 8 0.254 0.9799
Q26_2_ord_num To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme cold Exp1Type 8 0.241 0.9831 
anova_results %>%
  filter(!is.na(p.value)) %>%
  mutate(
    outcome_label = paste0(outcome, ": ", str_trunc(question, 65)),
    neg_log10_p = -log10(p.value)
  ) %>%
  ggplot(aes(x = fct_reorder(outcome_label, neg_log10_p), y = neg_log10_p)) +
  geom_col() +
  geom_hline(yintercept = -log10(.05), linetype = "dashed") +
  coord_flip() +
  facet_wrap(~ predictor) +
  labs(
    x = NULL,
    y = expression(-log[10](p)),
    title = "ANOVA evidence by outcome and experimental condition",
    subtitle = "Dashed line marks p = .05"
  ) +
  theme_minimal()

7 Multivariate logit models: Q32

The following binary logit models estimate whether respondents have each emergency item in Q32, controlling for age, gender, education, income, carer status, and urbanicity.

model_formula_rhs <- "age + gender + educ + income + carer + urbanicity"

q32_outcomes <- paste0(q32_cols, "_bin") %>% present_cols()

q32_logit_models <- map(q32_outcomes, function(outcome) {
  safe_glm(as.formula(paste(outcome, "~", model_formula_rhs)), dataframe)
})
names(q32_logit_models) <- q32_outcomes

q32_logit_results <- imap_dfr(
  q32_logit_models,
  ~ safe_tidy_model(.x, .y, exponentiate = TRUE)
) %>%
  filter(!str_detect(coalesce(term, ""), "\\|")) %>%
  mutate(
    question = map_chr(outcome, qlabel),
    odds_ratio = estimate,
    sig = p_stars(p.value)
  ) %>%
  dplyr::select(outcome, question, term, odds_ratio, conf.low, conf.high, p.value, sig, model_error)

q32_logit_results %>%
  mutate(across(c(odds_ratio, conf.low, conf.high, p.value), ~ round(.x, 3))) %>%
  knitr::kable(caption = "Q32 logit models: odds ratios with 95% confidence intervals.")
Q32 logit models: odds ratios with 95% confidence intervals.
outcome question term odds_ratio conf.low conf.high p.value sig model_error
Q32_1_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Essential medicines for household members who need them, lasting at least three days (Intercept) 0.817 0.379 1.775 0.606 NA
Q32_1_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Essential medicines for household members who need them, lasting at least three days age 1.016 1.006 1.026 0.002 ** NA
Q32_1_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Essential medicines for household members who need them, lasting at least three days genderFemale 1.488 1.106 2.006 0.009 ** NA
Q32_1_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Essential medicines for household members who need them, lasting at least three days educSecondary 1.454 0.865 2.386 0.147 NA
Q32_1_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Essential medicines for household members who need them, lasting at least three days educTertiary 2.080 1.197 3.550 0.008 ** NA
Q32_1_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Essential medicines for household members who need them, lasting at least three days incomeMiddle 1.344 0.967 1.861 0.077 † NA
Q32_1_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Essential medicines for household members who need them, lasting at least three days incomeHigh 2.574 1.050 7.757 0.059 † NA
Q32_1_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Essential medicines for household members who need them, lasting at least three days incomePrefer not to answer 1.603 0.856 3.204 0.158 NA
Q32_1_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Essential medicines for household members who need them, lasting at least three days carerCarer 1.247 0.924 1.683 0.149 NA
Q32_1_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Essential medicines for household members who need them, lasting at least three days urbanicitySuburban 0.705 0.504 0.990 0.042 * NA
Q32_1_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Essential medicines for household members who need them, lasting at least three days urbanicitySemi-rural 0.730 0.459 1.184 0.191 NA
Q32_1_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Essential medicines for household members who need them, lasting at least three days urbanicityRural 1.028 0.628 1.736 0.914 NA
Q32_2_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Baby food for baby/babies in your household lasting at least three days (Intercept) 0.574 0.267 1.222 0.152 NA
Q32_2_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Baby food for baby/babies in your household lasting at least three days age 0.977 0.967 0.986 0.000 *** NA
Q32_2_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Baby food for baby/babies in your household lasting at least three days genderFemale 0.892 0.671 1.184 0.429 NA
Q32_2_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Baby food for baby/babies in your household lasting at least three days educSecondary 0.727 0.443 1.214 0.215 NA
Q32_2_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Baby food for baby/babies in your household lasting at least three days educTertiary 0.595 0.351 1.021 0.056 † NA
Q32_2_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Baby food for baby/babies in your household lasting at least three days incomeMiddle 1.149 0.834 1.591 0.400 NA
Q32_2_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Baby food for baby/babies in your household lasting at least three days incomeHigh 1.422 0.684 2.869 0.333 NA
Q32_2_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Baby food for baby/babies in your household lasting at least three days incomePrefer not to answer 0.561 0.271 1.087 0.101 NA
Q32_2_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Baby food for baby/babies in your household lasting at least three days carerCarer 3.660 2.680 5.058 0.000 *** NA
Q32_2_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Baby food for baby/babies in your household lasting at least three days urbanicitySuburban 0.922 0.657 1.288 0.636 NA
Q32_2_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Baby food for baby/babies in your household lasting at least three days urbanicitySemi-rural 1.162 0.731 1.824 0.517 NA
Q32_2_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Baby food for baby/babies in your household lasting at least three days urbanicityRural 1.199 0.754 1.882 0.436 NA
Q32_3_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Pet food lasting at least three days (Intercept) 2.386 1.215 4.739 0.012 * NA
Q32_3_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Pet food lasting at least three days age 0.986 0.978 0.994 0.001 *** NA
Q32_3_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Pet food lasting at least three days genderFemale 1.158 0.907 1.480 0.240 NA
Q32_3_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Pet food lasting at least three days educSecondary 0.886 0.549 1.407 0.612 NA
Q32_3_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Pet food lasting at least three days educTertiary 0.685 0.417 1.107 0.127 NA
Q32_3_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Pet food lasting at least three days incomeMiddle 1.258 0.950 1.664 0.109 NA
Q32_3_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Pet food lasting at least three days incomeHigh 1.748 0.898 3.558 0.110 NA
Q32_3_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Pet food lasting at least three days incomePrefer not to answer 1.469 0.871 2.527 0.156 NA
Q32_3_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Pet food lasting at least three days carerCarer 1.212 0.947 1.551 0.127 NA
Q32_3_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Pet food lasting at least three days urbanicitySuburban 1.096 0.829 1.451 0.521 NA
Q32_3_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Pet food lasting at least three days urbanicitySemi-rural 1.265 0.841 1.926 0.264 NA
Q32_3_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Pet food lasting at least three days urbanicityRural 1.889 1.225 2.977 0.005 ** NA
Q32_4_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Spare glasses or contact lenses (Intercept) 0.591 0.308 1.130 0.113 NA
Q32_4_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Spare glasses or contact lenses age 1.007 1.000 1.015 0.064 † NA
Q32_4_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Spare glasses or contact lenses genderFemale 0.951 0.750 1.206 0.678 NA
Q32_4_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Spare glasses or contact lenses educSecondary 0.846 0.544 1.317 0.456 NA
Q32_4_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Spare glasses or contact lenses educTertiary 0.920 0.581 1.459 0.720 NA
Q32_4_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Spare glasses or contact lenses incomeMiddle 0.950 0.724 1.249 0.714 NA
Q32_4_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Spare glasses or contact lenses incomeHigh 1.899 1.003 3.672 0.051 † NA
Q32_4_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Spare glasses or contact lenses incomePrefer not to answer 1.005 0.607 1.652 0.985 NA
Q32_4_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Spare glasses or contact lenses carerCarer 1.177 0.925 1.499 0.185 NA
Q32_4_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Spare glasses or contact lenses urbanicitySuburban 0.885 0.670 1.166 0.386 NA
Q32_4_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Spare glasses or contact lenses urbanicitySemi-rural 0.995 0.668 1.474 0.979 NA
Q32_4_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Spare glasses or contact lenses urbanicityRural 1.134 0.763 1.683 0.532 NA 
q32_logit_results %>%
  filter(!is.na(odds_ratio), term != "(Intercept)") %>%
  mutate(outcome_label = paste0(outcome, ": ", str_trunc(question, 55))) %>%
  ggplot(aes(x = odds_ratio, y = fct_rev(term))) +
  geom_vline(xintercept = 1, linetype = "dashed") +
  geom_pointrange(aes(xmin = conf.low, xmax = conf.high)) +
  scale_x_log10() +
  facet_wrap(~ outcome_label) +
  labs(
    x = "Odds ratio, log scale",
    y = NULL,
    title = "Q32 logit model estimates"
  ) +
  theme_minimal()

8 Multivariate logit models: Q34

These binary logit models estimate whether each barrier to responding to a severe-weather warning applies to the respondent.

binary_outcomes <- paste0(q34_cols, "_bin") %>% present_cols()

binary_logit_models <- map(binary_outcomes, function(outcome) {
  safe_glm(as.formula(paste(outcome, "~", model_formula_rhs)), dataframe)
})
names(binary_logit_models) <- binary_outcomes

binary_logit_results <- imap_dfr(
  binary_logit_models,
  ~ safe_tidy_model(.x, .y, exponentiate = TRUE)
) %>%
  mutate(
    question = map_chr(outcome, qlabel),
    odds_ratio = estimate,
    sig = p_stars(p.value)
  ) %>%
  dplyr::select(outcome, question, term, odds_ratio, conf.low, conf.high, p.value, sig, model_error)

binary_logit_results %>%
  mutate(across(c(odds_ratio, conf.low, conf.high, p.value), ~ round(.x, 3))) %>%
  knitr::kable(caption = "Q34 logit models: odds ratios with 95% confidence intervals.")
Q34 logit models: odds ratios with 95% confidence intervals.
outcome question term odds_ratio conf.low conf.high p.value sig model_error
Q34_1_bin I cannot leave work early to avoid travelling (Intercept) 0.553 0.249 1.193 0.138 NA
Q34_1_bin I cannot leave work early to avoid travelling age 0.977 0.968 0.986 0.000 *** NA
Q34_1_bin I cannot leave work early to avoid travelling genderFemale 1.075 0.820 1.409 0.601 NA
Q34_1_bin I cannot leave work early to avoid travelling educSecondary 1.490 0.855 2.740 0.177 NA
Q34_1_bin I cannot leave work early to avoid travelling educTertiary 1.588 0.897 2.960 0.127 NA
Q34_1_bin I cannot leave work early to avoid travelling incomeMiddle 1.159 0.848 1.593 0.359 NA
Q34_1_bin I cannot leave work early to avoid travelling incomeHigh 1.481 0.728 2.909 0.263 NA
Q34_1_bin I cannot leave work early to avoid travelling incomePrefer not to answer 0.879 0.474 1.567 0.671 NA
Q34_1_bin I cannot leave work early to avoid travelling carerCarer 0.962 0.731 1.266 0.781 NA
Q34_1_bin I cannot leave work early to avoid travelling urbanicitySuburban 1.159 0.849 1.578 0.350 NA
Q34_1_bin I cannot leave work early to avoid travelling urbanicitySemi-rural 1.047 0.661 1.628 0.840 NA
Q34_1_bin I cannot leave work early to avoid travelling urbanicityRural 0.881 0.547 1.387 0.592 NA
Q34_2_bin I cannot work from home to avoid travelling (Intercept) 1.563 0.784 3.106 0.203 NA
Q34_2_bin I cannot work from home to avoid travelling age 0.978 0.969 0.986 0.000 *** NA
Q34_2_bin I cannot work from home to avoid travelling genderFemale 0.968 0.753 1.245 0.802 NA
Q34_2_bin I cannot work from home to avoid travelling educSecondary 0.956 0.600 1.548 0.851 NA
Q34_2_bin I cannot work from home to avoid travelling educTertiary 0.803 0.493 1.326 0.383 NA
Q34_2_bin I cannot work from home to avoid travelling incomeMiddle 1.205 0.904 1.612 0.205 NA
Q34_2_bin I cannot work from home to avoid travelling incomeHigh 0.963 0.472 1.882 0.914 NA
Q34_2_bin I cannot work from home to avoid travelling incomePrefer not to answer 0.708 0.394 1.230 0.232 NA
Q34_2_bin I cannot work from home to avoid travelling carerCarer 0.829 0.643 1.069 0.148 NA
Q34_2_bin I cannot work from home to avoid travelling urbanicitySuburban 1.030 0.769 1.376 0.844 NA
Q34_2_bin I cannot work from home to avoid travelling urbanicitySemi-rural 0.920 0.599 1.395 0.698 NA
Q34_2_bin I cannot work from home to avoid travelling urbanicityRural 1.123 0.740 1.691 0.580 NA
Q34_3_bin I cannot avoid doing outdoor tasks or physical exertion during extreme heat because of work responsibilities (Intercept) 0.323 0.128 0.786 0.014 * NA
Q34_3_bin I cannot avoid doing outdoor tasks or physical exertion during extreme heat because of work responsibilities age 0.979 0.968 0.990 0.000 *** NA
Q34_3_bin I cannot avoid doing outdoor tasks or physical exertion during extreme heat because of work responsibilities genderFemale 0.808 0.579 1.127 0.210 NA
Q34_3_bin I cannot avoid doing outdoor tasks or physical exertion during extreme heat because of work responsibilities educSecondary 1.273 0.704 2.454 0.445 NA
Q34_3_bin I cannot avoid doing outdoor tasks or physical exertion during extreme heat because of work responsibilities educTertiary 0.685 0.359 1.376 0.267 NA
Q34_3_bin I cannot avoid doing outdoor tasks or physical exertion during extreme heat because of work responsibilities incomeMiddle 0.945 0.656 1.370 0.764 NA
Q34_3_bin I cannot avoid doing outdoor tasks or physical exertion during extreme heat because of work responsibilities incomeHigh 0.663 0.217 1.658 0.419 NA
Q34_3_bin I cannot avoid doing outdoor tasks or physical exertion during extreme heat because of work responsibilities incomePrefer not to answer 0.576 0.242 1.221 0.177 NA
Q34_3_bin I cannot avoid doing outdoor tasks or physical exertion during extreme heat because of work responsibilities carerCarer 1.965 1.380 2.830 0.000 *** NA
Q34_3_bin I cannot avoid doing outdoor tasks or physical exertion during extreme heat because of work responsibilities urbanicitySuburban 1.429 0.979 2.078 0.063 † NA
Q34_3_bin I cannot avoid doing outdoor tasks or physical exertion during extreme heat because of work responsibilities urbanicitySemi-rural 0.891 0.493 1.540 0.690 NA
Q34_3_bin I cannot avoid doing outdoor tasks or physical exertion during extreme heat because of work responsibilities urbanicityRural 0.960 0.535 1.655 0.887 NA
Q34_4_bin I cannot avoid working in an exposed area during extreme heat or thunderstorms (Intercept) 0.629 0.260 1.484 0.295 NA
Q34_4_bin I cannot avoid working in an exposed area during extreme heat or thunderstorms age 0.975 0.964 0.987 0.000 *** NA
Q34_4_bin I cannot avoid working in an exposed area during extreme heat or thunderstorms genderFemale 0.719 0.513 1.003 0.053 † NA
Q34_4_bin I cannot avoid working in an exposed area during extreme heat or thunderstorms educSecondary 0.695 0.396 1.271 0.219 NA
Q34_4_bin I cannot avoid working in an exposed area during extreme heat or thunderstorms educTertiary 0.594 0.326 1.118 0.095 † NA
Q34_4_bin I cannot avoid working in an exposed area during extreme heat or thunderstorms incomeMiddle 1.087 0.741 1.611 0.674 NA
Q34_4_bin I cannot avoid working in an exposed area during extreme heat or thunderstorms incomeHigh 1.890 0.839 4.013 0.108 NA
Q34_4_bin I cannot avoid working in an exposed area during extreme heat or thunderstorms incomePrefer not to answer 1.391 0.691 2.664 0.334 NA
Q34_4_bin I cannot avoid working in an exposed area during extreme heat or thunderstorms carerCarer 1.298 0.922 1.840 0.138 NA
Q34_4_bin I cannot avoid working in an exposed area during extreme heat or thunderstorms urbanicitySuburban 0.987 0.656 1.468 0.949 NA
Q34_4_bin I cannot avoid working in an exposed area during extreme heat or thunderstorms urbanicitySemi-rural 1.725 1.038 2.808 0.031 * NA
Q34_4_bin I cannot avoid working in an exposed area during extreme heat or thunderstorms urbanicityRural 1.294 0.751 2.164 0.338 NA
Q34_5_bin Difficulty leaving my home in the event of a flood or wildfire warning because I depend on others for transport (Intercept) 0.305 0.119 0.758 0.012 * NA
Q34_5_bin Difficulty leaving my home in the event of a flood or wildfire warning because I depend on others for transport age 0.975 0.964 0.987 0.000 *** NA
Q34_5_bin Difficulty leaving my home in the event of a flood or wildfire warning because I depend on others for transport genderFemale 1.276 0.902 1.811 0.170 NA
Q34_5_bin Difficulty leaving my home in the event of a flood or wildfire warning because I depend on others for transport educSecondary 0.921 0.508 1.763 0.794 NA
Q34_5_bin Difficulty leaving my home in the event of a flood or wildfire warning because I depend on others for transport educTertiary 0.754 0.399 1.493 0.399 NA
Q34_5_bin Difficulty leaving my home in the event of a flood or wildfire warning because I depend on others for transport incomeMiddle 0.794 0.539 1.176 0.246 NA
Q34_5_bin Difficulty leaving my home in the event of a flood or wildfire warning because I depend on others for transport incomeHigh 1.341 0.562 2.942 0.483 NA
Q34_5_bin Difficulty leaving my home in the event of a flood or wildfire warning because I depend on others for transport incomePrefer not to answer 1.100 0.537 2.125 0.785 NA
Q34_5_bin Difficulty leaving my home in the event of a flood or wildfire warning because I depend on others for transport carerCarer 2.108 1.450 3.112 0.000 *** NA
Q34_5_bin Difficulty leaving my home in the event of a flood or wildfire warning because I depend on others for transport urbanicitySuburban 1.120 0.739 1.682 0.587 NA
Q34_5_bin Difficulty leaving my home in the event of a flood or wildfire warning because I depend on others for transport urbanicitySemi-rural 0.901 0.483 1.598 0.732 NA
Q34_5_bin Difficulty leaving my home in the event of a flood or wildfire warning because I depend on others for transport urbanicityRural 1.610 0.954 2.659 0.068 † NA
Q34_6_bin I cannot avoid travelling because of caring responsibilities for one or more vulnerable people (Intercept) 0.009 0.003 0.030 0.000 *** NA
Q34_6_bin I cannot avoid travelling because of caring responsibilities for one or more vulnerable people age 1.012 0.999 1.025 0.070 † NA
Q34_6_bin I cannot avoid travelling because of caring responsibilities for one or more vulnerable people genderFemale 0.883 0.607 1.283 0.515 NA
Q34_6_bin I cannot avoid travelling because of caring responsibilities for one or more vulnerable people educSecondary 1.919 0.893 4.776 0.122 NA
Q34_6_bin I cannot avoid travelling because of caring responsibilities for one or more vulnerable people educTertiary 2.054 0.928 5.219 0.098 † NA
Q34_6_bin I cannot avoid travelling because of caring responsibilities for one or more vulnerable people incomeMiddle 1.508 0.972 2.387 0.072 † NA
Q34_6_bin I cannot avoid travelling because of caring responsibilities for one or more vulnerable people incomeHigh 1.905 0.742 4.474 0.155 NA
Q34_6_bin I cannot avoid travelling because of caring responsibilities for one or more vulnerable people incomePrefer not to answer 1.062 0.410 2.431 0.893 NA
Q34_6_bin I cannot avoid travelling because of caring responsibilities for one or more vulnerable people carerCarer 5.574 3.484 9.337 0.000 *** NA
Q34_6_bin I cannot avoid travelling because of caring responsibilities for one or more vulnerable people urbanicitySuburban 0.919 0.581 1.430 0.711 NA
Q34_6_bin I cannot avoid travelling because of caring responsibilities for one or more vulnerable people urbanicitySemi-rural 1.773 1.013 3.024 0.039 * NA
Q34_6_bin I cannot avoid travelling because of caring responsibilities for one or more vulnerable people urbanicityRural 0.908 0.445 1.724 0.779 NA
Q34_7_bin Difficulty finding a cool place to stay during extreme heat (Intercept) 0.236 0.099 0.538 0.001 *** NA
Q34_7_bin Difficulty finding a cool place to stay during extreme heat age 0.999 0.989 1.009 0.829 NA
Q34_7_bin Difficulty finding a cool place to stay during extreme heat genderFemale 1.111 0.822 1.502 0.494 NA
Q34_7_bin Difficulty finding a cool place to stay during extreme heat educSecondary 1.280 0.724 2.400 0.416 NA
Q34_7_bin Difficulty finding a cool place to stay during extreme heat educTertiary 1.186 0.653 2.273 0.590 NA
Q34_7_bin Difficulty finding a cool place to stay during extreme heat incomeMiddle 0.852 0.608 1.200 0.356 NA
Q34_7_bin Difficulty finding a cool place to stay during extreme heat incomeHigh 0.707 0.277 1.583 0.430 NA
Q34_7_bin Difficulty finding a cool place to stay during extreme heat incomePrefer not to answer 0.839 0.428 1.554 0.592 NA
Q34_7_bin Difficulty finding a cool place to stay during extreme heat carerCarer 1.052 0.775 1.431 0.747 NA
Q34_7_bin Difficulty finding a cool place to stay during extreme heat urbanicitySuburban 0.766 0.535 1.086 0.139 NA
Q34_7_bin Difficulty finding a cool place to stay during extreme heat urbanicitySemi-rural 0.849 0.506 1.377 0.521 NA
Q34_7_bin Difficulty finding a cool place to stay during extreme heat urbanicityRural 0.575 0.320 0.979 0.051 † NA
Q34_8_bin Difficulty finding sheltered place to stay during heavy rain (Intercept) 0.119 0.046 0.294 0.000 *** NA
Q34_8_bin Difficulty finding sheltered place to stay during heavy rain age 1.004 0.994 1.015 0.419 NA
Q34_8_bin Difficulty finding sheltered place to stay during heavy rain genderFemale 1.007 0.729 1.393 0.964 NA
Q34_8_bin Difficulty finding sheltered place to stay during heavy rain educSecondary 1.289 0.694 2.592 0.446 NA
Q34_8_bin Difficulty finding sheltered place to stay during heavy rain educTertiary 1.399 0.734 2.866 0.330 NA
Q34_8_bin Difficulty finding sheltered place to stay during heavy rain incomeMiddle 0.915 0.630 1.338 0.642 NA
Q34_8_bin Difficulty finding sheltered place to stay during heavy rain incomeHigh 0.978 0.379 2.220 0.961 NA
Q34_8_bin Difficulty finding sheltered place to stay during heavy rain incomePrefer not to answer 1.813 0.987 3.244 0.049 * NA
Q34_8_bin Difficulty finding sheltered place to stay during heavy rain carerCarer 0.852 0.615 1.182 0.337 NA
Q34_8_bin Difficulty finding sheltered place to stay during heavy rain urbanicitySuburban 1.026 0.700 1.491 0.896 NA
Q34_8_bin Difficulty finding sheltered place to stay during heavy rain urbanicitySemi-rural 1.259 0.735 2.086 0.385 NA
Q34_8_bin Difficulty finding sheltered place to stay during heavy rain urbanicityRural 1.176 0.678 1.968 0.550 NA
Q34_9_bin I would not feel comfortable cancelling trips or travel to visit friends or family if avoiding unnecessary travel due to bad weather was recommended (Intercept) 0.047 0.008 0.199 0.000 *** NA
Q34_9_bin I would not feel comfortable cancelling trips or travel to visit friends or family if avoiding unnecessary travel due to bad weather was recommended age 0.976 0.959 0.993 0.005 ** NA
Q34_9_bin I would not feel comfortable cancelling trips or travel to visit friends or family if avoiding unnecessary travel due to bad weather was recommended genderFemale 0.959 0.591 1.556 0.866 NA
Q34_9_bin I would not feel comfortable cancelling trips or travel to visit friends or family if avoiding unnecessary travel due to bad weather was recommended educSecondary 1.963 0.680 8.315 0.274 NA
Q34_9_bin I would not feel comfortable cancelling trips or travel to visit friends or family if avoiding unnecessary travel due to bad weather was recommended educTertiary 2.176 0.730 9.388 0.217 NA
Q34_9_bin I would not feel comfortable cancelling trips or travel to visit friends or family if avoiding unnecessary travel due to bad weather was recommended incomeMiddle 1.154 0.665 2.060 0.618 NA
Q34_9_bin I would not feel comfortable cancelling trips or travel to visit friends or family if avoiding unnecessary travel due to bad weather was recommended incomeHigh 0.981 0.220 3.130 0.977 NA
Q34_9_bin I would not feel comfortable cancelling trips or travel to visit friends or family if avoiding unnecessary travel due to bad weather was recommended incomePrefer not to answer 1.050 0.336 2.738 0.926 NA
Q34_9_bin I would not feel comfortable cancelling trips or travel to visit friends or family if avoiding unnecessary travel due to bad weather was recommended carerCarer 2.313 1.364 4.087 0.003 ** NA
Q34_9_bin I would not feel comfortable cancelling trips or travel to visit friends or family if avoiding unnecessary travel due to bad weather was recommended urbanicitySuburban 1.215 0.683 2.122 0.498 NA
Q34_9_bin I would not feel comfortable cancelling trips or travel to visit friends or family if avoiding unnecessary travel due to bad weather was recommended urbanicitySemi-rural 1.193 0.518 2.498 0.657 NA
Q34_9_bin I would not feel comfortable cancelling trips or travel to visit friends or family if avoiding unnecessary travel due to bad weather was recommended urbanicityRural 1.348 0.604 2.784 0.439 NA
Q34_10_bin It is difficult for me or other household members to adopt the suggested measures because of mobility problems (Intercept) 0.070 0.024 0.190 0.000 *** NA
Q34_10_bin It is difficult for me or other household members to adopt the suggested measures because of mobility problems age 1.007 0.994 1.020 0.303 NA
Q34_10_bin It is difficult for me or other household members to adopt the suggested measures because of mobility problems genderFemale 1.050 0.712 1.549 0.806 NA
Q34_10_bin It is difficult for me or other household members to adopt the suggested measures because of mobility problems educSecondary 0.731 0.400 1.407 0.326 NA
Q34_10_bin It is difficult for me or other household members to adopt the suggested measures because of mobility problems educTertiary 0.665 0.346 1.337 0.235 NA
Q34_10_bin It is difficult for me or other household members to adopt the suggested measures because of mobility problems incomeMiddle 0.921 0.600 1.428 0.709 NA
Q34_10_bin It is difficult for me or other household members to adopt the suggested measures because of mobility problems incomeHigh 1.405 0.532 3.297 0.459 NA
Q34_10_bin It is difficult for me or other household members to adopt the suggested measures because of mobility problems incomePrefer not to answer 0.572 0.191 1.395 0.261 NA
Q34_10_bin It is difficult for me or other household members to adopt the suggested measures because of mobility problems carerCarer 2.552 1.662 4.019 0.000 *** NA
Q34_10_bin It is difficult for me or other household members to adopt the suggested measures because of mobility problems urbanicitySuburban 0.667 0.400 1.080 0.109 NA
Q34_10_bin It is difficult for me or other household members to adopt the suggested measures because of mobility problems urbanicitySemi-rural 1.007 0.521 1.833 0.981 NA
Q34_10_bin It is difficult for me or other household members to adopt the suggested measures because of mobility problems urbanicityRural 1.424 0.788 2.480 0.224 NA
Q34_11_bin I cannot or do not like to open windows at night to let cool air in (Intercept) 0.085 0.031 0.222 0.000 *** NA
Q34_11_bin I cannot or do not like to open windows at night to let cool air in age 1.004 0.993 1.016 0.454 NA
Q34_11_bin I cannot or do not like to open windows at night to let cool air in genderFemale 0.990 0.696 1.407 0.954 NA
Q34_11_bin I cannot or do not like to open windows at night to let cool air in educSecondary 1.256 0.671 2.541 0.499 NA
Q34_11_bin I cannot or do not like to open windows at night to let cool air in educTertiary 1.091 0.559 2.281 0.806 NA
Q34_11_bin I cannot or do not like to open windows at night to let cool air in incomeMiddle 0.853 0.575 1.274 0.432 NA
Q34_11_bin I cannot or do not like to open windows at night to let cool air in incomeHigh 0.608 0.175 1.625 0.370 NA
Q34_11_bin I cannot or do not like to open windows at night to let cool air in incomePrefer not to answer 1.504 0.766 2.824 0.218 NA
Q34_11_bin I cannot or do not like to open windows at night to let cool air in carerCarer 1.541 1.072 2.237 0.021 * NA
Q34_11_bin I cannot or do not like to open windows at night to let cool air in urbanicitySuburban 1.021 0.674 1.527 0.922 NA
Q34_11_bin I cannot or do not like to open windows at night to let cool air in urbanicitySemi-rural 0.584 0.275 1.123 0.130 NA
Q34_11_bin I cannot or do not like to open windows at night to let cool air in urbanicityRural 1.603 0.939 2.664 0.075 † NA
Q34_12_bin I do not have the physical strength to secure and protect all objects and equipment on my property, and have no one to help me do so (Intercept) 0.059 0.022 0.146 0.000 *** NA
Q34_12_bin I do not have the physical strength to secure and protect all objects and equipment on my property, and have no one to help me do so age 1.008 0.998 1.018 0.136 NA
Q34_12_bin I do not have the physical strength to secure and protect all objects and equipment on my property, and have no one to help me do so genderFemale 2.831 2.058 3.929 0.000 *** NA
Q34_12_bin I do not have the physical strength to secure and protect all objects and equipment on my property, and have no one to help me do so educSecondary 2.023 1.068 4.185 0.041 * NA
Q34_12_bin I do not have the physical strength to secure and protect all objects and equipment on my property, and have no one to help me do so educTertiary 1.757 0.904 3.708 0.114 NA
Q34_12_bin I do not have the physical strength to secure and protect all objects and equipment on my property, and have no one to help me do so incomeMiddle 0.856 0.607 1.213 0.377 NA
Q34_12_bin I do not have the physical strength to secure and protect all objects and equipment on my property, and have no one to help me do so incomeHigh 0.906 0.367 2.013 0.818 NA
Q34_12_bin I do not have the physical strength to secure and protect all objects and equipment on my property, and have no one to help me do so incomePrefer not to answer 0.811 0.416 1.506 0.521 NA
Q34_12_bin I do not have the physical strength to secure and protect all objects and equipment on my property, and have no one to help me do so carerCarer 0.859 0.631 1.172 0.337 NA
Q34_12_bin I do not have the physical strength to secure and protect all objects and equipment on my property, and have no one to help me do so urbanicitySuburban 1.117 0.786 1.579 0.534 NA
Q34_12_bin I do not have the physical strength to secure and protect all objects and equipment on my property, and have no one to help me do so urbanicitySemi-rural 1.150 0.691 1.863 0.580 NA
Q34_12_bin I do not have the physical strength to secure and protect all objects and equipment on my property, and have no one to help me do so urbanicityRural 0.557 0.303 0.970 0.048 * NA
Q34_13_bin Other response provided (Intercept) 0.003 0.000 0.021 0.000 *** NA
Q34_13_bin Other response provided age 1.048 1.025 1.074 0.000 *** NA
Q34_13_bin Other response provided genderFemale 0.992 0.506 1.947 0.982 NA
Q34_13_bin Other response provided educSecondary 1.304 0.406 5.898 0.688 NA
Q34_13_bin Other response provided educTertiary 1.540 0.463 7.128 0.523 NA
Q34_13_bin Other response provided incomeMiddle 1.054 0.479 2.470 0.899 NA
Q34_13_bin Other response provided incomeHigh 1.777 0.256 7.629 0.486 NA
Q34_13_bin Other response provided incomePrefer not to answer 2.093 0.609 6.400 0.209 NA
Q34_13_bin Other response provided carerCarer 0.363 0.169 0.734 0.006 ** NA
Q34_13_bin Other response provided urbanicitySuburban 0.873 0.382 1.882 0.737 NA
Q34_13_bin Other response provided urbanicitySemi-rural 0.793 0.181 2.444 0.718 NA
Q34_13_bin Other response provided urbanicityRural 1.797 0.619 4.587 0.244 NA 
binary_logit_results %>%
  filter(!is.na(p.value), p.value < .05, term != "(Intercept)") %>%
  arrange(p.value) %>%
  mutate(across(c(odds_ratio, conf.low, conf.high, p.value), ~ round(.x, 3))) %>%
  knitr::kable(caption = "Q34 predictors with p < .05.")
Q34 predictors with p < .05.
outcome question term odds_ratio conf.low conf.high p.value sig model_error
Q34_6_bin I cannot avoid travelling because of caring responsibilities for one or more vulnerable people carerCarer 5.574 3.484 9.337 0.000 *** NA
Q34_12_bin I do not have the physical strength to secure and protect all objects and equipment on my property, and have no one to help me do so genderFemale 2.831 2.058 3.929 0.000 *** NA
Q34_2_bin I cannot work from home to avoid travelling age 0.978 0.969 0.986 0.000 *** NA
Q34_1_bin I cannot leave work early to avoid travelling age 0.977 0.968 0.986 0.000 *** NA
Q34_4_bin I cannot avoid working in an exposed area during extreme heat or thunderstorms age 0.975 0.964 0.987 0.000 *** NA
Q34_10_bin It is difficult for me or other household members to adopt the suggested measures because of mobility problems carerCarer 2.552 1.662 4.019 0.000 *** NA
Q34_5_bin Difficulty leaving my home in the event of a flood or wildfire warning because I depend on others for transport age 0.975 0.964 0.987 0.000 *** NA
Q34_13_bin Other response provided age 1.048 1.025 1.074 0.000 *** NA
Q34_5_bin Difficulty leaving my home in the event of a flood or wildfire warning because I depend on others for transport carerCarer 2.108 1.450 3.112 0.000 *** NA
Q34_3_bin I cannot avoid doing outdoor tasks or physical exertion during extreme heat because of work responsibilities carerCarer 1.965 1.380 2.830 0.000 *** NA
Q34_3_bin I cannot avoid doing outdoor tasks or physical exertion during extreme heat because of work responsibilities age 0.979 0.968 0.990 0.000 *** NA
Q34_9_bin I would not feel comfortable cancelling trips or travel to visit friends or family if avoiding unnecessary travel due to bad weather was recommended carerCarer 2.313 1.364 4.087 0.003 ** NA
Q34_9_bin I would not feel comfortable cancelling trips or travel to visit friends or family if avoiding unnecessary travel due to bad weather was recommended age 0.976 0.959 0.993 0.005 ** NA
Q34_13_bin Other response provided carerCarer 0.363 0.169 0.734 0.006 ** NA
Q34_11_bin I cannot or do not like to open windows at night to let cool air in carerCarer 1.541 1.072 2.237 0.021 * NA
Q34_4_bin I cannot avoid working in an exposed area during extreme heat or thunderstorms urbanicitySemi-rural 1.725 1.038 2.808 0.031 * NA
Q34_6_bin I cannot avoid travelling because of caring responsibilities for one or more vulnerable people urbanicitySemi-rural 1.773 1.013 3.024 0.039 * NA
Q34_12_bin I do not have the physical strength to secure and protect all objects and equipment on my property, and have no one to help me do so educSecondary 2.023 1.068 4.185 0.041 * NA
Q34_12_bin I do not have the physical strength to secure and protect all objects and equipment on my property, and have no one to help me do so urbanicityRural 0.557 0.303 0.970 0.048 * NA
Q34_8_bin Difficulty finding sheltered place to stay during heavy rain incomePrefer not to answer 1.813 0.987 3.244 0.049 * NA

9 Ordered logit models: Q26

These ordered logit models estimate perceived changes in the frequency of weather and climate hazards. Following the original script, models are restricted to respondents who answered Yes to Q4.

q26_outcomes <- paste0(q26_cols, "_ord") %>% present_cols()

ordered_logit_models <- map(q26_outcomes, function(outcome) {
  safe_polr(
    as.formula(paste(outcome, "~", model_formula_rhs)),
    dataframe %>% filter(Q4 == "Yes")
  )
})
names(ordered_logit_models) <- q26_outcomes

ordered_logit_results <- imap_dfr(
  ordered_logit_models,
  ~ safe_tidy_model(.x, .y, exponentiate = TRUE)
) %>%
  filter(!str_detect(coalesce(term, ""), "\\|")) %>%
  mutate(
    question = map_chr(outcome, qlabel),
    odds_ratio = estimate,
    sig = p_stars(p.value)
  ) %>%
  dplyr::select(outcome, question, term, odds_ratio, conf.low, conf.high, p.value, sig, model_error)

ordered_logit_results %>%
  mutate(across(c(odds_ratio, conf.low, conf.high, p.value), ~ round(.x, 3))) %>%
  knitr::kable(caption = "Q26 ordered logit models: odds ratios with 95% confidence intervals.")
Q26 ordered logit models: odds ratios with 95% confidence intervals.
outcome question term odds_ratio conf.low conf.high p.value sig model_error
Q26_1_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme heat age 1.004 0.996 1.012 0.289 NA
Q26_1_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme heat genderFemale 1.325 1.040 1.691 0.023 * NA
Q26_1_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme heat educSecondary 1.249 0.791 1.972 0.340 NA
Q26_1_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme heat educTertiary 1.902 1.180 3.067 0.008 ** NA
Q26_1_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme heat incomeMiddle 1.351 1.019 1.792 0.037 * NA
Q26_1_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme heat incomeHigh 1.277 0.659 2.485 0.469 NA
Q26_1_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme heat incomePrefer not to answer 0.994 0.602 1.644 0.983 NA
Q26_1_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme heat carerCarer 0.778 0.609 0.993 0.044 * NA
Q26_1_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme heat urbanicitySuburban 0.778 0.588 1.027 0.077 † NA
Q26_1_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme heat urbanicitySemi-rural 0.974 0.646 1.471 0.901 NA
Q26_1_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme heat urbanicityRural 1.278 0.855 1.916 0.233 NA
Q26_2_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme cold age 1.005 0.997 1.012 0.228 NA
Q26_2_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme cold genderFemale 1.499 1.187 1.894 0.001 *** NA
Q26_2_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme cold educSecondary 0.932 0.605 1.433 0.749 NA
Q26_2_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme cold educTertiary 0.878 0.560 1.375 0.570 NA
Q26_2_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme cold incomeMiddle 0.901 0.689 1.178 0.445 NA
Q26_2_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme cold incomeHigh 0.737 0.401 1.360 0.327 NA
Q26_2_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme cold incomePrefer not to answer 0.748 0.459 1.221 0.244 NA
Q26_2_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme cold carerCarer 1.304 1.031 1.650 0.027 * NA
Q26_2_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme cold urbanicitySuburban 0.907 0.695 1.184 0.471 NA
Q26_2_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme cold urbanicitySemi-rural 0.893 0.609 1.314 0.565 NA
Q26_2_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme cold urbanicityRural 1.504 1.016 2.231 0.042 * NA
Q26_3_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Heavy rain age 1.007 0.999 1.015 0.082 † NA
Q26_3_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Heavy rain genderFemale 1.968 1.545 2.511 0.000 *** NA
Q26_3_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Heavy rain educSecondary 1.416 0.914 2.191 0.119 NA
Q26_3_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Heavy rain educTertiary 2.153 1.359 3.412 0.001 ** NA
Q26_3_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Heavy rain incomeMiddle 1.181 0.896 1.556 0.238 NA
Q26_3_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Heavy rain incomeHigh 1.395 0.732 2.682 0.314 NA
Q26_3_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Heavy rain incomePrefer not to answer 1.120 0.676 1.857 0.661 NA
Q26_3_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Heavy rain carerCarer 0.994 0.781 1.266 0.963 NA
Q26_3_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Heavy rain urbanicitySuburban 0.809 0.615 1.063 0.128 NA
Q26_3_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Heavy rain urbanicitySemi-rural 0.841 0.559 1.268 0.407 NA
Q26_3_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Heavy rain urbanicityRural 1.386 0.924 2.085 0.116 NA
Q26_4_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Strong wind age 1.007 0.999 1.015 0.075 † NA
Q26_4_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Strong wind genderFemale 1.539 1.212 1.958 0.000 *** NA
Q26_4_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Strong wind educSecondary 0.996 0.639 1.551 0.986 NA
Q26_4_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Strong wind educTertiary 1.456 0.916 2.313 0.111 NA
Q26_4_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Strong wind incomeMiddle 1.250 0.950 1.646 0.111 NA
Q26_4_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Strong wind incomeHigh 1.443 0.760 2.764 0.264 NA
Q26_4_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Strong wind incomePrefer not to answer 1.013 0.609 1.687 0.961 NA
Q26_4_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Strong wind carerCarer 1.142 0.899 1.452 0.277 NA
Q26_4_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Strong wind urbanicitySuburban 0.898 0.684 1.177 0.434 NA
Q26_4_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Strong wind urbanicitySemi-rural 0.916 0.612 1.374 0.671 NA
Q26_4_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Strong wind urbanicityRural 1.251 0.833 1.883 0.282 NA
Q26_5_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Snow age 1.013 1.005 1.021 0.001 ** NA
Q26_5_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Snow genderFemale 1.268 1.004 1.603 0.046 * NA
Q26_5_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Snow educSecondary 1.593 1.026 2.470 0.038 * NA
Q26_5_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Snow educTertiary 1.732 1.095 2.735 0.019 * NA
Q26_5_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Snow incomeMiddle 1.000 0.763 1.310 0.998 NA
Q26_5_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Snow incomeHigh 1.489 0.794 2.796 0.214 NA
Q26_5_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Snow incomePrefer not to answer 0.954 0.572 1.589 0.856 NA
Q26_5_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Snow carerCarer 0.879 0.696 1.110 0.278 NA
Q26_5_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Snow urbanicitySuburban 1.091 0.837 1.423 0.519 NA
Q26_5_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Snow urbanicitySemi-rural 1.064 0.714 1.587 0.760 NA
Q26_5_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Snow urbanicityRural 1.273 0.856 1.895 0.233 NA
Q26_6_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Ice age 1.012 1.004 1.020 0.002 ** NA
Q26_6_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Ice genderFemale 1.301 1.030 1.643 0.028 * NA
Q26_6_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Ice educSecondary 1.587 1.026 2.455 0.038 * NA
Q26_6_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Ice educTertiary 1.827 1.158 2.879 0.009 ** NA
Q26_6_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Ice incomeMiddle 0.967 0.736 1.271 0.811 NA
Q26_6_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Ice incomeHigh 1.434 0.747 2.752 0.278 NA
Q26_6_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Ice incomePrefer not to answer 0.757 0.457 1.254 0.280 NA
Q26_6_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Ice carerCarer 0.900 0.713 1.137 0.379 NA
Q26_6_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Ice urbanicitySuburban 1.108 0.850 1.445 0.447 NA
Q26_6_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Ice urbanicitySemi-rural 0.946 0.639 1.402 0.781 NA
Q26_6_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Ice urbanicityRural 1.785 1.196 2.667 0.005 ** NA
Q26_7_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Winter storms age 1.006 0.998 1.014 0.133 NA
Q26_7_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Winter storms genderFemale 1.710 1.350 2.167 0.000 *** NA
Q26_7_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Winter storms educSecondary 1.526 0.990 2.349 0.055 † NA
Q26_7_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Winter storms educTertiary 1.991 1.268 3.124 0.003 ** NA
Q26_7_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Winter storms incomeMiddle 0.992 0.756 1.300 0.952 NA
Q26_7_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Winter storms incomeHigh 0.968 0.512 1.837 0.920 NA
Q26_7_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Winter storms incomePrefer not to answer 1.045 0.627 1.743 0.866 NA
Q26_7_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Winter storms carerCarer 1.041 0.822 1.318 0.739 NA
Q26_7_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Winter storms urbanicitySuburban 0.897 0.685 1.174 0.427 NA
Q26_7_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Winter storms urbanicitySemi-rural 0.959 0.644 1.429 0.836 NA
Q26_7_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Winter storms urbanicityRural 1.486 0.999 2.218 0.052 † NA
Q26_8_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Summer storms age 1.012 1.004 1.020 0.003 ** NA
Q26_8_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Summer storms genderFemale 1.219 0.967 1.538 0.094 † NA
Q26_8_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Summer storms educSecondary 1.712 1.111 2.639 0.015 * NA
Q26_8_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Summer storms educTertiary 2.102 1.337 3.309 0.001 ** NA
Q26_8_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Summer storms incomeMiddle 1.050 0.802 1.376 0.721 NA
Q26_8_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Summer storms incomeHigh 0.907 0.477 1.729 0.765 NA
Q26_8_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Summer storms incomePrefer not to answer 0.809 0.495 1.321 0.397 NA
Q26_8_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Summer storms carerCarer 0.963 0.762 1.216 0.752 NA
Q26_8_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Summer storms urbanicitySuburban 0.987 0.755 1.289 0.922 NA
Q26_8_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Summer storms urbanicitySemi-rural 0.959 0.648 1.420 0.835 NA
Q26_8_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Summer storms urbanicityRural 1.570 1.060 2.328 0.025 * NA
Q26_9_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Wildfires age 1.009 1.001 1.017 0.029 * NA
Q26_9_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Wildfires genderFemale 1.217 0.955 1.552 0.113 NA
Q26_9_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Wildfires educSecondary 0.989 0.622 1.563 0.962 NA
Q26_9_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Wildfires educTertiary 1.422 0.876 2.298 0.152 NA
Q26_9_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Wildfires incomeMiddle 1.053 0.795 1.394 0.719 NA
Q26_9_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Wildfires incomeHigh 0.759 0.392 1.484 0.416 NA
Q26_9_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Wildfires incomePrefer not to answer 0.811 0.488 1.355 0.422 NA
Q26_9_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Wildfires carerCarer 0.979 0.766 1.251 0.866 NA
Q26_9_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Wildfires urbanicitySuburban 0.944 0.715 1.245 0.681 NA
Q26_9_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Wildfires urbanicitySemi-rural 1.359 0.903 2.060 0.144 NA
Q26_9_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Wildfires urbanicityRural 1.565 1.036 2.385 0.035 * NA
Q26_10_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Floods age 1.018 1.010 1.027 0.000 *** NA
Q26_10_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Floods genderFemale 1.673 1.318 2.126 0.000 *** NA
Q26_10_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Floods educSecondary 1.085 0.694 1.689 0.720 NA
Q26_10_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Floods educTertiary 1.502 0.943 2.388 0.085 † NA
Q26_10_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Floods incomeMiddle 1.225 0.931 1.610 0.147 NA
Q26_10_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Floods incomeHigh 1.148 0.599 2.218 0.679 NA
Q26_10_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Floods incomePrefer not to answer 1.512 0.903 2.546 0.117 NA
Q26_10_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Floods carerCarer 0.961 0.756 1.220 0.741 NA
Q26_10_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Floods urbanicitySuburban 1.052 0.801 1.384 0.714 NA
Q26_10_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Floods urbanicitySemi-rural 1.041 0.695 1.563 0.846 NA
Q26_10_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Floods urbanicityRural 1.280 0.861 1.909 0.224 NA
Q26_11_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Poor air quality age 1.022 1.014 1.030 0.000 *** NA
Q26_11_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Poor air quality genderFemale 1.809 1.429 2.293 0.000 *** NA
Q26_11_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Poor air quality educSecondary 1.375 0.875 2.160 0.167 NA
Q26_11_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Poor air quality educTertiary 1.526 0.951 2.448 0.080 † NA
Q26_11_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Poor air quality incomeMiddle 1.379 1.047 1.819 0.022 * NA
Q26_11_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Poor air quality incomeHigh 1.718 0.916 3.249 0.093 † NA
Q26_11_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Poor air quality incomePrefer not to answer 1.098 0.673 1.795 0.709 NA
Q26_11_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Poor air quality carerCarer 0.959 0.758 1.213 0.726 NA
Q26_11_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Poor air quality urbanicitySuburban 0.818 0.626 1.069 0.142 NA
Q26_11_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Poor air quality urbanicitySemi-rural 0.932 0.634 1.371 0.718 NA
Q26_11_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Poor air quality urbanicityRural 1.360 0.916 2.025 0.129 NA
Q26_12_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Landslides age 1.030 1.022 1.038 0.000 *** NA
Q26_12_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Landslides genderFemale 1.862 1.468 2.364 0.000 *** NA
Q26_12_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Landslides educSecondary 0.819 0.521 1.284 0.387 NA
Q26_12_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Landslides educTertiary 0.837 0.523 1.335 0.457 NA
Q26_12_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Landslides incomeMiddle 1.248 0.947 1.644 0.115 NA
Q26_12_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Landslides incomeHigh 1.257 0.670 2.373 0.478 NA
Q26_12_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Landslides incomePrefer not to answer 0.966 0.589 1.587 0.891 NA
Q26_12_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Landslides carerCarer 1.067 0.843 1.352 0.588 NA
Q26_12_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Landslides urbanicitySuburban 1.082 0.828 1.415 0.564 NA
Q26_12_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Landslides urbanicitySemi-rural 1.009 0.676 1.509 0.964 NA
Q26_12_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Landslides urbanicityRural 2.055 1.377 3.080 0.000 *** NA
Q26_13_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Sea swell age 1.021 1.012 1.029 0.000 *** NA
Q26_13_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Sea swell genderFemale 1.651 1.303 2.095 0.000 *** NA
Q26_13_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Sea swell educSecondary 1.131 0.728 1.757 0.583 NA
Q26_13_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Sea swell educTertiary 1.053 0.664 1.669 0.825 NA
Q26_13_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Sea swell incomeMiddle 1.164 0.886 1.529 0.275 NA
Q26_13_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Sea swell incomeHigh 1.353 0.704 2.620 0.366 NA
Q26_13_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Sea swell incomePrefer not to answer 1.049 0.641 1.719 0.848 NA
Q26_13_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Sea swell carerCarer 1.307 1.030 1.659 0.028 * NA
Q26_13_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Sea swell urbanicitySuburban 1.015 0.774 1.331 0.914 NA
Q26_13_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Sea swell urbanicitySemi-rural 0.953 0.642 1.415 0.810 NA
Q26_13_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Sea swell urbanicityRural 1.233 0.832 1.832 0.297 NA 
ordered_logit_results %>%
  filter(!is.na(odds_ratio), term != "(Intercept)") %>%
  mutate(outcome_label = paste0(outcome, ": ", str_trunc(question, 55))) %>%
  ggplot(aes(x = odds_ratio, y = fct_rev(term))) +
  geom_vline(xintercept = 1, linetype = "dashed") +
  geom_pointrange(aes(xmin = conf.low, xmax = conf.high)) +
  scale_x_log10() +
  facet_wrap(~ outcome_label) +
  labs(
    x = "Odds ratio, log scale",
    y = NULL,
    title = "Q26 ordered logit model estimates"
  ) +
  theme_minimal()

10 Compact model summary

compact_model_summary <- bind_rows(
  q32_logit_results %>% mutate(model_family = "Binary logit: Q32"),
  binary_logit_results %>% mutate(model_family = "Binary logit: Q34"),
  ordered_logit_results %>% mutate(model_family = "Ordered logit: Q26")
) %>%
  filter(!is.na(p.value), term != "(Intercept)") %>%
  group_by(model_family, outcome, question) %>%
  summarise(
    n_predictors_p_lt_05 = sum(p.value < .05, na.rm = TRUE),
    smallest_p = min(p.value, na.rm = TRUE),
    strongest_predictor = term[which.min(p.value)],
    strongest_predictor_or = odds_ratio[which.min(p.value)],
    .groups = "drop"
  ) %>%
  arrange(model_family, smallest_p)

compact_model_summary %>%
  mutate(
    smallest_p = round(smallest_p, 4),
    strongest_predictor_or = round(strongest_predictor_or, 3)
  ) %>%
  knitr::kable(caption = "Compact summary of strongest predictor per model.")
Compact summary of strongest predictor per model.
model_family outcome question n_predictors_p_lt_05 smallest_p strongest_predictor strongest_predictor_or
Binary logit: Q32 Q32_2_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Baby food for baby/babies in your household lasting at least three days 2 0.0000 carerCarer 3.660
Binary logit: Q32 Q32_3_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Pet food lasting at least three days 2 0.0009 age 0.986
Binary logit: Q32 Q32_1_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Essential medicines for household members who need them, lasting at least three days 4 0.0021 age 1.016
Binary logit: Q32 Q32_4_bin In an emergency causing interruption to utilities and other important services, do you have the following items at home? Please select “Not applicable” if this is not relevant to you or other members of your household. - Spare glasses or contact lenses 0 0.0513 incomeHigh 1.899
Binary logit: Q34 Q34_6_bin I cannot avoid travelling because of caring responsibilities for one or more vulnerable people 2 0.0000 carerCarer 5.574
Binary logit: Q34 Q34_12_bin I do not have the physical strength to secure and protect all objects and equipment on my property, and have no one to help me do so 3 0.0000 genderFemale 2.831
Binary logit: Q34 Q34_2_bin I cannot work from home to avoid travelling 1 0.0000 age 0.978
Binary logit: Q34 Q34_1_bin I cannot leave work early to avoid travelling 1 0.0000 age 0.977
Binary logit: Q34 Q34_4_bin I cannot avoid working in an exposed area during extreme heat or thunderstorms 2 0.0000 age 0.975
Binary logit: Q34 Q34_10_bin It is difficult for me or other household members to adopt the suggested measures because of mobility problems 1 0.0000 carerCarer 2.552
Binary logit: Q34 Q34_5_bin Difficulty leaving my home in the event of a flood or wildfire warning because I depend on others for transport 2 0.0000 age 0.975
Binary logit: Q34 Q34_13_bin Other response provided 2 0.0001 age 1.048
Binary logit: Q34 Q34_3_bin I cannot avoid doing outdoor tasks or physical exertion during extreme heat because of work responsibilities 2 0.0002 carerCarer 1.965
Binary logit: Q34 Q34_9_bin I would not feel comfortable cancelling trips or travel to visit friends or family if avoiding unnecessary travel due to bad weather was recommended 2 0.0026 carerCarer 2.313
Binary logit: Q34 Q34_11_bin I cannot or do not like to open windows at night to let cool air in 1 0.0209 carerCarer 1.541
Binary logit: Q34 Q34_8_bin Difficulty finding sheltered place to stay during heavy rain 1 0.0491 incomePrefer not to answer 1.813
Binary logit: Q34 Q34_7_bin Difficulty finding a cool place to stay during extreme heat 0 0.0510 urbanicityRural 0.575
Ordered logit: Q26 Q26_12_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Landslides 3 0.0000 age 1.030
Ordered logit: Q26 Q26_3_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Heavy rain 2 0.0000 genderFemale 1.968
Ordered logit: Q26 Q26_11_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Poor air quality 3 0.0000 age 1.022
Ordered logit: Q26 Q26_13_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Sea swell 3 0.0000 age 1.021
Ordered logit: Q26 Q26_10_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Floods 2 0.0000 age 1.018
Ordered logit: Q26 Q26_7_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Winter storms 2 0.0000 genderFemale 1.710
Ordered logit: Q26 Q26_4_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Strong wind 1 0.0004 genderFemale 1.539
Ordered logit: Q26 Q26_2_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme cold 3 0.0007 genderFemale 1.499
Ordered logit: Q26 Q26_5_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Snow 4 0.0012 age 1.013
Ordered logit: Q26 Q26_8_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Summer storms 4 0.0013 educTertiary 2.102
Ordered logit: Q26 Q26_6_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Ice 5 0.0024 age 1.012
Ordered logit: Q26 Q26_1_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Extreme heat 4 0.0083 educTertiary 1.902
Ordered logit: Q26 Q26_9_ord To what extent do you feel the following have become more common in Portugal over your lifetime? - Wildfires 2 0.0293 age 1.009

11 Objects saved in the knitting environment

The report creates the following main objects for follow-up inspection:

sessionInfo()
## R version 4.4.1 (2024-06-14 ucrt)
## Platform: x86_64-w64-mingw32/x64
## Running under: Windows 11 x64 (build 26200)
## 
## Matrix products: default
## 
## 
## locale:
## [1] LC_COLLATE=English_United Kingdom.utf8 
## [2] LC_CTYPE=English_United Kingdom.utf8   
## [3] LC_MONETARY=English_United Kingdom.utf8
## [4] LC_NUMERIC=C                           
## [5] LC_TIME=English_United Kingdom.utf8    
## 
## time zone: Europe/London
## tzcode source: internal
## 
## attached base packages:
## [1] stats     graphics  grDevices utils     datasets  methods   base     
## 
## other attached packages:
##  [1] tidyr_1.3.1   tibble_3.2.1  stringr_1.5.1 scales_1.3.0  readr_2.1.5  
##  [6] purrr_1.0.2   MASS_7.3-60.2 ggplot2_3.5.1 forcats_1.0.0 dplyr_1.1.4  
## [11] broom_1.0.6  
## 
## loaded via a namespace (and not attached):
##  [1] sass_0.4.9        utf8_1.2.4        generics_0.1.3    stringi_1.8.4    
##  [5] hms_1.1.3         digest_0.6.36     magrittr_2.0.3    evaluate_0.24.0  
##  [9] grid_4.4.1        fastmap_1.2.0     jsonlite_1.8.8    backports_1.5.0  
## [13] fansi_1.0.6       jquerylib_0.1.4   cli_3.6.6         rlang_1.2.0      
## [17] crayon_1.5.3      bit64_4.0.5       munsell_0.5.1     withr_3.0.0      
## [21] cachem_1.1.0      yaml_2.3.8        tools_4.4.1       parallel_4.4.1   
## [25] tzdb_0.4.0        colorspace_2.1-0  vctrs_0.6.5       R6_2.5.1         
## [29] lifecycle_1.0.4   bit_4.0.5         vroom_1.6.5       pkgconfig_2.0.3  
## [33] pillar_1.9.0      bslib_0.7.0       gtable_0.3.6      glue_1.8.1       
## [37] highr_0.11        xfun_0.52         tidyselect_1.2.1  rstudioapi_0.16.0
## [41] knitr_1.47        farver_2.1.2      htmltools_0.5.8.1 rmarkdown_2.27   
## [45] labeling_0.4.3    compiler_4.4.1