R/plot_q65_no_wrong_fishing_reserve.R

In Rare-Technology/HHS_Dashboard: PKG_TITLE

prep_q65_no_wrong_fishing_reserve <- function(.data){
  hhs_Q65 <- .data[,c("maa", "65_no_wrong_fishing_reserve")] %>%
    dplyr::filter(`65_no_wrong_fishing_reserve` %in% c(0:10))
  Q65_length <-
    tapply(hhs_Q65$`65_no_wrong_fishing_reserve`,
           hhs_Q65$maa,
           length)
  Q65_length <- as.vector(Q65_length)
  Q65_mean <-
    data.frame(avg = tapply(as.numeric(
      as.character(hhs_Q65$`65_no_wrong_fishing_reserve`)
    ),
    hhs_Q65$maa, mean) / 10)
  Q65_summary_bind <- cbind(N = Q65_length, round(Q65_mean, 3) * 100)
  Q65_summary <-
    rbind(Q65_summary_bind, "Mean ± SE" = c(
      sum(Q65_summary_bind$N),
      compute_summary_line(Q65_summary_bind$avg, 1)
    ))
  colnames(Q65_summary) <- c("N", "Proportion (%)")
  Q65_summary <-tibble::rownames_to_column(Q65_summary, "MA name")
  
  Q65 <- clean_plot_data(Q65_summary)
  Q65
}

plot_q65_no_wrong_fishing_reserve <- function(.data, ...){

  .data_plot <- prep_q65_no_wrong_fishing_reserve(.data)
  
  p <- plot_horiz_bar(
    .data_plot,
    title = "Proportion of community members who think is wrong to fish in the reserve"
  )

  result <- list(
    plot = p,
    data = .data_plot
  )
}