Rangos para temperaturas maximas.

tmmx

Equipo Data Science

DataIntelligence
29-12-2021 
library('DT')
library('tidyverse')
library('expss')

1 Introducción

Rankearemos las siguientes variables con la metodología que detallaremos más adelante:

  • tmmx: temperatura máxima.

2 El código comunal

2.0.1 Códigos comunales

codigos_comunales_actuales <- readxl::read_xls("CUT_2018_v04.xls")
names(codigos_comunales_actuales)[6] <- "Comuna"
codigos_comunales_actuales$Comuna  <- as.numeric(codigos_comunales_actuales$Comuna)

codigos_comunales_actuales
## # A tibble: 346 x 7
##    `Código Región` `Nombre Región` `Abreviatura Región` `Código Provincia`
##    <chr>           <chr>           <chr>                <chr>             
##  1 01              Tarapacá        TPCA                 011               
##  2 01              Tarapacá        TPCA                 011               
##  3 01              Tarapacá        TPCA                 014               
##  4 01              Tarapacá        TPCA                 014               
##  5 01              Tarapacá        TPCA                 014               
##  6 01              Tarapacá        TPCA                 014               
##  7 01              Tarapacá        TPCA                 014               
##  8 02              Antofagasta     ANTOF                021               
##  9 02              Antofagasta     ANTOF                021               
## 10 02              Antofagasta     ANTOF                021               
## # ... with 336 more rows, and 3 more variables: Nombre Provincia <chr>,
## #   Comuna <dbl>, Nombre Comuna <chr>
# datatable(codigos_comunales_actuales, extensions = 'Buttons', escape = FALSE, rownames = TRUE,
#           options = list(dom = 'Bfrtip',
#           buttons = list('colvis', list(extend = 'collection',
#           buttons = list(
#           list(extend='copy'),
#           list(extend='excel',
#             filename = 'tabla'),
#           list(extend='pdf',
#             filename= 'tabla')),
#           text = 'Download')), scrollX = TRUE))

Leamos los datos de temperaturas máximas imputadas (pero aún no corregidas por 0.1):

rangos GEE

Tenemos data para 21 años (2000-2020), por lo que debemos generar 20 intervalos de variaciones.

3 Temperaturas máximas (tmmx)

Mostraremos el proceso completo para Enero; para los demás, sólo las sd y mes y en otra los rangos y sus cotas.

3.1 Leemos la tabla imputada de temperaturas máximas

df_0 <- readRDS("TERRACLIMATE_MEDIA_tmmx_imp.rds")
names(df_0)[253] <- 'Comuna'
df_1 <- df_0[,c(253,1:252)]

3.2 Corregimos por 0.1

for (i in 2:253){
    df_1[,i] <- (df_1[,i]*.1)
}
options(digits=20)

4 Enero

minimas <- tibble()
df_min <- df_1[,c(1,seq(2, 253, 12))]
minimas <- rbind(minimas,df_min)
minimas <- minimas[order(minimas$Comuna),]
minimas$sd =  round(apply(minimas[,-1], 1, sd),3)
minimas$mean = rowMeans(minimas[,c(-1)])
names(minimas)[23] <- "sd_Enero"
names(minimas)[24] <- "mean_Enero"
minimas <- minimas[,-c(2:22,25)]
max <- max(minimas$sd_Enero)
min <- min(minimas$sd_Enero)
incremento <- (max-min)/5
dos <- min + incremento
tres <- dos  + incremento
cuatro <- tres  + incremento
quinto <- cuatro  + incremento
sexto <- quinto  + incremento + 0.001
rango <- c(1, 2, 3, 4, 5)
uno <- round(min, digits = 3)
dos <- round(dos, digits = 3)
primer_rango <- paste(uno, "-", dos)
dos <- round(dos, digits = 3)
tres <- round(tres, digits = 3)
segundo_rango <- paste(dos, "-", tres)
tres <- round(tres, digits = 3)
cuatro <- round(cuatro, digits = 3)
tercer_rango <- paste(tres, "-", cuatro)
cuatro <- round(cuatro, digits = 3)
quinto <- round(quinto, digits = 3)
cuarto_rango <- paste(cuatro, "-", quinto)
quinto <- round(quinto, digits = 3)
sexto <- round(sexto, digits = 3)
quinto_rango <- paste(quinto, "-", sexto)
rango <- c(1,2,3,4,5)
intervalos <- c(primer_rango, segundo_rango , tercer_rango , cuarto_rango , quinto_rango )
df <- data.frame(rango, intervalos)
receptaculo <- minimas
receptaculo$quintiles_rango_n <-
ifelse(receptaculo$sd_Enero >= uno & receptaculo$sd_Enero < dos , 1, ifelse(receptaculo$sd_Enero >= dos   & receptaculo$sd_Enero <  tres , 2, ifelse(receptaculo$sd_Enero >= tres   & receptaculo$sd_Enero < cuatro, 3, ifelse(receptaculo$sd_Enero >= cuatro    & receptaculo$sd_Enero < quinto , 4 , ifelse(receptaculo$sd_Enero >= quinto  & receptaculo$sd_Enero <  sexto,5,"")))))
receptaculo$quintiles_rango_n_valores <-
ifelse(receptaculo$quintiles_rango_n == 1 , df[1,2], ifelse(receptaculo$quintiles_rango_n == 2 , df[2,2], ifelse(receptaculo$quintiles_rango_n == 3, df[3,2], ifelse(receptaculo$quintiles_rango_n == 4, df[4,2] , ifelse(receptaculo$quintiles_rango_n == 5,df[5,2],"")))))

4.0.1 Enero

r_enero <- receptaculo
r_enero$Comuna <- as.double(r_enero$Comuna) 
r_enero$quintiles_rango_n <- as.double(r_enero$quintiles_rango_n) 
names(r_enero)[4] <- 'quintiles_rango_n_Enero'
names(r_enero)[5] <- 'quintiles_rango_n_valores_Enero'
r_Enero <- r_enero

5 Febrero

minimas <- tibble()
df_min <- df_1[,c(1,seq(3, 253, 12))]
minimas <- rbind(minimas,df_min)
minimas <- minimas[order(minimas$Comuna),]
minimas$sd =  round(apply(minimas[,-1], 1, sd),3)
minimas$mean = rowMeans(minimas[,c(-1)])
names(minimas)[23] <- "sd_Enero"
names(minimas)[24] <- "mean_Enero"
minimas <- minimas[,-c(2:22,25)]
max <- max(minimas$sd_Enero)
min <- min(minimas$sd_Enero)
incremento <- (max-min)/5
dos <- min + incremento
tres <- dos  + incremento
cuatro <- tres  + incremento
quinto <- cuatro  + incremento
sexto <- quinto  + incremento + 0.001
rango <- c(1, 2, 3, 4, 5)
uno <- round(min, digits = 3)
dos <- round(dos, digits = 3)
primer_rango <- paste(uno, "-", dos)
dos <- round(dos, digits = 3)
tres <- round(tres, digits = 3)
segundo_rango <- paste(dos, "-", tres)
tres <- round(tres, digits = 3)
cuatro <- round(cuatro, digits = 3)
tercer_rango <- paste(tres, "-", cuatro)
cuatro <- round(cuatro, digits = 3)
quinto <- round(quinto, digits = 3)
cuarto_rango <- paste(cuatro, "-", quinto)
quinto <- round(quinto, digits = 3)
sexto <- round(sexto, digits = 3)
quinto_rango <- paste(quinto, "-", sexto)
rango <- c(1,2,3,4,5)
intervalos <- c(primer_rango, segundo_rango , tercer_rango , cuarto_rango , quinto_rango )
df <- data.frame(rango, intervalos)
receptaculo <- minimas
receptaculo$quintiles_rango_n <-
ifelse(receptaculo$sd_Enero >= uno & receptaculo$sd_Enero < dos , 1, ifelse(receptaculo$sd_Enero >= dos   & receptaculo$sd_Enero <  tres , 2, ifelse(receptaculo$sd_Enero >= tres   & receptaculo$sd_Enero < cuatro, 3, ifelse(receptaculo$sd_Enero >= cuatro    & receptaculo$sd_Enero < quinto , 4 , ifelse(receptaculo$sd_Enero >= quinto  & receptaculo$sd_Enero <  sexto,5,"")))))
receptaculo$quintiles_rango_n_valores <-
ifelse(receptaculo$quintiles_rango_n == 1 , df[1,2], ifelse(receptaculo$quintiles_rango_n == 2 , df[2,2], ifelse(receptaculo$quintiles_rango_n == 3, df[3,2], ifelse(receptaculo$quintiles_rango_n == 4, df[4,2] , ifelse(receptaculo$quintiles_rango_n == 5,df[5,2],"")))))

