AMBIENT_TEMPERATURE
MARTIN SARMIENTO
2025-12-02
1 Configuración y Carga de Datos
##### UNIVERSIDAD CENTRAL DEL ECUADOR #####
#### AUTOR: MARTIN SARMIENTO ####
### CARRERA: INGENIERÍA EN PETRÓLEOS #####
#### VARIABLE TEMPERATURA AMBIENTE ####
## DATASET ##
setwd("~/R/AMBIENT_TEMPERATURE")
# 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, union2 Cálculo de Intervalos y Frecuencias
#Extraer variable
Variable <- na.omit(Datos$ambient_temperature)
N <- length(Variable)
# Cálculos básicos
min_val <- min(Variable)
max_val <- max(Variable)
Rango <- max_val - min_val
K <- floor(1 + 3.322 * log10(N))
Amplitud <- Rango / K
# Creación de Límites
lim_inf <- round(seq(from = min_val, to = max_val - Amplitud, by = Amplitud), 2)
# Ajuste para el último límite
lim_sup <- round(seq(from = min_val + Amplitud, to = max_val, by = Amplitud), 2)
lim_sup[K] <- round(max_val, 2)
# Marca de Clase
MC <- (lim_inf + lim_sup) / 2
# Conteo de Frecuencias (ni)
ni <- numeric(K)
for (i in 1:K) {
if (i < K) {
ni[i] <- length(subset(Variable, Variable >= lim_inf[i] & Variable < lim_sup[i]))
} else {
ni[i] <- length(subset(Variable, Variable >= lim_inf[i] & Variable <= lim_sup[i] + 0.001))
}
}
# Cálculos de Frecuencias
sum_ni <- sum(ni)
hi <- (ni / sum_ni) * 100
Ni_asc <- cumsum(ni)
Hi_asc <- cumsum(hi)
# Acumuladas Descendentes
Ni_desc <- rev(cumsum(rev(ni)))
Hi_desc <- rev(cumsum(rev(hi)))
# Construcción del Dataframe
TDF_Temperatura <- data.frame(
Li = lim_inf,
Ls = lim_sup,
MC = MC,
ni = ni,
hi = round(hi, 2),
Ni_asc = Ni_asc,
Ni_desc = Ni_desc,
Hi_asc = round(Hi_asc, 2),
Hi_desc = round(Hi_desc, 2)
)3 Tabla de Distribución de Frecuencias
#### Crear de fila de totales ####
totales <- c("TOTAL", "-", "-", sum(ni), sum(hi), "-", "-", "-", "-")
TDF_Temperatura_Char <- TDF_Temperatura %>% mutate(across(everything(), as.character))
TDF_Final <- rbind(TDF_Temperatura_Char, totales)
# TABLA DE FRECUENCIAS GT
tabla_gt <- TDF_Final %>%
gt() %>%
tab_header(
title = md("**TABLA DE DISTRIBUCIÓN DE FRECUENCIAS**"),
subtitle = "Variable: Temperatura Ambiente (°C)"
) %>%
tab_source_note(source_note = "Autor: Martin Sarmiento") %>%
cols_label(
Li = "Lim. Inf",
Ls = "Lim. Sup",
MC = "Marca Clase",
ni = "Frec. Abs (ni)",
hi = "Frec. Rel (%)",
Ni_asc = "Ni (Asc)",
Ni_desc = "Ni (Desc)",
Hi_asc = "Hi Asc (%)",
Hi_desc = "Hi Desc (%)"
) %>%
tab_options(
heading.title.font.size = px(16),
column_labels.background.color = "#f0f0f0",
table.font.size = px(14)
)
tabla_gt| TABLA DE DISTRIBUCIÓN DE FRECUENCIAS | ||||||||
| Variable: Temperatura Ambiente (°C) | ||||||||
| Lim. Inf | Lim. Sup | Marca Clase | Frec. Abs (ni) | Frec. Rel (%) | Ni (Asc) | Ni (Desc) | Hi Asc (%) | Hi Desc (%) |
|---|---|---|---|---|---|---|---|---|
| 0 | 2.32 | 1.16 | 116 | 1.62 | 116 | 7141 | 1.62 | 100 |
| 2.32 | 4.65 | 3.485 | 2 | 0.03 | 118 | 7025 | 1.65 | 98.38 |
| 4.65 | 6.97 | 5.81 | 17 | 0.24 | 135 | 7023 | 1.89 | 98.35 |
| 6.97 | 9.29 | 8.13 | 5 | 0.07 | 140 | 7006 | 1.96 | 98.11 |
| 9.29 | 11.62 | 10.455 | 7 | 0.1 | 147 | 7001 | 2.06 | 98.04 |
| 11.62 | 13.94 | 12.78 | 74 | 1.04 | 221 | 6994 | 3.09 | 97.94 |
| 13.94 | 16.26 | 15.1 | 145 | 2.03 | 366 | 6920 | 5.13 | 96.91 |
| 16.26 | 18.58 | 17.42 | 229 | 3.21 | 595 | 6775 | 8.33 | 94.87 |
| 18.58 | 20.91 | 19.745 | 364 | 5.1 | 959 | 6546 | 13.43 | 91.67 |
| 20.91 | 23.23 | 22.07 | 1027 | 14.38 | 1986 | 6182 | 27.81 | 86.57 |
| 23.23 | 25.55 | 24.39 | 1976 | 27.67 | 3962 | 5155 | 55.48 | 72.19 |
| 25.55 | 27.88 | 26.715 | 2232 | 31.26 | 6194 | 3179 | 86.74 | 44.52 |
| 27.88 | 30.2 | 29.04 | 947 | 13.26 | 7141 | 947 | 100 | 13.26 |
| TOTAL | - | - | 7141 | 100 | - | - | - | - |
| Autor: Martin Sarmiento | ||||||||
4 Gráfico 1 – Frecuencia Local
color_sutil <- "#EE5C42"
par(mar = c(8, 5, 4, 2))
barplot(TDF_Temperatura$ni,
names.arg = TDF_Temperatura$MC,
main = "Gráfica N°1: Frecuencia de Temperatura en Proyectos",
cex.main = 0.9,
xlab = "",
ylab = "Cantidad",
col = color_sutil,
space = 0,
las = 2,
cex.names = 0.7)
mtext("Temperatura (°C)", side = 1, line = 4)
5 Gráfico 2 – Frecuencia Global
color_grafico <- "#EE5C42"
barplot(TDF_Temperatura$ni,
main = "Gráfica N°2: Frecuencia Global de Temperatura",
xlab = "",
ylab = "Cantidad",
names.arg = TDF_Temperatura$MC,
col = color_sutil,
space = 0,
cex.main = 0.9,
cex.axis = 0.7,
cex.names = 0.7,
las = 2,
ylim = c(0, sum(TDF_Temperatura$ni)))
mtext("Temperatura (°C)", side = 1, line = 4)
6 Gráfico 3 – Porcentaje Local
color_grafico <- "#EE5C42"
par(mar = c(8, 5, 4, 2))
barplot(TDF_Temperatura$hi,
main = "Gráfica N°3: Porcentaje de Frecuencia de Temperatura",
xlab = "",
ylab = "Porcentaje (%)",
col = color_sutil,
space = 0,
names.arg = TDF_Temperatura$MC,
cex.main = 0.9,
cex.axis = 0.7,
cex.names = 0.7,
las = 2,
ylim = c(0, max(TDF_Temperatura$hi) * 1.1))
mtext("Temperatura (°C)", side = 1, line = 4)
7 Gráfico 4 – Porcentaje Global
color_grafico <- "#EE5C42"
par(mar = c(8, 5, 4, 2))
barplot(TDF_Temperatura$hi,
main = "Gráfica N°4: Porcentaje Global de Temperatura",
xlab = "",
ylab = "Porcentaje (%)",
col = color_sutil,
space = 0,
cex.main = 0.9,
names.arg = TDF_Temperatura$MC,
las = 2,
cex.names = 0.7,
ylim = c(0, 100))
mtext("Temperatura (°C)", side = 1, line = 4)
8 Gráfico 5 – Diagrama de Cajas (Boxplot)
par(mar = c(5, 5, 4, 2))
boxplot(Variable,
horizontal = TRUE,
col = color_sutil,
xlab = "Temperatura (°C)",
cex.main = 0.9,
main = "Gráfica N°5: Variabilidad de la Temperatura de los Proyectos")
9 Gráfico 6 – Ojivas de Frecuencia Acumulada
par(mar = c(5, 5, 4, 10), xpd = TRUE)
# Coordenadas
x_asc <- TDF_Temperatura$Ls
x_desc <- TDF_Temperatura$Li
y_asc <- TDF_Temperatura$Ni_asc
y_desc <- TDF_Temperatura$Ni_desc
# 1. Dibujar la Ascendente
plot(x_asc, y_asc,
type = "b",
main = "Gráfica N°6: Ojiva de Frecuencia Acumulada (Temperatura)",
cex.main = 0.9,
xlab = "Temperatura (°C)",
ylab = "Frecuencia acumulada",
col = "black",
pch = 19,
xlim = c(min(TDF_Temperatura$Li), max(x_asc)),
ylim = c(0, sum(ni)),
bty = "l"
)
# 2. Agregar la Descendente
lines(x_desc, y_desc, col = "#8B3626", type = "b", pch = 19)
grid()
legend("left",
legend = c("Ascendente", "Descendente"),
col = c("black", "#8B3626"),
lty = 1,
pch = 1,
cex = 0.6,
inset = c(0.05, 0.05),
bty = "n")
10 Indicadores Estadísticos
## INDICADORES DE TENDENCIA CENTRAL
# Media aritmética
media <- round(mean(Variable), 2)
# Moda
max_frecuencia <- max(TDF_Temperatura$ni)
moda_vals <- TDF_Temperatura$MC[TDF_Temperatura$ni == max_frecuencia]
moda_txt <- paste(round(moda_vals, 2), collapse = ", ")
# Mediana
mediana <- round(median(Variable), 2)
## INDICADORES DE DISPERSIÓN
# Varianza
varianza <- var(Variable)
# Desviación Estándar
sd_val <- sd(Variable)
# Coeficiente de Variación
cv <- round((sd_val / abs(media)) * 100, 2)
## INDICADORES DE FORMA
# Coeficiente de Asimetría
asimetria <- skewness(Variable, type = 2)
# Curtosis
curtosis <- kurtosis(Variable)
outliers <- boxplot.stats(Variable)$out
msg_atipicos <- if(length(outliers) > 0) "Presencia de valores atípicos" else "No hay presencia de valores atípicos"
tabla_indicadores <- data.frame(
"Variable" = c("Temperatura Ambiente (°C)"),
"Rango" = paste0("[", round(min(Variable), 2), "; ", round(max(Variable), 2), "]"),
"X" = c(media),
"Me" = c(round(mediana, 2)),
"Mo" = c(moda_txt),
"V" = c(round(varianza, 2)),
"Sd" = c(round(sd_val, 2)),
"Cv" = c(cv),
"As" = c(round(asimetria, 4)),
"K" = c(round(curtosis, 2)),
"Valores_Atipicos" = msg_atipicos
)
# Generar Tabla GT
tabla_conclusiones_gt <- tabla_indicadores %>%
gt() %>%
tab_header(
title = md("**CONCLUSIONES ESTADÍSTICAS**"),
subtitle = "Resumen de la variable Temperatura"
) %>%
tab_source_note(source_note = "Autor: Martin Sarmiento") %>%
cols_label(
Variable = "Variable",
Rango = "Rango",
X = "Media (X)",
Me = "Mediana (Me)",
Mo = "Moda (Mo)",
V = "Varianza (V)",
Sd = "Desv. Est. (Sd)",
Cv = "C.V. (%)",
As = "Asimetría (As)",
K = "Curtosis (K)",
Valores_Atipicos = "Valores Atípicos"
) %>%
tab_options(
heading.title.font.size = px(16),
column_labels.background.color = "#f0f0f0"
)
tabla_conclusiones_gt| CONCLUSIONES ESTADÍSTICAS | ||||||||||
| Resumen de la variable Temperatura | ||||||||||
| Variable | Rango | Media (X) | Mediana (Me) | Moda (Mo) | Varianza (V) | Desv. Est. (Sd) | C.V. (%) | Asimetría (As) | Curtosis (K) | Valores Atípicos |
|---|---|---|---|---|---|---|---|---|---|---|
| Temperatura Ambiente (°C) | [0; 30.2] | 24.14 | 25.2 | 26.72 | 21.23 | 4.61 | 19.09 | -2.8284 | 11.01 | Presencia de valores atípicos |
| Autor: Martin Sarmiento | ||||||||||