---
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)
```