5.0.1 Febrero

r_enero <- receptaculo
r_enero$Comuna <- as.double(r_enero$Comuna) 
r_enero$quintiles_rango_n <- as.double(r_enero$quintiles_rango_n) 
names(r_enero)[2] <- 'sd_Febrero'
names(r_enero)[3] <- 'mean_Febrero'
names(r_enero)[4] <- 'quintiles_rango_n_Febrero'
names(r_enero)[5] <- 'quintiles_rango_n_valores_Febrero'
r_febrero <- r_enero

6 Marzo

minimas <- tibble()
df_min <- df_1[,c(1,seq(4, 253, 12))]
minimas <- rbind(minimas,df_min)
minimas <- minimas[order(minimas$Comuna),]
minimas$sd =  round(apply(minimas[,-1], 1, sd),3)
minimas$mean = rowMeans(minimas[,c(-1)])
names(minimas)[23] <- "sd_Enero"
names(minimas)[24] <- "mean_Enero"
minimas <- minimas[,-c(2:22,25)]
max <- max(minimas$sd_Enero)
min <- min(minimas$sd_Enero)
incremento <- (max-min)/5
dos <- min + incremento
tres <- dos  + incremento
cuatro <- tres  + incremento
quinto <- cuatro  + incremento
sexto <- quinto  + incremento + 0.001
rango <- c(1, 2, 3, 4, 5)
uno <- round(min, digits = 3)
dos <- round(dos, digits = 3)
primer_rango <- paste(uno, "-", dos)
dos <- round(dos, digits = 3)
tres <- round(tres, digits = 3)
segundo_rango <- paste(dos, "-", tres)
tres <- round(tres, digits = 3)
cuatro <- round(cuatro, digits = 3)
tercer_rango <- paste(tres, "-", cuatro)
cuatro <- round(cuatro, digits = 3)
quinto <- round(quinto, digits = 3)
cuarto_rango <- paste(cuatro, "-", quinto)
quinto <- round(quinto, digits = 3)
sexto <- round(sexto, digits = 3)
quinto_rango <- paste(quinto, "-", sexto)
rango <- c(1,2,3,4,5)
intervalos <- c(primer_rango, segundo_rango , tercer_rango , cuarto_rango , quinto_rango )
df <- data.frame(rango, intervalos)
receptaculo <- minimas
receptaculo$quintiles_rango_n <-
ifelse(receptaculo$sd_Enero >= uno & receptaculo$sd_Enero < dos , 1, ifelse(receptaculo$sd_Enero >= dos   & receptaculo$sd_Enero <  tres , 2, ifelse(receptaculo$sd_Enero >= tres   & receptaculo$sd_Enero < cuatro, 3, ifelse(receptaculo$sd_Enero >= cuatro    & receptaculo$sd_Enero < quinto , 4 , ifelse(receptaculo$sd_Enero >= quinto  & receptaculo$sd_Enero <  sexto,5,"")))))
receptaculo$quintiles_rango_n_valores <-
ifelse(receptaculo$quintiles_rango_n == 1 , df[1,2], ifelse(receptaculo$quintiles_rango_n == 2 , df[2,2], ifelse(receptaculo$quintiles_rango_n == 3, df[3,2], ifelse(receptaculo$quintiles_rango_n == 4, df[4,2] , ifelse(receptaculo$quintiles_rango_n == 5,df[5,2],"")))))

6.0.1 Marzo

r_enero <- receptaculo
r_enero$Comuna <- as.double(r_enero$Comuna) 
r_enero$quintiles_rango_n <- as.double(r_enero$quintiles_rango_n) 
names(r_enero)[2] <- 'sd_Marzo'
names(r_enero)[3] <- 'mean_Marzo'
names(r_enero)[4] <- 'quintiles_rango_n_Marzo'
names(r_enero)[5] <- 'quintiles_rango_n_valores_Marzo'
r_marzo <- r_enero

7 Abril

minimas <- tibble()
df_min <- df_1[,c(1,seq(5, 253, 12))]
minimas <- rbind(minimas,df_min)
minimas <- minimas[order(minimas$Comuna),]
minimas$sd =  round(apply(minimas[,-1], 1, sd),3)
minimas$mean = rowMeans(minimas[,c(-1)])
names(minimas)[23] <- "sd_Enero"
names(minimas)[24] <- "mean_Enero"
minimas <- minimas[,-c(2:22,25)]
max <- max(minimas$sd_Enero)
min <- min(minimas$sd_Enero)
incremento <- (max-min)/5
dos <- min + incremento
tres <- dos  + incremento
cuatro <- tres  + incremento
quinto <- cuatro  + incremento
sexto <- quinto  + incremento + 0.001
rango <- c(1, 2, 3, 4, 5)
uno <- round(min, digits = 3)
dos <- round(dos, digits = 3)
primer_rango <- paste(uno, "-", dos)
dos <- round(dos, digits = 3)
tres <- round(tres, digits = 3)
segundo_rango <- paste(dos, "-", tres)
tres <- round(tres, digits = 3)
cuatro <- round(cuatro, digits = 3)
tercer_rango <- paste(tres, "-", cuatro)
cuatro <- round(cuatro, digits = 3)
quinto <- round(quinto, digits = 3)
cuarto_rango <- paste(cuatro, "-", quinto)
quinto <- round(quinto, digits = 3)
sexto <- round(sexto, digits = 3)
quinto_rango <- paste(quinto, "-", sexto)
rango <- c(1,2,3,4,5)
intervalos <- c(primer_rango, segundo_rango , tercer_rango , cuarto_rango , quinto_rango )
df <- data.frame(rango, intervalos)
receptaculo <- minimas
receptaculo$quintiles_rango_n <-
ifelse(receptaculo$sd_Enero >= uno & receptaculo$sd_Enero < dos , 1, ifelse(receptaculo$sd_Enero >= dos   & receptaculo$sd_Enero <  tres , 2, ifelse(receptaculo$sd_Enero >= tres   & receptaculo$sd_Enero < cuatro, 3, ifelse(receptaculo$sd_Enero >= cuatro    & receptaculo$sd_Enero < quinto , 4 , ifelse(receptaculo$sd_Enero >= quinto  & receptaculo$sd_Enero <  sexto,5,"")))))
receptaculo$quintiles_rango_n_valores <-
ifelse(receptaculo$quintiles_rango_n == 1 , df[1,2], ifelse(receptaculo$quintiles_rango_n == 2 , df[2,2], ifelse(receptaculo$quintiles_rango_n == 3, df[3,2], ifelse(receptaculo$quintiles_rango_n == 4, df[4,2] , ifelse(receptaculo$quintiles_rango_n == 5,df[5,2],"")))))

