Análisis inferencial de Depósitos Masivos de Sulfuro Volcánicos (Variable Presencia de Intrusión Post Mineralizacion)

Carga de Datos y Librerías

1. CARGA DE LIBRERIAS

# cargar librerías
library(countrycode)
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(knitr)

1.2. CARGA DE DATOS

#cargar datos
datos <- read.csv("C:\\Users\\joeja\\Desktop\\Proyecto Estadística\\Depositos_sulfuro.csv", 
                  header = TRUE, 
                  sep = ";", 
                  dec = ".")

Tabla de distribucion de probabilidad

2. TABLAS DE DISTRIBUCIÓN DE PROBABILIDAD

# Extraer variable
Intrusion_post <- datos$stringer

# Reemplazar NA, espacios, guiones, etc.
Intrusion_post[
  Intrusion_post == "" |
    Intrusion_post == " " |
    Intrusion_post == "-" |
    is.na(Intrusion_post)
] <- "Sin registro"

# Frecuencias y probabilidad
ni <- table(Intrusion_post)
hi <- round(prop.table(ni), 4)
P <- round(hi * 100, 2)

# Crear tabla base
tabla_finalintrusion <- data.frame(
  Intrusion_post = names(ni),
  ni = as.numeric(ni),
  hi = as.numeric(hi),
  P  = as.numeric(P)
)

# Fila TOTAL
fila_total <- data.frame(
  Intrusion_post = "TOTAL",
  ni = sum(tabla_finalintrusion$ni),
  hi = round(sum(tabla_finalintrusion$hi), ),
  P  = round(sum(tabla_finalintrusion$P), )
)

# Mostrar
tabla_finalintrusion <- rbind(tabla_finalintrusion, fila_total)

tabla_finalintrusion

##   Intrusion_post   ni     hi      P
## 1             no  431 0.3954  39.54
## 2   Sin registro  275 0.2523  25.23
## 3            yes  384 0.3523  35.23
## 4          TOTAL 1090 1.0000 100.00

2.1. Tabla Nº1-Distribución de Probabilidad de la presencia de Intrusiones Post Mineralización en Depósito masivos de sulfuros volcanicos

#Crear Tabla
tabla_intrusion_gt <- tabla_finalintrusion %>%
  gt() %>%
  tab_header(
    title = md("**Tabla N° 1**"),
    subtitle = md("Distribución de probabibilidad de la presencia de Intrusiones 
                  Post Mineralización en Depósitos Masivos de Sulfuro Volcánicos")
  ) %>%
  tab_source_note(
    source_note = md("Autor: Grupo 2")
  ) %>%
  tab_options(
    table.border.top.color = "black",
    table.border.bottom.color = "black",
    heading.border.bottom.color = "black",
    heading.border.bottom.width = px(2),
    column_labels.border.top.color = "black",
    column_labels.border.bottom.color = "black",
    column_labels.border.bottom.width = px(2),
    table_body.hlines.color = "gray",
    table_body.border.bottom.color = "black",
    row.striping.include_table_body = TRUE
  ) %>%
  tab_style(
    style = cell_text(weight = "bold"),
    locations = cells_body(rows = Intrusion_post == "Total")
  )

tabla_intrusion_gt

Tabla N° 1
Distribución de probabibilidad de la presencia de Intrusiones Post Mineralización en Depósitos Masivos de Sulfuro Volcánicos
Intrusion_post	ni	hi	P
no	431	0.3954	39.54
Sin registro	275	0.2523	25.23
yes	384	0.3523	35.23
TOTAL	1090	1.0000	100.00
Autor: Grupo 2

Gráficas de distribución de probabilidad

3. DIAGRAMAS DE BARRAS

# Extraer probabilidad (%) sin la fila TOTAL
P_global <- as.numeric(tabla_finalintrusion$P[1:(nrow(tabla_finalintrusion)-1)])

barplot(
  P_global,
  main = "Gráfica Nº1: Distribución de probabibilidad de la presencia de Intrusiones 
                  Post Mineralización en Depósitos Masivos de Sulfuro Volcánicos",
  cex.main = 0.7,
  xlab = "Presencia de Intrusión Post Mineralización",
  ylab = "Probabilidad (%)",
  col = "blue",
  names.arg = tabla_finalintrusion$Intrusion_post[1:(nrow(tabla_finalintrusion)-1)],
  cex.names = 0.9,
  ylim = c(0, 100),
  las = 1
)

Cáclulo de Probabilidad

# Eliminar fila TOTAL
tabla_sin_total <- tabla_finalintrusion[
  tabla_finalintrusion$Intrusion_post != "TOTAL",
]

# Extraer probabilidad de la categoría "no"
prob_no <- tabla_sin_total$P[
  tabla_sin_total$Intrusion_post == "no"
]

# 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é probabilidad existe de que un depósito\n",
    "masivo de sulfuros volcánicos NO presente\n",
    "intrusión post mineralización?\n\n",
    "Probabilidad = ", prob_no, " (%)",
    sep = ""
  ),
  cex = 1.4,
  col = "black",
  font = 2
)