Análisis de Depósitos Masivos de Sulfuro Volcánicos
Grupo 2
2025-12-02
1 CARGA DE DATOS
setwd("~/UNI/ESTADISTICA")
datos <- read.csv("Depositos_Sulfuro.csv", header = TRUE, sep = ";", dec = ".")Verificamos que nos lea correctamente los datos
pais<-datos$country
str(pais)## chr [1:1090] "Oman" "Japan" "Sweden" "United States" "Cyprus" "Finland" ...2 TABLA DE DISTRIBUCIÓN DE FRECUENCIA POR PAÍS
TablaPais <- as.data.frame(table(datos$country))
colnames(TablaPais) <- c("Pais", "ni")
TablaPais$hi <- round((TablaPais$ni / sum(TablaPais$ni)) * 100, 2)
total_pais <- data.frame(
Pais = "Total",
ni = sum(TablaPais$ni),
hi = format(round(sum(TablaPais$hi),), nsmall = 2)
)
TablaPaisFinal <- rbind(TablaPais, total_pais)
TablaPaisFinal## Pais ni hi
## 1 Argentina 2 0.18
## 2 Armenia 4 0.37
## 3 Australia 56 5.14
## 4 Bolivia 1 0.09
## 5 Brazil 3 0.28
## 6 Canada 317 29.08
## 7 Chile 2 0.18
## 8 China 37 3.39
## 9 Colombia 6 0.55
## 10 Cuba 13 1.19
## 11 Cyprus 18 1.65
## 12 Dominican Republic 3 0.28
## 13 Ecuador 2 0.18
## 14 Egypt 1 0.09
## 15 Eritrea 3 0.28
## 16 Fiji 3 0.28
## 17 Finland 12 1.1
## 18 France 4 0.37
## 19 Georgia 4 0.37
## 20 Great Britain 1 0.09
## 21 Guatemala 1 0.09
## 22 Guyana 2 0.18
## 23 India 2 0.18
## 24 Indonesia 4 0.37
## 25 Iran 2 0.18
## 26 Ireland 1 0.09
## 27 Japan 82 7.52
## 28 Kazakhstan 46 4.22
## 29 Mexico 16 1.47
## 30 Mongolia 1 0.09
## 31 Morocco 3 0.28
## 32 Norway 49 4.5
## 33 Oman 3 0.28
## 34 Pakistan 1 0.09
## 35 Peru 5 0.46
## 36 Philippines 19 1.74
## 37 Portugal 14 1.28
## 38 Russia 90 8.26
## 39 Saudi Arabia 21 1.93
## 40 Spain 61 5.6
## 41 Sweden 40 3.67
## 42 Turkey 26 2.39
## 43 Union of Myanmar 1 0.09
## 44 United States 100 9.17
## 45 Uzbekistan 5 0.46
## 46 Venezuela 3 0.28
## 47 Total 1090 100.002.1 Tabla N° 1 — Distribución por País de los depositos masivos de sulfuro volcanicos
tabla_pais_gt <- TablaPaisFinal %>%
gt() %>%
tab_header(
title = md("**Tabla N° 1**"),
subtitle = md("Distribución por País de los 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 = Pais == "Total")
)
tabla_pais_gt| Tabla N° 1 | ||
| Distribución por País de los Depósitos Masivos de Sulfuro Volcánicos |
||
| Pais | ni | hi |
|---|---|---|
| Argentina | 2 | 0.18 |
| Armenia | 4 | 0.37 |
| Australia | 56 | 5.14 |
| Bolivia | 1 | 0.09 |
| Brazil | 3 | 0.28 |
| Canada | 317 | 29.08 |
| Chile | 2 | 0.18 |
| China | 37 | 3.39 |
| Colombia | 6 | 0.55 |
| Cuba | 13 | 1.19 |
| Cyprus | 18 | 1.65 |
| Dominican Republic | 3 | 0.28 |
| Ecuador | 2 | 0.18 |
| Egypt | 1 | 0.09 |
| Eritrea | 3 | 0.28 |
| Fiji | 3 | 0.28 |
| Finland | 12 | 1.1 |
| France | 4 | 0.37 |
| Georgia | 4 | 0.37 |
| Great Britain | 1 | 0.09 |
| Guatemala | 1 | 0.09 |
| Guyana | 2 | 0.18 |
| India | 2 | 0.18 |
| Indonesia | 4 | 0.37 |
| Iran | 2 | 0.18 |
| Ireland | 1 | 0.09 |
| Japan | 82 | 7.52 |
| Kazakhstan | 46 | 4.22 |
| Mexico | 16 | 1.47 |
| Mongolia | 1 | 0.09 |
| Morocco | 3 | 0.28 |
| Norway | 49 | 4.5 |
| Oman | 3 | 0.28 |
| Pakistan | 1 | 0.09 |
| Peru | 5 | 0.46 |
| Philippines | 19 | 1.74 |
| Portugal | 14 | 1.28 |
| Russia | 90 | 8.26 |
| Saudi Arabia | 21 | 1.93 |
| Spain | 61 | 5.6 |
| Sweden | 40 | 3.67 |
| Turkey | 26 | 2.39 |
| Union of Myanmar | 1 | 0.09 |
| United States | 100 | 9.17 |
| Uzbekistan | 5 | 0.46 |
| Venezuela | 3 | 0.28 |
| Total | 1090 | 100.00 |
| Autor: Grupo 2 | ||
Debido a que se registró una gran cantidad de países, se decidió agrupar por continentes.
3 TABLA DE DISTRIBUCIÓN DE FRECUENCIA POR CONTINENTE
TablaPais$Continente <- countrycode(TablaPais$Pais, origin = "country.name", destination = "continent")
TablaContinente <- aggregate(ni ~ Continente, data = TablaPais, sum)
TablaContinente$hi <- round((TablaContinente$ni / sum(TablaContinente$ni)) * 100, 2)
total_continente <- data.frame(
Continente = "Total",
ni = sum(TablaContinente$ni),
hi = format(round(sum(TablaContinente$hi),), nsmall = 2)
)
TablaContinenteFinal <- rbind(TablaContinente, total_continente)
TablaContinenteFinal## Continente ni hi
## 1 Africa 7 0.64
## 2 Americas 476 43.67
## 3 Asia 276 25.32
## 4 Europe 272 24.95
## 5 Oceania 59 5.41
## 6 Total 1090 100.003.1 Tabla N° 2 — Distribución por continente de los depositos masivos de sulfuros volcanicos
tabla_continente_gt <- TablaContinenteFinal %>%
gt() %>%
tab_header(
title = md("**Tabla N° 2**"),
subtitle = md("Distribución por continente de Depósitos <br>
Masivos de Sulfuros 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 = Continente == "Total")
)
tabla_continente_gt| Tabla N° 2 | ||
| Distribución por continente de Depósitos Masivos de Sulfuros Volcánicos |
||
| Continente | ni | hi |
|---|---|---|
| Africa | 7 | 0.64 |
| Americas | 476 | 43.67 |
| Asia | 276 | 25.32 |
| Europe | 272 | 24.95 |
| Oceania | 59 | 5.41 |
| Total | 1090 | 100.00 |
| Autor: Grupo 2 | ||
4 DIAGRAMAS DE BARRAS Y DIAGRAMA CIRCULAR
4.1 Diagrama de barras ni local
barplot(
TablaContinenteFinal$ni[1:(nrow(TablaContinenteFinal)-1)],
main = "Gráfica Nº1: Distribución por Continente de Depósitos masivos de sulfuros volcanicos (Local)",
cex.main = 0.8,
xlab = "Continente",
ylab = "Cantidad (ni)",
col = "blue",
names.arg = TablaContinenteFinal$Continente[1:(nrow(TablaContinenteFinal)-1)],
ylim = c(0, max(TablaContinenteFinal$ni[1:(nrow(TablaContinenteFinal)-1)]))
)
…
4.2 Diagrama de barras ni global
barplot(
TablaContinenteFinal$ni[1:(nrow(TablaContinenteFinal)-1)],
main = "Gráfica Nº2: Distribución por Continente de Depósitos masivos de sulfuros volcanicos (Global)",
cex.main = 0.8,
xlab = "Continente",
ylab = "Cantidad (ni)",
col = "blue",
names.arg = TablaContinenteFinal$Continente[1:(nrow(TablaContinenteFinal)-1)],
ylim = c(0, sum(TablaContinenteFinal$ni[1:(nrow(TablaContinenteFinal)-1)]))
)
…
4.3 Diagrama de barras hi local
hi_local <- as.numeric(TablaContinenteFinal$hi[1:(nrow(TablaContinenteFinal)-1)])
continentes <- TablaContinenteFinal$Continente[1:(nrow(TablaContinenteFinal)-1)]
barplot(
hi_local,
main = "Gráfica Nº3:Distribución por Continente de Depósitos masivos de sulfuros volcanicos (Local)",
cex.main = 0.8,
xlab = "Continente",
ylab = "Porcentaje (%)",
col = "blue",
names.arg = continentes,
cex.names = 1,
ylim = c(0, max(hi_local))
)
…
4.4 Diagrama de barras hi global
hi_global <- as.numeric(TablaContinenteFinal$hi[1:(nrow(TablaContinenteFinal)-1)])
continentes <- TablaContinenteFinal$Continente[1:(nrow(TablaContinenteFinal)-1)]
barplot(
hi_global,
main = "Gráfica Nº4: Distribución por Continente de Depósitos masivos de sulfuros volcanicos (Global)",
cex.main = 0.8,
xlab = "Continente",
ylab = "Porcentaje (%)",
col = "blue",
names.arg = continentes,
cex.names = 1,
ylim = c(0, 100)
)
…
4.5 Diagrama circular
hi_continente <- as.numeric(TablaContinenteFinal$hi[1:(nrow(TablaContinenteFinal)-1)])
continentes <- TablaContinenteFinal$Continente[1:(nrow(TablaContinenteFinal)-1)]
Colores <- colorRampPalette(c("lightskyblue", "darkblue"))
etiquetas <- paste0(round(hi_continente), "%")
pie(
hi_continente,
radius = 0.8,
col = Colores(length(hi_continente)),
labels = etiquetas,
main = "Gráfica Nº5: Distribución por Continente de los depositos masivos de sulfuros volcanicos",
cex.main = 0.8
)
# LEYENDA
legend(
"bottomright",
title = "Continente",
legend = continentes,
fill = Colores(length(hi_continente)),
cex = 0.9
)