7.0.1 Abril

r_enero <- receptaculo
r_enero$Comuna <- as.double(r_enero$Comuna) 
r_enero$quintiles_rango_n <- as.double(r_enero$quintiles_rango_n) 
names(r_enero)[2] <- 'sd_Abril'
names(r_enero)[3] <- 'mean_Abril'
names(r_enero)[4] <- 'quintiles_rango_n_Abril'
names(r_enero)[5] <- 'quintiles_rango_n_valores_Abril'
r_abril <- r_enero

8 Mayo

minimas <- tibble()
df_min <- df_1[,c(1,seq(6, 253, 12))]
minimas <- rbind(minimas,df_min)
minimas <- minimas[order(minimas$Comuna),]
minimas$sd =  round(apply(minimas[,-1], 1, sd),3)
minimas$mean = rowMeans(minimas[,c(-1)])
names(minimas)[23] <- "sd_Enero"
names(minimas)[24] <- "mean_Enero"
minimas <- minimas[,-c(2:22,25)]
max <- max(minimas$sd_Enero)
min <- min(minimas$sd_Enero)
incremento <- (max-min)/5
dos <- min + incremento
tres <- dos  + incremento
cuatro <- tres  + incremento
quinto <- cuatro  + incremento
sexto <- quinto  + incremento + 0.001
rango <- c(1, 2, 3, 4, 5)
uno <- round(min, digits = 3)
dos <- round(dos, digits = 3)
primer_rango <- paste(uno, "-", dos)
dos <- round(dos, digits = 3)
tres <- round(tres, digits = 3)
segundo_rango <- paste(dos, "-", tres)
tres <- round(tres, digits = 3)
cuatro <- round(cuatro, digits = 3)
tercer_rango <- paste(tres, "-", cuatro)
cuatro <- round(cuatro, digits = 3)
quinto <- round(quinto, digits = 3)
cuarto_rango <- paste(cuatro, "-", quinto)
quinto <- round(quinto, digits = 3)
sexto <- round(sexto, digits = 3)
quinto_rango <- paste(quinto, "-", sexto)
rango <- c(1,2,3,4,5)
intervalos <- c(primer_rango, segundo_rango , tercer_rango , cuarto_rango , quinto_rango )
df <- data.frame(rango, intervalos)
receptaculo <- minimas
receptaculo$quintiles_rango_n <-
ifelse(receptaculo$sd_Enero >= uno & receptaculo$sd_Enero < dos , 1, ifelse(receptaculo$sd_Enero >= dos   & receptaculo$sd_Enero <  tres , 2, ifelse(receptaculo$sd_Enero >= tres   & receptaculo$sd_Enero < cuatro, 3, ifelse(receptaculo$sd_Enero >= cuatro    & receptaculo$sd_Enero < quinto , 4 , ifelse(receptaculo$sd_Enero >= quinto  & receptaculo$sd_Enero <  sexto,5,"")))))
receptaculo$quintiles_rango_n_valores <-
ifelse(receptaculo$quintiles_rango_n == 1 , df[1,2], ifelse(receptaculo$quintiles_rango_n == 2 , df[2,2], ifelse(receptaculo$quintiles_rango_n == 3, df[3,2], ifelse(receptaculo$quintiles_rango_n == 4, df[4,2] , ifelse(receptaculo$quintiles_rango_n == 5,df[5,2],"")))))

8.0.1 Mayo

r_enero <- receptaculo
r_enero$Comuna <- as.double(r_enero$Comuna) 
r_enero$quintiles_rango_n <- as.double(r_enero$quintiles_rango_n) 
names(r_enero)[2] <- 'sd_Mayo'
names(r_enero)[3] <- 'mean_Mayo'
names(r_enero)[4] <- 'quintiles_rango_n_Mayo'
names(r_enero)[5] <- 'quintiles_rango_n_valores_Mayo'
r_mayo <- r_enero

9 Junio

minimas <- tibble()
df_min <- df_1[,c(1,seq(7, 253, 12))]
minimas <- rbind(minimas,df_min)
minimas <- minimas[order(minimas$Comuna),]
minimas$sd =  round(apply(minimas[,-1], 1, sd),3)
minimas$mean = rowMeans(minimas[,c(-1)])
names(minimas)[23] <- "sd_Enero"
names(minimas)[24] <- "mean_Enero"
minimas <- minimas[,-c(2:22,25)]
max <- max(minimas$sd_Enero)
min <- min(minimas$sd_Enero)
incremento <- (max-min)/5
dos <- min + incremento
tres <- dos  + incremento
cuatro <- tres  + incremento
quinto <- cuatro  + incremento
sexto <- quinto  + incremento + 0.001
rango <- c(1, 2, 3, 4, 5)
uno <- round(min, digits = 3)
dos <- round(dos, digits = 3)
primer_rango <- paste(uno, "-", dos)
dos <- round(dos, digits = 3)
tres <- round(tres, digits = 3)
segundo_rango <- paste(dos, "-", tres)
tres <- round(tres, digits = 3)
cuatro <- round(cuatro, digits = 3)
tercer_rango <- paste(tres, "-", cuatro)
cuatro <- round(cuatro, digits = 3)
quinto <- round(quinto, digits = 3)
cuarto_rango <- paste(cuatro, "-", quinto)
quinto <- round(quinto, digits = 3)
sexto <- round(sexto, digits = 3)
quinto_rango <- paste(quinto, "-", sexto)
rango <- c(1,2,3,4,5)
intervalos <- c(primer_rango, segundo_rango , tercer_rango , cuarto_rango , quinto_rango )
df <- data.frame(rango, intervalos)
receptaculo <- minimas
receptaculo$quintiles_rango_n <-
ifelse(receptaculo$sd_Enero >= uno & receptaculo$sd_Enero < dos , 1, ifelse(receptaculo$sd_Enero >= dos   & receptaculo$sd_Enero <  tres , 2, ifelse(receptaculo$sd_Enero >= tres   & receptaculo$sd_Enero < cuatro, 3, ifelse(receptaculo$sd_Enero >= cuatro    & receptaculo$sd_Enero < quinto , 4 , ifelse(receptaculo$sd_Enero >= quinto  & receptaculo$sd_Enero <  sexto,5,"")))))
receptaculo$quintiles_rango_n_valores <-
ifelse(receptaculo$quintiles_rango_n == 1 , df[1,2], ifelse(receptaculo$quintiles_rango_n == 2 , df[2,2], ifelse(receptaculo$quintiles_rango_n == 3, df[3,2], ifelse(receptaculo$quintiles_rango_n == 4, df[4,2] , ifelse(receptaculo$quintiles_rango_n == 5,df[5,2],"")))))

9.0.1 Junio

r_enero <- receptaculo
r_enero$Comuna <- as.double(r_enero$Comuna) 
r_enero$quintiles_rango_n <- as.double(r_enero$quintiles_rango_n) 
names(r_enero)[2] <- 'sd_Junio'
names(r_enero)[3] <- 'mean_Junio'
names(r_enero)[4] <- 'quintiles_rango_n_Junio'
names(r_enero)[5] <- 'quintiles_rango_n_valores_Junio'
r_junio <- r_enero

