setwd("/cloud/project/PROYECTO/DATOS")
datos<-read.csv("Dataset Andres.csv",
header= TRUE, sep= ",", dec=".")
### ED DE VARIABLES CUANTITATIVAS CONTINUAS
## VARIABLE PROFUNDIDAD
prof_alcanzada<-datos$TOTAL_DEPTH
prof_alcanzada<-na.omit(prof_alcanzada)
# GDF
min(prof_alcanzada)
## [1] 1.37
max(prof_alcanzada)
## [1] 5000
R_profundidad<-max(prof_alcanzada)-min(prof_alcanzada)
k_profundidad<-1+(3.3*log10(length(prof_alcanzada)))
k_profundidad<-floor(k_profundidad)
A_profundidad<-R_profundidad/k_profundidad
HistoSturges_profundidad<-hist(prof_alcanzada,
main="Gráfica No.1:
Distribución de Profundidad Alcanzada en el pozo",
breaks=seq(min(prof_alcanzada),max(prof_alcanzada),A_profundidad),
xlab="Profundidad Alcanzada (m)",
ylab="Cantidad",
col="red")
HistoSturges_profundidad$counts
## [1] 11221 6748 1846 686 57 28 13 5 3 5 1 0
## [13] 0 0 1
Histograma_profundidad<-hist(prof_alcanzada,
main="Gráfica No.2:
Distribución de Profundidad Alcanzada en el pozo",
xlab="Profundidad Alcanzada (m)",
ylab="Cantidad",
col="red")
limites_profundidad<-Histograma_profundidad$breaks
liminf_profundidad<-limites_profundidad[1:15]
limsup_profundidad<-limites_profundidad[2:16]
MC_profundidad<-(liminf_profundidad+limsup_profundidad)/2
# TABLA PROFUNDIDAD
ni_profundidad<-c()
for(i in 1:k_profundidad) {
if(1==15)
ni_profundidad[i]<- length(subset(prof_alcanzada, prof_alcanzada>=liminf_profundidad[i] & prof_alcanzada<=limsup_profundidad[i]))
else
ni_profundidad[i]<- length(subset(prof_alcanzada, prof_alcanzada>=liminf_profundidad[i]& prof_alcanzada<limsup_profundidad[i]))
}
hi_profundidad<-(ni_profundidad/length(prof_alcanzada))*100
Niasc_profundidad<-cumsum(ni_profundidad)
Hiasc_profundidad<-cumsum(hi_profundidad)
Nidsc_profundidad<-rev(cumsum(rev(ni_profundidad)))
Hidsc_profundidad<-rev(cumsum(rev(hi_profundidad)))
Tabla_profundidad<-data.frame(liminf_profundidad,
limsup_profundidad,
MC_profundidad,
ni_profundidad,
round(hi_profundidad,2),
Niasc_profundidad,
round(Hiasc_profundidad,2),
Nidsc_profundidad,
round(Hidsc_profundidad,2))
Tabla_profundidad
## liminf_profundidad limsup_profundidad MC_profundidad ni_profundidad
## 1 0 500 250 15331
## 2 500 1000 750 4480
## 3 1000 1500 1250 732
## 4 1500 2000 1750 43
## 5 2000 2500 2250 17
## 6 2500 3000 2750 4
## 7 3000 3500 3250 6
## 8 3500 4000 3750 0
## 9 4000 4500 4250 0
## 10 4500 5000 4750 0
## 11 5000 NA NA 0
## 12 NA NA NA 0
## 13 NA NA NA 0
## 14 NA NA NA 0
## 15 NA NA NA 0
## round.hi_profundidad..2. Niasc_profundidad round.Hiasc_profundidad..2.
## 1 74.37 15331 74.37
## 2 21.73 19811 96.10
## 3 3.55 20543 99.66
## 4 0.21 20586 99.86
## 5 0.08 20603 99.95
## 6 0.02 20607 99.97
## 7 0.03 20613 100.00
## 8 0.00 20613 100.00
## 9 0.00 20613 100.00
## 10 0.00 20613 100.00
## 11 0.00 20613 100.00
## 12 0.00 20613 100.00
## 13 0.00 20613 100.00
## 14 0.00 20613 100.00
## 15 0.00 20613 100.00
## Nidsc_profundidad round.Hidsc_profundidad..2.
## 1 20613 100.00
## 2 5282 25.62
## 3 802 3.89
## 4 70 0.34
## 5 27 0.13
## 6 10 0.05
## 7 6 0.03
## 8 0 0.00
## 9 0 0.00
## 10 0 0.00
## 11 0 0.00
## 12 0 0.00
## 13 0 0.00
## 14 0 0.00
## 15 0 0.00
colnames(Tabla_profundidad)<-c("Xi-1","Xi+1","Xi (MC)","ni","hi","Ni ↓","Hi ↓","Ni ↑","Hi ↑")
# OJIVAS CONTINUAS PROFUNDIDAD ALCANZADA
plot(liminf_profundidad, Nidsc_profundidad, col="blue", type="b",
main="Gráfica No.3: Distribución de Profundidad Alcanzada en el Pozo",
ylab = "Cantidad",
xlab = "Profundidad Alcanzada (m)")
lines(limsup_profundidad, Niasc_profundidad, col="red", type="b")
# OJIVAS CONTINUAS PORCENTAJE PROFUNDIDAD ALCANZADA
plot(liminf_profundidad, Hidsc_profundidad,col="red", type="b",
main="Gráfica No.4: Distribución de Profundidad Alcanzada en el Pozo",
ylab = "Porcentaje",
xlab = "Profundidad Alcanzada (m)")
lines(limsup_profundidad, Hiasc_profundidad, col="blue", type="b")
# GDF CAJA Y BIGOTES
boxplot(prof_alcanzada,
horizontal=T,
col="darkgreen",
xlab= "Profundidad Alcanzada (m)",
main="Gráfica No. 5:
Distribución de Profundidad Alcanzada en el Pozo")
# CONCLUSIÓN
#El comportamiento de la profundidad alcanzada en el pozo fluctúan entre [1.37, 5000] y los valores se encuentran entorno a 301.75, con una desviación estándar de 273.6755 siendo un conjunto de valores muy heterogéneos donde los valores se agrupan en la parte izquierda de la variable y tiene un sesgo positivo con existencia de valores atípicos a partir de 1000.
## VARIABLE ELEVACIÓN
elevacion<-datos$GROUND_ELEVATION
elevacion<-na.omit(elevacion)
# GDF
min(elevacion)
## [1] 0.3
max(elevacion)
## [1] 694
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.1:
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] 33 185 75 6779 12218 1179 381 102 13 26 7 1
## [13] 0 1 1
Histograma_elevacion<-hist(elevacion, main="Gráfica No.2:
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:15]
limsup_elevacion<-limites_elevacion[2:16]
MC_elevacion<-(liminf_elevacion+limsup_elevacion)/2
# TABLA ELEVACIÓN
ni_elevacion<-c()
for(i in 1:k_elevacion) {
if(1==15)
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))
Tabla_elevacion
## liminf_elevacion limsup_elevacion MC_elevacion ni_elevacion
## 1 0 50 25 36
## 2 50 100 75 209
## 3 100 150 125 59
## 4 150 200 175 12764
## 5 200 250 225 7148
## 6 250 300 275 488
## 7 300 350 325 221
## 8 350 400 375 36
## 9 400 450 425 26
## 10 450 500 475 10
## 11 500 550 525 2
## 12 550 600 575 0
## 13 600 650 625 1
## 14 650 700 675 1
## 15 700 NA NA 0
## round.hi_elevacion..2. Niasc_elevacion round.Hiasc_elevacion..2.
## 1 0.17 36 0.17
## 2 1.00 245 1.17
## 3 0.28 304 1.45
## 4 60.78 13068 62.23
## 5 34.04 20216 96.26
## 6 2.32 20704 98.59
## 7 1.05 20925 99.64
## 8 0.17 20961 99.81
## 9 0.12 20987 99.93
## 10 0.05 20997 99.98
## 11 0.01 20999 99.99
## 12 0.00 20999 99.99
## 13 0.00 21000 100.00
## 14 0.00 21001 100.00
## 15 0.00 21001 100.00
## Nidsc_elevacion round.Hidsc_elevacion..2.
## 1 21001 100.00
## 2 20965 99.83
## 3 20756 98.83
## 4 20697 98.55
## 5 7933 37.77
## 6 785 3.74
## 7 297 1.41
## 8 76 0.36
## 9 40 0.19
## 10 14 0.07
## 11 4 0.02
## 12 2 0.01
## 13 2 0.01
## 14 1 0.00
## 15 0 0.00
colnames(Tabla_elevacion)<-c("Xi-1","Xi+1","Xi (MC)","ni","hi","Ni ↓","Hi ↓","Ni ↑","Hi ↑")
# OJIVAS CONTINUAS ELEVACIÓN
plot(limsup_elevacion, Niasc_elevacion, col="blue", type="b",
main="Gráfica No.3: 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.4: 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=T,
col="darkgreen",
xlab= "Profundidad Alcanzada (m)",
main="Gráfica No.5:
Distribución de Elevación en el Pozo")
# CONCLUSIÓN
#El comportamiento de la elevación del terreno fluctúan entre [0.3, 694] y los valores se encuentran entorno a 193.85, con una desviación estándar de 31.21711 siendo un conjunto de valores homogéneos donde los valores se agrupan en la parte izquierda de la variable y tiene un sesgo positivo con existencia de valores atípicos a partir de 250.
## VARIABLE VERTICALIDAD
prof_vertical<-datos$TRUE_VERTICAL_DEPTH
prof_vertical<-na.omit(prof_vertical)
# GDF
min(prof_vertical)
## [1] 1.37
max(prof_vertical)
## [1] 1828.8
R_vertical<-max(prof_vertical)-min(prof_vertical)
k_vertical<-1+(3.3*log10(length(prof_vertical)))
k_vertical<-floor(k_vertical)
A_vertical<-R_vertical/k_vertical
HistoSturges_vertical<-hist(prof_vertical, main="Gráfica No.1:
Distribución de Profundidad vertical en el pozo",
breaks = seq(min(prof_vertical),max(prof_vertical),A_vertical),
xlab="Profundidad vertical (m)",
ylab="Cantidad",
col="red")
HistoSturges_vertical$counts
## [1] 2017 5382 4533 3240 1957 1511 808 424 267 304 89 19 4 0 5
Histograma_vertical<-hist(prof_vertical,
main="Gráfica No.2:
Distribución de Profundidad vertical en el pozo",
xlab="Profundidad vertical (m)",
ylab="Cantidad",
col="red")
limites_vertical<-Histograma_vertical$breaks
liminf_vertical<-limites_vertical[1:15]
limsup_vertical<-limites_vertical[2:16]
MC_vertical<-(liminf_vertical+limsup_vertical)/2
# TABLA PROFUNDIDAD
ni_vertical<-c()
for(i in 1:k_vertical) {
if(1==15)
ni_vertical[i]<- length(subset(prof_vertical, prof_vertical>=liminf_vertical[i] & prof_vertical<=limsup_vertical[i]))
else
ni_vertical[i]<- length(subset(prof_vertical, prof_vertical>=liminf_vertical[i]& prof_vertical<limsup_vertical[i]))
}
hi_vertical<-(ni_vertical/length(prof_vertical))*100
Niasc_vertical<-cumsum(ni_vertical)
Hiasc_vertical<-cumsum(hi_vertical)
Nidsc_vertical<-rev(cumsum(rev(ni_vertical)))
Hidsc_vertical<-rev(cumsum(rev(hi_vertical)))
Tabla_vertical<-data.frame(liminf_vertical,
limsup_vertical,
MC_vertical,
ni_vertical,
round(hi_vertical,2),
Niasc_vertical,
round(Hiasc_vertical,2),
Nidsc_vertical,
round(Hidsc_vertical,2))
Tabla_vertical
## liminf_vertical limsup_vertical MC_vertical ni_vertical
## 1 0 100 50 784
## 2 100 200 150 4681
## 3 200 300 250 4772
## 4 300 400 350 2716
## 5 400 500 450 2357
## 6 500 600 550 1594
## 7 600 700 650 1414
## 8 700 800 750 855
## 9 800 900 850 521
## 10 900 1000 950 224
## 11 1000 1100 1050 225
## 12 1100 1200 1150 242
## 13 1200 1300 1250 133
## 14 1300 1400 1350 26
## 15 1400 1500 1450 11
## round.hi_vertical..2. Niasc_vertical round.Hiasc_vertical..2. Nidsc_vertical
## 1 3.81 784 3.81 20555
## 2 22.77 5465 26.58 19771
## 3 23.21 10237 49.79 15090
## 4 13.21 12953 63.00 10318
## 5 11.46 15310 74.46 7602
## 6 7.75 16904 82.22 5245
## 7 6.88 18318 89.10 3651
## 8 4.16 19173 93.25 2237
## 9 2.53 19694 95.79 1382
## 10 1.09 19918 96.88 861
## 11 1.09 20143 97.97 637
## 12 1.18 20385 99.15 412
## 13 0.65 20518 99.80 170
## 14 0.13 20544 99.92 37
## 15 0.05 20555 99.98 11
## round.Hidsc_vertical..2.
## 1 99.98
## 2 96.16
## 3 73.39
## 4 50.18
## 5 36.97
## 6 25.51
## 7 17.76
## 8 10.88
## 9 6.72
## 10 4.19
## 11 3.10
## 12 2.00
## 13 0.83
## 14 0.18
## 15 0.05
colnames(Tabla_vertical)<-c("Xi-1","Xi+1","Xi (MC)","ni","hi","Ni ↓","Hi ↓","Ni ↑","Hi ↑")
# OJIVAS CONTINUAS PROFUNDIDAD VERTICAL
plot(limsup_vertical, Niasc_vertical, col="blue", type="b",
main="Gráfica No.3: Distribución de Profundidad Vertical en el Pozo",
ylab = "Cantidad",
xlab = "Profundidad Vertical (m)")
lines(liminf_vertical, Nidsc_vertical, col="red", type="b")
# OJIVAS CONTINUAS PORCENTAJE PROFUNDIDAD ALCANZADA
plot(limsup_vertical, Hiasc_vertical, col="blue", type="b",
main="Gráfica No.4: Distribución de Profundidad Vertical en el Pozo",
ylab = "Porcentaje",
xlab = "Profundidad Vertical (m)")
lines(liminf_vertical, Hidsc_vertical, col="red", type="b")
# GDF CAJA Y BIGOTES
boxplot(prof_vertical,
horizontal=T,
col="darkgreen",
xlab= "Profundidad Vertical (m)",
main="Gráfica No.5:
Distribución de Profundidad Vertical en el Pozo")
# CONCLUSIÓN
#El comportamiento de la verticalidad fluctúan entre [1.37, 1828.8] y los valores se encuentran entorno a 301.75, con una desviación estándar de 249.0132 siendo un conjunto de valores heterogéneos donde los valores se agrupan en la parte izquierda de la variable y tiene un sesgo positivo con existencia de valores atípicos a partir de 1000.