# 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, union
library(ggplot2)
library(gt)
#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$postminint)
## chr [1:1090] "yes" "no" "" "" "no" "no" "yes" "yes" "" "yes" "yes" "" ...
Limpiar, corregir y estandarizar los valores
#Limpiar, corregir y estandarizar los valores de los datos
datos <- datos %>%
mutate(
postminint = trimws(postminint),
postminint = case_when(
postminint %in% c("yes") ~ "YES",
postminint %in% c("no") ~ "NO",
postminint == "" | is.na(postminint) ~ "Sin Registro",
TRUE ~ postminint
)
)
#cargar variable
datos$postminint <- trimws(datos$postminint)
IntrusionP <- datos$postminint
IntrusionP
## [1] "YES" "NO" "Sin Registro" "Sin Registro"
## [5] "NO" "NO" "YES" "YES"
## [9] "Sin Registro" "YES" "YES" "Sin Registro"
## [13] "YES" "Sin Registro" "NO" "YES"
## [17] "Sin Registro" "YES" "YES" "NO"
## [21] "Sin Registro" "YES" "Sin Registro" "YES"
## [25] "YES" "NO" "Sin Registro" "Sin Registro"
## [29] "Sin Registro" "YES" "YES" "YES"
## [33] "YES" "YES" "YES" "YES"
## [37] "Sin Registro" "YES" "YES" "Sin Registro"
## [41] "NO" "Sin Registro" "Sin Registro" "YES"
## [45] "YES" "NO" "YES" "NO"
## [49] "YES" "YES" "YES" "Sin Registro"
## [53] "Sin Registro" "Sin Registro" "Sin Registro" "Sin Registro"
## [57] "Sin Registro" "Sin Registro" "YES" "Sin Registro"
## [61] "YES" "Sin Registro" "YES" "YES"
## [65] "YES" "NO" "NO" "YES"
## [69] "NO" "Sin Registro" "Sin Registro" "YES"
## [73] "Sin Registro" "YES" "Sin Registro" "YES"
## [77] "YES" "YES" "Sin Registro" "Sin Registro"
## [81] "Sin Registro" "Sin Registro" "Sin Registro" "Sin Registro"
## [85] "YES" "YES" "YES" "YES"
## [89] "YES" "YES" "YES" "YES"
## [93] "NO" "NO" "YES" "YES"
## [97] "Sin Registro" "YES" "YES" "YES"
## [101] "YES" "YES" "NO" "Sin Registro"
## [105] "YES" "YES" "YES" "YES"
## [109] "Sin Registro" "YES" "YES" "YES"
## [113] "NO" "YES" "YES" "Sin Registro"
## [117] "YES" "YES" "Sin Registro" "YES"
## [121] "YES" "YES" "YES" "YES"
## [125] "YES" "Sin Registro" "YES" "Sin Registro"
## [129] "Sin Registro" "Sin Registro" "Sin Registro" "Sin Registro"
## [133] "YES" "YES" "YES" "Sin Registro"
## [137] "YES" "Sin Registro" "Sin Registro" "Sin Registro"
## [141] "YES" "YES" "YES" "YES"
## [145] "Sin Registro" "YES" "YES" "YES"
## [149] "Sin Registro" "YES" "Sin Registro" "YES"
## [153] "Sin Registro" "YES" "YES" "Sin Registro"
## [157] "YES" "Sin Registro" "YES" "YES"
## [161] "YES" "Sin Registro" "Sin Registro" "NO"
## [165] "Sin Registro" "Sin Registro" "YES" "YES"
## [169] "YES" "YES" "YES" "YES"
## [173] "YES" "Sin Registro" "YES" "Sin Registro"
## [177] "YES" "YES" "YES" "Sin Registro"
## [181] "NO" "NO" "Sin Registro" "YES"
## [185] "YES" "YES" "YES" "Sin Registro"
## [189] "YES" "YES" "YES" "Sin Registro"
## [193] "YES" "YES" "Sin Registro" "Sin Registro"
## [197] "YES" "Sin Registro" "YES" "YES"
## [201] "YES" "Sin Registro" "YES" "Sin Registro"
## [205] "YES" "Sin Registro" "YES" "NO"
## [209] "YES" "YES" "Sin Registro" "YES"
## [213] "NO" "Sin Registro" "Sin Registro" "Sin Registro"
## [217] "YES" "YES" "Sin Registro" "YES"
## [221] "NO" "YES" "Sin Registro" "Sin Registro"
## [225] "Sin Registro" "YES" "YES" "Sin Registro"
## [229] "Sin Registro" "YES" "YES" "YES"
## [233] "Sin Registro" "YES" "Sin Registro" "YES"
## [237] "YES" "YES" "YES" "YES"
## [241] "YES" "YES" "YES" "Sin Registro"
## [245] "YES" "NO" "NO" "NO"
## [249] "YES" "YES" "YES" "YES"
## [253] "Sin Registro" "YES" "YES" "YES"
## [257] "YES" "Sin Registro" "Sin Registro" "Sin Registro"
## [261] "YES" "YES" "YES" "NO"
## [265] "YES" "NO" "YES" "YES"
## [269] "Sin Registro" "YES" "Sin Registro" "YES"
## [273] "NO" "Sin Registro" "Sin Registro" "YES"
## [277] "YES" "Sin Registro" "Sin Registro" "YES"
## [281] "YES" "Sin Registro" "YES" "YES"
## [285] "NO" "Sin Registro" "NO" "YES"
## [289] "YES" "YES" "Sin Registro" "YES"
## [293] "Sin Registro" "YES" "YES" "Sin Registro"
## [297] "Sin Registro" "YES" "NO" "YES"
## [301] "YES" "YES" "NO" "YES"
## [305] "YES" "YES" "YES" "Sin Registro"
## [309] "NO" "Sin Registro" "YES" "Sin Registro"
## [313] "YES" "Sin Registro" "YES" "YES"
## [317] "YES" "YES" "NO" "YES"
## [321] "Sin Registro" "YES" "YES" "YES"
## [325] "YES" "Sin Registro" "YES" "Sin Registro"
## [329] "YES" "YES" "YES" "YES"
## [333] "Sin Registro" "YES" "YES" "YES"
## [337] "NO" "Sin Registro" "YES" "NO"
## [341] "NO" "Sin Registro" "NO" "YES"
## [345] "Sin Registro" "Sin Registro" "YES" "YES"
## [349] "Sin Registro" "YES" "YES" "Sin Registro"
## [353] "YES" "YES" "YES" "YES"
## [357] "YES" "YES" "NO" "YES"
## [361] "YES" "Sin Registro" "YES" "YES"
## [365] "Sin Registro" "Sin Registro" "YES" "YES"
## [369] "Sin Registro" "Sin Registro" "Sin Registro" "YES"
## [373] "NO" "NO" "YES" "Sin Registro"
## [377] "Sin Registro" "NO" "YES" "Sin Registro"
## [381] "Sin Registro" "Sin Registro" "YES" "Sin Registro"
## [385] "Sin Registro" "Sin Registro" "NO" "YES"
## [389] "YES" "YES" "NO" "YES"
## [393] "YES" "NO" "YES" "YES"
## [397] "Sin Registro" "Sin Registro" "YES" "NO"
## [401] "Sin Registro" "YES" "Sin Registro" "NO"
## [405] "YES" "NO" "YES" "Sin Registro"
## [409] "YES" "Sin Registro" "NO" "YES"
## [413] "YES" "Sin Registro" "YES" "YES"
## [417] "YES" "YES" "YES" "Sin Registro"
## [421] "YES" "Sin Registro" "Sin Registro" "YES"
## [425] "Sin Registro" "NO" "YES" "YES"
## [429] "YES" "YES" "YES" "YES"
## [433] "Sin Registro" "YES" "YES" "YES"
## [437] "YES" "YES" "Sin Registro" "Sin Registro"
## [441] "YES" "YES" "YES" "Sin Registro"
## [445] "Sin Registro" "Sin Registro" "YES" "YES"
## [449] "Sin Registro" "NO" "Sin Registro" "YES"
## [453] "Sin Registro" "NO" "YES" "Sin Registro"
## [457] "YES" "YES" "Sin Registro" "Sin Registro"
## [461] "Sin Registro" "Sin Registro" "Sin Registro" "YES"
## [465] "YES" "NO" "YES" "Sin Registro"
## [469] "YES" "NO" "YES" "YES"
## [473] "YES" "YES" "YES" "YES"
## [477] "YES" "YES" "YES" "YES"
## [481] "NO" "Sin Registro" "Sin Registro" "Sin Registro"
## [485] "Sin Registro" "YES" "Sin Registro" "YES"
## [489] "Sin Registro" "Sin Registro" "YES" "YES"
## [493] "Sin Registro" "Sin Registro" "Sin Registro" "Sin Registro"
## [497] "Sin Registro" "YES" "Sin Registro" "YES"
## [501] "Sin Registro" "Sin Registro" "Sin Registro" "Sin Registro"
## [505] "Sin Registro" "Sin Registro" "YES" "Sin Registro"
## [509] "Sin Registro" "Sin Registro" "YES" "Sin Registro"
## [513] "Sin Registro" "Sin Registro" "Sin Registro" "Sin Registro"
## [517] "Sin Registro" "NO" "Sin Registro" "Sin Registro"
## [521] "Sin Registro" "Sin Registro" "Sin Registro" "Sin Registro"
## [525] "YES" "Sin Registro" "Sin Registro" "Sin Registro"
## [529] "Sin Registro" "Sin Registro" "YES" "Sin Registro"
## [533] "YES" "NO" "YES" "YES"
## [537] "Sin Registro" "Sin Registro" "NO" "Sin Registro"
## [541] "Sin Registro" "NO" "NO" "NO"
## [545] "YES" "NO" "NO" "NO"
## [549] "NO" "NO" "NO" "YES"
## [553] "YES" "NO" "NO" "YES"
## [557] "NO" "NO" "NO" "YES"
## [561] "NO" "YES" "Sin Registro" "NO"
## [565] "Sin Registro" "Sin Registro" "Sin Registro" "NO"
## [569] "YES" "YES" "Sin Registro" "YES"
## [573] "YES" "YES" "YES" "YES"
## [577] "YES" "YES" "Sin Registro" "YES"
## [581] "YES" "Sin Registro" "Sin Registro" "YES"
## [585] "YES" "YES" "Sin Registro" "Sin Registro"
## [589] "Sin Registro" "Sin Registro" "Sin Registro" "Sin Registro"
## [593] "Sin Registro" "Sin Registro" "YES" "Sin Registro"
## [597] "Sin Registro" "Sin Registro" "YES" "YES"
## [601] "YES" "YES" "YES" "YES"
## [605] "YES" "NO" "YES" "YES"
## [609] "YES" "YES" "YES" "Sin Registro"
## [613] "YES" "Sin Registro" "NO" "Sin Registro"
## [617] "Sin Registro" "NO" "YES" "YES"
## [621] "YES" "Sin Registro" "Sin Registro" "Sin Registro"
## [625] "YES" "YES" "NO" "YES"
## [629] "YES" "Sin Registro" "YES" "YES"
## [633] "YES" "YES" "YES" "NO"
## [637] "NO" "YES" "NO" "YES"
## [641] "NO" "NO" "Sin Registro" "YES"
## [645] "YES" "YES" "NO" "YES"
## [649] "NO" "YES" "YES" "YES"
## [653] "YES" "NO" "NO" "NO"
## [657] "NO" "YES" "YES" "NO"
## [661] "Sin Registro" "Sin Registro" "NO" "NO"
## [665] "NO" "NO" "YES" "NO"
## [669] "NO" "YES" "YES" "NO"
## [673] "YES" "YES" "YES" "YES"
## [677] "Sin Registro" "YES" "Sin Registro" "Sin Registro"
## [681] "YES" "YES" "Sin Registro" "NO"
## [685] "YES" "YES" "NO" "NO"
## [689] "YES" "NO" "YES" "NO"
## [693] "Sin Registro" "YES" "NO" "YES"
## [697] "YES" "YES" "NO" "NO"
## [701] "Sin Registro" "NO" "YES" "YES"
## [705] "Sin Registro" "NO" "YES" "YES"
## [709] "YES" "YES" "YES" "YES"
## [713] "YES" "YES" "YES" "YES"
## [717] "YES" "YES" "YES" "Sin Registro"
## [721] "Sin Registro" "YES" "YES" "YES"
## [725] "Sin Registro" "YES" "YES" "YES"
## [729] "YES" "YES" "Sin Registro" "Sin Registro"
## [733] "YES" "YES" "YES" "YES"
## [737] "YES" "YES" "YES" "NO"
## [741] "Sin Registro" "YES" "YES" "Sin Registro"
## [745] "YES" "YES" "YES" "Sin Registro"
## [749] "YES" "YES" "YES" "YES"
## [753] "YES" "YES" "Sin Registro" "YES"
## [757] "YES" "YES" "YES" "NO"
## [761] "YES" "NO" "YES" "YES"
## [765] "YES" "YES" "YES" "YES"
## [769] "YES" "YES" "YES" "YES"
## [773] "NO" "Sin Registro" "YES" "YES"
## [777] "Sin Registro" "Sin Registro" "YES" "YES"
## [781] "YES" "Sin Registro" "YES" "YES"
## [785] "YES" "YES" "YES" "YES"
## [789] "YES" "YES" "NO" "YES"
## [793] "NO" "NO" "NO" "Sin Registro"
## [797] "YES" "YES" "YES" "YES"
## [801] "YES" "YES" "Sin Registro" "Sin Registro"
## [805] "Sin Registro" "Sin Registro" "Sin Registro" "Sin Registro"
## [809] "Sin Registro" "Sin Registro" "YES" "Sin Registro"
## [813] "YES" "YES" "YES" "Sin Registro"
## [817] "Sin Registro" "Sin Registro" "YES" "Sin Registro"
## [821] "Sin Registro" "YES" "YES" "Sin Registro"
## [825] "YES" "YES" "Sin Registro" "Sin Registro"
## [829] "Sin Registro" "Sin Registro" "Sin Registro" "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" "NO" "YES"
## [849] "YES" "YES" "Sin Registro" "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] "NO" "YES" "Sin Registro" "YES"
## [873] "Sin Registro" "YES" "YES" "YES"
## [877] "YES" "YES" "YES" "YES"
## [881] "YES" "Sin Registro" "Sin Registro" "YES"
## [885] "YES" "Sin Registro" "YES" "NO"
## [889] "YES" "Sin Registro" "Sin Registro" "YES"
## [893] "YES" "YES" "YES" "YES"
## [897] "YES" "Sin Registro" "YES" "YES"
## [901] "YES" "YES" "YES" "Sin Registro"
## [905] "YES" "YES" "YES" "Sin Registro"
## [909] "YES" "YES" "YES" "Sin Registro"
## [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" "Sin Registro" "YES" "YES"
## [949] "YES" "Sin Registro" "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" "YES" "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" "Sin Registro" "YES" "YES"
## [1029] "YES" "YES" "YES" "YES"
## [1033] "YES" "YES" "Sin Registro" "YES"
## [1037] "YES" "YES" "YES" "YES"
## [1041] "YES" "YES" "YES" "YES"
## [1045] "Sin Registro" "YES" "Sin Registro" "Sin Registro"
## [1049] "YES" "YES" "Sin Registro" "YES"
## [1053] "YES" "Sin Registro" "YES" "YES"
## [1057] "YES" "YES" "YES" "YES"
## [1061] "YES" "YES" "YES" "Sin Registro"
## [1065] "YES" "Sin Registro" "YES" "YES"
## [1069] "YES" "YES" "YES" "YES"
## [1073] "YES" "YES" "YES" "YES"
## [1077] "YES" "YES" "YES" "YES"
## [1081] "NO" "NO" "NO" "YES"
## [1085] "YES" "YES" "YES" "Sin Registro"
## [1089] "YES" "YES"
#Crear tabla de frecuencia y Convertir tabla en dataframe
TDFIntrusionP <- table(IntrusionP)
TDFIntrusionP <- as.data.frame(TDFIntrusionP)
TDFIntrusionP
## IntrusionP Freq
## 1 NO 122
## 2 Sin Registro 317
## 3 YES 651
#calcular frecuencia absolutas y simples
TDFIntrusionPFinal <- TDFIntrusionP %>%
group_by(IntrusionP) %>%
summarise(
ni = sum(Freq),
hi = round((sum(ni) / sum(TDFIntrusionP$Freq)) * 100, 2)
)
TDFIntrusionPFinal <- data.frame(TDFIntrusionPFinal)
TDFIntrusionPFinal
## IntrusionP ni hi
## 1 NO 122 11.19
## 2 Sin Registro 317 29.08
## 3 YES 651 59.72
#Agregar Total
Total_ni <- sum(TDFIntrusionPFinal$ni)
Total_hi <- 100
TDFIntrusionPFinalCompleto <- rbind(TDFIntrusionPFinal,
data.frame(IntrusionP="Total", ni = Total_ni, hi = Total_hi))
print(TDFIntrusionPFinalCompleto)
## IntrusionP ni hi
## 1 NO 122 11.19
## 2 Sin Registro 317 29.08
## 3 YES 651 59.72
## 4 Total 1090 100.00
#Crear Tabla
TablaIntrusionP <- TDFIntrusionPFinalCompleto %>%
gt(rowname_col = NULL) %>%
tab_header(
title = md("**Tabla Nº1**"),
subtitle = md("**Distribucion de las Intrusiones Post Mineralización**")
) %>%
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 = TDFIntrusionPFinalCompleto$Intrusion == "Total"
)
)
TablaIntrusionP
| Tabla Nº1 | ||
| Distribucion de las Intrusiones Post Mineralización | ||
| IntrusionP | ni | hi |
|---|---|---|
| NO | 122 | 11.19 |
| Sin Registro | 317 | 29.08 |
| YES | 651 | 59.72 |
| Total | 1090 | 100.00 |
| Autor: Grupo2 | ||
#Diagrama de barras ni (local)
barplot(TDFIntrusionPFinalCompleto$ni[1:(nrow(TDFIntrusionPFinalCompleto)-1)],
main = "Gráfica Nº1: Distribucion de las Intrusiones Post Mineralización en Depósito masivos de sulfuros volcanicos (Local)",
cex.main=0.7,
col = "gray",
xlab = "Presencia de Intrusiones Post Mineralización",
ylab = "Cantidad (ni)",
names.arg = TDFIntrusionPFinalCompleto$IntrusionP[1:(nrow(TDFIntrusionPFinalCompleto)-1)])
#Diagrama de barras ni (global)
barplot(
TDFIntrusionPFinalCompleto$ni[1:(nrow(TDFIntrusionPFinalCompleto)-1)],
main = "Gráfica Nº2: Distribucion de las Intrusiones Post Mineralización en Depósito masivos de sulfuros volcanicos (Global)",
cex.main = 0.7,
xlab = "Presencia de Intrusiones Post Mineralizacións",
ylab = "Cantidad (ni)",
col = "gray",
names.arg = TDFIntrusionPFinalCompleto$IntrusionP[1:(nrow(TDFIntrusionPFinalCompleto)-1)],
ylim = c(0, 1090)
)
#Diagrama de barras hi (local)
hi_local <- as.numeric(TDFIntrusionPFinalCompleto$hi[1:(nrow(TDFIntrusionPFinalCompleto)-1)])
barplot(
hi_local,
main = "Gráfica Nº3: Distribucion de las Intrusiones Post Mineralización en Depósito masivos de sulfuros volcanicos (Local)",
cex.main = 0.7,
xlab = "Presencia de Intrusiones Post Mineralización",
ylab = "Porcentaje (%)",
col = "gray",
names.arg = TDFIntrusionPFinalCompleto$IntrusionP[1:(nrow(TDFIntrusionPFinalCompleto)-1)] ,
cex.names = 1,
ylim = c(0, 60)
)
#Diagrama de barras hi (global)
hi_local <- as.numeric(TDFIntrusionPFinalCompleto$hi[1:(nrow(TDFIntrusionPFinalCompleto)-1)])
barplot(
hi_local,
main = "Gráfica Nº3: Distribucion de las Intrusiones Post Mineralización en Depósito masivos de sulfuros volcanicos (Local)",
cex.main = 0.7,
xlab = "Presencia de Intrusiones Post Mineralización",
ylab = "Porcentaje (%)",
col = "gray",
names.arg = TDFIntrusionPFinalCompleto$IntrusionP[1:(nrow(TDFIntrusionPFinalCompleto)-1)] ,
cex.names = 1,
ylim = c(0, 100)
)
#Diagrama circular
hi_instrusionp <- as.numeric(TDFIntrusionPFinalCompleto$hi[1:(nrow(TDFIntrusionPFinalCompleto)-1)])
Intrusiones <- TDFIntrusionPFinalCompleto$IntrusionP[1:(nrow(TDFIntrusionPFinalCompleto)-1)]
Colores <- colorRampPalette(c("lightskyblue", "darkblue"))(length(hi_instrusionp))
etiquetas <- paste0(hi_instrusionp, "%")
par(mfrow = c(1,2))
par(mar = c(2,2,4,2))
pie(
hi_instrusionp,
radius = 0.7,
col = Colores,
labels = etiquetas,
main = "Gráfica Nº5: Distribucion de las Intrusiones
Post Mineralización en Depósito
masivos de sulfuros volcanicos",
cex.main = 0.9
)
plot.new()
legend(
"center",
title = "Intrusiones Post Mineralización",
legend = Intrusiones,
fill = Colores,
cex = 0.7,
bg = "white",
box.lwd = 0.7
)