10 Julio

minimas <- tibble()
df_min <- df_1[,c(1,seq(8, 253, 12))]
minimas <- rbind(minimas,df_min)
minimas <- minimas[order(minimas$Comuna),]
minimas$sd =  round(apply(minimas[,-1], 1, sd),3)
minimas$mean = rowMeans(minimas[,c(-1)])
names(minimas)[23] <- "sd_Enero"
names(minimas)[24] <- "mean_Enero"
minimas <- minimas[,-c(2:22,25)]
max <- max(minimas$sd_Enero)
min <- min(minimas$sd_Enero)
incremento <- (max-min)/5
dos <- min + incremento
tres <- dos  + incremento
cuatro <- tres  + incremento
quinto <- cuatro  + incremento
sexto <- quinto  + incremento + 0.001
rango <- c(1, 2, 3, 4, 5)
uno <- round(min, digits = 3)
dos <- round(dos, digits = 3)
primer_rango <- paste(uno, "-", dos)
dos <- round(dos, digits = 3)
tres <- round(tres, digits = 3)
segundo_rango <- paste(dos, "-", tres)
tres <- round(tres, digits = 3)
cuatro <- round(cuatro, digits = 3)
tercer_rango <- paste(tres, "-", cuatro)
cuatro <- round(cuatro, digits = 3)
quinto <- round(quinto, digits = 3)
cuarto_rango <- paste(cuatro, "-", quinto)
quinto <- round(quinto, digits = 3)
sexto <- round(sexto, digits = 3)
quinto_rango <- paste(quinto, "-", sexto)
rango <- c(1,2,3,4,5)
intervalos <- c(primer_rango, segundo_rango , tercer_rango , cuarto_rango , quinto_rango )
df <- data.frame(rango, intervalos)
receptaculo <- minimas
receptaculo$quintiles_rango_n <-
ifelse(receptaculo$sd_Enero >= uno & receptaculo$sd_Enero < dos , 1, ifelse(receptaculo$sd_Enero >= dos   & receptaculo$sd_Enero <  tres , 2, ifelse(receptaculo$sd_Enero >= tres   & receptaculo$sd_Enero < cuatro, 3, ifelse(receptaculo$sd_Enero >= cuatro    & receptaculo$sd_Enero < quinto , 4 , ifelse(receptaculo$sd_Enero >= quinto  & receptaculo$sd_Enero <  sexto,5,"")))))
receptaculo$quintiles_rango_n_valores <-
ifelse(receptaculo$quintiles_rango_n == 1 , df[1,2], ifelse(receptaculo$quintiles_rango_n == 2 , df[2,2], ifelse(receptaculo$quintiles_rango_n == 3, df[3,2], ifelse(receptaculo$quintiles_rango_n == 4, df[4,2] , ifelse(receptaculo$quintiles_rango_n == 5,df[5,2],"")))))

10.0.1 Julio

r_enero <- receptaculo
r_enero$Comuna <- as.double(r_enero$Comuna) 
r_enero$quintiles_rango_n <- as.double(r_enero$quintiles_rango_n) 
names(r_enero)[2] <- 'sd_Julio'
names(r_enero)[3] <- 'mean_Julio'
names(r_enero)[4] <- 'quintiles_rango_n_Julio'
names(r_enero)[5] <- 'quintiles_rango_n_valores_Julio'
r_julio <- r_enero

11 Agosto

minimas <- tibble()
df_min <- df_1[,c(1,seq(9, 253, 12))]
minimas <- rbind(minimas,df_min)
minimas <- minimas[order(minimas$Comuna),]
minimas$sd =  round(apply(minimas[,-1], 1, sd),3)
minimas$mean = rowMeans(minimas[,c(-1)])
names(minimas)[23] <- "sd_Enero"
names(minimas)[24] <- "mean_Enero"
minimas <- minimas[,-c(2:22,25)]
max <- max(minimas$sd_Enero)
min <- min(minimas$sd_Enero)
incremento <- (max-min)/5
dos <- min + incremento
tres <- dos  + incremento
cuatro <- tres  + incremento
quinto <- cuatro  + incremento
sexto <- quinto  + incremento + 0.001
rango <- c(1, 2, 3, 4, 5)
uno <- round(min, digits = 3)
dos <- round(dos, digits = 3)
primer_rango <- paste(uno, "-", dos)
dos <- round(dos, digits = 3)
tres <- round(tres, digits = 3)
segundo_rango <- paste(dos, "-", tres)
tres <- round(tres, digits = 3)
cuatro <- round(cuatro, digits = 3)
tercer_rango <- paste(tres, "-", cuatro)
cuatro <- round(cuatro, digits = 3)
quinto <- round(quinto, digits = 3)
cuarto_rango <- paste(cuatro, "-", quinto)
quinto <- round(quinto, digits = 3)
sexto <- round(sexto, digits = 3)
quinto_rango <- paste(quinto, "-", sexto)
rango <- c(1,2,3,4,5)
intervalos <- c(primer_rango, segundo_rango , tercer_rango , cuarto_rango , quinto_rango )
df <- data.frame(rango, intervalos)
receptaculo <- minimas
receptaculo$quintiles_rango_n <-
ifelse(receptaculo$sd_Enero >= uno & receptaculo$sd_Enero < dos , 1, ifelse(receptaculo$sd_Enero >= dos   & receptaculo$sd_Enero <  tres , 2, ifelse(receptaculo$sd_Enero >= tres   & receptaculo$sd_Enero < cuatro, 3, ifelse(receptaculo$sd_Enero >= cuatro    & receptaculo$sd_Enero < quinto , 4 , ifelse(receptaculo$sd_Enero >= quinto  & receptaculo$sd_Enero <  sexto,5,"")))))
receptaculo$quintiles_rango_n_valores <-
ifelse(receptaculo$quintiles_rango_n == 1 , df[1,2], ifelse(receptaculo$quintiles_rango_n == 2 , df[2,2], ifelse(receptaculo$quintiles_rango_n == 3, df[3,2], ifelse(receptaculo$quintiles_rango_n == 4, df[4,2] , ifelse(receptaculo$quintiles_rango_n == 5,df[5,2],"")))))

11.0.1 Agosto

r_enero <- receptaculo
r_enero$Comuna <- as.double(r_enero$Comuna) 
r_enero$quintiles_rango_n <- as.double(r_enero$quintiles_rango_n) 
names(r_enero)[2] <- 'sd_Agosto'
names(r_enero)[3] <- 'mean_Agosto'
names(r_enero)[4] <- 'quintiles_rango_n_Agosto'
names(r_enero)[5] <- 'quintiles_rango_n_valores_Agosto'
r_agosto <- r_enero

12 Septiembre

