Exploracción y Producción de los Pozos Petrolíferos en Brasil
Anahi Macias
2026-01-25
Modelo de Probabilidad de los Bloques de Concesión
1 Cargar datos
Se importan los datos originales y se verifica su estructura para asegurar una lectura correcta.
setwd("C:/Users/majke/Downloads/Proyecto Estadistica/RMARKDOWN")
Datos <- read.csv("Pozos brasil 2.csv", header = TRUE, sep = ";", dec = ",", fileEncoding = "Latin1")
str(Datos)## 'data.frame': 29575 obs. of 59 variables:
## $ POCO : chr "7-RO-123HP-RJS" "1-BP-7-RJS" "7-ARGO-4H-ESS" "7-ARGO-5H-ESS" ...
## $ CADASTRO : num 7.43e+10 7.43e+10 3.43e+10 3.43e+10 3.43e+10 ...
## $ OPERADOR : chr "Petrobras" "BP Energy" "Shell Brasil" "Shell Brasil" ...
## $ POCO_OPERADOR : chr "7RO123HPRJS" "ANU" "7ARGO4HESS" "7ARGO5HESS" ...
## $ ESTADO : chr "RJ" "RJ" "ES" "ES" ...
## $ BACIA : chr "Campos" "Campos" "Campos" "Campos" ...
## $ BLOCO : chr "" "C-M-473" "" "" ...
## $ SIG_CAMPO : chr "RO " "" "ARGO " "ARGO " ...
## $ CAMPO : chr "RONCADOR" "" "ARGONAUTA" "ARGONAUTA" ...
## $ TERRA_MAR : chr "M" "M" "M" "M" ...
## $ POCO_POS_ANP : chr "S" "S" "S" "S" ...
## $ TIPO : chr "Explotatório" "Exploratório" "Explotatório" "Explotatório" ...
## $ CATEGORIA : chr "Desenvolvimento" "Pioneiro" "Desenvolvimento" "Desenvolvimento" ...
## $ RECLASSIFICACAO : chr "PRODUTOR COMERCIAL DE PETRÓLEO" "PORTADOR DE PETRÓLEO" "ABANDONADO POR OUTRAS RAZÕES" "INDEFINIDO" ...
## $ SITUACAO : chr "PRODUTOR" "ABANDONADO POR LOGÍSTICA EXPLORATÓRIA" "ABANDONADO DEFINITIVAMENTE" "PRODUTOR" ...
## $ INICIO : chr "03/03/2012" "14/03/2012" "12/04/2012" "14/04/2012" ...
## $ TERMINO : chr "10/03/2013" "20/09/2012" "19/09/2012" "02/05/2013" ...
## $ CONCLUSAO : chr "10/03/2013" "19/10/2012" "28/09/2012" "13/05/2013" ...
## $ TITULARIDADE : chr "Público" "Público" "Público" "Público" ...
## $ LATITUDE_BASE_4C : chr "-21:57:45,630" "-23:06:09,622" "-21:08:06,801" "-21:08:12,535" ...
## $ LONGITUDE_BASE_4C : chr "-39:44:01,130" "-40:00:41,295" "-39:46:53,226" "-39:46:49,296" ...
## $ LATITUDE_BASE_DD : num -22 -23.1 -21.1 -21.1 -21.1 ...
## $ LONGITUDE_BASE_DD : num -39.7 -40 -39.8 -39.8 -39.8 ...
## $ DATUM_HORIZONTAL : chr "SIRGAS2000" "SIRGAS2000" "SIRGAS2000" "SIRGAS2000" ...
## $ TIPO_DE_COORDENADA_DE_BASE: chr "Definitiva" "Definitiva" "Definitiva" "Definitiva" ...
## $ DIRECAO : chr "Horizontal" "Vertical" "Horizontal" "Horizontal" ...
## $ PROFUNDIDADE_VERTICAL_M : num -3145 6900 2937 2934 2953 ...
## $ PROFUNDIDADE_SONDADOR_M : num 4050 6925 3809 4575 4570 ...
## $ PROFUNDIDADE_MEDIDA_M : num 4050 6925 3809 4575 4570 ...
## $ REFERENCIA_DE_PROFUNDIDADE: chr "MR" "MR" "MR" "MR" ...
## $ MESA_ROTATIVA : num 24 25 24.2 24.2 24.2 ...
## $ COTA_ALTIMETRICA_M : num 0 0 0 0 0 ...
## $ LAMINA_D_AGUA_M : num 1827 2730 1706 1705 1654 ...
## $ DATUM_VERTICAL : chr "NM" "NM" "NM" "NM" ...
## $ UNIDADE_ESTRATIGRAFICA : chr "" "" "" "" ...
## $ GEOLOGIA_GRUPO_FINAL : chr "Campos" "Lagoa Feia" "Campos" "Campos" ...
## $ GEOLOGIA_FORMACAO_FINAL : chr "Carapebus" "Macabu" "Ubatuba" "Ubatuba" ...
## $ GEOLOGIA_MEMBRO_FINAL : chr "" "" "" "" ...
## $ CDPE : chr "Existe" "Existe" "Existe" "" ...
## $ AGP : chr "" "" "" "" ...
## $ PC : chr "" "Existe" "" "" ...
## $ PAG : chr "" "" "" "" ...
## $ PERFIS_CONVENCIONAIS : chr "" "Existe" "" "" ...
## $ DURANTE_PERFURACAO : chr "Existe" "Existe" "Existe" "Existe" ...
## $ PERFIS_DIGITAIS : chr "" "" "" "" ...
## $ PERFIS_PROCESSADOS : chr "" "" "" "" ...
## $ PERFIS_ESPECIAIS : chr "" "" "" "" ...
## $ AMOSTRA_LATERAL : chr "" "Existe" "" "" ...
## $ SISMICA : chr "" "Existe" "" "" ...
## $ TABELA_TEMPO_PROFUNDIDADE : chr "" "" "" "" ...
## $ DADOS_DIRECIONAIS : chr "Existe" "Existe" "Existe" "Existe" ...
## $ TESTE_A_CABO : chr "Existe" "Existe" "" "" ...
## $ TESTE_DE_FORMACAO : chr "" "" "" "" ...
## $ CANHONEIO : chr "" "" "" "" ...
## $ TESTEMUNHO : chr "" "" "" "" ...
## $ GEOQUIMICA : chr "" "" "" "" ...
## $ SIG_SONDA : chr "SS-49" "DS4" "NB2" "NB2" ...
## $ NOM_SONDA : chr "SEDCO 707" "Deep Ocean Clarion" "Bully 2" "Bully 2" ...
## $ DHA_ATUALIZACAO : chr "28/1/18 5:00" "28/1/18 5:00" "28/1/18 5:00" "28/1/18 5:00" ...2 Tabla de distribución de frecuencia
Se obtiene la frecuencia absoluta de los bloques de concesión según su denominación.
Bloque <- Datos$BLOCO
TDFBloque <- as.data.frame(table(Bloque, useNA = "ifany"))
TDFBloque## Bloque Freq
## 1 3253
## 2 A 28
## 3 A032 1
## 4 AB 172
## 5 ABL 95
## 6 AG 501
## 7 ALG 2
## 8 AM-T-62 1
## 9 AM-T-84 1
## 10 AM-T-85 2
## 11 AN 34
## 12 ANQ 6
## 13 APR 16
## 14 AR 348
## 15 ARB 13
## 16 ARG 798
## 17 ARI 7
## 18 ARL 3
## 19 ARR 5
## 20 ART 10
## 21 AT 33
## 22 ATS 21
## 23 B097 11
## 24 B104 1
## 25 BA 464
## 26 BA-1 3
## 27 BA-2 2
## 28 BA-3 8
## 29 BA-4 1
## 30 BAL 79
## 31 BALM-200 1
## 32 BAR 3
## 33 BB 18
## 34 BC-10 41
## 35 BC-100 3
## 36 BC-14 1
## 37 BC-15 6
## 38 BC-2 8
## 39 BC-20 10
## 40 BC-200 13
## 41 BC-200A 1
## 42 BC-3 1
## 43 BC-30 12
## 44 BC-4 1
## 45 BC-40 2
## 46 BC-400 12
## 47 BC-50 9
## 48 BC-500 3
## 49 BC-60 123
## 50 BC-600 8
## 51 BC-7 1
## 52 BC-8 1
## 53 BC-9 2
## 54 BCAM-1 2
## 55 BCAM-100 1
## 56 BCAM-2 2
## 57 BCAM-40 16
## 58 BCE-2 2
## 59 BCE-5 4
## 60 BCE-6 2
## 61 BCE-8 3
## 62 BCE-9 1
## 63 BCED 18
## 64 BCUM-1 4
## 65 BCUM-100 1
## 66 BCUM-20 2
## 67 BD 36
## 68 BE 29
## 69 BE-49 3
## 70 BEN 34
## 71 BES-100 83
## 72 BES-2 1
## 73 BES-200 4
## 74 BES-3 3
## 75 BFRD 7
## 76 BFZ-1 1
## 77 BFZ-2 3
## 78 BG 22
## 79 BI 68
## 80 BJ 45
## 81 BJ-30 1
## 82 BJZ 5
## 83 BL 3
## 84 BLV_R13AM 6
## 85 BM-BAR-1 2
## 86 BM-BAR-3 2
## 87 BM-C-10 1
## 88 BM-C-14 1
## 89 BM-C-25 3
## 90 BM-C-3 1
## 91 BM-C-4 2
## 92 BM-C-5 2
## 93 BM-C-6 1
## 94 BM-C-7 17
## 95 BM-C-8 52
## 96 BM-CAL-4 7
## 97 BM-CAL-5 2
## 98 BM-CAL-6 1
## 99 BM-CE-1 2
## 100 BM-CE-2 1
## 101 BM-ES-1 1
## 102 BM-ES-5 10
## 103 BM-ES-9 1
## 104 BM-FZA-1 1
## 105 BM-J-1 1
## 106 BM-J-2 2
## 107 BM-J-3 2
## 108 BM-PAMA-3 1
## 109 BM-POT-11 4
## 110 BM-POT-13 3
## 111 BM-S-10 2
## 112 BM-S-11 45
## 113 BM-S-12 4
## 114 BM-S-13 4
## 115 BM-S-14 3
## 116 BM-S-17 3
## 117 BM-S-21 1
## 118 BM-S-22 4
## 119 BM-S-24 7
## 120 BM-S-29 3
## 121 BM-S-3 9
## 122 BM-S-31 2
## 123 BM-S-4 9
## 124 BM-S-7 8
## 125 BM-S-8 12
## 126 BM-S-9 27
## 127 BM-SEAL-4 4
## 128 BMUC 1
## 129 BO 56
## 130 BP-1 1
## 131 BPAR-10 9
## 132 BPOT-10 10
## 133 BPOT-100 1
## 134 BPOT-100A 4
## 135 BPOT-2 2
## 136 BPOT-3 1
## 137 BPOT-4 21
## 138 BPOT-5 3
## 139 BPOT-6 9
## 140 BR 182
## 141 Brasil 2728
## 142 BREC-10 1
## 143 BREC-11 2
## 144 BREC-13 2
## 145 BREC-3 2
## 146 BREC-5 1
## 147 BREC-7 1
## 148 BREC-8 1
## 149 BREC-9 2
## 150 BREG 10
## 151 BRG 45
## 152 BS-1 11
## 153 BS-2 2
## 154 BS-3 5
## 155 BS-4 8
## 156 BS-400 36
## 157 BS-500 58
## 158 BSOL-4 4
## 159 BSOL-5 8
## 160 BSU 16
## 161 BT-ES-12 36
## 162 BT-ES-14 5
## 163 BT-ES-15 8
## 164 BT-POT-10 28
## 165 BT-POT-3 3
## 166 BT-POT-4 29
## 167 BT-POT-5 4
## 168 BT-POT-8 28
## 169 BT-POT-9 7
## 170 BT-REC-1 1
## 171 BT-REC-10 39
## 172 BT-REC-2 1
## 173 BT-REC-3 5
## 174 BT-REC-4 4
## 175 BT-REC-7 5
## 176 BT-REC-8 2
## 177 BT-SEAL-2 27
## 178 BT-SEAL-3 1
## 179 BT-SOL-1 6
## 180 BTUS-3 1
## 181 BTUS-5 1
## 182 BVS 69
## 183 C 252
## 184 C-M-101 7
## 185 C-M-146 1
## 186 C-M-151 1
## 187 C-M-202 3
## 188 C-M-231 3
## 189 C-M-299 1
## 190 C-M-333 2
## 191 C-M-401 11
## 192 C-M-403 1
## 193 C-M-466 19
## 194 C-M-471 4
## 195 C-M-473 5
## 196 C-M-498 1
## 197 C-M-499 12
## 198 C-M-529 4
## 199 C-M-530 1
## 200 C-M-535 2
## 201 C-M-539 9
## 202 C-M-560 8
## 203 C-M-591 1
## 204 C-M-592 40
## 205 C-M-61 5
## 206 C-M-620 5
## 207 C-M-621 7
## 208 C-M-622 3
## 209 C-M-95 3
## 210 C-M-96 1
## 211 C066 3
## 212 CA 25
## 213 CAB 7
## 214 CAC 27
## 215 CAL-M-314 1
## 216 CAL-M-60 2
## 217 CAM 1747
## 218 CAN 7
## 219 CAR 13
## 220 CAU 5
## 221 CB 53
## 222 CCN 11
## 223 CD 8
## 224 CER 28
## 225 CG 45
## 226 CGL 3
## 227 CH 63
## 228 CHA 1
## 229 CL 163
## 230 CM 120
## 231 CO 24
## 232 COG 4
## 233 CON 24
## 234 CP 1886
## 235 CR 35
## 236 CRJ 2
## 237 CRL 6
## 238 CRP 89
## 239 CRT 57
## 240 CRV 5
## 241 CS 229
## 242 CSF 3
## 243 CSM 38
## 244 CUP 1
## 245 CX 95
## 246 DA 2
## 247 DJ 239
## 248 DJM 643
## 249 DO 41
## 250 EM 4
## 251 EN 56
## 252 ENO 37
## 253 EO 2
## 254 EP 34
## 255 ES-M-413 3
## 256 ES-M-414 5
## 257 ES-M-416 1
## 258 ES-M-418 1
## 259 ES-M-466 3
## 260 ES-M-468 4
## 261 ES-M-472 2
## 262 ES-M-523 4
## 263 ES-M-525 9
## 264 ES-M-527 4
## 265 ES-M-529 1
## 266 ES-M-588 4
## 267 ES-M-590 2
## 268 ES-M-592 2
## 269 ES-M-594 3
## 270 ES-M-661 1
## 271 ES-M-663 1
## 272 ES-M-665 1
## 273 ES-M-737 2
## 274 ES-T-108 1
## 275 ES-T-125 1
## 276 ES-T-174 1
## 277 ES-T-188 1
## 278 ES-T-202 1
## 279 ES-T-215 1
## 280 ES-T-227 1
## 281 ES-T-291 1
## 282 ES-T-363 1
## 283 ES-T-364 15
## 284 ES-T-372 1
## 285 ES-T-373 9
## 286 ES-T-381 3
## 287 ES-T-382 8
## 288 ES-T-383 10
## 289 ES-T-390 3
## 290 ES-T-391 8
## 291 ES-T-392 5
## 292 ES-T-400 5
## 293 ES-T-401 1
## 294 ES-T-409 3
## 295 ES-T-410 1
## 296 ES-T-418 2
## 297 ES-T-419 2
## 298 ES-T-454 3
## 299 ES-T-466 5
## 300 ES-T-47 1
## 301 ES-T-486 11
## 302 ES-T-495 6
## 303 ES-T-496 5
## 304 ES-T-505 6
## 305 ES-T-506 3
## 306 ES-T-516 3
## 307 ES-T-527 2
## 308 ES-T-67 2
## 309 ES-T-68 1
## 310 ES-T-88 2
## 311 ESP 61
## 312 ET 1374
## 313 FA 16
## 314 FAL 129
## 315 FAP 1
## 316 FAV 146
## 317 FBE 73
## 318 FBL 27
## 319 FBM 300
## 320 FC 57
## 321 FCN 51
## 322 FGA 3
## 323 FGT 4
## 324 FI 117
## 325 FJ 3
## 326 FLORIM 2
## 327 FMQ 42
## 328 FMT 3
## 329 FO 13
## 330 FP 553
## 331 FPB 9
## 332 FQ 51
## 333 FR 54
## 334 FRANCO 12
## 335 FRB 8
## 336 FS 7
## 337 FSE 30
## 338 FSJ 54
## 339 FSL 57
## 340 FSP 1
## 341 FSR 62
## 342 FU 175
## 343 FVB 2
## 344 FZA-M-252 1
## 345 FZB 846
## 346 GA 92
## 347 GMR 93
## 348 GO 20
## 349 GP 43
## 350 GPN 10
## 351 GU 7
## 352 GUARA_SUL 2
## 353 I 77
## 354 IARA_ENT 5
## 355 IB 2
## 356 IC 11
## 357 ICA 12
## 358 IP 21
## 359 IR 15
## 360 JA 59
## 361 JD 17
## 362 JI 3
## 363 JR 8
## 364 JZ 6
## 365 LAR 14
## 366 LB 5
## 367 LD_R13AM 1
## 368 LE 2
## 369 LI 43
## 370 LIBRA 7
## 371 LM 19
## 372 LOR 80
## 373 LP 88
## 374 LPB 27
## 375 LPC 3
## 376 LPN 30
## 377 LPS 5
## 378 LPX 22
## 379 LS 69
## 380 LUC 86
## 381 LV 63
## 382 MA 102
## 383 MAG 79
## 384 MAN 5
## 385 MAO 2
## 386 MB 5
## 387 MDU 6
## 388 ME 10
## 389 MG 706
## 390 MGL 3
## 391 MGN 19
## 392 MJ 104
## 393 ML 48
## 394 MLH 18
## 395 MLL 81
## 396 MLS 219
## 397 MLZ 9
## 398 MO 127
## 399 MOR 24
## 400 MP 46
## 401 MRB 9
## 402 MRL 241
## 403 MUI 21
## 404 NEJ 3
## 405 NEN 5
## 406 NFC 56
## 407 NMR 12
## 408 NO 11
## 409 OC 4
## 410 P 7
## 411 PA 20
## 412 PARU 5
## 413 PC 41
## 414 PDR 5
## 415 PE 18
## 416 PER 10
## 417 PG 39
## 418 PIR 271
## 419 PJ 68
## 420 PJN 6
## 421 PM 80
## 422 PML 35
## 423 PN-T-102 4
## 424 PN-T-48 2
## 425 PN-T-49 10
## 426 PN-T-50 2
## 427 PN-T-66 1
## 428 PN-T-67 15
## 429 PN-T-68 65
## 430 PN-T-84 2
## 431 PN-T-85 4
## 432 PN-T-86 1
## 433 POT-M-665 1
## 434 POT-M-760 1
## 435 POT-M-855 2
## 436 POT-T-208 1
## 437 POT-T-210 1
## 438 POT-T-225 1
## 439 POT-T-240 1
## 440 POT-T-257 2
## 441 POT-T-298 1
## 442 POT-T-302 12
## 443 POT-T-321 1
## 444 POT-T-352 2
## 445 POT-T-353 1
## 446 POT-T-354 2
## 447 POT-T-355 2
## 448 POT-T-391 4
## 449 POT-T-392 1
## 450 POT-T-393 3
## 451 POT-T-394 2
## 452 POT-T-402 1
## 453 POT-T-403 1
## 454 POT-T-404 1
## 455 POT-T-432 5
## 456 POT-T-433 1
## 457 POT-T-434 1
## 458 POT-T-435 1
## 459 POT-T-440 2
## 460 POT-T-442 1
## 461 POT-T-445 3
## 462 POT-T-447 1
## 463 POT-T-476 1
## 464 POT-T-477 1
## 465 POT-T-478 4
## 466 POT-T-479 2
## 467 POT-T-480 3
## 468 POT-T-484 2
## 469 POT-T-485 1
## 470 POT-T-488 2
## 471 POT-T-512 4
## 472 POT-T-513 6
## 473 POT-T-514 8
## 474 POT-T-515 1
## 475 POT-T-520 4
## 476 POT-T-521 1
## 477 POT-T-523 2
## 478 POT-T-524 2
## 479 POT-T-525 1
## 480 POT-T-527 1
## 481 POT-T-528 1
## 482 POT-T-531 1
## 483 POT-T-556 2
## 484 POT-T-557 3
## 485 POT-T-558 6
## 486 POT-T-559 3
## 487 POT-T-560 2
## 488 POT-T-562 1
## 489 POT-T-563 2
## 490 POT-T-564 2
## 491 POT-T-565 1
## 492 POT-T-569 2
## 493 POT-T-573 1
## 494 POT-T-575 3
## 495 POT-T-600 1
## 496 POT-T-601 4
## 497 POT-T-602 2
## 498 POT-T-603 2
## 499 POT-T-605 1
## 500 POT-T-606 1
## 501 POT-T-607 1
## 502 POT-T-608 1
## 503 POT-T-609 3
## 504 POT-T-610 2
## 505 POT-T-612 2
## 506 POT-T-613 2
## 507 POT-T-614 2
## 508 POT-T-618 1
## 509 POT-T-650 1
## 510 POT-T-651 1
## 511 POT-T-655 1
## 512 POT-T-662 1
## 513 POT-T-696 1
## 514 POT-T-697 1
## 515 POT-T-699 2
## 516 POT-T-700 4
## 517 POT-T-704 2
## 518 POT-T-706 3
## 519 POT-T-743 3
## 520 POT-T-744 4
## 521 POT-T-745 1
## 522 POT-T-748 5
## 523 POT-T-749 3
## 524 POT-T-790 7
## 525 POT-T-791 2
## 526 POT-T-792 1
## 527 POT-T-794 2
## 528 POT-T-837 1
## 529 POT-T-881 1
## 530 POT-T-882 1
## 531 PRC-T-122 1
## 532 PSA 4
## 533 PU 25
## 534 PV 8
## 535 PV_R13AM 10
## 536 PX 12
## 537 QE 25
## 538 R035 8
## 539 R071 14
## 540 R093 10
## 541 R128 1
## 542 R134 1
## 543 RAL 1
## 544 RBS 6
## 545 RBU 130
## 546 RC 1
## 547 RCB 84
## 548 RD 7
## 549 RE 38
## 550 REC-T-102 1
## 551 REC-T-103 1
## 552 REC-T-106 1
## 553 REC-T-113 1
## 554 REC-T-126 1
## 555 REC-T-129 3
## 556 REC-T-130 3
## 557 REC-T-138 1
## 558 REC-T-139 1
## 559 REC-T-142 2
## 560 REC-T-153 2
## 561 REC-T-155 9
## 562 REC-T-163 1
## 563 REC-T-165 1
## 564 REC-T-166 1
## 565 REC-T-168 1
## 566 REC-T-178 1
## 567 REC-T-181 3
## 568 REC-T-182 2
## 569 REC-T-183 1
## 570 REC-T-192 4
## 571 REC-T-195 1
## 572 REC-T-196 1
## 573 REC-T-197 2
## 574 REC-T-206 1
## 575 REC-T-209 1
## 576 REC-T-210 3
## 577 REC-T-211 9
## 578 REC-T-220 1
## 579 REC-T-221 8
## 580 REC-T-225 1
## 581 REC-T-23 1
## 582 REC-T-235 6
## 583 REC-T-24 1
## 584 REC-T-31 2
## 585 REC-T-32 1
## 586 REC-T-39 2
## 587 REC-T-41 68
## 588 REC-T-42 2
## 589 REC-T-51 4
## 590 REC-T-52 1
## 591 REC-T-60 1
## 592 REC-T-67 1
## 593 REC-T-70 1
## 594 REC-T-71 1
## 595 REC-T-80 4
## 596 REC-T-91 1
## 597 REP 72
## 598 RFQ 51
## 599 RI 130
## 600 RIB 3
## 601 RIL 3
## 602 RIOP-T-20 2
## 603 RIOP-T-41 3
## 604 RIP 4
## 605 RJ 5
## 606 RM 15
## 607 RMO 8
## 608 RO 837
## 609 ROU 33
## 610 ROV 22
## 611 RP 42
## 612 RPJ 75
## 613 RPO 25
## 614 RPP 3
## 615 RPS 27
## 616 RS 27
## 617 RSE 1
## 618 RSM 57
## 619 RSP 9
## 620 RSU 4
## 621 RUC 81
## 622 RV 1
## 623 S-M-1035 1
## 624 S-M-1036 1
## 625 S-M-1100 1
## 626 S-M-1101 1
## 627 S-M-1102 2
## 628 S-M-1163 1
## 629 S-M-1165 4
## 630 S-M-1166 1
## 631 S-M-1226 1
## 632 S-M-1288 5
## 633 S-M-1289 3
## 634 S-M-1352 3
## 635 S-M-1354 1
## 636 S-M-1356 1
## 637 S-M-1358 1
## 638 S-M-1413 1
## 639 S-M-1476 1
## 640 S-M-166 1
## 641 S-M-172 3
## 642 S-M-225 3
## 643 S-M-226 3
## 644 S-M-239 1
## 645 S-M-268 4
## 646 S-M-270 5
## 647 S-M-314 4
## 648 S-M-322 2
## 649 S-M-330 1
## 650 S-M-405 2
## 651 S-M-415 2
## 652 S-M-417 2
## 653 S-M-506 1
## 654 S-M-508 3
## 655 S-M-518 3
## 656 S-M-557 1
## 657 S-M-613 1
## 658 S-M-623 2
## 659 S-M-672 2
## 660 S-M-673 2
## 661 S-M-674 1
## 662 S-M-675 1
## 663 S-M-731 2
## 664 S-M-789 1
## 665 S-M-790 1
## 666 S-M-792 2
## 667 S019 3
## 668 S107 6
## 669 SA 56
## 670 SAB 1
## 671 SC 21
## 672 SCE 4
## 673 SCR 187
## 674 SCX 10
## 675 SDS 16
## 676 SE 67
## 677 SEAL-1 2
## 678 SEAL-10 12
## 679 SEAL-100 23
## 680 SEAL-20 3
## 681 SEAL-30 3
## 682 SEAL-300 2
## 683 SEAL-M-347 1
## 684 SEAL-M-349 4
## 685 SEAL-M-424 6
## 686 SEAL-M-426 9
## 687 SEAL-M-495 6
## 688 SEAL-M-497 1
## 689 SEAL-M-499 4
## 690 SEAL-M-569 1
## 691 SEAL-T-112 1
## 692 SEAL-T-177 1
## 693 SEAL-T-187 1
## 694 SEAL-T-240 3
## 695 SEAL-T-252 2
## 696 SEAL-T-330 2
## 697 SEAL-T-356 2
## 698 SEAL-T-357 1
## 699 SEAL-T-367 3
## 700 SEAL-T-369 1
## 701 SEAL-T-390 1
## 702 SEAL-T-391 1
## 703 SEAL-T-412 4
## 704 SEAL-T-417 4
## 705 SEAL-T-418 1
## 706 SEAL-T-419 1
## 707 SEAL-T-420 2
## 708 SEAL-T-429 6
## 709 SEAL-T-430 2
## 710 SEAL-T-434 1
## 711 SEAL-T-455 1
## 712 SEAL-T-456 1
## 713 SEAL-T-460 2
## 714 SEAL-T-462 2
## 715 SER 31
## 716 SF 1
## 717 SF-T-104 1
## 718 SF-T-105 1
## 719 SF-T-106 2
## 720 SF-T-111 1
## 721 SF-T-112 2
## 722 SF-T-114 1
## 723 SF-T-115 1
## 724 SF-T-118 6
## 725 SF-T-121 1
## 726 SF-T-124 2
## 727 SF-T-125 3
## 728 SF-T-126 1
## 729 SF-T-128 1
## 730 SF-T-130 1
## 731 SF-T-131 2
## 732 SF-T-132 1
## 733 SF-T-133 1
## 734 SF-T-134 1
## 735 SF-T-137 1
## 736 SF-T-138 1
## 737 SF-T-139 1
## 738 SF-T-143 1
## 739 SF-T-81 1
## 740 SF-T-85 1
## 741 SF-T-94 1
## 742 SG 6
## 743 SI 53
## 744 SM 92
## 745 SMA 1
## 746 SMC 52
## 747 SMI 11
## 748 SOJ 2
## 749 SOL-T-150 1
## 750 SOL-T-168 1
## 751 SOL-T-169 2
## 752 SOL-T-170 2
## 753 SOL-T-171 10
## 754 SOL-T-172 1
## 755 SOL-T-191 1
## 756 SOL-T-192 2
## 757 SOL-T-193 2
## 758 SOL-T-194 2
## 759 SP 25
## 760 SUC 3
## 761 SUS 19
## 762 SV 1
## 763 SVM 8
## 764 SZ 513
## 765 TB 6
## 766 TG 19
## 767 TM 77
## 768 TQ 218
## 769 TR 14
## 770 TUPI_NE 4
## 771 TUPI_SUL 2
## 772 UB 137
## 773 UPN 39
## 774 VD 40
## 775 VL 36
## 776 VM 99
## 777 VQ 2
## 778 VR 2
## 779 VRG 8
## 780 XA 28
## 781 <NA> 68Agrupación
TDFBloque$Region <- ifelse(is.na(TDFBloque$Bloque), "Sin información",
ifelse(grepl("^(POT|BPOT|REC|BREC|SEAL|BA|BT|PN|CE|AR|AL|SE|AG|RC|RP)", TDFBloque$Bloque), "Nordeste",
ifelse(grepl("^(BC|BS|ES|C-M|S-M|BM|BES|TUPI|LIBRA|LULA|BUZIOS|IARA|SP|RJ|CM|CF|CR)", TDFBloque$Bloque), "Sudeste",
ifelse(grepl("^(AM|SOL|FZA|AC|SF|SO)", TDFBloque$Bloque), "Norte",
ifelse(TDFBloque$Bloque %in% c("Brasil", "CP", "CAM", "ET", "RO", "FZB", "ARG", "MG", "DJM", "FP", "SZ", "FBM", "PIR", "CS"), "Campos Maduros/Genérico", "Otros")))))
TDFBloque$Region[TDFBloque$Region == "Sin información"] <- "Otros"
head(TDFBloque$Region)## [1] "Otros" "Otros" "Otros" "Otros" "Otros" "Nordeste".
3 Cálculo de frecuencia
Se calculan las frecuencias absolutas y relativas por región.
TDFBloque$Freq <- as.numeric(as.character(TDFBloque$Freq))
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, unionTDFBloque_ <- Datos$TDFBloque
TDFBloque_ <- TDFBloque %>%
group_by(Region) %>%
summarise(
ni = sum(Freq),
hi = round(sum(Freq) / sum(TDFBloque$Freq)*100, 2))
TDFBloque_ <- data.frame(TDFBloque_)Agregar totales
TDFBloque_ <- TDFBloque_ [, c("Region", "ni", "hi")]
total_ni <- sum(TDFBloque_$ni)
total_hi <- 100
TDFBloque_1 <- rbind(TDFBloque_, data.frame( Region = "Total",
ni = total_ni,
hi = total_hi))
print(TDFBloque_1)## Region ni hi
## 1 Campos Maduros/Genérico 12633 42.72
## 2 Nordeste 3538 11.96
## 3 Norte 68 0.23
## 4 Otros 11654 39.40
## 5 Sudeste 1682 5.69
## 6 Total 29575 100.00library(gt)
gt(TDFBloque_1) %>%
tab_header(
title = md("**Cuadro No.1**"),
subtitle = "Distribución de los Bloques de Concesión por Ubicación") %>%
cols_align(align = "center", columns = everything()) %>%
tab_style(
style = list(cell_fill(color = "#2E4053"), cell_text(color = "white", weight = "bold")),
locations = cells_title()
) %>%
tab_style(
style = list(cell_fill(color = "#F2F3F4"), cell_text(weight = "bold", color = "#2E4053")),
locations = cells_column_labels()
) %>%
tab_options(
table.border.top.color = "#2E4053",
table.border.bottom.color = "#2E4053",
column_labels.border.bottom.color = "#2E4053",
data_row.padding = px(6))| Cuadro No.1 | ||
| Distribución de los Bloques de Concesión por Ubicación | ||
| Region | ni | hi |
|---|---|---|
| Campos Maduros/Genérico | 12633 | 42.72 |
| Nordeste | 3538 | 11.96 |
| Norte | 68 | 0.23 |
| Otros | 11654 | 39.40 |
| Sudeste | 1682 | 5.69 |
| Total | 29575 | 100.00 |
4 Gráficas
4.1 Histograma de frecuencia absoluta local
Se visualiza la cantidad de bloques por región.
TDFBloque_2 <- TDFBloque_1[TDFBloque_1$Region!= "Total", ]
par(mar = c(9, 4, 4, 2))
barplot(TDFBloque_2$ni,
main = "Gráfica N°1: Distribución en cantidad de los bloques de concesión por ubicación",
ylab = "Cantidad",
col = "#2E4053", names.arg = TDFBloque_2$Region,
las = 2, cex.names = 0.79, cex.axis = 0.8, cex.main = 1)
mtext("Región", side = 1, line = 7)
4.2 Histograma de frecuencia relativa local (hi)
Se muestra el porcentaje de bloques correspondiente a cada región.
par(mar = c(9, 4, 4, 2))
barplot(TDFBloque_2$hi,
main = "Gráfica N°2: Distribución en porcentaje de los bloques de concesión por ubicación",
ylab = "Porcentaje",
col = "#2E4053", names.arg = TDFBloque_2$Region,
las = 2, cex.names = 0.79, cex.axis = 0.8, cex.main = 1)
mtext("Región", side = 1, line = 7)
5 Modelo de probabilidad
Se construye la distribución de probabilidad empírica por región.
P_Bloque <- TDFBloque_2$ni / sum(TDFBloque_2$ni)
par(mar = c(9, 4, 4, 2))
barplot(P_Bloque,
main = "Gráfica N°3: Distribución de Probabilidad de los bloques de concesión por ubicación",
ylab = "Probabilidad",
col = "#2E4053", names.arg = TDFBloque_2$Region,
las = 2, cex.names = 0.79, cex.axis = 0.8, cex.main = 1)
mtext("Región", side = 1, line = 7)
5.1 Probabilidad
¿Cuál es la probabilidad de que un bloque de concesión seleccionado aleatoriamente pertenezca a la región Nordeste?
x <- round(
(TDFBloque_2$ni[TDFBloque_2$Region == "Nordeste"] /
sum(TDFBloque_2$ni)) * 100, 1)
print(paste("La probabilidad es de:", x, "%"))## [1] "La probabilidad es de: 12 %"La probabilidad calculada indica que cerca del 12 % de los bloques de concesión se ubican en la región Nordeste, lo que evidencia una presencia moderada frente al total nacional.
¿Cuántos bloques de concesión existen en la región Nordeste?
p <- x / 100
n <- 1000
cantidad_esperada <- n * p
Bloque_objetivo <- "Nordeste"
cat("Para n =", n, "bloques de concesión esperada en", Bloque_objetivo, "=", round(cantidad_esperada, 2), "\n")## Para n = 1000 bloques de concesión esperada en Nordeste = 120Con base en el modelo probabilístico, para un conjunto hipotético de 1000 bloques, se espera que aproximadamente 120 bloques pertenezcan a la región Nordeste, manteniendo la proporción observada en los datos reales.