R/plot_q51_fishers_permission.R

In Rare-Technology/HHS_Dashboard: PKG_TITLE

prep_q51_fishers_permission <- function(.data){
  hhs_Q51a <- .data[,c("maa", "51a_fishers_gear_not_permitted")] %>%
    dplyr::filter(`51a_fishers_gear_not_permitted` %in% c(0:10)) %>% 
    droplevels()
  ## Summary 
  Q51a_length <-
    tapply(hhs_Q51a$`51a_fishers_gear_not_permitted`,
           hhs_Q51a$maa,
           length)
  Q51a_length <- as.vector(Q51a_length)
  Q51a_mean <-
    as.data.frame(
      tapply(
        hhs_Q51a$`51a_fishers_gear_not_permitted`,
        hhs_Q51a$maa,
        mean
      ) / 10
    )
  Q51a_summary_bind <-
    cbind(N = Q51a_length, round(Q51a_mean, 3) * 100)
  colnames(Q51a_summary_bind) <- c("N", "Unapproved gear")
  Q51a_summary <- rbind(Q51a_summary_bind,
                        "Mean ± SE" = c(
                          sum(Q51a_summary_bind$N),
                          compute_summary_line(Q51a_summary_bind[[2]], 1)
                        ))
  Q51a_summary <-tibble::rownames_to_column(Q51a_summary, "MA name")
  
  ### Q51b Frequency of observed fishing in reserve ####
  hhs_Q51b <- .data[,c("maa", "51b_fishers_reserves")] %>%
    dplyr::filter(`51b_fishers_reserves` %in% c(0:10)) %>%
    droplevels()
  Q51b_length <-
    tapply(hhs_Q51b$`51b_fishers_reserves`,
           hhs_Q51b$maa,
           length)
  Q51b_length <- as.vector(Q51b_length)
  Q51b_mean <-
    data.frame(avg = tapply(
      hhs_Q51b$`51b_fishers_reserves`,
      hhs_Q51b$maa,
      mean
    ) / 10)
  Q51b_summary_bind <-
    cbind(N = Q51b_length, round(Q51b_mean, 3) * 100)
  Q51b_summary <-
    rbind(Q51b_summary_bind, "Mean ± SE" = c(
      sum(Q51b_summary_bind$N),
      compute_summary_line(Q51b_summary_bind$avg, 1)
    ))
  colnames(Q51b_summary) <- c("N", "Fishing in reserve")
  Q51b_summary <-tibble::rownames_to_column(Q51b_summary, "MA name")
  
  ### Q51c Frequency of observed unpermitted fishing in MA ####
  hhs_Q51c <- .data[,c("maa", "51c_fishers_ma_area")] %>%
    dplyr::filter(`51c_fishers_ma_area` %in% c(0:10)) %>%
    droplevels()
  Q51c_length <-
    tapply(hhs_Q51c$`51c_fishers_ma_area`,
           hhs_Q51c$maa,
           length)
  Q51c_length <- as.vector(Q51c_length)
  Q51c_mean <-
    data.frame(avg = tapply(
      hhs_Q51c$`51c_fishers_ma_area`,
      hhs_Q51c$maa,
      mean
    ) / 10)
  Q51c_summary_bind <-
    cbind(N = Q51c_length, round(Q51c_mean, 3) * 100)
  colnames(Q51c_summary_bind) <-
    c("N", "Fishing without permission")
  Q51c_summary <-
    rbind(Q51c_summary_bind, "Mean ± SE" = c(
      sum(Q51c_summary_bind$N),
      compute_summary_line(Q51c_summary_bind[[2]], 1)
    ))
  Q51c_summary <-tibble::rownames_to_column(Q51c_summary, "MA name")
  
  ### Combine Q51a, Q51b and Q51c ####
  Q51a_c_summary <-
    plyr::join_all(
      list(Q51a_summary, Q51b_summary, Q51c_summary),
      by = "MA name",
      type = "left"
    )
  
  #pivot table
  Q51a_c_summary_long <-
    as.data.frame(
      Q51a_c_summary %>% tidyr::pivot_longer(
        cols = c(
          "Unapproved gear",
          "Fishing in reserve",
          "Fishing without permission"
        ),
        names_to = "key",
        values_to = "Proportion (%)"
      )
    )
  #Fix Ns
  Q51a_c_summary_long$N <-
    as.data.frame(
      tidyr::pivot_longer(
        Q51a_c_summary,
        cols = c("N", "N", "N"),
        names_repair = "unique",
        names_to = "No",
        values_to = "N"
      )
    )$N
  
  Q51a_c <- clean_plot_data(Q51a_c_summary_long)
  Q51a_c
}

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

  .data_plot <- prep_q51_fishers_permission(.data)
 
  p <- plot_horiz_bar(
    .data_plot,
    title = "Perceived frequency of observing others violating regulations ",
    facet_var = key
  )
  
  result <- list(
    plot = p,
    data = .data_plot
  )
}