minimas <- tibble()
df_min <- df_1[,c(1,seq(10, 253, 12))]
minimas <- rbind(minimas,df_min)
minimas <- minimas[order(minimas$Comuna),]
minimas$sd =  round(apply(minimas[,-1], 1, sd),3)
minimas$mean = rowMeans(minimas[,c(-1)])
names(minimas)[23] <- "sd_Enero"
names(minimas)[24] <- "mean_Enero"
minimas <- minimas[,-c(2:22,25)]
max <- max(minimas$sd_Enero)
min <- min(minimas$sd_Enero)
incremento <- (max-min)/5
dos <- min + incremento
tres <- dos  + incremento
cuatro <- tres  + incremento
quinto <- cuatro  + incremento
sexto <- quinto  + incremento + 0.001
rango <- c(1, 2, 3, 4, 5)
uno <- round(min, digits = 3)
dos <- round(dos, digits = 3)
primer_rango <- paste(uno, "-", dos)
dos <- round(dos, digits = 3)
tres <- round(tres, digits = 3)
segundo_rango <- paste(dos, "-", tres)
tres <- round(tres, digits = 3)
cuatro <- round(cuatro, digits = 3)
tercer_rango <- paste(tres, "-", cuatro)
cuatro <- round(cuatro, digits = 3)
quinto <- round(quinto, digits = 3)
cuarto_rango <- paste(cuatro, "-", quinto)
quinto <- round(quinto, digits = 3)
sexto <- round(sexto, digits = 3)
quinto_rango <- paste(quinto, "-", sexto)
rango <- c(1,2,3,4,5)
intervalos <- c(primer_rango, segundo_rango , tercer_rango , cuarto_rango , quinto_rango )
df <- data.frame(rango, intervalos)
receptaculo <- minimas
receptaculo$quintiles_rango_n <-
ifelse(receptaculo$sd_Enero >= uno & receptaculo$sd_Enero < dos , 1, ifelse(receptaculo$sd_Enero >= dos   & receptaculo$sd_Enero <  tres , 2, ifelse(receptaculo$sd_Enero >= tres   & receptaculo$sd_Enero < cuatro, 3, ifelse(receptaculo$sd_Enero >= cuatro    & receptaculo$sd_Enero < quinto , 4 , ifelse(receptaculo$sd_Enero >= quinto  & receptaculo$sd_Enero <  sexto,5,"")))))
receptaculo$quintiles_rango_n_valores <-
ifelse(receptaculo$quintiles_rango_n == 1 , df[1,2], ifelse(receptaculo$quintiles_rango_n == 2 , df[2,2], ifelse(receptaculo$quintiles_rango_n == 3, df[3,2], ifelse(receptaculo$quintiles_rango_n == 4, df[4,2] , ifelse(receptaculo$quintiles_rango_n == 5,df[5,2],"")))))

12.0.1 Septiembre

r_enero <- receptaculo
r_enero$Comuna <- as.double(r_enero$Comuna) 
r_enero$quintiles_rango_n <- as.double(r_enero$quintiles_rango_n) 
names(r_enero)[2] <- 'sd_Septiembre'
names(r_enero)[3] <- 'mean_Septiembre'
names(r_enero)[4] <- 'quintiles_rango_n_Septiembre'
names(r_enero)[5] <- 'quintiles_rango_n_valores_Septiembre'
r_septiembre <- r_enero

13 Octubre

minimas <- tibble()
df_min <- df_1[,c(1,seq(11, 253, 12))]
minimas <- rbind(minimas,df_min)
minimas <- minimas[order(minimas$Comuna),]
minimas$sd =  round(apply(minimas[,-1], 1, sd),3)
minimas$mean = rowMeans(minimas[,c(-1)])
names(minimas)[23] <- "sd_Enero"
names(minimas)[24] <- "mean_Enero"
minimas <- minimas[,-c(2:22,25)]
max <- max(minimas$sd_Enero)
min <- min(minimas$sd_Enero)
incremento <- (max-min)/5
dos <- min + incremento
tres <- dos  + incremento
cuatro <- tres  + incremento
quinto <- cuatro  + incremento
sexto <- quinto  + incremento + 0.001
rango <- c(1, 2, 3, 4, 5)
uno <- round(min, digits = 3)
dos <- round(dos, digits = 3)
primer_rango <- paste(uno, "-", dos)
dos <- round(dos, digits = 3)
tres <- round(tres, digits = 3)
segundo_rango <- paste(dos, "-", tres)
tres <- round(tres, digits = 3)
cuatro <- round(cuatro, digits = 3)
tercer_rango <- paste(tres, "-", cuatro)
cuatro <- round(cuatro, digits = 3)
quinto <- round(quinto, digits = 3)
cuarto_rango <- paste(cuatro, "-", quinto)
quinto <- round(quinto, digits = 3)
sexto <- round(sexto, digits = 3)
quinto_rango <- paste(quinto, "-", sexto)
rango <- c(1,2,3,4,5)
intervalos <- c(primer_rango, segundo_rango , tercer_rango , cuarto_rango , quinto_rango )
df <- data.frame(rango, intervalos)
receptaculo <- minimas
receptaculo$quintiles_rango_n <-
ifelse(receptaculo$sd_Enero >= uno & receptaculo$sd_Enero < dos , 1, ifelse(receptaculo$sd_Enero >= dos   & receptaculo$sd_Enero <  tres , 2, ifelse(receptaculo$sd_Enero >= tres   & receptaculo$sd_Enero < cuatro, 3, ifelse(receptaculo$sd_Enero >= cuatro    & receptaculo$sd_Enero < quinto , 4 , ifelse(receptaculo$sd_Enero >= quinto  & receptaculo$sd_Enero <  sexto,5,"")))))
receptaculo$quintiles_rango_n_valores <-
ifelse(receptaculo$quintiles_rango_n == 1 , df[1,2], ifelse(receptaculo$quintiles_rango_n == 2 , df[2,2], ifelse(receptaculo$quintiles_rango_n == 3, df[3,2], ifelse(receptaculo$quintiles_rango_n == 4, df[4,2] , ifelse(receptaculo$quintiles_rango_n == 5,df[5,2],"")))))

13.0.1 Octubre

r_enero <- receptaculo
r_enero$Comuna <- as.double(r_enero$Comuna) 
r_enero$quintiles_rango_n <- as.double(r_enero$quintiles_rango_n) 
names(r_enero)[2] <- 'sd_Octubre'
names(r_enero)[3] <- 'mean_Octubre'
names(r_enero)[4] <- 'quintiles_rango_n_Octubre'
names(r_enero)[5] <- 'quintiles_rango_n_valores_Octubre'
r_octubre <- r_enero

14 Noviembre

minimas <- tibble()
df_min <- df_1[,c(1,seq(12, 253, 12))]
minimas <- rbind(minimas,df_min)
minimas <- minimas[order(minimas$Comuna),]
minimas$sd =  round(apply(minimas[,-1], 1, sd),3)
minimas$mean = rowMeans(minimas[,c(-1)])
names(minimas)[23] <- "sd_Enero"
names(minimas)[24] <- "mean_Enero"
minimas <- minimas[,-c(2:22,25)]
max <- max(minimas$sd_Enero)
min <- min(minimas$sd_Enero)
incremento <- (max-min)/5
dos <- min + incremento
tres <- dos  + incremento
cuatro <- tres  + incremento
quinto <- cuatro  + incremento
sexto <- quinto  + incremento + 0.001
rango <- c(1, 2, 3, 4, 5)
uno <- round(min, digits = 3)
dos <- round(dos, digits = 3)
primer_rango <- paste(uno, "-", dos)
dos <- round(dos, digits = 3)
tres <- round(tres, digits = 3)
segundo_rango <- paste(dos, "-", tres)
tres <- round(tres, digits = 3)
cuatro <- round(cuatro, digits = 3)
tercer_rango <- paste(tres, "-", cuatro)
cuatro <- round(cuatro, digits = 3)
quinto <- round(quinto, digits = 3)
cuarto_rango <- paste(cuatro, "-", quinto)
quinto <- round(quinto, digits = 3)
sexto <- round(sexto, digits = 3)
quinto_rango <- paste(quinto, "-", sexto)
rango <- c(1,2,3,4,5)
intervalos <- c(primer_rango, segundo_rango , tercer_rango , cuarto_rango , quinto_rango )
df <- data.frame(rango, intervalos)
receptaculo <- minimas
receptaculo$quintiles_rango_n <-
ifelse(receptaculo$sd_Enero >= uno & receptaculo$sd_Enero < dos , 1, ifelse(receptaculo$sd_Enero >= dos   & receptaculo$sd_Enero <  tres , 2, ifelse(receptaculo$sd_Enero >= tres   & receptaculo$sd_Enero < cuatro, 3, ifelse(receptaculo$sd_Enero >= cuatro    & receptaculo$sd_Enero < quinto , 4 , ifelse(receptaculo$sd_Enero >= quinto  & receptaculo$sd_Enero <  sexto,5,"")))))
receptaculo$quintiles_rango_n_valores <-
ifelse(receptaculo$quintiles_rango_n == 1 , df[1,2], ifelse(receptaculo$quintiles_rango_n == 2 , df[2,2], ifelse(receptaculo$quintiles_rango_n == 3, df[3,2], ifelse(receptaculo$quintiles_rango_n == 4, df[4,2] , ifelse(receptaculo$quintiles_rango_n == 5,df[5,2],"")))))

