--- title: "Análisis de desigualdades" author: "Carlos Castañeda-Orjuela" date: "abril 2016" output: html_document --- http://www.cdc.gov/cancer/npcr/tools/registryplus/lp.htm ```{r} setwd("~/Documents/Observatorio Nacional de Salud/Analisis/7 informe") library(readxl) library(RGoogleDocs) library(ggplot2) library(stringr) library(xlsx) # download.file("https://docs.google.com/spreadsheets/d/1962-OC70X16qq8wcocRRERmW6dgTQHaluCzF_no1RN4/edit?ts=56e9c6c5#gid=964737524", destfile = "parametros.xlsx", method = "curl") parametros <- read_excel("Indicadores 7Informe.xlsx", na = "ND")[1:34, ] var <- read_excel("Indicadores 7Informe.xlsx", sheet = 2, skip = 1)[1:41, ] parametros <- as.data.frame(parametros) parametros$itd_2013 <- as.numeric(parametros$itd_2013) adjust_perc <- function(parameter, expansor){ return(parameter / expansor) } adjust_par <- function(parameter, expansor){ return(parameter * expansor) } for(i in 3:41){ parametros[i] <- adjust_perc(parametros[i], var$multiplicador[i - 2]) parametros[i] <- adjust_par(parametros[i], var$multiplicadorMostrar[i - 2]) } nacional <- parametros[34, ] parametros <- parametros[1:33, ] ``` Parámetros en el contexto - Antioquia: Naranja - Bogotá: Aguamarina - Chocó: Verde - La Guajira: Morado - Tolima: Rojo - Valle del Cauca: Verde oscuro ```{r} # Función para generar gráfico por departamento plot_data <- function(parameter, dpto){ vector <- parametros[parameter] dptoData <- vector[parametros$departamento == dpto, ] colData <- nacional[parameter] boxplot(vector, horizontal = TRUE, col = "gray", outline = FALSE, staplewex = 1, width = 0.6, boxwex = 1.1, main = c(dpto, parameter), range = 0) abline(v = colData, col = "red", lwd = 3) points(dptoData, 1, col = "orange", pch = 19, cex = 2) } # Función para generar gráfico por indicador, varios departamentos en el mismo gráfico plot_multiData <- function(parameter, dptos, colores){ vector <- parametros[parameter] dptoData <- vector[parametros$departamento %in% (dptos), ] colData <- nacional[parameter] nombre <- var$indicador[var$`nombre corto` == parameter] boxplot(vector, horizontal = TRUE, col = "gray", outline = TRUE, staplewex = 1, width = 0.6, boxwex = 1.1, main = nombre, range = 0) abline(v = colData, col = "red", lwd = 3) points(dptoData, rep(1, 6), col = colores, pch = 19, cex = 2) } # for(dp in c("ANTIOQUIA", "BOGOTA", "CHOCO", "LA GUAJIRA", "TOLIMA", "VALLE DEL CAUCA")){ # for(i in 3:41){ # plot_data(names(parametros)[i], dp) # } # } for(i in 3:41){ plot_multiData(names(parametros)[i], c("ANTIOQUIA", "BOGOTA", "CHOCO", "LA GUAJIRA", "TOLIMA", "VALLE DEL CAUCA"), c("orange","lightblue", "green", "purple", "red", "darkgreen")) } ``` ```{r} dataexport <- rbind(parametros, nacional) write.xlsx(dataexport, "indicadores.xlsx", row.names = FALSE) ```