```
#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)
```