knitr::opts_chunk$set(message = FALSE, warning = FALSE)

Data quakes

Utilizaremos la infromación del dataset quakes

library(datasets)
head(quakes)
##      lat   long depth mag stations
## 1 -20.42 181.62   562 4.8       41
## 2 -20.62 181.03   650 4.2       15
## 3 -26.00 184.10    42 5.4       43
## 4 -17.97 181.66   626 4.1       19
## 5 -20.42 181.96   649 4.0       11
## 6 -19.68 184.31   195 4.0       12

Posición de sismos

library(leaflet)
leaflet() %>% 
  addTiles() %>% 
  addMarkers(lng = as.numeric(quakes[,2]),
             lat = as.numeric(quakes[,1]),
             clusterOptions = markerClusterOptions())

Clásificación de sismos por nivel de magnitud

quakes$Class <- cut(quakes$mag,breaks = c (0,4,5,6,7,10))

# Definición de paleta de colores
pal <- colorFactor(palette = c("green", "darkgreen", 'blue', "darkblue", "red"), 
                   levels = c("(0,4]", "(4,5]", "(5,6]", "(6,7]", "(7,10]"))
# Elaborando el mapa
library(dplyr)
quakes %>% leaflet() %>% 
  addTiles() %>% 
  addCircleMarkers(lng = as.numeric(quakes[,2]), lat =as.numeric(quakes[,1]), radius = 1,
                   color = ~pal(Class)) %>% 
  addLegend(position = "bottomright",
            pal = pal, 
            values = c("(0,4]", "(4,5]", "(5,6]", "(6,7]", "(7,10]"))