#---------------------CARGA DE DATOS---------------------#

# 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 = ".")

#mostrar 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

#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
#Comprobación del tamaño
sum(TDFIntrusionPFinal$ni)
## [1] 1090
sum(TDFIntrusionPFinal$hi)
## [1] 99.99
#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°**"),
    subtitle = md("**Tabla de Distribucion de Frecuencias Absolutas y Relativas Simples
                  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°
Tabla de Distribucion de Frecuencias Absolutas y Relativas Simples 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