【原】用R語言制作交互式圖表和地圖

拓端數(shù)據(jù) 2020-05-19

展開全文

原文鏈接：http:///?p=8032

可以直接從R / RStudio制作在線交互式圖表和地圖。

配置

啟動(dòng)RStudio，創(chuàng)建一個(gè)新的RScript,，然后將工作目錄設(shè)置為下載的數(shù)據(jù)文件夾。

使ggplot2圖表成為交互式Plotly圖表

制作交互式點(diǎn)線圖

以下代碼將安裝并加載程序包（該程序還將自動(dòng)加載ggplot2）,，加載readr和dplyr,，然后加載我們之前使用的食品券數(shù)據(jù)。

# install and load plotly, load readr and dplyrinstall.packages("plotly")library(plotly)library(readr)library(dplyr)

將其轉(zhuǎn)變?yōu)镻lotly圖,。

# load datafood_stamps <- read_csv("food_stamps.csv")

# dot-and-line chartfood_stamps_chart <- ggplot(food_stamps, aes(x = year, y = participants)) + xlab("Year") + ylab("Participants (millions)") + theme_minimal(base_size = 14, base_family = "Georgia") geom_point() + geom_line()

plot(food_stamps_chart)

如前所述,，這會(huì)將ggplot2圖表保存在您的環(huán)境中。以下代碼將其轉(zhuǎn)換為“繪圖”圖表,，該圖表應(yīng)顯示在Viewer右下角的標(biāo)簽中：

當(dāng)您將鼠標(biāo)懸停在默認(rèn)圖表上時(shí),，一些控件會(huì)顯示在右上角。

要重新格式化工具提示,，我們需要同時(shí)修改ggplot2和plotly代碼

最終layout功能中的代碼代碼為工具提示設(shè)置樣式,，使用白色背景并更改字體系列，使其與圖表的其余部分保持一致,。將代碼插入頁面并按如下所示進(jìn)行編輯：

<div class="container">

<iframe width="100%" height="450" frameborder="0" scrolling="no" src="food_stamps_interactive.html"></iframe>

</div> 

結(jié)果應(yīng)該是：

制作疾病和民主散點(diǎn)圖的交互式版本

此代碼創(chuàng)建具有線性趨勢(shì)線的基本散點(diǎn)圖,，而沒有映射到點(diǎn)income_group的顏色：

結(jié)果應(yīng)該是：

以下代碼使用定性的ColorBrewer調(diào)色板為這些點(diǎn)創(chuàng)建此圖表的版本，以由世界銀行收入組為其著色,。

這是靜態(tài)版本：

這是交互式的：

請(qǐng)注意,，交互式ggplot版本未繼承scale_color_brewer對(duì)靜態(tài)ggplot圖表中的圖例中的項(xiàng)目進(jìn)行排序的代碼,。

但是我們可以通過首先運(yùn)行以下代碼來解決此問題：

這將轉(zhuǎn)換income_group為分類變量或factor，然后按其類別或級(jí)別的順序?qū)?shù)據(jù)框進(jìn)行排序,。

現(xiàn)在,，像以前一樣運(yùn)行圖表代碼將可以固定交互式圖表中圖例中項(xiàng)目的順序。

制作交互式版本的加州幼兒園免疫接種熱圖

這是另一個(gè)示例,。

結(jié)果應(yīng)該是：

練習(xí)制作其他交互式圖表

在課堂上,，在時(shí)間允許的情況下，我們將使用plotly 創(chuàng)建這些ggplot2圖表的交互式版本：

使用Leaflet制作地震風(fēng)險(xiǎn)圖和地震圖

制作交互式在線地圖的最廣泛使用的JavaScript庫,。

# install and load leaflet and rdgalinstall.packages("leaflet")install.packages("rgdal")library(leaflet)library(rgdal)

首先,，讓我們看看如何制作以伯克利為中心的基本Leaflet地圖：

該地圖應(yīng)顯示在Viewer：

該leaflet函數(shù)創(chuàng)建一個(gè) 地圖。

該addProviderTiles函數(shù)使用Leaflet Providers插件將各種圖塊添加到地圖,。

現(xiàn)在seismic使用rgdal中的readOGR函數(shù)加載從shapefile 開始的地震地圖所需的數(shù)據(jù),。

提到的兩個(gè)分別seismic指向文件夾和其中的shapefile。

現(xiàn)在你應(yīng)該在你的環(huán)境中調(diào)用的對(duì)象seismic是一個(gè)SpatialPolygonsDataFrame,。

我們還可以直接從美國地質(zhì)調(diào)查局地震API中加載地震數(shù)據(jù),，如有關(guān)查找和下載數(shù)據(jù)的課程說明中所述：

使用該網(wǎng)址，我們加載了自1960年初以來6級(jí)及以上的地震,，地震發(fā)生在美國大陸地理中心的6,000公里半徑內(nèi),。

讓我們看一下seismic數(shù)據(jù)摘要：

# view summary of seismic_risk datasummary(seismic)

定義破壞性地震的年度風(fēng)險(xiǎn)的數(shù)據(jù)是可變的ValueRange。但是,，此合并變量的類別順序不正確,。要更正此問題，我們應(yīng)將變量從文本轉(zhuǎn)換為factor,，或levels按正確的順序進(jìn)行分類,。

現(xiàn)在，類別應(yīng)按正確的順序排列：

請(qǐng)注意,，要運(yùn)行dplyr代碼以mutate

接下來,，我們將地震風(fēng)險(xiǎn)數(shù)據(jù)加載到地圖中：

現(xiàn)在，您應(yīng)該leaflet在環(huán)境中看到一個(gè)類型的對(duì)象,。

現(xiàn)在,，我們只需兩行代碼就可以創(chuàng)建具有兩層的地圖：

結(jié)果應(yīng)該是：

該函數(shù)colorFactor將命名的ColorBrewer調(diào)色板分配給類別變量。

該函數(shù)addPolygons將多邊形添加到地圖：不stroke = FALSE給它們輪廓,；fillOpacity = 0.7使它們稍微透明,；color = ~pal(ValueRange))使用調(diào)色板根據(jù)ValueRange數(shù)據(jù)中的值為多邊形著色。