[設定所需的函式庫(libraries)以及載入資料]

setwd("/media/hsusir/DATA/Rdata Practice/04Spatial/dengue fever")
dengue <- read.csv("./dengue-20151107-utf8.csv",  fileEncoding="utf8")

[Part 1].Data-distributed

1-1.敘述性統計

str(dengue)
## 'data.frame':    22258 obs. of  6 variables:
##  $ 確診日  : Factor w/ 141 levels "2015/01/06","2015/01/19",..: 1 2 3 4 5 6 7 8 9 10 ...
##  $ 區別    : Factor w/ 41 levels "七股區","下營區",..: 38 27 7 7 7 7 7 16 7 7 ...
##  $ 里別    : Factor w/ 530 levels "一甲里","七甲里",..: 103 262 77 77 77 323 323 219 253 323 ...
##  $ 道路名稱: Factor w/ 1489 levels "","30鄰五福街",..: 878 1310 296 296 296 183 183 1414 1357 183 ...
##  $ 緯度座標: num  23 23 23 23 23 ...
##  $ 經度座標: num  120 120 120 120 120 ...
summary(dengue)
##         確診日           區別           里別          道路名稱    
##  2015/09/24:  712   北區   :3529   正覺里 :  332   金華路 :  690  
##  2015/09/18:  642   中西區 :3455   永祥里 :  319   公園路 :  671  
##  2015/09/21:  604   永康區 :2632   大豐里 :  299   西門路 :  627  
##  2015/09/22:  599   北 區 :2179   勝利里 :  272   文賢路 :  430  
##  2015/09/23:  584   安南區 :1831   五王里 :  269   長榮路 :  423  
##  2015/09/10:  553   南區   :1824   成德里 :  266   府安路 :  343  
##  (Other)   :18564   (Other):6808   (Other):20501   (Other):19074  
##     緯度座標        經度座標    
##  Min.   :22.63   Min.   :118.8  
##  1st Qu.:22.99   1st Qu.:120.2  
##  Median :23.00   Median :120.2  
##  Mean   :23.01   Mean   :120.2  
##  3rd Qu.:23.01   3rd Qu.:120.2  
##  Max.   :24.95   Max.   :121.5  
## 

1-2.畫在地圖上來看

install.packages("ggmap")
install.packages("mapproj")
library(ggmap)
library(mapproj)

map <- get_map(location = "Taiwan", zoom = 7, language = "zh-TW", maptype = "roadmap")
ggmap(map, darken = c(0.5, "white")) +geom_point(aes(x = 經度座標, y = 緯度座標),
  color = "red", data = dengue)