#==============================ENCABEZADO=======================================
# TEMA: ED Variables Continuas - PRODUCCION TOTAL DE GAS 2023
# AUTOR: GRUPO 4
# FECHA: 18-12-2025
#==============================CARGA DE DATOS===================================
library(gt)
library(dplyr)
##
## 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
library(moments)
# Funcion para calcular la moda
get_mode <- function(v) {
uniqv <- unique(v)
uniqv[which.max(tabulate(match(v, uniqv)))]
}
# Configuracion de directorio y carga
setwd("C:/Users/HP/Documents/PROYECTO ESTADISTICA/RStudio")
datos <- read.csv("tablap.csv", header = TRUE, dec = ",", sep = ";")
#================================================================================
# VARIABLE: PRODUCCION TOTAL DE GAS HACIA 2023
#================================================================================
# Limpieza: Conversion a numerico y eliminacion de NAs
gas_production <- as.numeric(gsub(",", ".", as.character(datos$Total.gas.production.by.2023)))
gas_production <- na.omit(gas_production)
# CALCULOS REGLA DE STURGES
R <- max(gas_production) - min(gas_production)
k <- floor(1 + (3.3 * log10(length(gas_production))))
A <- R / k
# Definicion de intervalos
liminf <- seq(from = min(gas_production), by = A, length.out = k)
limsup <- liminf + A
limsup[k] <- max(gas_production)
MC <- (liminf + limsup) / 2
# Calculo de frecuencias absolutas (ni)
ni <- c()
for (i in 1:k) {
if (i == k)
ni[i] <- length(subset(gas_production, gas_production >= liminf[i] & gas_production <= limsup[i]))
else
ni[i] <- length(subset(gas_production, gas_production >= liminf[i] & gas_production < limsup[i]))
}
# Calculo de frecuencias relativas y acumuladas
hi <- (ni / length(gas_production)) * 100
Niasc <- cumsum(ni)
Nidsc <- rev(cumsum(rev(ni)))
Hiasc <- cumsum(hi)
Hidsc <- rev(cumsum(rev(hi)))
#--- DATA FRAME PARA LA TABLA ---
TDFcu_gas_prod <- data.frame(
N = 1:k,
liminf = round(liminf, 2),
limsup = round(limsup, 2),
MC = round(MC, 2),
ni = ni,
hi_perc = round(hi, 2),
Niasc = Niasc,
Nidsc = Nidsc,
Hiasc_perc = round(Hiasc, 2),
Hidsc_perc = round(Hidsc, 2)
)
# Fila de totales
fila_total <- data.frame(
N = NA, liminf = NA, limsup = NA, MC = NA,
ni = sum(ni), hi_perc = 100,
Niasc = NA, Nidsc = NA, Hiasc_perc = NA, Hidsc_perc = NA
)
TDF_final_gt <- rbind(TDFcu_gas_prod, fila_total)
#============================== TABLA GT PROFESIONAL ==============================
tabla_gt_gas <- gt(TDF_final_gt) %>%
tab_header(
title = "Distribucion de Produccion Total de Gas",
) %>%
cols_label(
N = "Item",
liminf = "Lim. Inf.",
limsup = "Lim. Sup.",
MC = "Marca Clase",
ni = "ni",
hi_perc = "hi (%)",
Niasc = "Ni Asc.",
Nidsc = "Ni Desc.",
Hiasc_perc = "Hi Asc. (%)",
Hidsc_perc = "Hi Desc. (%)"
) %>%
sub_missing(columns = everything(), missing_text = "---") %>%
tab_style(
style = list(cell_fill(color = "lightcyan"), cell_text(weight = "bold")),
locations = cells_body(rows = is.na(N))
)
print(tabla_gt_gas)
## <div id="vslhcxfxpw" style="padding-left:0px;padding-right:0px;padding-top:10px;padding-bottom:10px;overflow-x:auto;overflow-y:auto;width:auto;height:auto;">
## <style>#vslhcxfxpw table {
## font-family: system-ui, 'Segoe UI', Roboto, Helvetica, Arial, sans-serif, 'Apple Color Emoji', 'Segoe UI Emoji', 'Segoe UI Symbol', 'Noto Color Emoji';
## -webkit-font-smoothing: antialiased;
## -moz-osx-font-smoothing: grayscale;
## }
##
## #vslhcxfxpw thead, #vslhcxfxpw tbody, #vslhcxfxpw tfoot, #vslhcxfxpw tr, #vslhcxfxpw td, #vslhcxfxpw th {
## border-style: none;
## }
##
## #vslhcxfxpw p {
## margin: 0;
## padding: 0;
## }
##
## #vslhcxfxpw .gt_table {
## display: table;
## border-collapse: collapse;
## line-height: normal;
## margin-left: auto;
## margin-right: auto;
## color: #333333;
## font-size: 16px;
## font-weight: normal;
## font-style: normal;
## background-color: #FFFFFF;
## width: auto;
## border-top-style: solid;
## border-top-width: 2px;
## border-top-color: #A8A8A8;
## border-right-style: none;
## border-right-width: 2px;
## border-right-color: #D3D3D3;
## border-bottom-style: solid;
## border-bottom-width: 2px;
## border-bottom-color: #A8A8A8;
## border-left-style: none;
## border-left-width: 2px;
## border-left-color: #D3D3D3;
## }
##
## #vslhcxfxpw .gt_caption {
## padding-top: 4px;
## padding-bottom: 4px;
## }
##
## #vslhcxfxpw .gt_title {
## color: #333333;
## font-size: 125%;
## font-weight: initial;
## padding-top: 4px;
## padding-bottom: 4px;
## padding-left: 5px;
## padding-right: 5px;
## border-bottom-color: #FFFFFF;
## border-bottom-width: 0;
## }
##
## #vslhcxfxpw .gt_subtitle {
## color: #333333;
## font-size: 85%;
## font-weight: initial;
## padding-top: 3px;
## padding-bottom: 5px;
## padding-left: 5px;
## padding-right: 5px;
## border-top-color: #FFFFFF;
## border-top-width: 0;
## }
##
## #vslhcxfxpw .gt_heading {
## background-color: #FFFFFF;
## text-align: center;
## border-bottom-color: #FFFFFF;
## border-left-style: none;
## border-left-width: 1px;
## border-left-color: #D3D3D3;
## border-right-style: none;
## border-right-width: 1px;
## border-right-color: #D3D3D3;
## }
##
## #vslhcxfxpw .gt_bottom_border {
## border-bottom-style: solid;
## border-bottom-width: 2px;
## border-bottom-color: #D3D3D3;
## }
##
## #vslhcxfxpw .gt_col_headings {
## border-top-style: solid;
## border-top-width: 2px;
## border-top-color: #D3D3D3;
## border-bottom-style: solid;
## border-bottom-width: 2px;
## border-bottom-color: #D3D3D3;
## border-left-style: none;
## border-left-width: 1px;
## border-left-color: #D3D3D3;
## border-right-style: none;
## border-right-width: 1px;
## border-right-color: #D3D3D3;
## }
##
## #vslhcxfxpw .gt_col_heading {
## color: #333333;
## background-color: #FFFFFF;
## font-size: 100%;
## font-weight: normal;
## text-transform: inherit;
## border-left-style: none;
## border-left-width: 1px;
## border-left-color: #D3D3D3;
## border-right-style: none;
## border-right-width: 1px;
## border-right-color: #D3D3D3;
## vertical-align: bottom;
## padding-top: 5px;
## padding-bottom: 6px;
## padding-left: 5px;
## padding-right: 5px;
## overflow-x: hidden;
## }
##
## #vslhcxfxpw .gt_column_spanner_outer {
## color: #333333;
## background-color: #FFFFFF;
## font-size: 100%;
## font-weight: normal;
## text-transform: inherit;
## padding-top: 0;
## padding-bottom: 0;
## padding-left: 4px;
## padding-right: 4px;
## }
##
## #vslhcxfxpw .gt_column_spanner_outer:first-child {
## padding-left: 0;
## }
##
## #vslhcxfxpw .gt_column_spanner_outer:last-child {
## padding-right: 0;
## }
##
## #vslhcxfxpw .gt_column_spanner {
## border-bottom-style: solid;
## border-bottom-width: 2px;
## border-bottom-color: #D3D3D3;
## vertical-align: bottom;
## padding-top: 5px;
## padding-bottom: 5px;
## overflow-x: hidden;
## display: inline-block;
## width: 100%;
## }
##
## #vslhcxfxpw .gt_spanner_row {
## border-bottom-style: hidden;
## }
##
## #vslhcxfxpw .gt_group_heading {
## padding-top: 8px;
## padding-bottom: 8px;
## padding-left: 5px;
## padding-right: 5px;
## color: #333333;
## background-color: #FFFFFF;
## font-size: 100%;
## font-weight: initial;
## text-transform: inherit;
## border-top-style: solid;
## border-top-width: 2px;
## border-top-color: #D3D3D3;
## border-bottom-style: solid;
## border-bottom-width: 2px;
## border-bottom-color: #D3D3D3;
## border-left-style: none;
## border-left-width: 1px;
## border-left-color: #D3D3D3;
## border-right-style: none;
## border-right-width: 1px;
## border-right-color: #D3D3D3;
## vertical-align: middle;
## text-align: left;
## }
##
## #vslhcxfxpw .gt_empty_group_heading {
## padding: 0.5px;
## color: #333333;
## background-color: #FFFFFF;
## font-size: 100%;
## font-weight: initial;
## border-top-style: solid;
## border-top-width: 2px;
## border-top-color: #D3D3D3;
## border-bottom-style: solid;
## border-bottom-width: 2px;
## border-bottom-color: #D3D3D3;
## vertical-align: middle;
## }
##
## #vslhcxfxpw .gt_from_md > :first-child {
## margin-top: 0;
## }
##
## #vslhcxfxpw .gt_from_md > :last-child {
## margin-bottom: 0;
## }
##
## #vslhcxfxpw .gt_row {
## padding-top: 8px;
## padding-bottom: 8px;
## padding-left: 5px;
## padding-right: 5px;
## margin: 10px;
## border-top-style: solid;
## border-top-width: 1px;
## border-top-color: #D3D3D3;
## border-left-style: none;
## border-left-width: 1px;
## border-left-color: #D3D3D3;
## border-right-style: none;
## border-right-width: 1px;
## border-right-color: #D3D3D3;
## vertical-align: middle;
## overflow-x: hidden;
## }
##
## #vslhcxfxpw .gt_stub {
## color: #333333;
## background-color: #FFFFFF;
## font-size: 100%;
## font-weight: initial;
## text-transform: inherit;
## border-right-style: solid;
## border-right-width: 2px;
## border-right-color: #D3D3D3;
## padding-left: 5px;
## padding-right: 5px;
## }
##
## #vslhcxfxpw .gt_stub_row_group {
## color: #333333;
## background-color: #FFFFFF;
## font-size: 100%;
## font-weight: initial;
## text-transform: inherit;
## border-right-style: solid;
## border-right-width: 2px;
## border-right-color: #D3D3D3;
## padding-left: 5px;
## padding-right: 5px;
## vertical-align: top;
## }
##
## #vslhcxfxpw .gt_row_group_first td {
## border-top-width: 2px;
## }
##
## #vslhcxfxpw .gt_row_group_first th {
## border-top-width: 2px;
## }
##
## #vslhcxfxpw .gt_summary_row {
## color: #333333;
## background-color: #FFFFFF;
## text-transform: inherit;
## padding-top: 8px;
## padding-bottom: 8px;
## padding-left: 5px;
## padding-right: 5px;
## }
##
## #vslhcxfxpw .gt_first_summary_row {
## border-top-style: solid;
## border-top-color: #D3D3D3;
## }
##
## #vslhcxfxpw .gt_first_summary_row.thick {
## border-top-width: 2px;
## }
##
## #vslhcxfxpw .gt_last_summary_row {
## padding-top: 8px;
## padding-bottom: 8px;
## padding-left: 5px;
## padding-right: 5px;
## border-bottom-style: solid;
## border-bottom-width: 2px;
## border-bottom-color: #D3D3D3;
## }
##
## #vslhcxfxpw .gt_grand_summary_row {
## color: #333333;
## background-color: #FFFFFF;
## text-transform: inherit;
## padding-top: 8px;
## padding-bottom: 8px;
## padding-left: 5px;
## padding-right: 5px;
## }
##
## #vslhcxfxpw .gt_first_grand_summary_row {
## padding-top: 8px;
## padding-bottom: 8px;
## padding-left: 5px;
## padding-right: 5px;
## border-top-style: double;
## border-top-width: 6px;
## border-top-color: #D3D3D3;
## }
##
## #vslhcxfxpw .gt_last_grand_summary_row_top {
## padding-top: 8px;
## padding-bottom: 8px;
## padding-left: 5px;
## padding-right: 5px;
## border-bottom-style: double;
## border-bottom-width: 6px;
## border-bottom-color: #D3D3D3;
## }
##
## #vslhcxfxpw .gt_striped {
## background-color: rgba(128, 128, 128, 0.05);
## }
##
## #vslhcxfxpw .gt_table_body {
## border-top-style: solid;
## border-top-width: 2px;
## border-top-color: #D3D3D3;
## border-bottom-style: solid;
## border-bottom-width: 2px;
## border-bottom-color: #D3D3D3;
## }
##
## #vslhcxfxpw .gt_footnotes {
## color: #333333;
## background-color: #FFFFFF;
## border-bottom-style: none;
## border-bottom-width: 2px;
## border-bottom-color: #D3D3D3;
## border-left-style: none;
## border-left-width: 2px;
## border-left-color: #D3D3D3;
## border-right-style: none;
## border-right-width: 2px;
## border-right-color: #D3D3D3;
## }
##
## #vslhcxfxpw .gt_footnote {
## margin: 0px;
## font-size: 90%;
## padding-top: 4px;
## padding-bottom: 4px;
## padding-left: 5px;
## padding-right: 5px;
## }
##
## #vslhcxfxpw .gt_sourcenotes {
## color: #333333;
## background-color: #FFFFFF;
## border-bottom-style: none;
## border-bottom-width: 2px;
## border-bottom-color: #D3D3D3;
## border-left-style: none;
## border-left-width: 2px;
## border-left-color: #D3D3D3;
## border-right-style: none;
## border-right-width: 2px;
## border-right-color: #D3D3D3;
## }
##
## #vslhcxfxpw .gt_sourcenote {
## font-size: 90%;
## padding-top: 4px;
## padding-bottom: 4px;
## padding-left: 5px;
## padding-right: 5px;
## }
##
## #vslhcxfxpw .gt_left {
## text-align: left;
## }
##
## #vslhcxfxpw .gt_center {
## text-align: center;
## }
##
## #vslhcxfxpw .gt_right {
## text-align: right;
## font-variant-numeric: tabular-nums;
## }
##
## #vslhcxfxpw .gt_font_normal {
## font-weight: normal;
## }
##
## #vslhcxfxpw .gt_font_bold {
## font-weight: bold;
## }
##
## #vslhcxfxpw .gt_font_italic {
## font-style: italic;
## }
##
## #vslhcxfxpw .gt_super {
## font-size: 65%;
## }
##
## #vslhcxfxpw .gt_footnote_marks {
## font-size: 75%;
## vertical-align: 0.4em;
## position: initial;
## }
##
## #vslhcxfxpw .gt_asterisk {
## font-size: 100%;
## vertical-align: 0;
## }
##
## #vslhcxfxpw .gt_indent_1 {
## text-indent: 5px;
## }
##
## #vslhcxfxpw .gt_indent_2 {
## text-indent: 10px;
## }
##
## #vslhcxfxpw .gt_indent_3 {
## text-indent: 15px;
## }
##
## #vslhcxfxpw .gt_indent_4 {
## text-indent: 20px;
## }
##
## #vslhcxfxpw .gt_indent_5 {
## text-indent: 25px;
## }
##
## #vslhcxfxpw .katex-display {
## display: inline-flex !important;
## margin-bottom: 0.75em !important;
## }
##
## #vslhcxfxpw div.Reactable > div.rt-table > div.rt-thead > div.rt-tr.rt-tr-group-header > div.rt-th-group:after {
## height: 0px !important;
## }
## </style>
## <table class="gt_table" data-quarto-disable-processing="false" data-quarto-bootstrap="false">
## <thead>
## <tr class="gt_heading">
## <td colspan="10" class="gt_heading gt_title gt_font_normal gt_bottom_border" style>Distribucion de Produccion Total de Gas</td>
## </tr>
##
## <tr class="gt_col_headings">
## <th class="gt_col_heading gt_columns_bottom_border gt_right" rowspan="1" colspan="1" scope="col" id="N">Item</th>
## <th class="gt_col_heading gt_columns_bottom_border gt_right" rowspan="1" colspan="1" scope="col" id="liminf">Lim. Inf.</th>
## <th class="gt_col_heading gt_columns_bottom_border gt_right" rowspan="1" colspan="1" scope="col" id="limsup">Lim. Sup.</th>
## <th class="gt_col_heading gt_columns_bottom_border gt_right" rowspan="1" colspan="1" scope="col" id="MC">Marca Clase</th>
## <th class="gt_col_heading gt_columns_bottom_border gt_right" rowspan="1" colspan="1" scope="col" id="ni">ni</th>
## <th class="gt_col_heading gt_columns_bottom_border gt_right" rowspan="1" colspan="1" scope="col" id="hi_perc">hi (%)</th>
## <th class="gt_col_heading gt_columns_bottom_border gt_right" rowspan="1" colspan="1" scope="col" id="Niasc">Ni Asc.</th>
## <th class="gt_col_heading gt_columns_bottom_border gt_right" rowspan="1" colspan="1" scope="col" id="Nidsc">Ni Desc.</th>
## <th class="gt_col_heading gt_columns_bottom_border gt_right" rowspan="1" colspan="1" scope="col" id="Hiasc_perc">Hi Asc. (%)</th>
## <th class="gt_col_heading gt_columns_bottom_border gt_right" rowspan="1" colspan="1" scope="col" id="Hidsc_perc">Hi Desc. (%)</th>
## </tr>
## </thead>
## <tbody class="gt_table_body">
## <tr><td headers="N" class="gt_row gt_right">1</td>
## <td headers="liminf" class="gt_row gt_right">25946</td>
## <td headers="limsup" class="gt_row gt_right">5024793</td>
## <td headers="MC" class="gt_row gt_right">2525369</td>
## <td headers="ni" class="gt_row gt_right">11815</td>
## <td headers="hi_perc" class="gt_row gt_right">94.06</td>
## <td headers="Niasc" class="gt_row gt_right">11815</td>
## <td headers="Nidsc" class="gt_row gt_right">12561</td>
## <td headers="Hiasc_perc" class="gt_row gt_right">94.06</td>
## <td headers="Hidsc_perc" class="gt_row gt_right">100.00</td></tr>
## <tr><td headers="N" class="gt_row gt_right">2</td>
## <td headers="liminf" class="gt_row gt_right">5024793</td>
## <td headers="limsup" class="gt_row gt_right">10023639</td>
## <td headers="MC" class="gt_row gt_right">7524216</td>
## <td headers="ni" class="gt_row gt_right">552</td>
## <td headers="hi_perc" class="gt_row gt_right">4.39</td>
## <td headers="Niasc" class="gt_row gt_right">12367</td>
## <td headers="Nidsc" class="gt_row gt_right">746</td>
## <td headers="Hiasc_perc" class="gt_row gt_right">98.46</td>
## <td headers="Hidsc_perc" class="gt_row gt_right">5.94</td></tr>
## <tr><td headers="N" class="gt_row gt_right">3</td>
## <td headers="liminf" class="gt_row gt_right">10023639</td>
## <td headers="limsup" class="gt_row gt_right">15022486</td>
## <td headers="MC" class="gt_row gt_right">12523063</td>
## <td headers="ni" class="gt_row gt_right">107</td>
## <td headers="hi_perc" class="gt_row gt_right">0.85</td>
## <td headers="Niasc" class="gt_row gt_right">12474</td>
## <td headers="Nidsc" class="gt_row gt_right">194</td>
## <td headers="Hiasc_perc" class="gt_row gt_right">99.31</td>
## <td headers="Hidsc_perc" class="gt_row gt_right">1.54</td></tr>
## <tr><td headers="N" class="gt_row gt_right">4</td>
## <td headers="liminf" class="gt_row gt_right">15022486</td>
## <td headers="limsup" class="gt_row gt_right">20021333</td>
## <td headers="MC" class="gt_row gt_right">17521910</td>
## <td headers="ni" class="gt_row gt_right">38</td>
## <td headers="hi_perc" class="gt_row gt_right">0.30</td>
## <td headers="Niasc" class="gt_row gt_right">12512</td>
## <td headers="Nidsc" class="gt_row gt_right">87</td>
## <td headers="Hiasc_perc" class="gt_row gt_right">99.61</td>
## <td headers="Hidsc_perc" class="gt_row gt_right">0.69</td></tr>
## <tr><td headers="N" class="gt_row gt_right">5</td>
## <td headers="liminf" class="gt_row gt_right">20021333</td>
## <td headers="limsup" class="gt_row gt_right">25020180</td>
## <td headers="MC" class="gt_row gt_right">22520756</td>
## <td headers="ni" class="gt_row gt_right">23</td>
## <td headers="hi_perc" class="gt_row gt_right">0.18</td>
## <td headers="Niasc" class="gt_row gt_right">12535</td>
## <td headers="Nidsc" class="gt_row gt_right">49</td>
## <td headers="Hiasc_perc" class="gt_row gt_right">99.79</td>
## <td headers="Hidsc_perc" class="gt_row gt_right">0.39</td></tr>
## <tr><td headers="N" class="gt_row gt_right">6</td>
## <td headers="liminf" class="gt_row gt_right">25020180</td>
## <td headers="limsup" class="gt_row gt_right">30019026</td>
## <td headers="MC" class="gt_row gt_right">27519603</td>
## <td headers="ni" class="gt_row gt_right">13</td>
## <td headers="hi_perc" class="gt_row gt_right">0.10</td>
## <td headers="Niasc" class="gt_row gt_right">12548</td>
## <td headers="Nidsc" class="gt_row gt_right">26</td>
## <td headers="Hiasc_perc" class="gt_row gt_right">99.90</td>
## <td headers="Hidsc_perc" class="gt_row gt_right">0.21</td></tr>
## <tr><td headers="N" class="gt_row gt_right">7</td>
## <td headers="liminf" class="gt_row gt_right">30019026</td>
## <td headers="limsup" class="gt_row gt_right">35017873</td>
## <td headers="MC" class="gt_row gt_right">32518450</td>
## <td headers="ni" class="gt_row gt_right">4</td>
## <td headers="hi_perc" class="gt_row gt_right">0.03</td>
## <td headers="Niasc" class="gt_row gt_right">12552</td>
## <td headers="Nidsc" class="gt_row gt_right">13</td>
## <td headers="Hiasc_perc" class="gt_row gt_right">99.93</td>
## <td headers="Hidsc_perc" class="gt_row gt_right">0.10</td></tr>
## <tr><td headers="N" class="gt_row gt_right">8</td>
## <td headers="liminf" class="gt_row gt_right">35017873</td>
## <td headers="limsup" class="gt_row gt_right">40016720</td>
## <td headers="MC" class="gt_row gt_right">37517296</td>
## <td headers="ni" class="gt_row gt_right">1</td>
## <td headers="hi_perc" class="gt_row gt_right">0.01</td>
## <td headers="Niasc" class="gt_row gt_right">12553</td>
## <td headers="Nidsc" class="gt_row gt_right">9</td>
## <td headers="Hiasc_perc" class="gt_row gt_right">99.94</td>
## <td headers="Hidsc_perc" class="gt_row gt_right">0.07</td></tr>
## <tr><td headers="N" class="gt_row gt_right">9</td>
## <td headers="liminf" class="gt_row gt_right">40016720</td>
## <td headers="limsup" class="gt_row gt_right">45015566</td>
## <td headers="MC" class="gt_row gt_right">42516143</td>
## <td headers="ni" class="gt_row gt_right">1</td>
## <td headers="hi_perc" class="gt_row gt_right">0.01</td>
## <td headers="Niasc" class="gt_row gt_right">12554</td>
## <td headers="Nidsc" class="gt_row gt_right">8</td>
## <td headers="Hiasc_perc" class="gt_row gt_right">99.94</td>
## <td headers="Hidsc_perc" class="gt_row gt_right">0.06</td></tr>
## <tr><td headers="N" class="gt_row gt_right">10</td>
## <td headers="liminf" class="gt_row gt_right">45015566</td>
## <td headers="limsup" class="gt_row gt_right">50014413</td>
## <td headers="MC" class="gt_row gt_right">47514990</td>
## <td headers="ni" class="gt_row gt_right">1</td>
## <td headers="hi_perc" class="gt_row gt_right">0.01</td>
## <td headers="Niasc" class="gt_row gt_right">12555</td>
## <td headers="Nidsc" class="gt_row gt_right">7</td>
## <td headers="Hiasc_perc" class="gt_row gt_right">99.95</td>
## <td headers="Hidsc_perc" class="gt_row gt_right">0.06</td></tr>
## <tr><td headers="N" class="gt_row gt_right">11</td>
## <td headers="liminf" class="gt_row gt_right">50014413</td>
## <td headers="limsup" class="gt_row gt_right">55013260</td>
## <td headers="MC" class="gt_row gt_right">52513837</td>
## <td headers="ni" class="gt_row gt_right">3</td>
## <td headers="hi_perc" class="gt_row gt_right">0.02</td>
## <td headers="Niasc" class="gt_row gt_right">12558</td>
## <td headers="Nidsc" class="gt_row gt_right">6</td>
## <td headers="Hiasc_perc" class="gt_row gt_right">99.98</td>
## <td headers="Hidsc_perc" class="gt_row gt_right">0.05</td></tr>
## <tr><td headers="N" class="gt_row gt_right">12</td>
## <td headers="liminf" class="gt_row gt_right">55013260</td>
## <td headers="limsup" class="gt_row gt_right">60012107</td>
## <td headers="MC" class="gt_row gt_right">57512683</td>
## <td headers="ni" class="gt_row gt_right">2</td>
## <td headers="hi_perc" class="gt_row gt_right">0.02</td>
## <td headers="Niasc" class="gt_row gt_right">12560</td>
## <td headers="Nidsc" class="gt_row gt_right">3</td>
## <td headers="Hiasc_perc" class="gt_row gt_right">99.99</td>
## <td headers="Hidsc_perc" class="gt_row gt_right">0.02</td></tr>
## <tr><td headers="N" class="gt_row gt_right">13</td>
## <td headers="liminf" class="gt_row gt_right">60012107</td>
## <td headers="limsup" class="gt_row gt_right">65010953</td>
## <td headers="MC" class="gt_row gt_right">62511530</td>
## <td headers="ni" class="gt_row gt_right">0</td>
## <td headers="hi_perc" class="gt_row gt_right">0.00</td>
## <td headers="Niasc" class="gt_row gt_right">12560</td>
## <td headers="Nidsc" class="gt_row gt_right">1</td>
## <td headers="Hiasc_perc" class="gt_row gt_right">99.99</td>
## <td headers="Hidsc_perc" class="gt_row gt_right">0.01</td></tr>
## <tr><td headers="N" class="gt_row gt_right">14</td>
## <td headers="liminf" class="gt_row gt_right">65010953</td>
## <td headers="limsup" class="gt_row gt_right">70009800</td>
## <td headers="MC" class="gt_row gt_right">67510377</td>
## <td headers="ni" class="gt_row gt_right">1</td>
## <td headers="hi_perc" class="gt_row gt_right">0.01</td>
## <td headers="Niasc" class="gt_row gt_right">12561</td>
## <td headers="Nidsc" class="gt_row gt_right">1</td>
## <td headers="Hiasc_perc" class="gt_row gt_right">100.00</td>
## <td headers="Hidsc_perc" class="gt_row gt_right">0.01</td></tr>
## <tr><td headers="N" class="gt_row gt_right" style="background-color: #E0FFFF; font-weight: bold;">—</td>
## <td headers="liminf" class="gt_row gt_right" style="background-color: #E0FFFF; font-weight: bold;">—</td>
## <td headers="limsup" class="gt_row gt_right" style="background-color: #E0FFFF; font-weight: bold;">—</td>
## <td headers="MC" class="gt_row gt_right" style="background-color: #E0FFFF; font-weight: bold;">—</td>
## <td headers="ni" class="gt_row gt_right" style="background-color: #E0FFFF; font-weight: bold;">12561</td>
## <td headers="hi_perc" class="gt_row gt_right" style="background-color: #E0FFFF; font-weight: bold;">100.00</td>
## <td headers="Niasc" class="gt_row gt_right" style="background-color: #E0FFFF; font-weight: bold;">—</td>
## <td headers="Nidsc" class="gt_row gt_right" style="background-color: #E0FFFF; font-weight: bold;">—</td>
## <td headers="Hiasc_perc" class="gt_row gt_right" style="background-color: #E0FFFF; font-weight: bold;">—</td>
## <td headers="Hidsc_perc" class="gt_row gt_right" style="background-color: #E0FFFF; font-weight: bold;">—</td></tr>
## </tbody>
##
## </table>
## </div>
#============================== GRAFICOS ==============================
colores_gas <- heat.colors(k)
# 1. HISTOGRAMA
hist(gas_production,
main = "Distribucion de Produccion Total de Gas",
breaks = seq(min(gas_production), max(gas_production), A),
xlab = "PRODUCCION DE GAS (M3)", ylab = "FRECUENCIA (ni)",
col = colores_gas, border = "white")
# 2. OJIVAS
par(mar = c(5, 4, 4, 6), xpd = TRUE)
x_asc <- c(min(TDFcu_gas_prod$liminf), TDFcu_gas_prod$limsup)
y_asc <- c(0, TDFcu_gas_prod$Niasc)
x_desc <- c(TDFcu_gas_prod$liminf, max(TDFcu_gas_prod$limsup))
y_desc <- c(TDFcu_gas_prod$Nidsc, 0)
plot(x_asc, y_asc, type = "o", pch = 19, lwd = 2, col = "darkgreen",
main = "Distribucion de Produccion Total de Gas",
xlab = "PRODUCCION DE GAS", ylab = "CANTIDAD DE SITIOS")
lines(x_desc, y_desc, type = "o", pch = 18, lwd = 2, col = "darkred")
legend("topright", inset = c(-0.05, 0), legend = c("Asc.", "Desc."),
col = c("darkgreen", "darkred"), lty = 1, pch = c(19, 18), bty = "n")
# 3. BOXPLOT
boxplot(gas_production, horizontal = TRUE, col = "darkcyan",
main = "Distribucion de Produccion Total de Gas",
xlab = "PRODUCCION DE GAS (M3)")
#============================== INDICADORES ESTADISTICOS ==============================
mean_g <- mean(gas_production)
sd_g <- sd(gas_production)
indicadores_g <- data.frame(
Indicador = c("Moda", "Mediana", "Media", "Desv. Estandar", "Varianza", "Coef. Variacion (%)", "Asimetria", "Curtosis"),
Valor = c(
round(get_mode(gas_production), 2),
round(median(gas_production), 2),
round(mean_g, 2),
round(sd_g, 2),
round(var(gas_production), 2),
round((sd_g / mean_g) * 100, 2),
round(skewness(gas_production), 2),
round(kurtosis(gas_production), 2)
)
)
print("TABLA DE INDICADORES: PRODUCCION TOTAL DE GAS")
## [1] "TABLA DE INDICADORES: PRODUCCION TOTAL DE GAS"
print(indicadores_g)
## Indicador Valor
## 1 Moda 1.725565e+06
## 2 Mediana 1.001777e+06
## 3 Media 1.740013e+06
## 4 Desv. Estandar 2.784581e+06
## 5 Varianza 7.753889e+12
## 6 Coef. Variacion (%) 1.600300e+02
## 7 Asimetria 7.960000e+00
## 8 Curtosis 1.154800e+02