1. Carga de Datos 2.Tabla de distribución de frecuencia 3.Graficas de distribución de frecuencia

1 . CARGA DE DATOS

setwd("~/UNI/ESTADISTICA")
datos <- read.csv("Depositos_Sulfuro.csv", header = TRUE, sep = ";", dec = ".")

Verificamos que nos lea correctamente los datos

VETAS_S <-datos$stringer
str(VETAS_S)
##  chr [1:1090] "yes" "no" "yes" "yes" "yes" "no" "no" "yes" "" "no" "no" ...

2 . TABLA DE DISTRIBUCIÓN DE FRECUENCIA POR PRESENCIA DE VETAS SECUNDARIAS

#1 Extraer variable
Vetas_Secundarias <- datos$stringer

#2 Reemplazar NA, vacíos, espacios, guiones, etc.
Vetas_Secundarias[Vetas_Secundarias == "" |
                    Vetas_Secundarias == " " |
                    Vetas_Secundarias == "-" |
                    is.na(Vetas_Secundarias)] <- "Sin registro"

#3 Frecuencia absoluta
ni <- table(Vetas_Secundarias)

#4 Frecuencia relativa %
hi <- round(prop.table(ni) * 100,2)

#5 Crear tabla base
tabla_finalvetas <- data.frame(
  Vetas_Secundarias = names(ni),
  ni = as.numeric(ni),
  hi = as.numeric(hi)
)

#6 Crear fila total
fila_total <- data.frame(
  Vetas_Secundarias = "TOTAL",
  ni = sum(tabla_finalvetas$ni),
  hi = sum(tabla_finalvetas$hi)
)

#7 Unir tabla
tabla_finalvetas <- rbind(tabla_finalvetas, fila_total)

#8  Mostrar
tabla_finalvetas
##   Vetas_Secundarias   ni     hi
## 1                no  431  39.54
## 2      Sin registro  275  25.23
## 3               yes  384  35.23
## 4             TOTAL 1090 100.00

3 . Tabla N° 1 — Distribución por presencia de vetas secundarias de depositos masivos de sulfuro volcanicos

tabla_vetas_gt <- tabla_finalvetas %>%
  gt() %>%
  tab_header(
    title = md("**Tabla N° 1**"),
    subtitle = md("Distribución por presencia de vetas secundarias 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 = Vetas_Secundarias == "Total")
  )

tabla_vetas_gt
Tabla N° 1
Distribución por presencia de vetas secundarias en Depósitos
Masivos de Sulfuro Volcánicos
Vetas_Secundarias ni hi
no 431 39.54
Sin registro 275 25.23
yes 384 35.23
TOTAL 1090 100.00
Autor: Grupo 2

4 . DIAGRAMAS DE BARRAS Y DIAGRAMA CIRCULAR

4.1 . Diagrama de barras ni local

barplot(tabla_finalvetas$ni[1:(nrow(tabla_finalvetas)-1)],
        main = "Gráfica Nº1: Distribución por presencia de vetas secundarias en Depósito masivos de sulfuros volcanicos  (Local)",
        cex.main=0.7,
        col = "gray",
        xlab = "Presencia de vetas secundarias",
        ylab = "Cantidad (ni)",
        names.arg = tabla_finalvetas$Vetas_Secundarias[1:(nrow(tabla_finalvetas)-1)])

…

4.2 . Diagrama de barras ni global

barplot(
  tabla_finalvetas$ni[1:(nrow(tabla_finalvetas)-1)],
  main = "Gráfica Nº2: Distribución por presencia de vetas secundarias en Depósitos masivos de sulfuros volcanicos (Global)",
  cex.main = 0.7,
  xlab = "Presencia de vetas secundarias",
  ylab = "Cantidad (ni)",
  col = "gray",
  names.arg = tabla_finalvetas$Vetas_Secundarias[1:(nrow(tabla_finalvetas)-1)],
  ylim = c(0, 1090)
)

…

4.3 . Diagrama de barras hi local

hi_local <- as.numeric(tabla_finalvetas$hi[1:(nrow(tabla_finalvetas)-1)])
barplot(
  hi_local,
  main = "Gráfica Nº3: Distribución por presencia de vetas secundarias en Depósitos masivos de sulfuros volcanicos (Local)",
  cex.main = 0.7,
  xlab = "Presencia de vetas secundarias",
  ylab = "Porcentaje (%)",
  col = "gray",
  names.arg = tabla_finalvetas$Vetas_Secundarias[1:(nrow(tabla_finalvetas)-1)] ,
  cex.names = 1,
  ylim = c(0, 40)
)

…

4.4 . Diagrama de barras hi global

hi_global <- as.numeric(tabla_finalvetas$hi[1:(nrow(tabla_finalvetas)-1)])
barplot(
  hi_global,
  main = "Gráfica Nº4: Distribución por presencia de vetas secundarias en Depósitos masivos de sulfuros volcanicos (Global)",
  cex.main = 0.7,
  xlab = "Presencia de vetas secundarias",
  ylab = "Porcentaje (%)",
  col = "blue",
  names.arg = tabla_finalvetas$Vetas_Secundarias[1:(nrow(tabla_finalvetas)-1)],
  cex.names = 1,
  ylim = c(0, 100)
)

…

4.5 . Diagrama circular

hi_vetas <- as.numeric(tabla_finalvetas$hi[1:(nrow(tabla_finalvetas)-1)])
vetas <- tabla_finalvetas$Vetas_Secundarias[1:(nrow(tabla_finalvetas)-1)]

Colores <- colorRampPalette(c("lightskyblue", "darkblue"))(length(hi_vetas))

etiquetas <- paste0(hi_vetas, "%")

par(mfrow = c(1,2))     
par(mar = c(2,2,4,2))   

pie(
  hi_vetas,
  radius = 0.7,
  col = Colores,
  labels = etiquetas,
  main = "Gráfica Nº5: Distribución por presencia de vetas
  secundarias en depósitos masivos de sulfuros 
  volcánicos",
  cex.main = 0.9
)

plot.new()  
legend(
  "center",
  title = "Presencia de Vetas Secundarias",
  legend = vetas,
  fill = Colores,
  cex = 0.7,
  bg = "white",
  box.lwd = 0.7
)