# cargar 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(ggplot2)
library(gt)
library(knitr)
#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"
#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] 1090
sum(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.00
#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 | ||
#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)
)
#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
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
)
#MODA VARIABLE MINERAL
# Asegurar que ERA sea vector
DEF_limpia <- as.character(Deformacion)
DEF_limpia <- DEF_limpia [!is.na(DEF_limpia)]
# Tabla de frecuencias
tabla_DEF <- table(DEF_limpia)
# Moda
moda_def <- names(tabla_DEF)[which.max(as.numeric(tabla_DEF))]
moda_def
## [1] "YES"
#MODA
Variable<-c("Deformación Post Mineralización")
TablaIndicadores<-data.frame(Variable,moda_def)
colnames(TablaIndicadores)<-c("Variable","Moda")
kable(TablaIndicadores, format = "markdown", caption = "Tabla N°2. Indicadores Estadíticos de la variable Deformación Post Mineralización en depósitos masivos de sulfuros volcanicos")
| Variable | Moda |
|---|---|
| Deformación Post Mineralización | YES |
La variable Deformación Post Mineralización indica que la moda es “YES”, lo que significa que la presencia de deformación post mineralización es la condición más frecuente en los depósitos masivos de sulfuros volcánicos analizados. Dado que se trata de una variable cualitativa nominal, la moda es el estadístico más representativo, permitiendo identificar una clara predominancia de eventos deformacionales posteriores a la mineralización dentro del conjunto de datos.