CARGA DE DATOS
#Carga de datos
setwd("~/UNI/ESTADISTICA")
datos <- read.csv("Depositos_Sulfuro.csv", header = TRUE, sep = ";", dec = ".")
datos2 <- read.csv2("Clasificacion depage.csv", header = TRUE, sep = ";", dec = ".")
CARGA DE LIBRERIAS
#Carga de librerias
library(gt)
library(dplyr)
library(knitr)
TABLA DE DISTRIBUCION DE FRECUENCIA POR EDAD GEOLOGICA
# Extraer variable
edad <- datos$depage
# Frecuencia absoluta
ni <- table(edad)
# Frecuencia relativa en %
hi <- prop.table(ni) * 100
# Tabla principal
tabla_final <- data.frame(
edad = names(ni),
ni = as.numeric(ni),
hi = as.numeric(hi)
)
Fila total de las sumas de ni y hi
# Fila TOTAL
fila_total <- data.frame(
edad = "TOTAL",
ni = sum(tabla_final$ni),
hi = sum(tabla_final$hi)
)
# Unir
tabla_final <- rbind(tabla_final, fila_total)
tabla_final
## edad ni hi
## 1 Upper Devonian (378.3\xb139 - 364\xb115) 1 0.09174312
## 2 Upper Triassic (217\xb128, Re-Os) 1 0.09174312
## 3 Miocene-Pliocene (14\xb12 - 2.2) 1 0.09174312
## 4 Archean 114 10.45871560
## 5 Archean-Mid Proterozoic 1 0.09174312
## 6 Archean-Proterozoic 1 0.09174312
## 7 Cambrian 31 2.84403670
## 8 Cambrian-Ordovician 13 1.19266055
## 9 Cambrian-Silurian 3 0.27522936
## 10 Cambrian(?) 2 0.18348624
## 11 Carboniferous 35 3.21100917
## 12 Cretaceous 27 2.47706422
## 13 Cretaceous-Eocene 1 0.09174312
## 14 Cretaceous or Paleocene-Eocene 1 0.09174312
## 15 Devonian 27 2.47706422
## 16 Devonian-Carboniferous 6 0.55045872
## 17 Devonian-Mississippian 7 0.64220183
## 18 Devonian-Pennsylvanian 1 0.09174312
## 19 Devonian-Permian 1 0.09174312
## 20 Early- Middle Devonian (Emsian-Eifelian) 1 0.09174312
## 21 Early-Middle Devonian 2 0.18348624
## 22 Early-Middle Devonian (Emsian-Eifelian) 1 0.09174312
## 23 Early-Middle Devonian(Emsian-Eifelian) 1 0.09174312
## 24 Early Cambrian 9 0.82568807
## 25 Early Carboniferous 1 0.09174312
## 26 Early Carboniferous (Visean-Namurian) 4 0.36697248
## 27 Early Carboniferous (Visean-Namurian?) 1 0.09174312
## 28 Early Cretaceous 6 0.55045872
## 29 Early Devonian 5 0.45871560
## 30 Early Devonian (Emsian) 8 0.73394495
## 31 Early Ordovician 1 0.09174312
## 32 Early Proterozoic 1 0.09174312
## 33 Early Silurian 1 0.09174312
## 34 Early Silurian (Llandovery) 2 0.18348624
## 35 Early Silurian(Llandovery) 1 0.09174312
## 36 Early to Middle Cretaceous 1 0.09174312
## 37 Eocene 8 0.73394495
## 38 Eocene-Middle Miocene 1 0.09174312
## 39 Eocene-Pliocene 8 0.73394495
## 40 Jurassic 22 2.01834862
## 41 Jurassic-Early Cretaceous 1 0.09174312
## 42 Jurassic-Lower Cretaceous 4 0.36697248
## 43 Jurassic or Cretaceous 1 0.09174312
## 44 Late-Middle Triassic 1 0.09174312
## 45 Late Ordovician 1 0.09174312
## 46 Late Archean 27 2.47706422
## 47 Late Cambrian-Early Ordovician 1 0.09174312
## 48 Late Carboniferous 2 0.18348624
## 49 Late Cretaceous 1 0.09174312
## 50 Late Cretaceous (Turonian) 4 0.36697248
## 51 Late Devonian 3 0.27522936
## 52 Late Devonian (Frasnian) 3 0.27522936
## 53 Late Jurassic 1 0.09174312
## 54 Late Jurassic-Lower Cretaceous 4 0.36697248
## 55 Late Ordovician 6 0.55045872
## 56 Late Permian 2 0.18348624
## 57 Late Permian-Early Triassic 1 0.09174312
## 58 Late Precambrian 1 0.09174312
## 59 Late Proterozoic 1 0.09174312
## 60 Late Silurian 1 0.09174312
## 61 Late Silurian-Early Devonian 2 0.18348624
## 62 Late Silurian-Lower Devonian 1 0.09174312
## 63 Late Triassic 1 0.09174312
## 64 Late Triassic-Late Jurassic 1 0.09174312
## 65 Lower-Middle Carboniferous 1 0.09174312
## 66 Lower-Middle Devonian (Emsian-early Eifelian) 1 0.09174312
## 67 Lower-Middle Devonian (Emsian-Eifelian) 2 0.18348624
## 68 Lower-Middle Jurassic 1 0.09174312
## 69 Lower Cretaceous 14 1.28440367
## 70 Lower Devonian 1 0.09174312
## 71 Lower Devonian (407, U-Pb) 1 0.09174312
## 72 Lower Devonian (Emsian) 2 0.18348624
## 73 Lower Devonian (late Emsian) 2 0.18348624
## 74 Lower Jurassic 2 0.18348624
## 75 Lower or mid-Cretaceous 1 0.09174312
## 76 Lower Ordovician 7 0.64220183
## 77 Lower Paleocene-Upper Eocene 2 0.18348624
## 78 Lower Permian 1 0.09174312
## 79 Lower Permian-Late Carboniferous (292-283) 1 0.09174312
## 80 Lower Proterozoic 16 1.46788991
## 81 Lower Tertiary 5 0.45871560
## 82 Meso-Neoproterozoic 1 0.09174312
## 83 Mesoproterozoic 3 0.27522936
## 84 Mesoproterozoic (1216-1026, Sm-Nd) 1 0.09174312
## 85 Mesoproterozoic? (Riphean?) 1 0.09174312
## 86 Mesozoic 2 0.18348624
## 87 Mid Cretaceous 1 0.09174312
## 88 Middle-Late Devonian (Givetian-Frasnian) 1 0.09174312
## 89 Middle-Upper Devonian 3 0.27522936
## 90 Middle-Upper Devonian (Givetian-Frasnian) 6 0.55045872
## 91 Middle-Upper Jurassic 1 0.09174312
## 92 Middle Cambrian 5 0.45871560
## 93 Middle Devonian 23 2.11009174
## 94 Middle Devonian (Eifelian-Givetian) 12 1.10091743
## 95 Middle Devonian (Eifelian to early Givetian) 1 0.09174312
## 96 Middle Devonian (Eifelian) 13 1.19266055
## 97 Middle Devonian (Emsian-Eifelian) 1 0.09174312
## 98 Middle Devonian (Givetian) 4 0.36697248
## 99 Middle Devonian (Late Givetian) 1 0.09174312
## 100 Middle Devonian (lower Givetian) 1 0.09174312
## 101 Middle Devonian(Givetian) 2 0.18348624
## 102 Middle Jurassic (late Bajocian) 3 0.27522936
## 103 Middle Miocene 1 0.09174312
## 104 Middle Ordovician 2 0.18348624
## 105 Middle Proterozoic 2 0.18348624
## 106 Middle Proterozoic? 1 0.09174312
## 107 Miocene 39 3.57798165
## 108 Miocene? 1 0.09174312
## 109 Mississippian 4 0.36697248
## 110 Neoarchean-Paleoproterozoic (2.8-2.6 Ga) 1 0.09174312
## 111 Neoproterozoic 3 0.27522936
## 112 Neoproterozoic (976.4-802.3) 1 0.09174312
## 113 Neoproterozoic (Early Paleozoic?) 1 0.09174312
## 114 Neoproterozoic, Sinian 1 0.09174312
## 115 Oligocene 1 0.09174312
## 116 Ordovician 108 9.90825688
## 117 Ordovician-Silurian 5 0.45871560
## 118 Ordovician to Precambrian 1 0.09174312
## 119 Orodivician 1 0.09174312
## 120 Paleoproterozoic 2 0.18348624
## 121 Paleoproterozoic (1700-2000) 2 0.18348624
## 122 Paleozoic 29 2.66055046
## 123 Paleozoic-Mid Mesozoic 1 0.09174312
## 124 Paleozoic-Triassic 1 0.09174312
## 125 Pennsylvanian-Permian 1 0.09174312
## 126 Permian 9 0.82568807
## 127 Permian-Cretaceous 1 0.09174312
## 128 Permian-Triassic 2 0.18348624
## 129 Pliocene 1 0.09174312
## 130 pre-Eocene 1 0.09174312
## 131 Pre-Tertiary 1 0.09174312
## 132 Precambrian 2 0.18348624
## 133 Precambrian-Late Cambrian 2 0.18348624
## 134 Precambrian Z (570-800 Ma) 1 0.09174312
## 135 Proterozoic 154 14.12844037
## 136 Proterozoic or Ordovician 2 0.18348624
## 137 Silurian 39 3.57798165
## 138 Silurian-Devonian 2 0.18348624
## 139 Silurian to Carboniferous 1 0.09174312
## 140 Sin registro 3 0.27522936
## 141 Tertiary 1 0.09174312
## 142 Triassic 5 0.45871560
## 143 Triassic-Jurassic 1 0.09174312
## 144 Triassic or Jurassic 1 0.09174312
## 145 Upper-Devonian-Lower Carboniferous 1 0.09174312
## 146 Upper Carboniferous 1 0.09174312
## 147 Upper Cretaceous 29 2.66055046
## 148 Upper Devonian 2 0.18348624
## 149 Upper Devonian- Lower Carboniferous 1 0.09174312
## 150 Upper Devonian-Lower Carboniferou 2 0.18348624
## 151 Upper Devonian-Lower Carboniferous 27 2.47706422
## 152 Upper Devonian (Famennian) 2 0.18348624
## 153 Upper Devonian (Frasnian) 1 0.09174312
## 154 Upper Jurassic 1 0.09174312
## 155 Upper Jurassic-Lower Cretaceous 4 0.36697248
## 156 Upper Proterozoic 2 0.18348624
## 157 Upper Silurian 1 0.09174312
## 158 Upper Triassic 1 0.09174312
## 159 Upper Devonian (373\xb115) 1 0.09174312
## 160 TOTAL 1090 100.00000000
Debido a que la tabla presenta numerosos registros de edades geológicas, se decidió agruparlos por era geológica, pasando de una variable cualitativa nominal a una variable cualitativa ordinal
TABLA DE DISTRIBUCION DE FRECUENCIA AGRUPADA
# Extraer variable
Era <- datos2$Classificacion
# Corrección desde los datos (por si viene mal escrito)
Era <- gsub("precámbico", "precambrico", Era, ignore.case = TRUE)
Era <- gsub("Precámbico", "precambrico", Era, ignore.case = TRUE)
# Frecuencia absoluta
ni <- table(Era)
# Frecuencia relativa porcentual
hi <- round(prop.table(ni) * 100, 2)
# Crear tabla
tabla_era <- data.frame(
Era = names(ni),
ni = as.numeric(ni),
hi = as.numeric(hi)
)
# Corrección definitiva dentro de la tabla
tabla_era$Era <- gsub("precámbico", "Precambrico", tabla_era$Era, ignore.case = TRUE)
tabla_era$Era <- gsub("Precámbico", "Precambrico", tabla_era$Era, ignore.case = TRUE)
Fila total de las sumas de ni y hi
# Fila TOTAL
fila_total2 <- data.frame(
Era = "TOTAL",
ni = sum(tabla_era$ni),
hi = round(sum(tabla_era$hi), )
)
# Unir tablas
tabla_era <- rbind(tabla_era, fila_total2)
tabla_era
## Era ni hi
## 1 Cenozoico 73 6.70
## 2 Mesozoico 149 13.67
## 3 Paleozoico 522 47.89
## 4 Precámbrico 343 31.47
## 5 Sin registro 3 0.28
## 6 TOTAL 1090 100.00
TABLA DE DISTRIBUCION DE FRECUENCIA AGRUPADA FINAL
tabla_era_gt <- tabla_era %>%
gt() %>%
tab_header(
title = md("**Tabla N° 1**"),
subtitle = md("Distribución de frecuencia simples por era geologica en <br>
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
)
tabla_era_gt
| Tabla N° 1 | ||
| Distribución de frecuencia simples por era geologica en depósitos Masivos de Sulfuro Volcánicos |
||
| Era | ni | hi |
|---|---|---|
| Cenozoico | 73 | 6.70 |
| Mesozoico | 149 | 13.67 |
| Paleozoico | 522 | 47.89 |
| Precámbrico | 343 | 31.47 |
| Sin registro | 3 | 0.28 |
| TOTAL | 1090 | 100.00 |
| Autor: Grupo 2 | ||
Diagrama de barras de frecuencia absoluta local
barplot(
tabla_era$ni[1:(nrow(tabla_era)-1)],
main = "Gráfica Nº1: Distribución de frecuencia absoluta local de era
geologica de los depositos masivos de sulfuros volcanicos ",
cex.main = 0.8,
col = "gray",
xlab = "Era geologica",
ylab = "Cantidad (ni)",
names.arg = tabla_era$Era[1:(nrow(tabla_era)-1)]
)
Diagrama de barras de frecuencia absoluta global
barplot(
tabla_era$ni[1:(nrow(tabla_era)-1)],
main = "Gráfica Nº2: Distribución de frecuencia absoluta global de era
geologica de los depositos masivos de sulfuros volcanicos",
cex.main = 0.8,
xlab = "Era geologica",
ylab = "Cantidad (ni)",
col = "gray",
names.arg = tabla_era$Era[1:(nrow(tabla_era)-1)],
ylim = c(0, 1090)
)
Diagrama de barras de frecuencia relativa local
hi_local <- as.numeric(tabla_era$hi[1:(nrow(tabla_era)-1)])
barplot(
hi_local,
main = "Gráfica Nº3: Distribución de frecuencia relativa local de era
geologica de los depositos masivos de sulfuros volcanicos",
cex.main = 0.8,
xlab = "Era geologica",
ylab = "Porcentaje (%)",
col = "gray",
names.arg = tabla_era$Era[1:(nrow(tabla_era)-1)]
)
Diagrama de barras de frecuencia relativa global
hi_global <- as.numeric(tabla_era$hi[1:(nrow(tabla_era)-1)])
barplot(
hi_global,
main = "Gráfica Nº4: Distribución de frecuencia relativa global de era
geologica de los depositos masivos de sulfuros volcanicos",
cex.main = 0.8,
xlab = "Era geologica",
ylab = "Porcentaje (%)",
col = "gray",
names.arg = tabla_era$Era[1:(nrow(tabla_era)-1)],
ylim = c(0, 100)
)
DIAGRAMA CIRCULAR
hi_era <- as.numeric(tabla_era$hi[1:(nrow(tabla_era)-1)])
eras <- tabla_era$Era[1:(nrow(tabla_era)-1)]
Colores <- colorRampPalette(c("lightskyblue", "darkblue"))(length(hi_era))
etiquetas <- paste0(hi_era, "%")
par(mfrow = c(1,2))
par(mar = c(2,2,4,2))
pie(
hi_era,
radius = 0.8,
col = Colores,
labels = etiquetas,
main = "Gráfica Nº5: Distribución de frecuencias
relativas de era geologica de los depositos masivos
de sulfuros volcanicos",
cex.main = 0.8
)
plot.new()
legend(
"center",
title = "Era geologica",
legend = eras,
fill = Colores,
cex = 0.7,
bg = "white",
box.lwd = 0.7
)
Indicadores Estadisticos
POSICION MODA
# Asegurar que ERA sea vector
ERA_limpia <- as.character(Era)
ERA_limpia <- ERA_limpia[!is.na(ERA_limpia)]
# Tabla de frecuencias
tabla_era <- table(ERA_limpia)
# Moda
moda_era <- names(tabla_era)[which.max(as.numeric(tabla_era))]
moda_era
## [1] "Paleozoico"
TABLA DE INDICADORES ESTADISTICOS
Variable<-c("Era Geologica")
TablaIndicadores<-data.frame(Variable,moda_era)
colnames(TablaIndicadores)<-c("Variable","Moda")
kable(TablaIndicadores, format = "markdown", caption = "Tabla N°3. Indicadores estadíticos de la variable era geologica de los depositos masivos de sulfuros volcanicos")
| Variable | Moda |
|---|---|
| Era Geologica | Paleozoico |
CONCLUSÍON
La variable Era Geológica presenta registros que fluctúan entre el Precámbrico y el Cenozoico, concentrándose principalmente en la era Paleozoica. Esta mayor presencia resulta beneficiosa para el análisis minero, ya que dicho periodo geológico está asociado a la formación de numerosos yacimientos minerales, lo que facilita la evaluación del potencial geológico y la orientación de las actividades de exploración.