14.0.1 Noviembre

r_enero <- receptaculo
r_enero$Comuna <- as.double(r_enero$Comuna) 
r_enero$quintiles_rango_n <- as.double(r_enero$quintiles_rango_n) 
names(r_enero)[2] <- 'sd_Noviembre'
names(r_enero)[3] <- 'mean_Noviembre'
names(r_enero)[4] <- 'quintiles_rango_n_Noviembre'
names(r_enero)[5] <- 'quintiles_rango_n_valores_Noviembre'
r_noviembre <- r_enero

15 Diciembre

minimas <- tibble()
df_min <- df_1[,c(1,seq(13, 253, 12))]
minimas <- rbind(minimas,df_min)
minimas <- minimas[order(minimas$Comuna),]
minimas$sd =  round(apply(minimas[,-1], 1, sd),3)
minimas$mean = rowMeans(minimas[,c(-1)])
names(minimas)[23] <- "sd_Enero"
names(minimas)[24] <- "mean_Enero"
minimas <- minimas[,-c(2:22,25)]
max <- max(minimas$sd_Enero)
min <- min(minimas$sd_Enero)
incremento <- (max-min)/5
dos <- min + incremento
tres <- dos  + incremento
cuatro <- tres  + incremento
quinto <- cuatro  + incremento
sexto <- quinto  + incremento + 0.001
rango <- c(1, 2, 3, 4, 5)
uno <- round(min, digits = 3)
dos <- round(dos, digits = 3)
primer_rango <- paste(uno, "-", dos)
dos <- round(dos, digits = 3)
tres <- round(tres, digits = 3)
segundo_rango <- paste(dos, "-", tres)
tres <- round(tres, digits = 3)
cuatro <- round(cuatro, digits = 3)
tercer_rango <- paste(tres, "-", cuatro)
cuatro <- round(cuatro, digits = 3)
quinto <- round(quinto, digits = 3)
cuarto_rango <- paste(cuatro, "-", quinto)
quinto <- round(quinto, digits = 3)
sexto <- round(sexto, digits = 3)
quinto_rango <- paste(quinto, "-", sexto)
rango <- c(1,2,3,4,5)
intervalos <- c(primer_rango, segundo_rango , tercer_rango , cuarto_rango , quinto_rango )
df <- data.frame(rango, intervalos)
receptaculo <- minimas
receptaculo$quintiles_rango_n <-
ifelse(receptaculo$sd_Enero >= uno & receptaculo$sd_Enero < dos , 1, ifelse(receptaculo$sd_Enero >= dos   & receptaculo$sd_Enero <  tres , 2, ifelse(receptaculo$sd_Enero >= tres   & receptaculo$sd_Enero < cuatro, 3, ifelse(receptaculo$sd_Enero >= cuatro    & receptaculo$sd_Enero < quinto , 4 , ifelse(receptaculo$sd_Enero >= quinto  & receptaculo$sd_Enero <  sexto,5,"")))))
receptaculo$quintiles_rango_n_valores <-
ifelse(receptaculo$quintiles_rango_n == 1 , df[1,2], ifelse(receptaculo$quintiles_rango_n == 2 , df[2,2], ifelse(receptaculo$quintiles_rango_n == 3, df[3,2], ifelse(receptaculo$quintiles_rango_n == 4, df[4,2] , ifelse(receptaculo$quintiles_rango_n == 5,df[5,2],"")))))

15.0.1 Diciembre

r_enero <- receptaculo
r_enero$Comuna <- as.double(r_enero$Comuna) 
r_enero$quintiles_rango_n <- as.double(r_enero$quintiles_rango_n) 
names(r_enero)[2] <- 'sd_Diciembre'
names(r_enero)[3] <- 'mean_Diciembre'
names(r_enero)[4] <- 'quintiles_rango_n_Diciembre'
names(r_enero)[5] <- 'quintiles_rango_n_valores_Diciembre'
r_diciembre <- r_enero

16 Hay que apilar horizontalmente

datos_tmmn = cbind(r_Enero, r_febrero, r_marzo, r_abril, r_mayo, r_junio, r_julio, r_agosto,  r_septiembre, r_octubre, r_noviembre, r_diciembre)

# datatable(datos_tmmn, extensions = 'Buttons', escape = FALSE, rownames = TRUE,
#           options = list(dom = 'Bfrtip',
#           buttons = list('colvis', list(extend = 'collection',
#           buttons = list(
#           list(extend='copy'),
#           list(extend='excel',
#             filename = 'tabla'),
#           list(extend='pdf',
#             filename= 'tabla')),
#           text = 'Download')), scrollX = TRUE))

17 sumamos los rangos:

