Read earthquake dataset into R and map it with animation

#1
Dear All,
I would like to read Turkey's eartquake data set into R and then plot with map with the animation according to the date(earlist to soonest
) that earthquake happened. Can you help me please?
Below is what I got so far but can't read the data.
Thank you.
Code:
install.packages('XML')
library(XML)
url <- "<http://www.koeri.boun.edu.tr/scripts/lst5.asp>"
data_df <- readHTMLTable(url)
 
#2
It isn't a table but fixed width data embedded in the page.

You can read it in using the following:

Code:
library(rvest)

widths <- c(11,10,10,10,14,5,5,5,48,100)

dat <- "http://www.koeri.boun.edu.tr/scripts/lst5.asp" %>%
   read_html %>%
   html_nodes("pre") %>%
   html_text %>%
   textConnection %>%
   read.fwf(widths = widths, stringsAsFactors = FALSE) %>%
   setNames(nm = .[6,]) %>%
   tail(-7) %>%
   head(-2)
 
#3
Thank you for your response jamesmartinn. I also would like to plot the map of the earthqakes happened. Below is what I got so far. But my question below:

Code:
#Loading necessary packages
library(rvest)
library(ggmap)
library(ggplot2)
library(dplyr)

##reading data set
widths <- c(11,10,10,10,14,5,5,5,48,100)
dat <- "http://www.koeri.boun.edu.tr/scripts/lst5.asp" %>%
  read_html %>%
  html_nodes("pre") %>%
  html_text %>%
  textConnection %>%
  read.fwf(widths = widths, stringsAsFactors = FALSE) %>%
  setNames(nm = .[6,]) %>%
  tail(-7) %>%
  head(-2)

##converting to numeric vector
dat$"Boylam(E)" <- as.numeric(as.character(dat$"Boylam(E)"))
dat$"Enlem(N)" <- as.numeric(as.character(dat$"Enlem(N)"))
dat<- dat[, -c(3,4)]

##adding range to ML
dat$range = case_when(dat$ML>0 & dat$ML<=1 ~ 1, #0<ML<=1
                      dat$ML>1 & dat$ML<=2 ~ 2,#1<ML<=2
                      dat$ML>2 & dat$ML<=3 ~ 3,#2<ML<=3
                      dat$ML>3 & dat$ML<=4 ~ 4,#3<ML<=4
                      dat$ML>4 & dat$ML<=5 ~ 5,#4<ML<=5
                      dat$ML>5 & dat$ML<=6 ~ 6,#5<ML<=6
                      dat$ML>6 & dat$ML<=7 ~ 7)#6<ML<=7


##first ploting Country
map <- get_map("Turkey", zoom = 5, maptype = "roadmap")
ggmap(map) +geom_point(data = dat,aes(x = "Boylam(E)", y ="Enlem(N)"),size = range ,alpha = .5)
 
Last edited:
#4
Dear jamesmartinn,
When I run the top code, I got error message below. Can you help with that?
Thank you.
Code:
 ggmap(map) +geom_point(data = dat,aes(x = "Boylam(E)", y ="Enlem(N)"),size = range ,alpha = .5)
Error: Discrete value supplied to continuous scale