##### UNIVERSIDAD CENTRAL DEL ECUADOR #####
#### AUTOR: MARTIN SARMIENTO ####
### CARRERA: INGENIERÍA EN PETRÓLEOS #####
#### VARIABLE TIPO DE CURVATURA ####
## DATASET ##
setwd("~/R/CURVATURE_TYPE")
# Cargar dataset
Datos <- read.csv("DataSet_prov.csv", sep = ";", fileEncoding = "latin1")
# Estructura de los datos
str(Datos)## 'data.frame': 5075 obs. of 30 variables:
## $ FID_ : int 0 2 3 4 5 6 10 11 12 13 ...
## $ OBJECTID : int 127 129 130 131 132 133 137 138 139 140 ...
## $ code : chr "00127-ARG-P" "00129-ARG-G" "00130-ARG-P" "00131-ARG-P" ...
## $ plant_name : chr "Aconcagua solar farm" "Altiplano 200 Solar Power Plant" "Altiplano 200 Solar Power Plant" "Anchoris solar farm" ...
## $ country : chr "Argentina" "Argentina" "Argentina" "Argentina" ...
## $ operational_status : chr "announced" "operating" "operating" "construction" ...
## $ longitude : chr "-68,8713" "-66,895798" "-66,926102" "-68,915001" ...
## $ latitude : chr "-32,998501" "-24,1392" "-24,073999" "-33,330101" ...
## $ elevation : int 929 4000 4000 937 865 858 570 1612 665 3989 ...
## $ area : chr "250,337006" "4397290" "5774,399902" "645,163025" ...
## $ size : chr "Pequeña" "Grande" "Pequeña" "Pequeña" ...
## $ slope : chr "0,574179" "1,60257" "6,24265" "0,902748" ...
## $ slope_type : chr "Plano o casi plano" "Plano o casi plano" "Moderado" "Plano o casi plano" ...
## $ curvature : chr "0,000795" "-0,002781" "-0,043699" "0,002781" ...
## $ curvature_type : chr "Superficies planas o intermedias" "Superficies planas o intermedias" "Superficies cóncavas / Valles" "Superficies planas o intermedias" ...
## $ aspect : chr "55,124672" "188,707367" "270,913513" "108,434952" ...
## $ aspect_type : chr "Northeast" "South" "West" "East" ...
## $ dist_to_road : chr "127,2827045" "56014,95403" "52696,78572" "335,9280031" ...
## $ ambient_temperature : chr "12,6" "6,8" "6,8" "13,1" ...
## $ ghi : chr "6,11" "8,012" "7,878" "6,119" ...
## $ humidity : chr "53,74" "53,74" "53,74" "53,74" ...
## $ wind_speed : chr "3,7789" "7,02062" "8,32836" "3,87037" ...
## $ wind_direction : chr "55,099998" "55,099998" "55,099998" "55,099998" ...
## $ dt_wind : chr "Northeast" "Northeast" "Northeast" "Northeast" ...
## $ solar_aptitude : chr "0,746197" "0,8" "0,726996" "0,595309" ...
## $ solar_aptitude_rounded: int 7 8 7 6 7 7 7 8 7 8 ...
## $ solar_aptittude_class : chr "Alta" "Alta" "Alta" "Media" ...
## $ capacity : chr "25" "101" "107" "180" ...
## $ optimal_tilt : int 31 26 26 31 33 30 31 29 31 27 ...
## $ pv_potential : chr "4,983" "6,389" "6,392" "4,969" ...
##
## 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 la variable
Curvature <- Datos$curvature_type
# EDA variable nominal
TDF_curvature <- table(Curvature)
tabla_curvature <- as.data.frame(TDF_curvature)
hi <- tabla_curvature$Freq/sum(tabla_curvature$Freq)
hi_porc <- hi*100
sum(hi_porc)## [1] 100
### Agrupación ###
tabla_CURVATURE$grupo <- case_when(
grepl("Valles", tabla_CURVATURE$TipoCurvatura, ignore.case = TRUE) ~ "Superficies cóncavas / Valles",
grepl("convexas", tabla_CURVATURE$TipoCurvatura, ignore.case = TRUE) ~ "Superficies convexas / Crestas",
grepl("planas", tabla_CURVATURE$TipoCurvatura, ignore.case = TRUE) ~ "Superficies planas o intermedias",
TRUE ~ "Otros")
tabla_resumen <- tabla_CURVATURE %>%
group_by(grupo) %>%
summarise(
Frecuencia = sum(Freq),
Porcentaje = sum(hi_porc)) %>%
arrange(desc(Frecuencia))
# Renombramos columnas
colnames(tabla_resumen) <- c("Curvatura","ni","hi (%)")
# Tabla Intermedia GT
tabla_resumen_gt <- tabla_resumen %>%
gt() %>%
tab_header(
title = md("**Tabla N°1 de Agrupación por Tipo de Curvatura de las Plantas Solares**")) %>%
tab_source_note(source_note = "Autor: Martin Sarmiento") %>%
cols_label(
Curvatura = "Tipo de Curvatura",
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 N°1 de Agrupación por Tipo de Curvatura de las Plantas Solares | ||
| Tipo de Curvatura | Frecuencia (ni) | Porcentaje (hi%) |
|---|---|---|
| Superficies planas o intermedias | 4269 | 84.12 |
| Superficies cóncavas / Valles | 405 | 7.98 |
| Superficies convexas / Crestas | 401 | 7.90 |
| Autor: Martin Sarmiento | ||
#### Crear fila de totales ####
totales <- c("TOTAL", sum(tabla_resumen$ni), sum(tabla_resumen$`hi (%)`))
tabla_Curvature_Final <- rbind(tabla_resumen, totales)
# Convertir a números para GT
tabla_Curvature_Final$ni <- as.numeric(tabla_Curvature_Final$ni)
tabla_Curvature_Final$`hi (%)` <- as.numeric(tabla_Curvature_Final$`hi (%)`)
# TABLA 1
tabla_final_gt <- tabla_Curvature_Final %>%
gt() %>%
tab_header(
title = md("**Tabla N°2 de Distribución de Frecuencias de Tipo de Curvatura de las Plantas Solares**")) %>%
tab_source_note(source_note = "Autor: Martin Sarmiento") %>%
cols_label(
Curvatura = "Tipo de Curvatura",
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 N°2 de Distribución de Frecuencias de Tipo de Curvatura de las Plantas Solares | ||
| Tipo de Curvatura | Frecuencia (ni) | Porcentaje (hi%) |
|---|---|---|
| Superficies planas o intermedias | 4269 | 84.12 |
| Superficies cóncavas / Valles | 405 | 7.98 |
| Superficies convexas / Crestas | 401 | 7.90 |
| TOTAL | 5075 | 100.00 |
| Autor: Martin Sarmiento | ||
par(mar = c(13, 4, 4, 2))
barplot(tabla_resumen$ni,
main="",
xlab = "",
ylab = "Cantidad",
col = "#EEE685",
names.arg=tabla_resumen$Curvatura,
cex.names = 0.7,
las = 2)
mtext("Tipo de Curvatura", side = 1, line = 11)
mtext("Gráfica N°1: Distribución de Cantidad de Plantas Solares por Tipo de Curvatura",
side = 3,
line = 2,
adj = 0.5,
cex = 1,
font = 2)par(mar = c(13, 4, 4, 2))
barplot(tabla_resumen$ni,
main="",
xlab = "",
ylab = "Cantidad",
col = "#EEE685",
ylim = c(0,8000),
names.arg=tabla_resumen$Curvatura,
cex.names = 0.7,
las = 2)
mtext("Tipo de Curvatura", side = 1, line = 11)
mtext("Gráfica N°2: Distribución de Cantidad de Plantas Solares por Tipo de Curvatura",
side = 3,
line = 2,
adj = 0.5,
cex = 1,
font = 2)par(mar = c(13, 4, 4, 2))
barplot(tabla_resumen$`hi (%)`,
main="",
xlab = "",
ylab = "Porcentaje %",
col = "#EEE685",
names.arg=tabla_resumen$Curvatura,
cex.names = 0.7,
las = 2)
mtext("Tipo de Curvatura", side = 1, line = 11)
mtext("Gráfica N°3: Distribución Porcentual de las Plantas Solares por Tipo de Curvatura",
side = 3,
line = 2,
adj = 0.5,
cex = 1,
font = 2)par(mar = c(13, 4, 4, 2))
barplot(tabla_resumen$`hi (%)`,
main="",
xlab = "",
ylab = "Porcentaje %",
col = "#EEE685",
ylim = c(0,100),
names.arg=tabla_resumen$Curvatura,
cex.names = 0.7,
las = 2)
mtext("Tipo de Curvatura", side = 1, line = 11)
mtext("Gráfica N°4: Distribución Porcentual de las Plantas Solares por Tipo de Curvatura",
side = 3,
line = 2,
adj = 0.5,
cex = 1,
font = 2)par(mar = c(4, 4, 4, 4), xpd = TRUE)
colores <- c("#FFFACD","#F0E68C", "#8B864E")
colores <- rep(colores, length.out = length(tabla_resumen$Curvatura))
pie(tabla_resumen$`hi (%)`,
main = "",
radius = 0.9,
labels = paste0(round(tabla_resumen$`hi (%)`, 1), "%"),
col = colores,
cex = 0.7)
mtext("Gráfica N°5: Distribución Porcentual de las Plantas Solares por Tipo de Curvatura",
side = 3,
line = 2,
adj = 0.5,
cex = 1,
font = 2)
legend(x = 1.1, y = 1.1,
legend = tabla_resumen$Curvatura,
fill = colores,
cex = 0.6,
title = "Tipo de Curvatura",
bty = "n")# Cálculo de la Moda
moda_curvatura <- tabla_resumen$Curvatura[which.max(tabla_resumen$ni)]
# Tabla de Indicadores
tabla_indicadores <- data.frame(
"Variable" = "Tipo de Curvatura",
"Rango" = "Cóncavo, Convexo, Plano",
"Media (X)" = "-",
"Mediana (Me)" = "-",
"Moda (Mo)" = moda_curvatura,
"Varianza (V)" = "-",
"Desv. Est. (Sd)" = "-",
"C.V. (%)" = "-",
"Asimetría (As)" = "-",
"Curtosis (K)" = "-",
check.names = FALSE)
# Generar Tabla de Indicadores
tabla_conclusiones_gt <- tabla_indicadores %>%
gt() %>%
tab_header(title = md("**Tabla N°3 de Conclusiones de Tipo de Curvatura de las Plantas Solares**")) %>%
tab_source_note(source_note = "Autor: Martin Sarmiento") %>%
tab_options(column_labels.background.color = "#F0F0F0")
tabla_conclusiones_gt| Tabla N°3 de Conclusiones de Tipo de Curvatura de las Plantas Solares | |||||||||
| Variable | Rango | Media (X) | Mediana (Me) | Moda (Mo) | Varianza (V) | Desv. Est. (Sd) | C.V. (%) | Asimetría (As) | Curtosis (K) |
|---|---|---|---|---|---|---|---|---|---|
| Tipo de Curvatura | Cóncavo, Convexo, Plano | - | - | Superficies planas o intermedias | - | - | - | - | - |
| Autor: Martin Sarmiento | |||||||||
La variable “Tipo de Curvatura”, presenta como valor más frecuente es Superficies planas o intermedias, con una participación destacada en la muestra.