inst/shiny_tipsae/server_files/plots/results_raneff_dens_spat.R

In tipsae: Tools for Handling Indices and Proportions in Small Area Estimation

### Requested inputs ------

# none

### Creating plot ------

plot_dens_spat <- shiny::reactive({
  data_reff_s <- data.frame(reff = scale(res_model()$summary$raneff$spatial[, "mean"]))
  xlim_reff <- range(c(data_reff_s$reff, 2, -2))


  ggplot2::ggplot(data_reff_s, ggplot2::aes(x = reff)) +
    ggplot2::geom_function(fun = dnorm,  ggplot2::aes(color = "Standard normal")) +
    ggplot2::stat_density(ggplot2::aes(color = "Scaled random effects"),
                          geom = "line", position = "identity") +
    ggplot2::ylab("Density") + ggplot2::xlim(xlim_reff) +
    ggplot2::theme_bw(base_size = 15) + ggplot2::labs(color = "") +
    ggplot2::scale_color_manual(values = c(
      "Scaled random effects" = "black",
      "Standard normal" = "grey")) +
    ggplot2::xlab("Spatial random effect") +
    ggplot2::theme(aspect.ratio = 2/3)
  })

### Output: plot and save -----

output$dens_spat <- shiny::renderPlot({
  plot_dens_spat()
}, bg = "transparent")


output$download_dens_spat <- shiny::downloadHandler(
  filename = 'tipsae_dens_spat.RData',
  content = function(file) {
    tipsae_dens_spat <- plot_dens_spat()
    save(tipsae_dens_spat, file = file)
  }
)

output$save_pdf_dens_spat = shiny::downloadHandler(
  filename = "tipsae_dens_spat.pdf",
  content = function(file) {
    ggplot2::ggsave(file, plot = plot_dens_spat(), device = "pdf")
  }
)