Análisis de Depósitos Masivos de Sulfuro Volcánicos (Variable Deformación Post Mineralizacion)
Grupo 2
2025-12-09
Carga de Datos
1. CARGA DE LIBRERIAS
# cargar librerías
library(readxl)
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, unionlibrary(ggplot2)
library(gt)1.2. CARGA DE DATOS
#cargar datos
datos <- read.csv("C:\\Users\\joeja\\Desktop\\Proyecto Estadística\\Depositos_sulfuro.csv",
header = TRUE,
sep = ";",
dec = ".")Verificamos los datos de la variable a trabajar
str(datos$postmindef)## chr [1:1090] "yes" "no" "" "yes" "yes" "yes" "yes" "yes" "" "yes" "yes" "" ...Limpiar, corregir y estandarizar los valores
#Limpiar, corregir y estandarizar los valores de los datos
datos <- datos %>%
mutate(
postmindef = trimws(postmindef),
postmindef = case_when(
postmindef %in% c("yee", "Yes", "yes") ~ "YES",
postmindef %in% c("no", "No", "NO") ~ "NO",
postmindef == "" | is.na(postmindef) ~ "Sin Registro",
TRUE ~ postmindef
)
)
#cargar variable
datos$postmindef <- trimws(datos$postmindef)
Deformacion <- datos$postmindef
Deformacion## [1] "YES" "NO" "Sin Registro" "YES"
## [5] "YES" "YES" "YES" "YES"
## [9] "Sin Registro" "YES" "YES" "Sin Registro"
## [13] "YES" "Sin Registro" "YES" "YES"
## [17] "YES" "YES" "YES" "Sin Registro"
## [21] "YES" "YES" "Sin Registro" "YES"
## [25] "YES" "YES" "YES" "Sin Registro"
## [29] "Sin Registro" "YES" "YES" "YES"
## [33] "Sin Registro" "YES" "NO" "YES"
## [37] "YES" "NO" "YES" "Sin Registro"
## [41] "YES" "Sin Registro" "Sin Registro" "YES"
## [45] "YES" "YES" "YES" "YES"
## [49] "YES" "YES" "YES" "YES"
## [53] "YES" "YES" "YES" "YES"
## [57] "YES" "Sin Registro" "YES" "YES"
## [61] "YES" "Sin Registro" "YES" "YES"
## [65] "YES" "YES" "YES" "YES"
## [69] "YES" "YES" "YES" "YES"
## [73] "YES" "YES" "Sin Registro" "YES"
## [77] "YES" "YES" "YES" "YES"
## [81] "YES" "YES" "YES" "YES"
## [85] "YES" "YES" "YES" "YES"
## [89] "YES" "YES" "YES" "YES"
## [93] "YES" "YES" "YES" "YES"
## [97] "YES" "YES" "YES" "Sin Registro"
## [101] "YES" "YES" "Sin Registro" "YES"
## [105] "YES" "YES" "YES" "YES"
## [109] "Sin Registro" "YES" "YES" "YES"
## [113] "YES" "YES" "YES" "Sin Registro"
## [117] "YES" "YES" "YES" "YES"
## [121] "YES" "YES" "YES" "YES"
## [125] "YES" "YES" "YES" "YES"
## [129] "YES" "NO" "Sin Registro" "YES"
## [133] "YES" "YES" "YES" "Sin Registro"
## [137] "YES" "YES" "Sin Registro" "YES"
## [141] "YES" "YES" "YES" "YES"
## [145] "Sin Registro" "YES" "YES" "YES"
## [149] "Sin Registro" "YES" "YES" "NO"
## [153] "YES" "YES" "YES" "Sin Registro"
## [157] "YES" "YES" "YES" "YES"
## [161] "Sin Registro" "YES" "Sin Registro" "YES"
## [165] "YES" "YES" "YES" "YES"
## [169] "YES" "Sin Registro" "NO" "NO"
## [173] "Sin Registro" "YES" "YES" "Sin Registro"
## [177] "YES" "YES" "YES" "YES"
## [181] "YES" "YES" "YES" "YES"
## [185] "Sin Registro" "YES" "YES" "Sin Registro"
## [189] "YES" "YES" "YES" "YES"
## [193] "YES" "YES" "YES" "Sin Registro"
## [197] "YES" "Sin Registro" "NO" "YES"
## [201] "YES" "YES" "YES" "YES"
## [205] "YES" "YES" "YES" "YES"
## [209] "YES" "YES" "Sin Registro" "YES"
## [213] "NO" "Sin Registro" "Sin Registro" "Sin Registro"
## [217] "YES" "YES" "YES" "YES"
## [221] "YES" "YES" "YES" "YES"
## [225] "YES" "YES" "YES" "Sin Registro"
## [229] "Sin Registro" "Sin Registro" "NO" "NO"
## [233] "Sin Registro" "YES" "Sin Registro" "NO"
## [237] "NO" "YES" "YES" "NO"
## [241] "YES" "YES" "YES" "YES"
## [245] "YES" "YES" "YES" "YES"
## [249] "YES" "YES" "YES" "YES"
## [253] "Sin Registro" "YES" "YES" "YES"
## [257] "YES" "YES" "Sin Registro" "Sin Registro"
## [261] "YES" "YES" "YES" "YES"
## [265] "YES" "YES" "YES" "YES"
## [269] "Sin Registro" "YES" "Sin Registro" "YES"
## [273] "YES" "YES" "Sin Registro" "YES"
## [277] "YES" "Sin Registro" "YES" "Sin Registro"
## [281] "YES" "Sin Registro" "YES" "YES"
## [285] "YES" "YES" "YES" "YES"
## [289] "YES" "YES" "Sin Registro" "Sin Registro"
## [293] "Sin Registro" "YES" "YES" "YES"
## [297] "YES" "YES" "YES" "Sin Registro"
## [301] "NO" "YES" "NO" "YES"
## [305] "YES" "YES" "YES" "Sin Registro"
## [309] "YES" "YES" "YES" "YES"
## [313] "YES" "YES" "YES" "YES"
## [317] "YES" "YES" "YES" "YES"
## [321] "YES" "YES" "YES" "YES"
## [325] "YES" "Sin Registro" "YES" "Sin Registro"
## [329] "NO" "YES" "YES" "YES"
## [333] "Sin Registro" "YES" "YES" "YES"
## [337] "YES" "YES" "YES" "YES"
## [341] "NO" "Sin Registro" "YES" "YES"
## [345] "YES" "Sin Registro" "YES" "Sin Registro"
## [349] "YES" "YES" "YES" "Sin Registro"
## [353] "YES" "YES" "YES" "YES"
## [357] "YES" "YES" "Sin Registro" "YES"
## [361] "YES" "Sin Registro" "YES" "YES"
## [365] "Sin Registro" "Sin Registro" "YES" "YES"
## [369] "YES" "YES" "Sin Registro" "YES"
## [373] "YES" "YES" "YES" "YES"
## [377] "YES" "YES" "YES" "YES"
## [381] "Sin Registro" "Sin Registro" "YES" "YES"
## [385] "YES" "Sin Registro" "YES" "YES"
## [389] "YES" "YES" "YES" "YES"
## [393] "YES" "YES" "YES" "YES"
## [397] "Sin Registro" "YES" "YES" "YES"
## [401] "Sin Registro" "YES" "YES" "Sin Registro"
## [405] "YES" "YES" "YES" "Sin Registro"
## [409] "YES" "YES" "YES" "YES"
## [413] "YES" "Sin Registro" "YES" "YES"
## [417] "YES" "YES" "NO" "Sin Registro"
## [421] "YES" "YES" "Sin Registro" "Sin Registro"
## [425] "Sin Registro" "NO" "YES" "YES"
## [429] "YES" "NO" "YES" "YES"
## [433] "YES" "YES" "YES" "YES"
## [437] "YES" "YES" "Sin Registro" "Sin Registro"
## [441] "YES" "Sin Registro" "Sin Registro" "Sin Registro"
## [445] "Sin Registro" "Sin Registro" "YES" "YES"
## [449] "Sin Registro" "YES" "YES" "YES"
## [453] "Sin Registro" "NO" "YES" "YES"
## [457] "YES" "YES" "Sin Registro" "YES"
## [461] "Sin Registro" "YES" "Sin Registro" "YES"
## [465] "YES" "YES" "YES" "YES"
## [469] "YES" "YES" "YES" "YES"
## [473] "YES" "YES" "YES" "YES"
## [477] "YES" "YES" "NO" "NO"
## [481] "YES" "YES" "Sin Registro" "YES"
## [485] "Sin Registro" "YES" "Sin Registro" "YES"
## [489] "YES" "YES" "YES" "YES"
## [493] "Sin Registro" "YES" "YES" "Sin Registro"
## [497] "YES" "YES" "YES" "YES"
## [501] "YES" "YES" "YES" "YES"
## [505] "YES" "YES" "YES" "YES"
## [509] "YES" "YES" "YES" "Sin Registro"
## [513] "YES" "YES" "YES" "YES"
## [517] "YES" "YES" "YES" "YES"
## [521] "YES" "Sin Registro" "YES" "YES"
## [525] "YES" "YES" "YES" "YES"
## [529] "YES" "YES" "YES" "Sin Registro"
## [533] "YES" "YES" "YES" "YES"
## [537] "Sin Registro" "Sin Registro" "Sin Registro" "Sin Registro"
## [541] "Sin Registro" "Sin Registro" "YES" "YES"
## [545] "YES" "YES" "YES" "YES"
## [549] "YES" "YES" "YES" "YES"
## [553] "YES" "YES" "YES" "YES"
## [557] "YES" "YES" "YES" "YES"
## [561] "YES" "YES" "YES" "YES"
## [565] "YES" "YES" "YES" "YES"
## [569] "YES" "YES" "YES" "YES"
## [573] "YES" "YES" "YES" "YES"
## [577] "YES" "YES" "Sin Registro" "YES"
## [581] "YES" "Sin Registro" "Sin Registro" "Sin Registro"
## [585] "Sin Registro" "YES" "Sin Registro" "Sin Registro"
## [589] "Sin Registro" "YES" "Sin Registro" "YES"
## [593] "YES" "YES" "YES" "Sin Registro"
## [597] "YES" "Sin Registro" "YES" "YES"
## [601] "YES" "YES" "YES" "YES"
## [605] "YES" "YES" "YES" "YES"
## [609] "YES" "Sin Registro" "Sin Registro" "Sin Registro"
## [613] "NO" "YES" "YES" "Sin Registro"
## [617] "Sin Registro" "YES" "YES" "Sin Registro"
## [621] "Sin Registro" "Sin Registro" "Sin Registro" "Sin Registro"
## [625] "YES" "YES" "YES" "YES"
## [629] "YES" "Sin Registro" "YES" "YES"
## [633] "YES" "YES" "YES" "Sin Registro"
## [637] "Sin Registro" "YES" "YES" "YES"
## [641] "YES" "YES" "YES" "YES"
## [645] "YES" "YES" "YES" "YES"
## [649] "YES" "YES" "YES" "YES"
## [653] "YES" "YES" "YES" "YES"
## [657] "YES" "YES" "YES" "YES"
## [661] "YES" "Sin Registro" "YES" "YES"
## [665] "YES" "YES" "YES" "YES"
## [669] "YES" "YES" "YES" "YES"
## [673] "YES" "YES" "NO" "YES"
## [677] "Sin Registro" "YES" "YES" "YES"
## [681] "YES" "YES" "Sin Registro" "YES"
## [685] "YES" "YES" "YES" "YES"
## [689] "YES" "YES" "YES" "YES"
## [693] "YES" "YES" "YES" "YES"
## [697] "YES" "YES" "YES" "YES"
## [701] "Sin Registro" "YES" "YES" "Sin Registro"
## [705] "YES" "YES" "YES" "NO"
## [709] "YES" "YES" "YES" "YES"
## [713] "YES" "YES" "NO" "YES"
## [717] "YES" "YES" "NO" "YES"
## [721] "YES" "YES" "YES" "YES"
## [725] "YES" "YES" "NO" "NO"
## [729] "YES" "YES" "YES" "YES"
## [733] "YES" "YES" "YES" "YES"
## [737] "YES" "YES" "YES" "YES"
## [741] "Sin Registro" "YES" "Sin Registro" "Sin Registro"
## [745] "YES" "YES" "YES" "Sin Registro"
## [749] "YES" "YES" "YES" "YES"
## [753] "YES" "YES" "Sin Registro" "Sin Registro"
## [757] "YES" "YES" "YES" "YES"
## [761] "YES" "YES" "YES" "YES"
## [765] "YES" "YES" "YES" "YES"
## [769] "YES" "YES" "YES" "YES"
## [773] "YES" "YES" "YES" "YES"
## [777] "YES" "YES" "YES" "YES"
## [781] "YES" "YES" "YES" "YES"
## [785] "YES" "YES" "YES" "YES"
## [789] "NO" "YES" "YES" "YES"
## [793] "YES" "YES" "YES" "YES"
## [797] "YES" "YES" "YES" "YES"
## [801] "YES" "YES" "YES" "Sin Registro"
## [805] "Sin Registro" "Sin Registro" "Sin Registro" "Sin Registro"
## [809] "Sin Registro" "YES" "YES" "Sin Registro"
## [813] "YES" "YES" "YES" "YES"
## [817] "Sin Registro" "Sin Registro" "YES" "Sin Registro"
## [821] "Sin Registro" "Sin Registro" "Sin Registro" "YES"
## [825] "YES" "YES" "Sin Registro" "Sin Registro"
## [829] "YES" "Sin Registro" "YES" "YES"
## [833] "Sin Registro" "YES" "YES" "Sin Registro"
## [837] "YES" "YES" "YES" "Sin Registro"
## [841] "Sin Registro" "Sin Registro" "Sin Registro" "Sin Registro"
## [845] "Sin Registro" "YES" "YES" "YES"
## [849] "YES" "YES" "YES" "Sin Registro"
## [853] "YES" "Sin Registro" "Sin Registro" "Sin Registro"
## [857] "Sin Registro" "Sin Registro" "Sin Registro" "Sin Registro"
## [861] "Sin Registro" "YES" "YES" "YES"
## [865] "YES" "YES" "Sin Registro" "Sin Registro"
## [869] "YES" "YES" "YES" "YES"
## [873] "Sin Registro" "YES" "YES" "YES"
## [877] "YES" "YES" "YES" "YES"
## [881] "YES" "YES" "YES" "YES"
## [885] "YES" "YES" "YES" "Sin Registro"
## [889] "Sin Registro" "YES" "YES" "YES"
## [893] "YES" "YES" "YES" "YES"
## [897] "YES" "YES" "YES" "YES"
## [901] "YES" "YES" "YES" "Sin Registro"
## [905] "YES" "YES" "YES" "YES"
## [909] "YES" "YES" "YES" "YES"
## [913] "YES" "YES" "Sin Registro" "YES"
## [917] "YES" "YES" "YES" "YES"
## [921] "YES" "YES" "YES" "YES"
## [925] "YES" "YES" "YES" "YES"
## [929] "YES" "YES" "YES" "YES"
## [933] "YES" "YES" "YES" "YES"
## [937] "YES" "YES" "YES" "YES"
## [941] "YES" "YES" "YES" "YES"
## [945] "YES" "YES" "YES" "YES"
## [949] "YES" "YES" "YES" "YES"
## [953] "YES" "YES" "YES" "YES"
## [957] "YES" "YES" "YES" "YES"
## [961] "YES" "YES" "YES" "YES"
## [965] "YES" "YES" "YES" "YES"
## [969] "YES" "YES" "YES" "YES"
## [973] "YES" "YES" "YES" "YES"
## [977] "YES" "YES" "YES" "YES"
## [981] "YES" "YES" "YES" "YES"
## [985] "Sin Registro" "Sin Registro" "YES" "YES"
## [989] "YES" "YES" "YES" "YES"
## [993] "YES" "YES" "YES" "YES"
## [997] "YES" "YES" "YES" "YES"
## [1001] "YES" "YES" "YES" "Sin Registro"
## [1005] "YES" "YES" "YES" "YES"
## [1009] "YES" "YES" "YES" "YES"
## [1013] "YES" "YES" "YES" "YES"
## [1017] "YES" "YES" "YES" "YES"
## [1021] "YES" "YES" "YES" "YES"
## [1025] "Sin Registro" "YES" "YES" "YES"
## [1029] "YES" "YES" "YES" "YES"
## [1033] "YES" "YES" "Sin Registro" "YES"
## [1037] "YES" "YES" "YES" "YES"
## [1041] "YES" "YES" "YES" "YES"
## [1045] "YES" "YES" "YES" "YES"
## [1049] "YES" "YES" "Sin Registro" "YES"
## [1053] "YES" "YES" "YES" "YES"
## [1057] "YES" "YES" "YES" "YES"
## [1061] "YES" "YES" "YES" "YES"
## [1065] "YES" "YES" "YES" "YES"
## [1069] "YES" "YES" "YES" "YES"
## [1073] "YES" "YES" "YES" "YES"
## [1077] "YES" "YES" "YES" "YES"
## [1081] "YES" "YES" "YES" "YES"
## [1085] "YES" "NO" "YES" "YES"
## [1089] "YES" "YES"Tabla de distribucion de frecuencia
2. TABLAS DE DISTRIBUCIÓN DE FRECUENCIA
#Crear tabla de frecuencia y Convertir tabla en dataframe
TDFDeformacion <- table(Deformacion)
TDFDeformacion <- as.data.frame(TDFDeformacion)
TDFDeformacion## Deformacion Freq
## 1 NO 33
## 2 Sin Registro 188
## 3 YES 869#calcular frecuencia absolutas y simples
TDFDeformacionFinal <- TDFDeformacion %>%
group_by(Deformacion) %>%
summarise(
ni = sum(Freq),
hi = round((sum(ni) / sum(TDFDeformacion$Freq)) * 100, 2)
)
TDFDeformacionFinal <- data.frame(TDFDeformacionFinal)
TDFDeformacionFinal## Deformacion ni hi
## 1 NO 33 3.03
## 2 Sin Registro 188 17.25
## 3 YES 869 79.72#Comprobación del tamaño
sum(TDFDeformacionFinal$ni)## [1] 1090sum(TDFDeformacionFinal$hi)## [1] 100#Agregar Total
Total_ni <- sum(TDFDeformacionFinal$ni)
Total_hi <- sum(TDFDeformacionFinal$hi)
TDFDeformacionFinalCompleto <- rbind(TDFDeformacionFinal,
data.frame(Deformacion="Total", ni = Total_ni, hi = Total_hi))
print(TDFDeformacionFinalCompleto)## Deformacion ni hi
## 1 NO 33 3.03
## 2 Sin Registro 188 17.25
## 3 YES 869 79.72
## 4 Total 1090 100.002.1. Tabla Nº1-Distribución de Deformaciones Post Mineralización en Depósito masivos de sulfuros volcanicos
#Crear Tabla
TablaDeformacion <- TDFDeformacionFinalCompleto %>%
gt(rowname_col = NULL) %>%
tab_header(
title = md("**Tabla Nº1**"),
subtitle = md("**Distribución de Deformaciones Post Mineralización
en Depósito masivos de sulfuros volcanicos**")
) %>%
tab_source_note(
source_note = md("Autor: Grupo2")
) %>%
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 = TDFDeformacionFinalCompleto$Deformacion == "Total"
)
)
TablaDeformacion| Tabla Nº1 | ||
| Distribución de Deformaciones Post Mineralización en Depósito masivos de sulfuros volcanicos | ||
| Deformacion | ni | hi |
|---|---|---|
| NO | 33 | 3.03 |
| Sin Registro | 188 | 17.25 |
| YES | 869 | 79.72 |
| Total | 1090 | 100.00 |
| Autor: Grupo2 | ||
Diagramas de barras
3. DIAGRAMAS DE BARRAS
3.1. DIAGRAMA DE BARRAS ni (LOCAL/GLOBAL)
#Diagrama de barras ni (local)
barplot(TDFDeformacionFinalCompleto$ni[1:(nrow(TDFDeformacionFinalCompleto)-1)],
main = "Gráfica Nº1: Distribucion de Deformaciones Post Mineralización
en Depósito masivos de sulfuros volcanicos (Local)",
cex.main=0.7,
col = "gray",
xlab = "Presencia de Deformaciones Post Mineralización",
ylab = "Cantidad (ni)",
names.arg = TDFDeformacionFinalCompleto$Deformacion[1:(nrow(TDFDeformacionFinalCompleto)-1)])
#Diagrama de barras ni (global)
barplot(
TDFDeformacionFinalCompleto$ni[1:(nrow(TDFDeformacionFinalCompleto)-1)],
main = "Gráfica Nº2: Distribucion de Deformaciones Post Mineralización en Depósito masivos de sulfuros volcanicos (Global)",
cex.main = 0.7,
xlab = "Presencia Deformaciones Post Mineralización",
ylab = "Cantidad (ni)",
col = "gray",
names.arg = TDFDeformacionFinalCompleto$Deformacion[1:(nrow(TDFDeformacionFinalCompleto)-1)],
ylim = c(0, 1090)
)
3.2. DIAGRAMA DE BARRAS hi (LOCAL/GLOBAL)
#Diagrama de barras hi (local)
hi_local <- as.numeric(TDFDeformacionFinalCompleto$hi[1:(nrow(TDFDeformacionFinalCompleto)-1)])
barplot(
hi_local,
main = "Gráfica Nº3: Distribucion de Deformaciones Post Mineralización en Depósito masivos de sulfuros volcanicos (Local)",
cex.main = 0.7,
xlab = "Presencia de Deformaciones Post Mineralización",
ylab = "Porcentaje (%)",
col = "gray",
names.arg = TDFDeformacionFinalCompleto$Deformacion[1:(nrow(TDFDeformacionFinalCompleto)-1)] ,
cex.names = 1,
ylim = c(0, 80)
)
#Diagrama de barras hi (global)
hi_local <- as.numeric(TDFDeformacionFinalCompleto$hi[1:(nrow(TDFDeformacionFinalCompleto)-1)])
barplot(
hi_local,
main = "Gráfica Nº3: Distribucion de Deformaciones Post Mineralización en Depósito masivos de sulfuros volcanicos (Local)",
cex.main = 0.7,
xlab = "Presencia de Deformaciones Post Mineralización",
ylab = "Porcentaje (%)",
col = "gray",
names.arg = TDFDeformacionFinalCompleto$Deformacion[1:(nrow(TDFDeformacionFinalCompleto)-1)] ,
cex.names = 1,
ylim = c(0, 100)
)
Diagrama circular
4. DIAGRAMA CIRCULAR
#Diagrama circular
hi_deformacion <- as.numeric(TDFDeformacionFinalCompleto$hi[1:(nrow(TDFDeformacionFinalCompleto)-1)])
deformacion <- TDFDeformacionFinalCompleto$Deformacion[1:(nrow(TDFDeformacionFinalCompleto)-1)]
Colores <- colorRampPalette(c("lightskyblue", "darkblue"))(length(hi_deformacion))
etiquetas <- paste0(hi_deformacion, "%")
par(mfrow = c(1,2))
par(mar = c(2,2,4,2))
pie(
hi_deformacion,
radius = 0.7,
col = Colores,
labels = etiquetas,
main = "Gráfica Nº5: Distribucion de Deformaciones
Post Mineralización
en Depósito masivos de sulfuros volcanicos",
cex.main = 0.9
)
plot.new()
legend(
"center",
title = "Deformaciones Post Mineralización",
legend = deformacion,
fill = Colores,
cex = 0.7,
bg = "white",
box.lwd = 0.7
)