drinking_age_count <- sum(data$drinking_age == 1, na.rm = TRUE) non_drinking_age_count <- sum(data$drinking_age == 0, na.rm = TRUE) cat("Drinking Age Counts\n") cat ...