setwd("~/UNI/ESTADISTICA")
datos <- read.csv("mineral_clasificado.csv", header = TRUE, sep = ";", dec = ".")
Verificamos que nos lea correctamente los datos
TIPO_MINERALES <-datos$mineral
str(TIPO_MINERALES)
## chr [1:13045] "arsenopyrite" "chalcopyrite" "chlorite" "sphalerite" ...
#1 Extraer la variable
mineral<-datos$mineral
#2 crear y visualizar la tabla de frecuencia
TDFmineral<-table(mineral)
Se agrupo por el tipo de mineral (metalico o no metalico) debido ya que se registraron 528 minerales, y la tabla de distribucion de frecuencia va hacer muy extensa
#1 extraer la variable
Tipo_mineral<-datos$tipo
#2 crear y visualizar la tabla de frecuencia
TDF_tipomineral<-table(Tipo_mineral)
#3 Reemplazar NA, vacíos, espacios, guiones, etc.
Tipo_mineral[Tipo_mineral == "" |
Tipo_mineral == " " |
Tipo_mineral == "-" |
is.na(Tipo_mineral)] <- "Sin registro"
#4 Frecuencia absoluta
ni <- table(Tipo_mineral)
#5 Frecuencia relativa %
hi <- round(prop.table(ni) * 100,2)
#6 Crear tabla base
tabla_final_tipo_mineral <- data.frame(
Tipo_mineral = names(ni),
ni = as.numeric(ni),
hi = as.numeric(hi)
)
#7 Crear fila total
fila_total <- data.frame(
Tipo_mineral = "TOTAL",
ni = sum(tabla_final_tipo_mineral$ni),
hi = sum(tabla_final_tipo_mineral$hi)
)
#8 Unir tabla
tabla_final_tipo_mineral<- rbind(tabla_final_tipo_mineral, fila_total)
#9 Mostrar tabla
tabla_final_tipo_mineral
## Tipo_mineral ni hi
## 1 Metálico 7471 57.27
## 2 No metálico 5574 42.73
## 3 TOTAL 13045 100.00
tabla_mineral_gt <- tabla_final_tipo_mineral %>%
gt() %>%
tab_header(
title = md("**Tabla N° 1**"),
subtitle = md("Distribución por tipo de mineral metalico o no metalico en Depósitos <br>
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 = Tipo_mineral == "Total")
)
tabla_mineral_gt
| Tabla N° 1 | ||
| Distribución por tipo de mineral metalico o no metalico en Depósitos Masivos de Sulfuro Volcánicos |
||
| Tipo_mineral | ni | hi |
|---|---|---|
| Metálico | 7471 | 57.27 |
| No metálico | 5574 | 42.73 |
| TOTAL | 13045 | 100.00 |
| Autor: Grupo 2 | ||
barplot(tabla_final_tipo_mineral$ni[1:(nrow(tabla_final_tipo_mineral)-1)],
main = "Gráfica Nº1: Distribución por tipo de mineral en Depósito masivos de
sulfuros volcanicos (Local)",
cex.main=0.8,
col = "gray",
xlab = "Tipo de mineral",
ylab = "Cantidad (ni)",
names.arg = tabla_final_tipo_mineral$Tipo_mineral[1:(nrow(tabla_final_tipo_mineral)-1)]
)
…
barplot(
tabla_final_tipo_mineral$ni[1:(nrow(tabla_final_tipo_mineral)-1)],
main = "Gráfica Nº2: Distribución por tipo de mineral en Depósito masivos de
sulfuros volcanicos (Global)",
cex.main = 0.8,
xlab = "Tipo de mineral",
ylab = "Cantidad (ni)",
col = "gray",
names.arg = tabla_final_tipo_mineral$Tipo_mineral[1:(nrow(tabla_final_tipo_mineral)-1)],
ylim = c(0, 13045)
)
…
hi_local <- as.numeric(tabla_final_tipo_mineral$hi[1:(nrow(tabla_final_tipo_mineral)-1)])
barplot(
hi_local,
main = "Gráfica Nº3: Distribución por tipo de mineralen Depósito masivos de
sulfuros volcanicos (Local)",
cex.main = 0.8,
xlab = "Tipo de mineral",
ylab = "Porcentaje (%)",
col = "gray",
names.arg = tabla_final_tipo_mineral$Tipo_mineral[1:(nrow(tabla_final_tipo_mineral)-1)] ,
cex.names = 1,
)
…
hi_global <- as.numeric(tabla_final_tipo_mineral$hi[1:(nrow(tabla_final_tipo_mineral)-1)])
barplot(
hi_global,
main = "Gráfica Nº4: Distribución por tipo de mineral en Depósito masivos de
sulfuros volcanicos (Global)",
cex.main = 0.8,
xlab = "Tipo de mineral",
ylab = "Porcentaje (%)",
col = "gray",
names.arg = tabla_final_tipo_mineral$Tipo_mineral[1:(nrow(tabla_final_tipo_mineral)-1)],
cex.names = 1,
ylim = c(0, 100)
)
…
hi_mineral<- as.numeric(tabla_final_tipo_mineral$hi[1:(nrow(tabla_final_tipo_mineral)-1)])
vetas <- tabla_final_tipo_mineral$Tipo_mineral[1:(nrow(tabla_final_tipo_mineral)-1)]
Colores <- colorRampPalette(c("lightskyblue", "darkblue"))(length(hi_mineral))
etiquetas <- paste0(hi_mineral, "%")
par(mfrow = c(1,2))
par(mar = c(2,2,4,2))
pie(
hi_mineral,
radius = 0.9,
col = Colores,
labels = etiquetas,
main = "Gráfica Nº5: Distribución por tipo de
mineral en Depósito masivos de
sulfuros volcanicos",
cex.main = 1.2
)
plot.new()
legend(
"center",
title = "Tipo de mineral",
legend = vetas,
fill = Colores,
cex = 0.9,
bg = "white",
box.lwd = 0.9
)