library(leaflet)
Mgraf <- readxl::read_excel("C:/Users/MELANY/OneDrive/Escritorio/Exmod/pregun/baseuno.xlsx","Hoja1")
Dat <- colorBin("magma", domain = Mgraf$Fe,bins = 4)

leaflet() %>% 
  addProviderTiles('OpenTopoMap') %>%
  
  addCircleMarkers(data = Mgraf,
                   lng = Mgraf$Longitud,
                   lat = Mgraf$Latitud,
                   color = ~Dat(Mgraf$Fe),
                   label = ~Mgraf$Fe) %>% 
  addLegend('topright',pal = Dat, values = Mgraf$Fe, title = 'parámetro Fe')