CARGA DE DATOS Y LIBRERÍAS

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(gt)
library(knitr)
library(moments)

datos <- read.csv("~/UNI/estadistica/2026/dataset_geologico_limpio_80.csv")
cat("Dimensiones:", nrow(datos), "observaciones:", ncol(datos), "variables\n\n")
## Dimensiones: 27784 observaciones: 58 variables
# Extraer variable
dispM_raw <- datos$DEVICE
print("Primeros valores de DEVICE:")
## [1] "Primeros valores de DEVICE:"
print(head(dispM_raw))
## [1] "DIGBY DREDGE"         "1 METER PLANKTON NET" "DIGBY DREDGE"        
## [4] "DIGBY DREDGE"         "DIGBY DREDGE"         "DIGBY DREDGE"
dispM <- na.omit(dispM_raw)
summary(dispM)
##    Length     Class      Mode 
##     27783 character character
n <- length(dispM)
print(paste("Número de observaciones válidas (n):", n))
## [1] "Número de observaciones válidas (n): 27783"

TABLA DE DISTRIBUCIÓN DE PROBABILIDAD

Tabla de Distribucion de Probabilidad

# Frecuencia absoluta
ni <- table(dispM)

# Frecuencia relativa (hi)
hi <- prop.table(ni) * 100

# Probabilidad (p)
p <- prop.table(ni)

# Crear tabla con ni, hi y p
TablaProb_DEVICE <- data.frame(
  Device = names(ni),
  ni = as.numeric(ni),
  hi = round(as.numeric(hi), 6),
  p = round(as.numeric(p), 6)
)

# Agregar fila TOTAL
TablaProb_DEVICE[nrow(TablaProb_DEVICE) + 1, ] <- 
  c("TOTAL", sum(ni), sum(hi), sum(p))

