##### UNIVERSIDAD CENTRAL DEL ECUADOR #####
#### AUTOR: Liss Murillo ####
### CARRERA: INGENIERÍA EN PETROLEOS #####
###Grupo 4###
# 1. CARGAR DATOS
library(readxl)
datos <- read_excel("r-graficas/Producción Campo Sacha.csv.xlsx")
View(datos)
str(datos)
## tibble [8,344 × 31] (S3: tbl_df/tbl/data.frame)
## $ mes : chr [1:8344] "Ene" "Ene" "Ene" "Ene" ...
## $ día : num [1:8344] 1 1 1 1 1 1 1 1 1 1 ...
## $ Pozo : chr [1:8344] "SACHA-001A" "SACHA-019A" "SACHA-052B" "SACHA-083A" ...
## $ Campo : chr [1:8344] "SACHA" "SACHA" "SACHA" "SACHA" ...
## $ Reservorio : chr [1:8344] "U" "U" "U INFERIOR" "HOLLIN INFERIOR" ...
## $ Bpd : num [1:8344] NA 53 249 139 186 136 NA 456 161 164 ...
## $ Bppd_BH : num [1:8344] 159 NA NA NA NA NA 155 NA NA NA ...
## $ Bfpd_BE : num [1:8344] NA 534 346 1158 1163 ...
## $ Bfpd_BH : num [1:8344] 695 NA NA NA NA NA 441 NA NA NA ...
## $ Bapd_BE : num [1:8344] NA 481 97 1019 977 ...
## $ Bapd_BH : num [1:8344] 536 NA NA NA NA NA 286 NA NA NA ...
## $ Bsw_BE : num [1:8344] NA 90.1 28 88 84 ...
## $ Bsw_BH : num [1:8344] 77.1 NA NA NA NA ...
## $ Api_BE : num [1:8344] NA 26.7 27.8 27.7 24 20.5 NA 28.5 29.9 26.3 ...
## $ Api_BH : num [1:8344] 27.8 NA NA NA NA NA 23.2 NA NA NA ...
## $ Gas_BE : num [1:8344] NA 10.76 50.55 1.11 27.9 ...
## $ Gas_BH : num [1:8344] 32.3 NA NA NA NA ...
## $ Salinidad_BE : num [1:8344] NA 15920 30227 1600 13000 ...
## $ Salinidad_BH : num [1:8344] 10800 NA NA NA NA NA 3800 NA NA NA ...
## $ Rgl_BE : num [1:8344] NA 20.15 146.1 0.96 23.99 ...
## $ Rgl_BH : num [1:8344] 46.5 NA NA NA NA ...
## $ Gor_BE : num [1:8344] NA 203.02 203.01 7.99 150 ...
## $ Gor_BH : num [1:8344] 203 NA NA NA NA ...
## $ Horas_BE : num [1:8344] NA 4 5 4 4 10 NA 4 10 10 ...
## $ Horas_BH : num [1:8344] 4 NA NA NA NA NA 4 NA NA NA ...
## $ Bomba_BE : chr [1:8344] NA "SF-320|SF-320|SF-900|SFGH2500/520/180/9259" "RC 1000|RC 1000|RC 1000/300/120/9250" "P23/68/30/7000" ...
## $ Bomba_BH : chr [1:8344] "JET 12K/0//0" NA NA NA ...
## $ Frecuencia Operaciones: num [1:8344] NA 65 62 46 59 52 NA 58.5 57 54 ...
## $ Voltaje : num [1:8344] NA 479 457 364 440 452 NA 475 455 439 ...
## $ Amperaje : num [1:8344] NA 29 35 14 59 30 NA 23 35 34 ...
## $ Presión Intake : num [1:8344] NA 484 406 0 345 162 NA 546 338 0 ...
# 2. EXTRAER Gas_BE
if(!"Gas_BE" %in% colnames(datos)) {
col_gas <- grep("gas_be|gas be|gasbe|gas", colnames(datos), ignore.case = TRUE, value = TRUE)[1]
datos$Gas_BE <- datos[[col_gas]]
}
# Extraer variable
Gas_BE <- datos$Gas_BE
# Convertir a numérico y eliminar NAs
Gas_BE <- as.numeric(Gas_BE)
ValorNulo <- is.na(Gas_BE)
Gas_BE <- na.omit(Gas_BE)
# 3. CALCULAR INTERVALOS (STURGES)
n <- length(Gas_BE)
k <- 1 + (3.3 * log10(n))
k <- floor(k)
min_val <- min(Gas_BE)
max_val <- max(Gas_BE)
R <- max_val - min_val
A <- R / k
A <- round(A, 2)
min_val <- round(min_val, 2)
max_val <- round(max_val, 2)
# 4. CREAR LÍMITES DE INTERVALOS
Li <- round(seq(from = min_val, to = max_val - A, by = A), 4)
Ls <- round(seq(from = min_val + A, to = max_val, by = A), 4)
MC <- round((Li + Ls) / 2, 2)
# 5. FRECUENCIAS
ni <- numeric(length(Li))
for (i in 1:length(Li)) {
if (i < length(Li)) {
ni[i] <- sum(Gas_BE >= Li[i] & Gas_BE < Ls[i])
} else {
ni[i] <- sum(Gas_BE >= Li[i] & Gas_BE <= max_val)
}
}
# 6. CALCULAR TODAS LAS FRECUENCIAS
hi <- ni / sum(ni) * 100
Niasc <- cumsum(ni)
Nidsc <- rev(cumsum(rev(ni)))
Hiasc <- cumsum(hi)
Hidsc <- rev(cumsum(rev(hi)))
TDF_Gas_BE_es <- data.frame(
Li = Li,
Ls = Ls,
MC = MC,
"ni (FA)" = ni,
"hi (FR)" = round(hi, 4),
"Ni (FAAa)" = Niasc,
"Hi (FRAa)" = round(Hiasc, 4),
"Ni (FAAd)" = Nidsc,
"Hi (FRAd)" = round(Hidsc, 4)
)
# 7. MOSTRAR TABLA DE DISTRIBUCIÓN
cat("\n" , rep("=", 80), "\n", sep = "")
##
## ================================================================================
cat("TABLA 1: DISTRIBUCIÓN DE FRECUENCIAS - Gas_BE\n")
## TABLA 1: DISTRIBUCIÓN DE FRECUENCIAS - Gas_BE
cat(rep("=", 80), "\n\n", sep = "")
## ================================================================================
print(TDF_Gas_BE_es)
## Li Ls MC ni..FA. hi..FR. Ni..FAAa. Hi..FRAa. Ni..FAAd.
## 1 0.02 95.02 47.52 5531 71.7846 5531 71.7846 7705
## 2 95.02 190.02 142.52 1053 13.6665 6584 85.4510 2174
## 3 190.02 285.02 237.52 461 5.9831 7045 91.4341 1121
## 4 285.02 380.02 332.52 318 4.1272 7363 95.5613 660
## 5 380.02 475.02 427.52 217 2.8164 7580 98.3777 342
## 6 475.02 570.02 522.52 44 0.5711 7624 98.9487 125
## 7 570.02 665.02 617.52 10 0.1298 7634 99.0785 81
## 8 665.02 760.02 712.52 31 0.4023 7665 99.4809 71
## 9 760.02 855.02 807.52 26 0.3374 7691 99.8183 40
## 10 855.02 950.02 902.52 0 0.0000 7691 99.8183 14
## 11 950.02 1045.02 997.52 0 0.0000 7691 99.8183 14
## 12 1045.02 1140.02 1092.52 14 0.1817 7705 100.0000 14
## Hi..FRAd.
## 1 100.0000
## 2 28.2154
## 3 14.5490
## 4 8.5659
## 5 4.4387
## 6 1.6223
## 7 1.0513
## 8 0.9215
## 9 0.5191
## 10 0.1817
## 11 0.1817
## 12 0.1817
# =============================================================================
# GRÁFICAS CORTAS: LOCAL vs GLOBAL
# =============================================================================
cat("\n" , rep("=", 60), "\n", sep = "")
##
## ============================================================
cat("GRÁFICAS: LOCAL vs GLOBAL\n")
## GRÁFICAS: LOCAL vs GLOBAL
cat(rep("=", 60), "\n")
## = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =
# 1. PREPARAR
if("Pozo" %in% colnames(datos)) {
pozos <- unique(datos$Pozo)
local_pozos <- pozos[1:min(9, length(pozos))]
Gas_local <- na.omit(as.numeric(datos$Gas_BE[datos$Pozo %in% local_pozos]))
Gas_global <- Gas_BE
cat("LOCAL:", length(local_pozos), "pozos,", length(Gas_local), "datos\n")
cat("GLOBAL:", length(Gas_global), "datos\n\n")
# 2. CREAR 2 GRÁFICAS
par(mfrow = c(1, 2))
# Gráfica LOCAL
hist(Gas_local, breaks = 10, col = rgb(0.2,0.4,0.8,0.6), border = "white",
main = paste("LOCAL\n", length(local_pozos), " pozos", sep = ""),
xlab = "Gas", ylab = "Frecuencia", las = 1)
# Gráfica GLOBAL
hist(Gas_global, breaks = k, col = rgb(0.8,0.4,0.2,0.6), border = "white",
main = "GLOBAL\nTodos los pozos",
xlab = "Gas", ylab = "Frecuencia", las = 1)
par(mfrow = c(1, 1))
} else {
cat("SOLO GLOBAL (no hay columna Pozo)\n\n")
hist(Gas_BE, breaks = k, col = "lightblue", border = "white",
main = "Gas_BE (GLOBAL)", xlab = "Gas", ylab = "Frecuencia", las = 1)
}
## LOCAL: 9 pozos, 300 datos
## GLOBAL: 7705 datos
# =============================================================================
# GRÁFICA OJIVA ASCENDENTE
# =============================================================================
cat("\n" , rep("=", 60), "\n", sep = "")
##
## ============================================================
cat("OJIVA ASCENDENTE - Gas_BE\n")
## OJIVA ASCENDENTE - Gas_BE
cat(rep("=", 60), "\n\n", sep = "")
## ============================================================
# Crear puntos para la ojiva
puntos_x <- c(Ls)
puntos_y <- cumsum(ni)
# 1. Configurar márgenes para espacio a la derecha
par(mar = c(5, 4, 4, 8))
# 2. Crear ojiva
plot(puntos_x, puntos_y,
type = "b",
pch = 16,
col = "#A0F0ED",
lwd = 2,
main = "OJIVA ASCENDENTE\nGas_BE",
xlab = "Límite Superior del Intervalo",
ylab = "Frecuencia Acumulada",
las = 1,
xlim = c(min_val, max_val),
ylim = c(0, n * 1.05))
# 3. Grid
grid(col = "gray80", lty = "dotted")
# 4. Eje derecho con porcentajes
porcentajes <- c(0, 25, 50, 75, 100)
axis(4,
at = porcentajes/100 * n,
labels = paste0(porcentajes, "%"),
las = 1,
col.axis = "darkred")
mtext("Porcentaje Acumulado", side = 4, line = 3, col = "darkred")
# 5. Puntos importantes (50% y 75%)
# 50%
indice_50 <- which.max(Hiasc >= 50)
if(length(indice_50) > 0) {
x50 <- puntos_x[indice_50 + 1]
y50 <- puntos_y[indice_50 + 1]
points(x50, y50, col = "red", pch = 17, cex = 1.3)
segments(x50, 0, x50, y50, col = "red", lty = 2)
segments(min_val, y50, x50, y50, col = "red", lty = 2)
}
# 75%
indice_75 <- which.max(Hiasc >= 75)
if(length(indice_75) > 0) {
x75 <- puntos_x[indice_75 + 1]
y75 <- puntos_y[indice_75 + 1]
points(x75, y75, col = "darkgreen", pch = 15, cex = 1.3)
segments(x75, 0, x75, y75, col = "darkgreen", lty = 2)
segments(min_val, y75, x75, y75, col = "darkgreen", lty = 2)
}
# =============================================================================
# DIAGRAMA DE CAJA
# =============================================================================
cat("\n" , rep("=", 60), "\n", sep = "")
##
## ============================================================
cat("4.2 Diagrama de caja\n")
## 4.2 Diagrama de caja
cat(rep("=", 60), "\n\n", sep = "")
## ============================================================
# Exactamente como el ejemplo
boxplot(Gas_BE,
horizontal = TRUE,
col = "#5AD9D2",
main = "Distribución de Gas (Gas_BE)",
xlab = "Gas")
cat("ESTADÍSTICAS:\n")
## ESTADÍSTICAS:
s <- boxplot.stats(Gas_BE)$stats
cat("• Mín:", round(s[1], 2), "\n")
## • Mín: 0.02
cat("• Q1: ", round(s[2], 2), "\n")
## • Q1: 26.73
cat("• Med:", round(s[3], 2), "\n")
## • Med: 52.17
cat("• Q3: ", round(s[4], 2), "\n")
## • Q3: 109.73
cat("• Máx:", round(s[5], 2), "\n")
## • Máx: 234
# =============================================================================
# 5. INDICADORES ESTADÍSTICOS - Gas_BE
# =============================================================================
cat("\n" , rep("=", 70), "\n", sep = "")
##
## ======================================================================
cat("5. INDICADORES ESTADÍSTICOS\n")
## 5. INDICADORES ESTADÍSTICOS
cat(rep("=", 70), "\n\n", sep = "")
## ======================================================================
# Cargar librería para asimetría y curtosis (e1071 tiene skewness y kurtosis)
if(!require(e1071)) {
install.packages("e1071")
library(e1071)
}
## Cargando paquete requerido: e1071
# =============================================================================
# 5.1 MEDIDAS DE POSICIÓN
# =============================================================================
cat("5.1 MEDIDAS DE POSICIÓN:\n")
## 5.1 MEDIDAS DE POSICIÓN:
cat("------------------------\n")
## ------------------------
# Asegurar que Gas_BE esté definido y sea numérico
Gas_BE <- as.numeric(Gas_BE)
Gas_BE <- na.omit(Gas_BE)
# Tamaño de muestra
n <- length(Gas_BE)
# Media aritmética
x <- mean(Gas_BE)
cat("Media aritmética (x̄):", round(x, 4), "Gas\n")
## Media aritmética (x̄): 98.3445 Gas
# Mediana (de los datos originales)
Me <- median(Gas_BE)
cat("Mediana (Me): ", round(Me, 4), "Gas\n")
## Mediana (Me): 52.17 Gas
# Mediana de las marcas de clase (opcional, si MC existe)
if(exists("MC")) {
Me_MC <- median(MC)
cat("Mediana de MC: ", round(Me_MC, 4), "Gas\n")
} else {
Me_MC <- NA
}
## Mediana de MC: 570.02 Gas
# Mínimo y máximo
ri <- min(Gas_BE)
rs <- max(Gas_BE)
cat("Mínimo (ri): ", round(ri, 4), "Gas\n")
## Mínimo (ri): 0.02 Gas
cat("Máximo (rs): ", round(rs, 4), "Gas\n")
## Máximo (rs): 1235 Gas
# Moda (aproximada)
calcular_moda <- function(x) {
freq <- table(round(x, 1)) # Redondear a 1 decimal
moda_val <- as.numeric(names(freq)[which.max(freq)])
return(moda_val)
}
Mo <- calcular_moda(Gas_BE)
cat("Moda aproximada (Mo): ", round(Mo, 4), "Gas\n")
## Moda aproximada (Mo): 1.1 Gas
# Cuartiles
Q1 <- quantile(Gas_BE, 0.25)
Q2 <- quantile(Gas_BE, 0.50) # Igual que mediana
Q3 <- quantile(Gas_BE, 0.75)
cat("\nCuartiles:\n")
##
## Cuartiles:
cat(" Q1 (25%): ", round(Q1, 4), "Gas\n")
## Q1 (25%): 26.73 Gas
cat(" Q2 (50% - Mediana): ", round(Q2, 4), "Gas\n")
## Q2 (50% - Mediana): 52.17 Gas
cat(" Q3 (75%): ", round(Q3, 4), "Gas\n")
## Q3 (75%): 109.73 Gas
# Percentiles
percentiles <- quantile(Gas_BE, probs = c(0.10, 0.25, 0.50, 0.75, 0.90, 0.95, 0.99))
cat("\nPercentiles seleccionados:\n")
##
## Percentiles seleccionados:
for(i in 1:length(percentiles)) {
cat(" P", names(percentiles)[i], ": ",
round(percentiles[i], 4), "Gas\n", sep = "")
}
## P10%: 9.328Gas
## P25%: 26.73Gas
## P50%: 52.17Gas
## P75%: 109.73Gas
## P90%: 251.718Gas
## P95%: 356.88Gas
## P99%: 576.04Gas
# =============================================================================
# 5.2 MEDIDAS DE DISPERSIÓN
# =============================================================================
cat("\n\n5.2 MEDIDAS DE DISPERSIÓN:\n")
##
##
## 5.2 MEDIDAS DE DISPERSIÓN:
cat("---------------------------\n")
## ---------------------------
# Rango
R <- rs - ri
cat("Rango (R): ", round(R, 4), "Gas\n")
## Rango (R): 1234.98 Gas
# Rango intercuartílico
RIQ <- Q3 - Q1
cat("Rango intercuartílico:", round(RIQ, 4), "Gas\n")
## Rango intercuartílico: 83 Gas
# Varianza
varianza <- var(Gas_BE)
cat("Varianza (s²): ", round(varianza, 4), "\n")
## Varianza (s²): 16194.58
# Desviación estándar
s <- sd(Gas_BE)
cat("Desviación estándar (s):", round(s, 4), "Gas\n")
## Desviación estándar (s): 127.2579 Gas
# Coeficiente de variación
CV <- (s / x) * 100
cat("Coeficiente de variación (CV):", round(CV, 2), "%\n")
## Coeficiente de variación (CV): 129.4 %
# Error estándar de la media
EE <- s / sqrt(n)
cat("Error estándar de la media: ", round(EE, 4), "Gas\n")
## Error estándar de la media: 1.4498 Gas
# =============================================================================
# 5.3 MEDIDAS DE FORMA
# =============================================================================
cat("\n\n5.3 MEDIDAS DE FORMA:\n")
##
##
## 5.3 MEDIDAS DE FORMA:
cat("---------------------\n")
## ---------------------
# Coeficiente de asimetría (skewness)
As <- skewness(Gas_BE)
cat("Coeficiente de asimetría (As):", round(As, 4), "\n")
## Coeficiente de asimetría (As): 3.0686
# Interpretación de asimetría
if(abs(As) < 0.5) {
cat(" Interpretación: Distribución aproximadamente simétrica\n")
} else if(As > 0) {
cat(" Interpretación: Asimetría positiva (cola a la derecha)\n")
} else {
cat(" Interpretación: Asimetría negativa (cola a la izquierda)\n")
}
## Interpretación: Asimetría positiva (cola a la derecha)
# Coeficiente de curtosis (kurtosis)
k <- kurtosis(Gas_BE)
cat("\nCoeficiente de curtosis (k): ", round(k, 4), "\n")
##
## Coeficiente de curtosis (k): 13.8933
# Curtosis exceso (comparada con normal = 3)
k_exceso <- k - 3
cat("Curtosis exceso (k-3): ", round(k_exceso, 4), "\n")
## Curtosis exceso (k-3): 10.8933
# Interpretación de curtosis
if(abs(k_exceso) < 0.5) {
cat(" Interpretación: Mesocúrtica (similar a normal)\n")
} else if(k_exceso > 0) {
cat(" Interpretación: Leptocúrtica (más picuda que normal)\n")
} else {
cat(" Interpretación: Platicúrtica (más aplanada que normal)\n")
}
## Interpretación: Leptocúrtica (más picuda que normal)
# =============================================================================
# 5.4 RESUMEN EN TABLA
# =============================================================================
cat("\n\n" , rep("-", 70), "\n", sep = "")
##
##
## ----------------------------------------------------------------------
cat("RESUMEN DE INDICADORES ESTADÍSTICOS - Gas_BE\n")
## RESUMEN DE INDICADORES ESTADÍSTICOS - Gas_BE
cat(rep("-", 70), "\n\n", sep = "")
## ----------------------------------------------------------------------
# Crear tabla resumen
resumen_estadisticas <- data.frame(
"Categoría" = c(
rep("POSICIÓN", 8),
rep("DISPERSIÓN", 6),
rep("FORMA", 3)
),
"Indicador" = c(
"n (tamaño muestra)",
"Media aritmética",
"Mediana",
"Moda aproximada",
"Mínimo",
"Máximo",
"Q1 (25%)",
"Q3 (75%)",
"Rango",
"Rango intercuartílico",
"Varianza",
"Desviación estándar",
"Coeficiente variación",
"Error estándar",
"Asimetría",
"Curtosis",
"Curtosis exceso"
),
"Valor" = c(
n,
round(x, 4),
round(Me, 4),
round(Mo, 4),
round(ri, 4),
round(rs, 4),
round(Q1, 4),
round(Q3, 4),
round(R, 4),
round(RIQ, 4),
round(varianza, 4),
round(s, 4),
paste0(round(CV, 2), "%"),
round(EE, 4),
round(As, 4),
round(k, 4),
round(k_exceso, 4)
),
"Unidad" = c(
"observaciones",
"Gas disuelto (Gas_BE)",
"Gas disuelto (Gas_BE)",
"Gas disuelto (Gas_BE)",
"Gas disuelto (Gas_BE)",
"Gas disuelto (Gas_BE)",
"Gas disuelto (Gas_BE)",
"Gas disuelto (Gas_BE)",
"Gas_BE",
"Gas_BE",
"(Gas_BE)²",
"Gas_BE",
"%",
"Gas_BE",
"adimensional",
"adimensional",
"adimensional"
)
)
print(resumen_estadisticas)
## Categoría Indicador Valor Unidad
## 1 POSICIÓN n (tamaño muestra) 7705 observaciones
## 2 POSICIÓN Media aritmética 98.3445 Gas disuelto (Gas_BE)
## 3 POSICIÓN Mediana 52.17 Gas disuelto (Gas_BE)
## 4 POSICIÓN Moda aproximada 1.1 Gas disuelto (Gas_BE)
## 5 POSICIÓN Mínimo 0.02 Gas disuelto (Gas_BE)
## 6 POSICIÓN Máximo 1235 Gas disuelto (Gas_BE)
## 7 POSICIÓN Q1 (25%) 26.73 Gas disuelto (Gas_BE)
## 8 POSICIÓN Q3 (75%) 109.73 Gas disuelto (Gas_BE)
## 9 DISPERSIÓN Rango 1234.98 Gas_BE
## 10 DISPERSIÓN Rango intercuartílico 83 Gas_BE
## 11 DISPERSIÓN Varianza 16194.5765 (Gas_BE)²
## 12 DISPERSIÓN Desviación estándar 127.2579 Gas_BE
## 13 DISPERSIÓN Coeficiente variación 129.4% %
## 14 DISPERSIÓN Error estándar 1.4498 Gas_BE
## 15 FORMA Asimetría 3.0686 adimensional
## 16 FORMA Curtosis 13.8933 adimensional
## 17 FORMA Curtosis exceso 10.8933 adimensional
# =============================================================================
# 6. TABLA RESUMEN FINAL
# =============================================================================
cat("\n" , rep("=", 70), "\n", sep = "")
##
## ======================================================================
cat("6. TABLA RESUMEN - Gas_BE\n")
## 6. TABLA RESUMEN - Gas_BE
cat(rep("=", 70), "\n\n", sep = "")
## ======================================================================
cat("La tabla anterior (Resumen de Indicadores Estadísticos) es la tabla resumen de Gas_BE.\n")
## La tabla anterior (Resumen de Indicadores Estadísticos) es la tabla resumen de Gas_BE.
cat("Contiene todas las medidas calculadas organizadas por categoría.\n")
## Contiene todas las medidas calculadas organizadas por categoría.
# =============================================================================
# 7. CONCLUSIÓN
# =============================================================================
cat("ANÁLISIS DEL GAS DISUELTO (Gas_BE) EN EL CAMPO SACHA:\n\n")
## ANÁLISIS DEL GAS DISUELTO (Gas_BE) EN EL CAMPO SACHA:
cat("1. El análisis de", n, "observaciones de Gas_BE muestra:\n")
## 1. El análisis de 7705 observaciones de Gas_BE muestra:
cat(" • Promedio:", round(x, 2), " (Nivel de gas ")
## • Promedio: 98.34 (Nivel de gas
# Clasificación simple según gas disuelto — puedes ajustarla si deseas
if(x < 50) {
cat("bajo")
} else if(x < 150) {
cat("medio")
} else if(x < 300) {
cat("alto")
} else {
cat("muy alto")
}
## medio
cat(")\n")
## )
cat(" • Variabilidad:", round(CV, 1), "% (")
## • Variabilidad: 129.4 % (
if(CV > 30) {
cat("alta")
} else if(CV > 15) {
cat("moderada")
} else {
cat("baja")
}
## alta
cat(")\n")
## )
cat(" • Distribución: ")
## • Distribución:
if(abs(As) < 0.5) {
cat("aproximadamente simétrica")
} else if(As > 0) {
cat("sesgo positivo (valores altos más frecuentes)")
} else {
cat("sesgo negativo (valores bajos más frecuentes)")
}
## sesgo positivo (valores altos más frecuentes)
cat("\n\n")
cat("2. El análisis estadístico completo proporciona una caracterización\n")
## 2. El análisis estadístico completo proporciona una caracterización
cat(" detallada del comportamiento del gas disuelto en el campo Sacha.\n")
## detallada del comportamiento del gas disuelto en el campo Sacha.
cat("\n✓ Análisis estadístico completado exitosamente para Gas_BE\n")
##
## ✓ Análisis estadístico completado exitosamente para Gas_BE