datos_tmmn$suma_de_rangos <- datos_tmmn$quintiles_rango_n_Enero + datos_tmmn$quintiles_rango_n_Febrero + datos_tmmn$quintiles_rango_n_Marzo + datos_tmmn$quintiles_rango_n_Abril + datos_tmmn$quintiles_rango_n_Mayo +
datos_tmmn$quintiles_rango_n_Junio + datos_tmmn$quintiles_rango_n_Julio + datos_tmmn$quintiles_rango_n_Agosto +
datos_tmmn$quintiles_rango_n_Septiembre + datos_tmmn$quintiles_rango_n_Octubre + datos_tmmn$quintiles_rango_n_Noviembre +
datos_tmmn$quintiles_rango_n_Diciembre 
# datatable(datos_tmmn, extensions = 'Buttons', escape = FALSE, rownames = TRUE,
#           options = list(dom = 'Bfrtip',
#           buttons = list('colvis', list(extend = 'collection',
#           buttons = list(
#           list(extend='copy'),
#           list(extend='excel',
#             filename = 'tabla'),
#           list(extend='pdf',
#             filename= 'tabla')),
#           text = 'Download')), scrollX = TRUE))
max <- max(datos_tmmn$suma_de_rangos)
min <- min(datos_tmmn$suma_de_rangos)
incremento <- (max-min)/5
dos <- min + incremento
tres <- dos  + incremento
cuatro <- tres  + incremento
quinto <- cuatro  + incremento
sexto <- quinto  + incremento + 0.001
rango <- c(1, 2, 3, 4, 5)
uno <- round(min, digits = 3)
dos <- round(dos, digits = 3)
primer_rango <- paste(uno, "-", dos)
dos <- round(dos, digits = 3)
tres <- round(tres, digits = 3)
segundo_rango <- paste(dos, "-", tres)
tres <- round(tres, digits = 3)
cuatro <- round(cuatro, digits = 3)
tercer_rango <- paste(tres, "-", cuatro)
cuatro <- round(cuatro, digits = 3)
quinto <- round(quinto, digits = 3)
cuarto_rango <- paste(cuatro, "-", quinto)
quinto <- round(quinto, digits = 3)
sexto <- round(sexto, digits = 3)
quinto_rango <- paste(quinto, "-", sexto)
rango <- c(1,2,3,4,5)
intervalos <- c(primer_rango, segundo_rango , tercer_rango , cuarto_rango , quinto_rango )
df <- data.frame(rango, intervalos)
df
##   rango    intervalos
## 1     1     12 - 19.8
## 2     2   19.8 - 27.6
## 3     3   27.6 - 35.4
## 4     4   35.4 - 43.2
## 5     5 43.2 - 51.001
receptaculo <- datos_tmmn
receptaculo$quintiles_rango_n_suma <-
ifelse(receptaculo$suma_de_rangos >= uno & receptaculo$suma_de_rangos < dos , 1, ifelse(receptaculo$suma_de_rangos >= dos   & receptaculo$suma_de_rangos <  tres , 2, ifelse(receptaculo$suma_de_rangos >= tres   & receptaculo$suma_de_rangos < cuatro, 3, ifelse(receptaculo$suma_de_rangos >= cuatro    & receptaculo$suma_de_rangos < quinto , 4 , ifelse(receptaculo$suma_de_rangos >= quinto  & receptaculo$suma_de_rangos <  sexto,5,"")))))
receptaculo$quintiles_rango_n_valores_suma <-
ifelse(receptaculo$quintiles_rango_n_suma == 1 , df[1,2], ifelse(receptaculo$quintiles_rango_n_suma == 2 , df[2,2], ifelse(receptaculo$quintiles_rango_n_suma == 3, df[3,2], ifelse(receptaculo$quintiles_rango_n_suma == 4, df[4,2] , ifelse(receptaculo$quintiles_rango_n_suma == 5,df[5,2],"")))))
names(receptaculo)[6] <- "comuna"
names(receptaculo)[11] <- "comuna"
names(receptaculo)[16] <- "comuna"
names(receptaculo)[21] <- "comuna"
names(receptaculo)[26] <- "comuna"
names(receptaculo)[31] <- "comuna"
names(receptaculo)[36] <- "comuna"
names(receptaculo)[41] <- "comuna"
names(receptaculo)[46] <- "comuna"
names(receptaculo)[51] <- "comuna"
names(receptaculo)[56] <- "comuna"
aaa <- merge(receptaculo, y = codigos_comunales_actuales)
## Warning in merge.data.frame(receptaculo, y = codigos_comunales_actuales): column
## names 'comuna', 'comuna', 'comuna', 'comuna', 'comuna', 'comuna', 'comuna',
## 'comuna', 'comuna', 'comuna' are duplicated in the result

aaa es nuestra tabla final provisional. Guademosla como rds

saveRDS(aaa, "tmmx.rds")
bbb <- readRDS("tmmx.rds")

Sumamos las sd de cada mes:

bbb$promedios_sd_parte_tmmx_del_ivacc <- (bbb[, 2] +  bbb[, 7] + bbb[, 12] +  bbb[, 17] + bbb[, 22] +  bbb[, 27] + bbb[, 32] +  bbb[, 37] +
    bbb[, 42] +  bbb[, 47] + bbb[, 52] +  bbb[, 57])/12

Anadimos columna de tmmx_promedios_sd estandarizada:

maxi <- max(bbb$promedios_sd_parte_tmmx_del_ivacc)
mini <- min(bbb$promedios_sd_parte_tmmx_del_ivacc)
bbb$promedios_sd_estandarizado <- (bbb$promedios_sd_parte_tmmx_del_ivacc - mini)/(maxi - mini)
bbb$tmmn_quintiles_rango_n_suma_leyenda <- ifelse(bbb$quintiles_rango_n_suma == 1 , "muy baja vulnerabilidad - tmmx", ifelse(bbb$quintiles_rango_n_suma == 2 ,"baja vulnerabilidad - tmmx", ifelse(bbb$quintiles_rango_n_suma == 3, "moderada vulnerabilidad - tmmx", ifelse(bbb$quintiles_rango_n_suma == 4, "alta vulnerabilidad - tmmx" , ifelse(bbb$quintiles_rango_n_suma == 5,"muy alta vulnerabilidad - tmmx","")))))

en la columna 61 hay que agregar arteriscos

bbb$x <-  "*"

Ordenamos las columnas:

en 61 debe ir la 73

bbb1 <- bbb[, c(  1,  2 , 3 , 4 , 5 , 6 , 7 , 8 , 9 ,10, 11, 12, 13, 14, 15 ,16 ,17, 18 ,19 ,20 ,21 ,22 ,23 ,24 ,25 ,26 ,27 ,28 ,29 ,30 ,31 ,32 ,33 ,34 ,35 ,36 ,37 ,38 ,39 ,40 ,41 ,42 ,43, 44 ,45 ,46 ,47,48 ,49 ,50, 51 ,52 ,53 ,54, 55 ,56 ,57 ,58 ,59 ,60 ,73 ,62 ,63, 64 ,65 ,66 ,67, 68, 69 ,70 ,71, 72,61)]

en 62 debe ir la 73

bbb2 <- bbb1[, c(  1,  2 , 3 , 4 , 5 , 6 , 7 , 8 , 9 ,10, 11, 12, 13, 14, 15 ,16 ,17, 18 ,19 ,20 ,21 ,22 ,23 ,24 ,25 ,26 ,27 ,28 ,29 ,30 ,31 ,32 ,33 ,34 ,35 ,36 ,37 ,38 ,39 ,40 ,41 ,42 ,43, 44 ,45 ,46 ,47,48 ,49 ,50, 51 ,52 ,53 ,54, 55 ,56 ,57 ,58 ,59 ,60 ,61 ,73 ,63, 64 ,65 ,66 ,67, 68, 69 ,70 ,71, 72,62)]

en 63 debe ir la 73

bbb3 <- bbb2[, c(  1,  2 , 3 , 4 , 5 , 6 , 7 , 8 , 9 ,10, 11, 12, 13, 14, 15 ,16 ,17, 18 ,19 ,20 ,21 ,22 ,23 ,24 ,25 ,26 ,27 ,28 ,29 ,30 ,31 ,32 ,33 ,34 ,35 ,36 ,37 ,38 ,39 ,40 ,41 ,42 ,43, 44 ,45 ,46 ,47,48 ,49 ,50, 51 ,52 ,53 ,54, 55 ,56 ,57 ,58 ,59 ,60 ,61 ,62 ,73, 64 ,65 ,66 ,67, 68, 69 ,70 ,71, 72,63)]

en 64 y 65 deben ir 72 y 73

bbb4 <- bbb3[, c(  1,  2 , 3 , 4 , 5 , 6 , 7 , 8 , 9 ,10, 11, 12, 13, 14, 15 ,16 ,17, 18 ,19 ,20 ,21 ,22 ,23 ,24 ,25 ,26 ,27 ,28 ,29 ,30 ,31 ,32 ,33 ,34 ,35 ,36 ,37 ,38 ,39 ,40 ,41 ,42 ,43, 44 ,45 ,46 ,47,48 ,49 ,50, 51 ,52 ,53 ,54, 55 ,56 ,57 ,58 ,59 ,60 ,61 ,62 ,63, 72 ,73 ,66 ,67, 68, 69 ,70 ,71, 64,65)]