TablaProb_DEVICE
##                                                   Device    ni        hi
## 1                                                          297  1.068999
## 2                                         1-METER DREDGE    58  0.208761
## 3                                   1 METER PLANKTON NET     2  0.007199
## 4                                 3/4 METER PLANKTON NET    26  0.093582
## 5                               35MM STEREO CAMERA+DREDG     1  0.003599
## 6                                          ANCHOR DREDGE    19  0.068387
## 7                                          ANCHOR FLUKES     1  0.003599
## 8                                          ANCHOR SAMPLE     1  0.003599
## 9                                       BIG SEABOSS GRAB   302  1.086996
## 10                                           BLAKE TRAWL     1  0.003599
## 11                                       BOOMERANG CORER    15   0.05399
## 12                               BOT DRIFTER SALINOMETER     2  0.007199
## 13                                        BOTTOM DRIFTER     1  0.003599
## 14                                    BOTTOM DRIFTER BUG     2  0.007199
## 15                                        BOTTOM SKIMMER     5  0.017997
## 16                                              BOX CORE   262  0.943023
## 17                                           BOX SAMPLER     3  0.010798
## 18                              CAM GRAB W/O CAM+PLANK T     2  0.007199
## 19                                           CAMERA SLED    10  0.035993
## 20                               CAMPBELL GRAB W/ CAMERA   544  1.958032
## 21                              CAMPBELL GRAB W/O CAMERA   522  1.878847
## 22                                             CHAIN BAG     2  0.007199
## 23                                      CHAIN BAG DREDGE   238  0.856639
## 24                                 CHAIN BAG+PIPE DREDGE    45   0.16197
## 25                              CHAIN BAG+SM PIPE DREDGE    63  0.226757
## 26                                          CHAIN DREDGE     5  0.017997
## 27                               CORING TUBE IN MECH ARM     1  0.003599
## 28                                                 Craib   100  0.359932
## 29                                           CRAIB CORER   125  0.449915
## 30                                           D-L SNAPPER     1  0.003599
## 31                                  D-L SNAPPER + DREDGE     1  0.003599
## 32                                D-L SNAPPER + GRAV COR     4  0.014397
## 33                                          DIETZ-LAFOND     7  0.025195
## 34                                     DIETZ-LAFOND GRAB    10  0.035993
## 35                                  DIETZ-LAFOND SNAPPER   242  0.871036
## 36                                          DIGBY DREDGE    26  0.093582
## 37                              DIGBY DREDGE W/ CNVS BAG     2  0.007199
## 38                                                 DIVER     5  0.017997
## 39                                           DIVER GRABS    44   0.15837
## 40                                       DIVER PUSH CORE    12  0.043192
## 41                                           diver, hand    27  0.097182
## 42                                       DIVERS, BY HAND     2  0.007199
## 43                                           DRAG DREDGE     1  0.003599
## 44                                                DREDGE     5  0.017997
## 45                                            DRILL CORE     9  0.032394
## 46                                             DRILL RIG     1  0.003599
## 47                                          DRILLED CORE    20  0.071986
## 48                                        Drilled Sample    20  0.071986
## 49                                        DRILLED SAMPLE   455  1.637692
## 50                                       EDGERTON CAMERA    11  0.039593
## 51                                           EPA VANVEEN    74   0.26635
## 52                                      GOS ANCHOR CHAIN     1  0.003599
## 53                                                  Grab    24  0.086384
## 54                                                  GRAB   181  0.651478
## 55                                          GRAB SAMPLER   967  3.480546
## 56                                        GRAPPLING HOOK     1  0.003599
## 57                                          Gravity Core    12  0.043192
## 58                                          GRAVITY CORE    82  0.295145
## 59                                         GRAVITY CORER  1302  4.686319
## 60                                                  hand    34  0.122377
## 61                                           Hand Sample    76  0.273549
## 62                                           HAND SAMPLE   320  1.151783
## 63                                      HYDROSTATIC CORE   192   0.69107
## 64                                     Hydrostatic CORER    10  0.035993
## 65                                     HYDROSTATIC CORER  1999  7.195047
## 66                                              IKU GRAB     2  0.007199
## 67                    InterOcean SEABOSS (Osprey Config)   350  1.259763
## 68                                       JAPANESE DREDGE     3  0.010798
## 69                              M-VV GRAB + D-LF SNAPPER     7  0.025195
## 70                              M-VV GRAB + GRAVITY CORE     2  0.007199
## 71                                 MECH ARM + CORE DRILL     1  0.003599
## 72                                     MECH ARM W/ CORER     8  0.028795
## 73                                   MECH ARM W/ PRY BAR     6  0.021596
## 74                                        MECHANICAL ARM    28  0.100781
## 75                                           metal scoop     3  0.010798
## 76                                           METAL SCOOP    15   0.05399
## 77  MicroSEABOSS (a.k.a SEABOSSpro, a.k.a. GoProSEABOSS)    87  0.313141
## 78                                MID-SIZED SEABOSS GRAB    66  0.237555
## 79                                          MINI-SEABOSS    41  0.147572
## 80                                          MINI SEABOSS   342  1.230969
## 81                                    MINIATURE VAN VEEN    69  0.248353
## 82                             MINIATURE VAN VEEN W/O C.     1  0.003599
## 83                                MODIFIED PIERCE DREDGE    21  0.075586
## 84                                     NATURALIST DREDGE   233  0.838642
## 85                                        NISKIN SAMPLER     1  0.003599
## 86                                                 OTHER   453  1.630493
## 87                                            PAUL CORER     1  0.003599
## 88                                           PIPE DREDGE   158  0.568693
## 89                              PIPE DREDGE AND D-L SNAP     1  0.003599
## 90                                  PIPE DREDGE, 12 INCH    19  0.068387
## 91                                           PISTON CORE   130  0.467912
## 92                                          PISTON CORER  1238  4.455962
## 93                            PISTON SAMPLER (HYDRAULIC)     1  0.003599
## 94                                            PONAR GRAB     6  0.021596
## 95                                         PONAR SAMPLER    74   0.26635
## 96                                           Pulse Auger    84  0.302343
## 97                                           PULSE AUGER    80  0.287946
## 98                                            PUSH CORER    11  0.039593
## 99                              RAKED INTO BASKET W/ ARM     1  0.003599
## 100                                             RING NET    42  0.151172
## 101                                          ROCK DREDGE     4  0.014397
## 102                                        ROCKER DREDGE     5  0.017997
## 103                                         ROTARY CORER    20  0.071986
## 104                                       ROTOSONIC CORE   218  0.784652
## 105                                      ROTOSONIC CORER   242  0.871036
## 106                                        ROV PUSH CORE    35  0.125976
## 107                             SALINOMETER,VANDORN BOTS     1  0.003599
## 108                                       SCALLOP DREDGE    66  0.237555
## 109                                                scoop     6  0.021596
## 110                                           SCOOP FISH     4  0.014397
## 111                                              SEABOSS   127  0.457114
## 112                                        SEABOSS (BIG)   246  0.885434
## 113                                        SEDIMENT TRAP   572  2.058813
## 114                                          SHIPEK GRAB   581  2.091207
## 115                                            SLOW CORE     5  0.017997
## 116                                   SMALL SEABOSS GRAB   311   1.11939
## 117                                       SMALL VAN VEEN    10  0.035993
## 118                             SMITH-MAC+MIN VV, CAMERA     1  0.003599
## 119                                 SMITH-MACINTYRE GRAB   244  0.878235
## 120                            SMITH-MACINTYRE W/ CAMERA   104   0.37433
## 121                           SMITH-MACINTYRE W/O CAMERA   367  1.320952
## 122                             SMITH/MAC AND NISKIN SPL     2  0.007199
## 123                                 SUBMERSIBLE ARM GRAB   496  1.785264
## 124                                   Submersible Sample    17  0.061188
## 125                                   SUBMERSIBLE SAMPLE   176  0.633481
## 126                                        TEFLON SHOVEL    14  0.050391
## 127                                                TRAWL     3  0.010798
## 128                                         TRIGGER CORE     7  0.025195
## 129                                               TROWEL    14  0.050391
## 130                                              UNKNOWN   274  0.986215
## 131                                      USGS Slow Corer    20  0.071986
## 132                                    VAN VEEN AND CORE     5  0.017997
## 133                                        Van Veen grab     7  0.025195
## 134                                        Van Veen Grab  1308  4.707915
## 135                                        VAN VEEN GRAB  7835 28.200698
## 136                                              VANVEEN     6  0.021596
## 137                                            Vibracore    15   0.05399
## 138                                            VIBRACORE  1657  5.964079
## 139                                           VIBRACORER    38  0.136774
## 140                                   VISUAL OBSERVATION    87  0.313141
## 141                                VV GRAB + D-L SNAPPER    14  0.050391
## 142                                      WATER PUMP ONLY     1  0.003599
## 143                             WHOI GIANT GRAVITY CORER   198  0.712666
## 144                                                TOTAL 27783       100
##            p
## 1    0.01069
## 2   0.002088
## 3    7.2e-05
## 4   0.000936
## 5    3.6e-05
## 6   0.000684
## 7    3.6e-05
## 8    3.6e-05
## 9    0.01087
## 10   3.6e-05
## 11   0.00054
## 12   7.2e-05
## 13   3.6e-05
## 14   7.2e-05
## 15   0.00018
## 16   0.00943
## 17  0.000108
## 18   7.2e-05
## 19   0.00036
## 20   0.01958
## 21  0.018788
## 22   7.2e-05
## 23  0.008566
## 24   0.00162
## 25  0.002268
## 26   0.00018
## 27   3.6e-05
## 28  0.003599
## 29  0.004499
## 30   3.6e-05
## 31   3.6e-05
## 32  0.000144
## 33  0.000252
## 34   0.00036
## 35   0.00871
## 36  0.000936
## 37   7.2e-05
## 38   0.00018
## 39  0.001584
## 40  0.000432
## 41  0.000972
## 42   7.2e-05
## 43   3.6e-05
## 44   0.00018
## 45  0.000324
## 46   3.6e-05
## 47   0.00072
## 48   0.00072
## 49  0.016377
## 50  0.000396
## 51  0.002663
## 52   3.6e-05
## 53  0.000864
## 54  0.006515
## 55  0.034805
## 56   3.6e-05
## 57  0.000432
## 58  0.002951
## 59  0.046863
## 60  0.001224
## 61  0.002735
## 62  0.011518
## 63  0.006911
## 64   0.00036
## 65   0.07195
## 66   7.2e-05
## 67  0.012598
## 68  0.000108
## 69  0.000252
## 70   7.2e-05
## 71   3.6e-05
## 72  0.000288
## 73  0.000216
## 74  0.001008
## 75  0.000108
## 76   0.00054
## 77  0.003131
## 78  0.002376
## 79  0.001476
## 80   0.01231
## 81  0.002484
## 82   3.6e-05
## 83  0.000756
## 84  0.008386
## 85   3.6e-05
## 86  0.016305
## 87   3.6e-05
## 88  0.005687
## 89   3.6e-05
## 90  0.000684
## 91  0.004679
## 92   0.04456
## 93   3.6e-05
## 94  0.000216
## 95  0.002663
## 96  0.003023
## 97  0.002879
## 98  0.000396
## 99   3.6e-05
## 100 0.001512
## 101 0.000144
## 102  0.00018
## 103  0.00072
## 104 0.007847
## 105  0.00871
## 106  0.00126
## 107  3.6e-05
## 108 0.002376
## 109 0.000216
## 110 0.000144
## 111 0.004571
## 112 0.008854
## 113 0.020588
## 114 0.020912
## 115  0.00018
## 116 0.011194
## 117  0.00036
## 118  3.6e-05
## 119 0.008782
## 120 0.003743
## 121  0.01321
## 122  7.2e-05
## 123 0.017853
## 124 0.000612
## 125 0.006335
## 126 0.000504
## 127 0.000108
## 128 0.000252
## 129 0.000504
## 130 0.009862
## 131  0.00072
## 132  0.00018
## 133 0.000252
## 134 0.047079
## 135 0.282007
## 136 0.000216
## 137  0.00054
## 138 0.059641
## 139 0.001368
## 140 0.003131
## 141 0.000504
## 142  3.6e-05
## 143 0.007127
## 144        1
TablaProb_DEVICE %>%
  gt() %>%
  tab_header(
    title = md("**Tabla Nº1**"),
    subtitle = md("**Distribución de probabilidad simple de la variable DEVICE 
                  en depósitos marinos**")
  ) %>%
  tab_source_note(
    source_note = md("**Autor: Grupo 3**")
  ) %>%
  cols_label(
    Device = "Categoría",
    ni = "Frecuencia (ni)",
    hi = "Frecuencia relativa hi (%)",
    p = "Probabilidad p (%)"
  ) %>%
  tab_options(
    table.border.top.color = "black",
    table.border.bottom.color = "black",
    table.border.top.style = "solid",
    table.border.bottom.style = "solid",
    column_labels.border.top.color = "black",
    column_labels.border.bottom.color = "black",
    column_labels.border.bottom.width = px(2),
    row.striping.include_table_body = TRUE,
    heading.border.bottom.color = "black",
    heading.border.bottom.width = px(2),
    table_body.hlines.color = "gray",
    table_body.border.bottom.color = "black"
  ) %>%
  tab_style(
    style = cell_text(weight = "bold"),
    locations = cells_body(
      rows = Device == "TOTAL"
    )
  )
Tabla Nº1
Distribución de probabilidad simple de la variable DEVICE en depósitos marinos
Categoría Frecuencia (ni) Frecuencia relativa hi (%) Probabilidad p (%)
297 1.068999 0.01069
1-METER DREDGE 58 0.208761 0.002088
1 METER PLANKTON NET 2 0.007199 7.2e-05
3/4 METER PLANKTON NET 26 0.093582 0.000936
35MM STEREO CAMERA+DREDG 1 0.003599 3.6e-05
ANCHOR DREDGE 19 0.068387 0.000684
ANCHOR FLUKES 1 0.003599 3.6e-05
ANCHOR SAMPLE 1 0.003599 3.6e-05
BIG SEABOSS GRAB 302 1.086996 0.01087
BLAKE TRAWL 1 0.003599 3.6e-05
BOOMERANG CORER 15 0.05399 0.00054
BOT DRIFTER SALINOMETER 2 0.007199 7.2e-05
BOTTOM DRIFTER 1 0.003599 3.6e-05
BOTTOM DRIFTER BUG 2 0.007199 7.2e-05
BOTTOM SKIMMER 5 0.017997 0.00018
BOX CORE 262 0.943023 0.00943
BOX SAMPLER 3 0.010798 0.000108
CAM GRAB W/O CAM+PLANK T 2 0.007199 7.2e-05
CAMERA SLED 10 0.035993 0.00036
CAMPBELL GRAB W/ CAMERA 544 1.958032 0.01958
CAMPBELL GRAB W/O CAMERA 522 1.878847 0.018788
CHAIN BAG 2 0.007199 7.2e-05
CHAIN BAG DREDGE 238 0.856639 0.008566
CHAIN BAG+PIPE DREDGE 45 0.16197 0.00162
CHAIN BAG+SM PIPE DREDGE 63 0.226757 0.002268
CHAIN DREDGE 5 0.017997 0.00018
CORING TUBE IN MECH ARM 1 0.003599 3.6e-05
Craib 100 0.359932 0.003599
CRAIB CORER 125 0.449915 0.004499
D-L SNAPPER 1 0.003599 3.6e-05
D-L SNAPPER + DREDGE 1 0.003599 3.6e-05
D-L SNAPPER + GRAV COR 4 0.014397 0.000144
DIETZ-LAFOND 7 0.025195 0.000252
DIETZ-LAFOND GRAB 10 0.035993 0.00036
DIETZ-LAFOND SNAPPER 242 0.871036 0.00871
DIGBY DREDGE 26 0.093582 0.000936
DIGBY DREDGE W/ CNVS BAG 2 0.007199 7.2e-05
DIVER 5 0.017997 0.00018
DIVER GRABS 44 0.15837 0.001584
DIVER PUSH CORE 12 0.043192 0.000432
diver, hand 27 0.097182 0.000972
DIVERS, BY HAND 2 0.007199 7.2e-05
DRAG DREDGE 1 0.003599 3.6e-05
DREDGE 5 0.017997 0.00018
DRILL CORE 9 0.032394 0.000324
DRILL RIG 1 0.003599 3.6e-05
DRILLED CORE 20 0.071986 0.00072
Drilled Sample 20 0.071986 0.00072
DRILLED SAMPLE 455 1.637692 0.016377
EDGERTON CAMERA 11 0.039593 0.000396
EPA VANVEEN 74 0.26635 0.002663
GOS ANCHOR CHAIN 1 0.003599 3.6e-05
Grab 24 0.086384 0.000864
GRAB 181 0.651478 0.006515
GRAB SAMPLER 967 3.480546 0.034805
GRAPPLING HOOK 1 0.003599 3.6e-05
Gravity Core 12 0.043192 0.000432
GRAVITY CORE 82 0.295145 0.002951
GRAVITY CORER 1302 4.686319 0.046863
hand 34 0.122377 0.001224
Hand Sample 76 0.273549 0.002735
HAND SAMPLE 320 1.151783 0.011518
HYDROSTATIC CORE 192 0.69107 0.006911
Hydrostatic CORER 10 0.035993 0.00036
HYDROSTATIC CORER 1999 7.195047 0.07195
IKU GRAB 2 0.007199 7.2e-05
InterOcean SEABOSS (Osprey Config) 350 1.259763 0.012598
JAPANESE DREDGE 3 0.010798 0.000108
M-VV GRAB + D-LF SNAPPER 7 0.025195 0.000252
M-VV GRAB + GRAVITY CORE 2 0.007199 7.2e-05
MECH ARM + CORE DRILL 1 0.003599 3.6e-05
MECH ARM W/ CORER 8 0.028795 0.000288
MECH ARM W/ PRY BAR 6 0.021596 0.000216
MECHANICAL ARM 28 0.100781 0.001008
metal scoop 3 0.010798 0.000108
METAL SCOOP 15 0.05399 0.00054
MicroSEABOSS (a.k.a SEABOSSpro, a.k.a. GoProSEABOSS) 87 0.313141 0.003131
MID-SIZED SEABOSS GRAB 66 0.237555 0.002376
MINI-SEABOSS 41 0.147572 0.001476
MINI SEABOSS 342 1.230969 0.01231
MINIATURE VAN VEEN 69 0.248353 0.002484
MINIATURE VAN VEEN W/O C. 1 0.003599 3.6e-05
MODIFIED PIERCE DREDGE 21 0.075586 0.000756
NATURALIST DREDGE 233 0.838642 0.008386
NISKIN SAMPLER 1 0.003599 3.6e-05
OTHER 453 1.630493 0.016305
PAUL CORER 1 0.003599 3.6e-05
PIPE DREDGE 158 0.568693 0.005687
PIPE DREDGE AND D-L SNAP 1 0.003599 3.6e-05
PIPE DREDGE, 12 INCH 19 0.068387 0.000684
PISTON CORE 130 0.467912 0.004679
PISTON CORER 1238 4.455962 0.04456
PISTON SAMPLER (HYDRAULIC) 1 0.003599 3.6e-05
PONAR GRAB 6 0.021596 0.000216
PONAR SAMPLER 74 0.26635 0.002663
Pulse Auger 84 0.302343 0.003023
PULSE AUGER 80 0.287946 0.002879
PUSH CORER 11 0.039593 0.000396
RAKED INTO BASKET W/ ARM 1 0.003599 3.6e-05
RING NET 42 0.151172 0.001512
ROCK DREDGE 4 0.014397 0.000144
ROCKER DREDGE 5 0.017997 0.00018
ROTARY CORER 20 0.071986 0.00072
ROTOSONIC CORE 218 0.784652 0.007847
ROTOSONIC CORER 242 0.871036 0.00871
ROV PUSH CORE 35 0.125976 0.00126
SALINOMETER,VANDORN BOTS 1 0.003599 3.6e-05
SCALLOP DREDGE 66 0.237555 0.002376
scoop 6 0.021596 0.000216
SCOOP FISH 4 0.014397 0.000144
SEABOSS 127 0.457114 0.004571
SEABOSS (BIG) 246 0.885434 0.008854
SEDIMENT TRAP 572 2.058813 0.020588
SHIPEK GRAB 581 2.091207 0.020912
SLOW CORE 5 0.017997 0.00018
SMALL SEABOSS GRAB 311 1.11939 0.011194
SMALL VAN VEEN 10 0.035993 0.00036
SMITH-MAC+MIN VV, CAMERA 1 0.003599 3.6e-05
SMITH-MACINTYRE GRAB 244 0.878235 0.008782
SMITH-MACINTYRE W/ CAMERA 104 0.37433 0.003743
SMITH-MACINTYRE W/O CAMERA 367 1.320952 0.01321
SMITH/MAC AND NISKIN SPL 2 0.007199 7.2e-05
SUBMERSIBLE ARM GRAB 496 1.785264 0.017853
Submersible Sample 17 0.061188 0.000612
SUBMERSIBLE SAMPLE 176 0.633481 0.006335
TEFLON SHOVEL 14 0.050391 0.000504
TRAWL 3 0.010798 0.000108
TRIGGER CORE 7 0.025195 0.000252
TROWEL 14 0.050391 0.000504
UNKNOWN 274 0.986215 0.009862
USGS Slow Corer 20 0.071986 0.00072
VAN VEEN AND CORE 5 0.017997 0.00018
Van Veen grab 7 0.025195 0.000252
Van Veen Grab 1308 4.707915 0.047079
VAN VEEN GRAB 7835 28.200698 0.282007
VANVEEN 6 0.021596 0.000216
Vibracore 15 0.05399 0.00054
VIBRACORE 1657 5.964079 0.059641
VIBRACORER 38 0.136774 0.001368
VISUAL OBSERVATION 87 0.313141 0.003131
VV GRAB + D-L SNAPPER 14 0.050391 0.000504
WATER PUMP ONLY 1 0.003599 3.6e-05
WHOI GIANT GRAVITY CORER 198 0.712666 0.007127
TOTAL 27783 100 1
Autor: Grupo 3

TABLA DE DISTRIBUCIÓN DE PROBABILIDAD AGRUPADA

Tabla de Distribucion de Probabilidad por Agrupación Debido a que la variable presenta numerosos tipos de dispositivos, se decidió agruparlos por categorías de muestreo

cat("TABLA DE DISTRIBUCIÓN DE PROBABILIDAD AGRUPADA\n")
## TABLA DE DISTRIBUCIÓN DE PROBABILIDAD AGRUPADA
device_grupo <- ifelse(grepl("grab|van veen|seaboss|ponar|smith|craib", dispM, ignore.case = TRUE),
                       "Pinzas de Sedimento",
                       ifelse(grepl("core|corer|vibracore|piston|gravity|rotosonic", dispM, ignore.case = TRUE),
                              "Extracción de Núcleos",
                              ifelse(grepl("dredge", dispM, ignore.case = TRUE),
                                     "Dragado",
                                     ifelse(grepl("net|trawl", dispM, ignore.case = TRUE),
                                            "Redes",
                                            ifelse(grepl("drill|auger", dispM, ignore.case = TRUE),
                                                   "Perforación",
                                                   ifelse(grepl("submersible|rov|diver", dispM, ignore.case = TRUE),
                                                          "Sumergible",
                                                          ifelse(grepl("camera|visual", dispM, ignore.case = TRUE),
                                                                 "Observación/Cámara",
                                                                 ifelse(grepl("hand|scoop|shovel|trowel", dispM, ignore.case = TRUE),
                                                                        "Muestreo Manual",
                                                                        ifelse(grepl("pump|niskin|salinometer", dispM, ignore.case = TRUE),
                                                                               "Muestreo Químico",
                                                                               "Otros")))))))))

# Frecuencia absoluta
ni_device <- table(device_grupo)

# Frecuencia relativa (hi)
hi_device <- prop.table(ni_device) * 100

# Probabilidad (p)
p_device <- prop.table(ni_device)

# Tabla agrupada con ni, hi y p
tabla_device <- data.frame(
  Device = names(ni_device),
  ni = as.numeric(ni_device),
  hi = round(as.numeric(hi_device), 6),
  P = round(as.numeric(p_device), 6)
)

# Agregar fila TOTAL
tabla_device[nrow(tabla_device) + 1, ] <- c(
  "TOTAL",
  sum(tabla_device$ni),
  sum(tabla_device$hi),
  sum(tabla_device$P)
)

tabla_device
##                   Device    ni        hi        P
## 1                Dragado   973  3.502142 0.035021
## 2  Extracción de Núcleos  7760 27.930749 0.279307
## 3        Muestreo Manual   486  1.749271 0.017493
## 4       Muestreo Químico     5  0.017997  0.00018
## 5     Observación/Cámara   109  0.392326 0.003923
## 6                  Otros  1983  7.137458 0.071375
## 7            Perforación   640  2.303567 0.023036
## 8    Pinzas de Sedimento 15526 55.883094 0.558831
## 9                  Redes    74   0.26635 0.002663
## 10            Sumergible   227  0.817046  0.00817
## 11                 TOTAL 27783       100 0.999999
tabla_device %>%
  gt() %>%
  tab_header(
    title = md("**Tabla Nº2**"),
    subtitle = md("**Distribución de probabilidad agrupada de la variable DEVICE**")
  ) %>%
  tab_source_note(
    source_note = md("**Autor: Grupo 3**")
  ) %>%
  cols_label(
    Device = "Tipo de dispositivo",
    ni = "Frecuencia (ni)",
    hi = "Frecuencia relativa hi (%)",
    P = "Probabilidad P (%)"
  ) %>%
  tab_options(
    table.border.top.color = "black",
    table.border.bottom.color = "black",
    table.border.top.style = "solid",
    table.border.bottom.style = "solid",
    column_labels.border.top.color = "black",
    column_labels.border.bottom.color = "black",
    column_labels.border.bottom.width = px(2),
    row.striping.include_table_body = TRUE,
    heading.border.bottom.color = "black",
    heading.border.bottom.width = px(2),
    table_body.hlines.color = "gray",
    table_body.border.bottom.color = "black"
  ) %>%
  tab_style(
    style = cell_text(weight = "bold"),
    locations = cells_body(
      rows = Device == "TOTAL"
    )
  )
Tabla Nº2
Distribución de probabilidad agrupada de la variable DEVICE
Tipo de dispositivo Frecuencia (ni) Frecuencia relativa hi (%) Probabilidad P (%)
Dragado 973 3.502142 0.035021
Extracción de Núcleos 7760 27.930749 0.279307
Muestreo Manual 486 1.749271 0.017493
Muestreo Químico 5 0.017997 0.00018
Observación/Cámara 109 0.392326 0.003923
Otros 1983 7.137458 0.071375
Perforación 640 2.303567 0.023036
Pinzas de Sedimento 15526 55.883094 0.558831
Redes 74 0.26635 0.002663
Sumergible 227 0.817046 0.00817
TOTAL 27783 100 0.999999
Autor: Grupo 3

GRÁFICAS DE DISTRIBUCIÓN DE PROBABILIDAD

# Extraer probabilidades agrupadas en porcentaje, sin incluir TOTAL
tabla_device$P <- as.numeric(tabla_device$P)

P_global <- as.numeric(
  tabla_device$P[1:(nrow(tabla_device) - 1)]
) * 100

dispositivos <- tabla_device$Device[
  1:(nrow(tabla_device) - 1)
]

# Diagrama de barras
par(mfrow = c(1,1))
par(mar = c(10,4,4,2))

barplot(
  P_global,
  main = "Gráfica Nº1: Distribución de probabilidad de los dispositivos de muestreo en depósitos marinos",
  xlab = "Tipo de dispositivo",
  ylab = "Probabilidad (%)",
  col = "blue",
  names.arg = dispositivos,
  cex.names = 0.8,
  las = 2,
  ylim = c(0, max(P_global) * 1.15)
)

abline(h = 0, lwd = 1)

CALCULOS DE PROBABILIDAD

# Identificar el dispositivo con mayor probabilidad (excluye TOTAL)
tabla_sin_total <- tabla_device[tabla_device$Device != "TOTAL", ]

dispositivo_mayor <- tabla_sin_total$Device[
  which.max(as.numeric(tabla_sin_total$P))
]

prob_mayor <- tabla_sin_total$P[
  which.max(as.numeric(tabla_sin_total$P))
]

# Convertir a numérico
prob_mayor <- as.numeric(prob_mayor)

# Gráfico de texto explicativo
plot(1, type = "n", axes = FALSE, xlab = "", ylab = "")

text(
  x = 1, y = 1,
  labels = paste(
    "Cálculo de probabilidad\n(Estimación general)\n\n",
    "¿Qué tipo de dispositivo es más probable\n",
    "que sea utilizado en la recolección de\n",
    "depósitos marinos?\n\n",
    "R: ", dispositivo_mayor, "\n",
    "Probabilidad = ", round(prob_mayor*100, 2), " %",
    sep = ""
  ),
  cex = 1.4,
  col = "black",
  font = 2
)

CONCLUSIONES

El dispositivo más probable de ser utilizado es Pinzas de Sedimento, con una probabilidad del 55.88%, lo que indica que más de la mitad de las observaciones corresponden a este tipo de instrumento.