Elevacion—continua.R

setwd("/cloud/project/datos")
datos <- read.csv("Petroleo_Ontaro.csv", header=T, dec=".", sep=";")


# VARIABLE ELEVACIÓN

elevacion <- datos$GROUND_ELEVATION
elevacion <- na.omit(elevacion)

min_elevacion <- min(elevacion)
max_elevacion <- max(elevacion)
R_elevacion <- max_elevacion - min_elevacion

k_elevacion <- 1 + (3.3 * log10(length(elevacion)))
k_elevacion <- floor(k_elevacion)
A_elevacion <- R_elevacion / k_elevacion

HistoSturges_elevacion <- hist(elevacion, 
                               main = "Gráfica No.48: Distribución de Elevación en el pozo", 
                               breaks = seq(min_elevacion, max_elevacion, A_elevacion),
                               xlab = "Elevación (m)", 
                               ylab = "Cantidad", 
                               col = "red")

HistoSturges_elevacion$counts

##  [1]    37   222    88  7695 14491  1481   420   104    16    27    14     4
## [13]     0     1     1

Histograma_elevacion <- hist(elevacion, 
                             main = "Gráfica No. : Distribución de Elevación en el pozo", 
                             xlab = "Elevación (m)", 
                             ylab = "Cantidad", 
                             col = "red")

limites_elevacion <- Histograma_elevacion$breaks
liminf_elevacion <- limites_elevacion[1:(k_elevacion)]
limsup_elevacion <- limites_elevacion[2:(k_elevacion+1)]
MC_elevacion <- (liminf_elevacion + limsup_elevacion) / 2



# TABLA ELEVACIÓN


ni_elevacion <- c()
for(i in 1:k_elevacion) {
  if(i == k_elevacion) {
    ni_elevacion[i] <- length(subset(elevacion, elevacion >= liminf_elevacion[i] & elevacion <= limsup_elevacion[i]))
  } else {
    ni_elevacion[i] <- length(subset(elevacion, elevacion >= liminf_elevacion[i] & elevacion < limsup_elevacion[i]))
  }
}

hi_elevacion <- (ni_elevacion / length(elevacion)) * 100

Niasc_elevacion <- cumsum(ni_elevacion)
Hiasc_elevacion <- cumsum(hi_elevacion)
Nidsc_elevacion <- rev(cumsum(rev(ni_elevacion)))
Hidsc_elevacion <- rev(cumsum(rev(hi_elevacion)))

Tabla_elevacion <- data.frame(liminf_elevacion, limsup_elevacion, MC_elevacion, ni_elevacion, 
                              round(hi_elevacion, 2), Niasc_elevacion, 
                              round(Hiasc_elevacion, 2), Nidsc_elevacion, 
                              round(Hidsc_elevacion, 2))

colnames(Tabla_elevacion) <- c("Xi-1", "Xi+1", "Xi (MC)", "ni", "hi", "Ni ↓", "Hi ↓", "Ni ↑", "Hi ↑")

Tabla_elevacion

##    Xi-1 Xi+1 Xi (MC)    ni    hi  Ni ↓   Hi ↓  Ni ↑   Hi ↑
## 1     0   50      25    49  0.20    49   0.20 24601 100.00
## 2    50  100      75   239  0.97   288   1.17 24552  99.80
## 3   100  150     125    70  0.28   358   1.46 24313  98.83
## 4   150  200     175 14479 58.86 14837  60.31 24243  98.54
## 5   200  250     225  8893 36.15 23730  96.46  9764  39.69
## 6   250  300     275   545  2.22 24275  98.67   871   3.54
## 7   300  350     325   235  0.96 24510  99.63   326   1.33
## 8   350  400     375    39  0.16 24549  99.79    91   0.37
## 9   400  450     425    27  0.11 24576  99.90    52   0.21
## 10  450  500     475    18  0.07 24594  99.97    25   0.10
## 11  500  550     525     5  0.02 24599  99.99     7   0.03
## 12  550  600     575     0  0.00 24599  99.99     2   0.01
## 13  600  650     625     1  0.00 24600 100.00     2   0.01
## 14  650  700     675     1  0.00 24601 100.00     1   0.00
## 15  700   NA      NA     0  0.00 24601 100.00     0   0.00

# OJIVAS CONTINUAS ELEVACIÓN  

plot(limsup_elevacion, Niasc_elevacion, col = "blue", type = "b", 
     main = "Gráfica No.1: Distribución de Elevación en el Pozo", 
     ylab = "Cantidad", xlab = "Elevación (m)")
lines(liminf_elevacion, Nidsc_elevacion, col = "red", type = "b")

# OJIVAS CONTINUAS PORCENTAJE ELEVACIÓN


plot(limsup_elevacion, Hiasc_elevacion, col = "blue", type = "b", 
     main = "Gráfica No.2: Distribución de Elevación en el Pozo", 
     ylab = "Porcentaje", xlab = "Elevación (m)")
lines(liminf_elevacion, Hidsc_elevacion, col = "red", type = "b")

# GDF CAJA Y BIGOTES


boxplot(elevacion, 
        horizontal = TRUE, 
        col = "darkgreen", 
        xlab = "Profundidad Alcanzada (m)", 
        main = "Gráfica No. 3: Distribución de Elevación en el Pozo")