en 66 y 67 deben ir 70 y 71

bbb5 <- bbb4[, c(  1,  2 , 3 , 4 , 5 , 6 , 7 , 8 , 9 ,10, 11, 12, 13, 14, 15 ,16 ,17, 18 ,19 ,20 ,21 ,22 ,23 ,24 ,25 ,26 ,27 ,28 ,29 ,30 ,31 ,32 ,33 ,34 ,35 ,36 ,37 ,38 ,39 ,40 ,41 ,42 ,43, 44 ,45 ,46 ,47,48 ,49 ,50, 51 ,52 ,53 ,54, 55 ,56 ,57 ,58 ,59 ,60 ,61 ,62 ,63, 64 ,65 ,70 ,71, 68, 69 ,66 ,67,72,73)]

Asignemosle correctos nombres a las columnas desde el *:

names(bbb5)[62] <- "suma de los quintiles de sd por mes"
names(bbb5)[63] <- "quintiles_sobre_la_suma_de_los_quintiles_de_sd_por_mes"
names(bbb5)[64] <- "categoría"
names(bbb5)[65] <- "rango"
names(bbb5)[66] <- "promedios_sd.parte_del_ivacc"
names(bbb5)[67] <- "promedios_sd_estandarizados"

debemos asignarle el postfijo tmmx

colnames(bbb5)[1:73] <- str_c( colnames(bbb5)[1:73], "_tmmx" )
# datatable(bbb5, extensions = 'Buttons', escape = FALSE, rownames = TRUE,
#           options = list(dom = 'Bfrtip',
#           buttons = list('colvis', list(extend = 'collection',
#           buttons = list(
#           list(extend='copy'),
#           list(extend='excel',
#             filename = 'tmmx.xlsx'),
#           list(extend='pdf',
#             filename= 'tabla')),
#           text = 'Download')), scrollX = TRUE))

Debemos establecer 5 rangos numericos para la columna promedios_sd.parte_del_ivacc_tmmx

max <- max(bbb5$promedios_sd.parte_del_ivacc_tmmx)
min <- min(bbb5$promedios_sd.parte_del_ivacc_tmmx)
incremento <- (max-min)/5
dos <- min + incremento
tres <- dos  + incremento
cuatro <- tres  + incremento
quinto <- cuatro  + incremento
sexto <- quinto  + incremento + 0.001
rango <- c(1, 2, 3, 4, 5)
uno <- round(min, digits = 3)
dos <- round(dos, digits = 3)
primer_rango <- paste(uno, "-", dos)
dos <- round(dos, digits = 3)
tres <- round(tres, digits = 3)
segundo_rango <- paste(dos, "-", tres)
tres <- round(tres, digits = 3)
cuatro <- round(cuatro, digits = 3)
tercer_rango <- paste(tres, "-", cuatro)
cuatro <- round(cuatro, digits = 3)
quinto <- round(quinto, digits = 3)
cuarto_rango <- paste(cuatro, "-", quinto)
quinto <- round(quinto, digits = 3)
sexto <- round(sexto, digits = 3)
quinto_rango <- paste(quinto, "-", sexto)
rango <- c(1,2,3,4,5)
intervalos <- c(primer_rango, segundo_rango , tercer_rango , cuarto_rango , quinto_rango )
df <- data.frame(rango, intervalos)

debemos asociar promedios_sd.parte_del_ivacc_tmmx con valor del rango:

receptaculo <- bbb5
receptaculo$valor_del_rango_de_pro_sd <-
ifelse(receptaculo$promedios_sd.parte_del_ivacc_tmmx >= 0.516 & receptaculo$promedios_sd.parte_del_ivacc_tmmx < 0.672 , 1, ifelse(receptaculo$promedios_sd.parte_del_ivacc_tmmx >= 0.672   & receptaculo$promedios_sd.parte_del_ivacc_tmmx <  0.828 , 2, ifelse(receptaculo$promedios_sd.parte_del_ivacc_tmmx >= 0.828   & receptaculo$promedios_sd.parte_del_ivacc_tmmx < 0.984, 3, ifelse(receptaculo$promedios_sd.parte_del_ivacc_tmmx >= 0.984    & receptaculo$promedios_sd.parte_del_ivacc_tmmx < 1.14 , 4 , ifelse(receptaculo$promedios_sd.parte_del_ivacc_tmmx >= 1.14  & receptaculo$promedios_sd.parte_del_ivacc_tmmx <  1.297,5,"")))))
receptaculo$rango <-
ifelse(receptaculo$valor_del_rango_de_pro_sd == 1 , df[1,2], ifelse(receptaculo$valor_del_rango_de_pro_sd == 2 , df[2,2], ifelse(receptaculo$valor_del_rango_de_pro_sd == 3, df[3,2], ifelse(receptaculo$valor_del_rango_de_pro_sd == 4, df[4,2] , ifelse(receptaculo$valor_del_rango_de_pro_sd == 5,df[5,2],"")))))
receptaculo$nombre_del_rango <- ifelse(receptaculo$quintiles_sobre_la_suma_de_los_quintiles_de_sd_por_mes_tmmx == 1 , "muy baja vulnerabilidad - tmmx", ifelse(receptaculo$quintiles_sobre_la_suma_de_los_quintiles_de_sd_por_mes_tmmx == 2 ,"baja vulnerabilidad - tmmx", ifelse(receptaculo$quintiles_sobre_la_suma_de_los_quintiles_de_sd_por_mes_tmmx == 3, "moderada vulnerabilidad - tmmx", ifelse(receptaculo$quintiles_sobre_la_suma_de_los_quintiles_de_sd_por_mes_tmmx == 4, "alta vulnerabilidad - tmmx" , ifelse(receptaculo$quintiles_sobre_la_suma_de_los_quintiles_de_sd_por_mes_tmmx == 5,"muy alta vulnerabilidad - tmmx","")))))

Ordenamos las columnas:

receptaculo_final <- receptaculo[, c(  1,  2 , 3 , 4 , 5 , 6 , 7 , 8 , 9 ,10, 11, 12, 13, 14, 15 ,16 ,17, 18 ,19 ,20 ,21 ,22 ,23 ,24 ,25 ,26 ,27 ,28 ,29 ,30 ,31 ,32 ,33 ,34 ,35 ,36 ,37 ,38 ,39 ,40 ,41 ,42 ,43, 44 ,45 ,46 ,47,48 ,49 ,50, 51 ,52 ,53 ,54, 55 ,56 ,57 ,58 ,59 ,60 ,61 ,62 ,63, 64 ,65 ,66 ,67, 75, 74 ,76 ,71,72,73,69,68,70)]
datatable(receptaculo_final, extensions = 'Buttons', escape = FALSE, rownames = TRUE,
          options = list(dom = 'Bfrtip',
          buttons = list('colvis', list(extend = 'collection',
          buttons = list(
          list(extend='copy'),
          list(extend='excel',
            filename = 'tmmx_z'),
          list(extend='pdf',
            filename= 'tabla')),
          text = 'Download')), scrollX = TRUE))

La fórmula aplicada al conjunto de variables IVACC se hace sobre los promedios_sd.parte_del_ivacc_tmmx para todas las variables!!!


Se deben pintar:

promedios_sd.parte_del_ivacc_tmmx y promedios_sd_estandarizados_tmmx