##### UNIVERSIDAD CENTRAL DEL ECUADOR #####
#### AUTOR: MARTIN SARMIENTO ####
### CARRERA: INGENIERÍA EN PETRÓLEOS #####
#### VARIABLE ESTADO OPERATIVO ####
## DATASET ##
setwd("~/R/OPERATIONAL_STATUS")
# Cargar dataset
Datos <- read.csv("DataSet_.csv", sep = ";", fileEncoding = "latin1")
# Estructura de los datos
str(Datos)## 'data.frame': 7142 obs. of 26 variables:
## $ fid : int 1 2 3 4 5 6 7 8 9 10 ...
## $ objectid : int 127 128 129 130 131 132 133 134 135 136 ...
## $ code : chr "Arg-00001" "Arg-00002" "Arg-00003" "Arg-00004" ...
## $ country : chr "Argentina" "Argentina" "Argentina" "Argentina" ...
## $ plant_name : chr "Aconcagua solar farm" "Aconcagua solar farm" "Altiplano 200 Solar Power Plant" "Altiplano 200 Solar Power Plant" ...
## $ operational_status : chr "announced" "announced" "operating" "operating" ...
## $ longitude : num -68.9 -68.9 -66.9 -66.9 -68.9 ...
## $ latitude : num -33 -33 -24.1 -24.1 -33.3 ...
## $ elevation : int 929 929 4000 4000 937 865 858 858 858 858 ...
## $ area : num 0 0 4397290 5774 0 ...
## $ slope : num 0.574 0.574 1.603 6.243 0.903 ...
## $ slope_type : chr "Plano o casi plano" "Plano o casi plano" "Plano o casi plano" "Moderado" ...
## $ curvature : num 0.000795 0.000795 -0.002781 -0.043699 0.002781 ...
## $ curvature_type : chr "Superficies planas o intermedias" "Superficies planas o intermedias" "Superficies planas o intermedias" "Superficies cóncavas / Valles" ...
## $ aspect : num 55.1 55.1 188.7 270.9 108.4 ...
## $ aspect_type : chr "Northeast" "Northeast" "South" "West" ...
## $ ghi : num 6.11 6.11 8.01 7.88 6.12 ...
## $ solar_aptitude : num 0.746 0.746 0.8 0.727 0.595 ...
## $ solar_aptittude_class: chr "Alta" "Alta" "Alta" "Alta" ...
## $ humidity : num 0 0 53.7 53.7 0 ...
## $ wind_speed : num 3.78 3.78 7.02 8.33 3.87 ...
## $ wind_direction : num 0 0 55.1 55.1 0 ...
## $ ambient_temperature : num 12.6 12.6 6.8 6.8 13.1 ...
## $ optimal_tilt : int 31 31 26 26 31 33 30 30 30 30 ...
## $ peak_power_per_hour : num 4.98 4.98 6.39 6.39 4.97 ...
## $ total_power : num 25 66.2 101 107 180 ...
##
## Adjuntando el paquete: 'dplyr'
## The following objects are masked from 'package:stats':
##
## filter, lag
## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, union
#Extraer variable
Status <- Datos$operational_status
# EDAvariable nominal
TDF_status <- table(Status)
tabla_status <- as.data.frame(TDF_status)
hi <- tabla_status$Freq/sum(tabla_status$Freq)
hi_porc <- hi*100
sum(hi_porc)## [1] 100
### Agrupación ###
tabla_resumen <- tabla_STATUS %>%
mutate(grupo = case_when(
grepl("cancelled - inferred 4 y", Status, ignore.case = TRUE) ~ "Cancelado (Inferred 4y)",
grepl("shelved - inferred 2 y", Status, ignore.case = TRUE) ~ "Archivado (Inferred 2y)",
grepl("operating", Status, ignore.case = TRUE) ~ "Operativo",
grepl("announced", Status, ignore.case = TRUE) ~ "Anunciado",
grepl("pre-construction", Status, ignore.case = TRUE) ~ "Pre-Construcción",
grepl("construction", Status, ignore.case = TRUE) ~ "Construcción",
grepl("mothballed", Status, ignore.case = TRUE) ~ "Inactivo (Mothballed)",
grepl("shelved", Status, ignore.case = TRUE) ~ "Archivado (Shelved)",
grepl("cancelled", Status, ignore.case = TRUE) ~ "Cancelado",
TRUE ~ "Otros"))
# Definimos el orden
orden_especifico <- c(
"Operativo",
"Anunciado",
"Construcción",
"Pre-Construcción",
"Inactivo (Mothballed)",
"Archivado (Shelved)",
"Archivado (Inferred 2y)",
"Cancelado",
"Cancelado (Inferred 4y)",
"Otros")
tabla_resumen <- tabla_resumen %>%
mutate(grupo = factor(grupo, levels = orden_especifico)) %>%
group_by(grupo) %>%
summarise(
Frecuencia = sum(Freq),
Porcentaje = sum(hi_porc)) %>%
arrange(grupo)
# Renombramos columnas
colnames(tabla_resumen) <- c("Estado","ni","hi (%)")
# Tabla Intermedia GT
tabla_resumen_gt <- tabla_resumen %>%
gt() %>%
tab_header(
title = md("**TABLA DE AGRUPACIÓN**"),
subtitle = "Frecuencias por Estado Operativo (Orden Específico)") %>%
tab_source_note(source_note = "Autor: Martin Sarmiento") %>%
cols_label(
Estado = "Estado Operativo",
ni = "Frecuencia (ni)",
`hi (%)` = "Porcentaje (hi%)") %>%
fmt_number(columns = c(`hi (%)`), decimals = 2) %>%
tab_options(
heading.title.font.size = px(16),
heading.subtitle.font.size = px(14),
column_labels.background.color = "#f0f0f0")
# Mostramos la tabla
tabla_resumen_gt| TABLA DE AGRUPACIÓN | ||
| Frecuencias por Estado Operativo (Orden Específico) | ||
| Estado Operativo | Frecuencia (ni) | Porcentaje (hi%) |
|---|---|---|
| Operativo | 3035 | 42.50 |
| Anunciado | 93 | 1.30 |
| Construcción | 161 | 2.25 |
| Pre-Construcción | 3412 | 47.77 |
| Inactivo (Mothballed) | 2 | 0.03 |
| Archivado (Shelved) | 37 | 0.52 |
| Archivado (Inferred 2y) | 105 | 1.47 |
| Cancelado | 267 | 3.74 |
| Cancelado (Inferred 4y) | 29 | 0.41 |
| Otros | 1 | 0.01 |
| Autor: Martin Sarmiento | ||
#### Crear de fila de totales ####
tabla_resumen$Estado <- as.character(tabla_resumen$Estado)
totales <- c("TOTAL", sum(tabla_resumen$ni, na.rm=TRUE), sum(tabla_resumen$`hi (%)`, na.rm=TRUE))
tabla_Status_Final <- rbind(tabla_resumen, totales)
# Convertir a números para GT
tabla_Status_Final$ni <- as.numeric(tabla_Status_Final$ni)
tabla_Status_Final$`hi (%)` <- as.numeric(tabla_Status_Final$`hi (%)`)
# TABLA 1
tabla_final_gt <- tabla_Status_Final %>%
gt() %>%
tab_header(
title = md("**TABLA 1**"),
subtitle = "Distribución Detallada del Estado Operativo") %>%
tab_source_note(source_note = "Autor: Martin Sarmiento") %>%
cols_label(
Estado = "Estado Operativo",
ni = "Frecuencia (ni)",
`hi (%)` = "Porcentaje (hi%)") %>%
fmt_number(columns = c(`hi (%)`), decimals = 2) %>%
tab_options(
heading.title.font.size = px(16),
heading.subtitle.font.size = px(14),
column_labels.background.color = "#f0f0f0")
tabla_final_gt| TABLA 1 | ||
| Distribución Detallada del Estado Operativo | ||
| Estado Operativo | Frecuencia (ni) | Porcentaje (hi%) |
|---|---|---|
| Operativo | 3035 | 42.50 |
| Anunciado | 93 | 1.30 |
| Construcción | 161 | 2.25 |
| Pre-Construcción | 3412 | 47.77 |
| Inactivo (Mothballed) | 2 | 0.03 |
| Archivado (Shelved) | 37 | 0.52 |
| Archivado (Inferred 2y) | 105 | 1.47 |
| Cancelado | 267 | 3.74 |
| Cancelado (Inferred 4y) | 29 | 0.41 |
| Otros | 1 | 0.01 |
| TOTAL | 7142 | 100.00 |
| Autor: Martin Sarmiento | ||
# Diagrama de barrras local ni
par(mar = c(12, 4, 4, 2))
barplot(tabla_resumen$ni,
main="Gráfica N°1: Distribución del estado operativo",
xlab = "",
ylab = "Cantidad",
col = "#D8BFD8",
names.arg=tabla_resumen$Estado,
cex.names = 0.6,
las = 2)
mtext("Estado Operativo", side = 1, line = 10)# Diagrama de barrras global ni
par(mar = c(12, 4, 4, 2))
barplot(tabla_resumen$ni,
main="Gráfica N°2: Distribución Global",
xlab = "",
ylab = "Cantidad",
col ="#D8BFD8",
ylim = c(0,8000),
names.arg=tabla_resumen$Estado,
cex.names = 0.6,
las = 2)
mtext("Estado Operativo", side = 1, line = 10)# Diagrama de barrras local hi
par(mar = c(12, 4, 4, 2))
barplot(tabla_resumen$`hi (%)`,
main="Gráfica N°3: Distribución porcentual",
xlab = "",
ylab = "Porcentaje %",
col = "#D8BFD8",
names.arg=tabla_resumen$Estado,
cex.names = 0.6,
las = 2)
mtext("Estado Operativo", side = 1, line = 10)# Diagrama de barrras global hi
par(mar = c(12, 4, 4, 2))
barplot(tabla_resumen$`hi (%)`,
main="Gráfica N°4: Distribución porcentual global",
xlab = "",
ylab = "Porcentaje %",
col = "#D8BFD8",
ylim = c(0,100),
names.arg=tabla_resumen$Estado,
cex.names = 0.6,
las = 2)
mtext("Estado Operativo", side = 1, line = 10)# Diagrma circular
par(mar = c(5, 6, 4, 10), xpd = TRUE)
labels_formato <- paste0(round(tabla_resumen$`hi (%)`, 2), "%")
tonos_base <- c("#4A148C", "#7B1FA2", "#AB47BC", "#CE93D8", "#F3E5F5")
paleta_violeta <- colorRampPalette(tonos_base)
colores_finales <- paleta_violeta(length(tabla_resumen$Estado))
pie (tabla_resumen$`hi (%)`,
main = "Gráfica N°5: Distribución porcentual del estado operativo",
radius = 1,
labels = paste0(round(tabla_resumen$`hi (%)`,2),"%"),
col = colores_finales,
cex = 0.6,)
# Leyenda del diagrama circular
legend(x = 1.25, y = 1,
legend = tabla_resumen$Estado,
fill = colores_finales,
cex = 0.5,
title = "Estado Operativo")# Tabla de conclusiones
tabla_indicadores <- data.frame(
"Variable" = "Estado Operativo",
"Rango" = "Ciclo de vida del proyecto (Ordinal)",
"X" = "-", "Me" = "-",
"Mo" = "Anunciado / Operativo",
"V" = "-", "Sd" = "-", "Cv" = "-", "As" = "-", "K" = "-",
"Valores_Atipicos" = "-")
tabla_conclusiones_gt <- tabla_indicadores %>%
gt() %>%
tab_header(
title = md("**CONCLUSIONES**"),
subtitle = "Análisis de la variable Estado Operativo") %>%
tab_source_note(source_note = "Autor: Martin Sarmiento") %>%
cols_label(Valores_Atipicos = "Valores Atípicos") %>%
tab_options(column_labels.background.color = "#f0f0f0")
tabla_conclusiones_gt| CONCLUSIONES | ||||||||||
| Análisis de la variable Estado Operativo | ||||||||||
| Variable | Rango | X | Me | Mo | V | Sd | Cv | As | K | Valores Atípicos |
|---|---|---|---|---|---|---|---|---|---|---|
| Estado Operativo | Ciclo de vida del proyecto (Ordinal) | - | - | Anunciado / Operativo | - | - | - | - | - | - |
| Autor: Martin Sarmiento | ||||||||||