為什么要學(xué)習(xí)R語言

作者：佚名 2017-04-05 18:10:05

作為開源軟件的R能夠迅速發(fā)展，很大程度上取決于其活躍的社區(qū)。學(xué)習(xí)R，很大程度上也是學(xué)習(xí)各種R包的使用。截止目前(2017年2月25日)，CRAN(Comprehensive R Archive Network)上已經(jīng)有10162個(gè)可以獲取的R擴(kuò)展包，內(nèi)容涉及各行各業(yè)，可以適用于各種復(fù)雜的統(tǒng)計(jì)。各地的CRAN鏡像都是R網(wǎng)站的備份文件，內(nèi)容完全一樣，你可以選擇離自己最近的去訪問。

R的誕生

1992年，肉絲 (Ross Ihaka)和蘿卜特 (Robert Gentleman)兩個(gè)人在S語言(貝爾實(shí)驗(yàn)室開發(fā)的一種統(tǒng)計(jì)用編程語言)的基礎(chǔ)上開始構(gòu)思一種新的用于統(tǒng)計(jì)學(xué)分析的開源語言，直到1995年第一個(gè)版本正式發(fā)布(和各位年齡相仿)。因?yàn)樗麄兠值牡谝粋€(gè)字母都是R，所以這門語言就被叫做R。這兩個(gè)人都是統(tǒng)計(jì)學(xué)教授出身，再加上R語言的生父S語言，所以 R語言在統(tǒng)計(jì)學(xué)方面有著純正的血統(tǒng) !

如果你平時(shí)的工作會(huì)涉及到統(tǒng)計(jì)學(xué)，那么接觸R語言實(shí)在是太正常不過了。

另外，關(guān)于R語言的開發(fā)者，看名字(Ross和Robert)部分人以為是伉儷，其實(shí)就是兩個(gè)大老爺們。如圖所示

[[187484]]

R的發(fā)展

R的特長

在R官網(wǎng)有這樣幾句介紹

R provides a wide variety of statistical (linear and nonlinear modelling, classical statistical tests, time-series analysis, classification, clustering, …) and graphical techniques, and is highly extensible.

One of R’s strengths is the ease with which well-designed publication-quality plots can be produced, including mathematical symbols and formulae where needed.

因?yàn)镽語言本身為統(tǒng)計(jì)而生，所以你能想到的所有統(tǒng)計(jì)相關(guān)的工作，R都可以非常簡潔的用幾行命令(甚至1行命令)幫你完成。
R高度的可擴(kuò)展性正是體現(xiàn)在它那1萬多個(gè)包上，你想做的幾乎所有事情都可以用現(xiàn)有的R包來輔助完成(當(dāng)然，有些工作即便能完成但也不適合)。
R另一個(gè)殺手锏就是其強(qiáng)大的繪圖功能，正如上面的英文介紹所言，R可以畫圖，畫各種各樣的圖，畫各種各樣高逼格的圖，畫各種各樣高逼格可以直接出版的圖。
完善的統(tǒng)計(jì)學(xué)功能再加上強(qiáng)大的繪圖功能，就是你學(xué)習(xí)的最大理由。
對于生物相關(guān)的工作者而言，他們還有一個(gè)巨大的福利就是Bioconductor，這里面的一千多個(gè)R包都是用來解決生物(信息)問題的。

R應(yīng)用示例

在這一部分，僅僅是給展示幾個(gè)用R可以輕松完成的相對有趣的工作。

安裝對應(yīng)包后應(yīng)該可以直接運(yùn)行

示例1 ggplot2畫圖

#第一次使用，需要安裝相應(yīng)的包  
#以后只需要調(diào)用即可  
#install.packages("ggplot2")  
library("ggplot2")  
theta <- seq(0,24*pi, len=2000)  
radius <- exp(cos(theta)) - 2*cos(4*theta) + sin(theta/12)^5  
dd <- data.frame(x=radius*sin(theta), y=radius*cos(theta))  
ggplot(dd, aes(x, y))+geom_path()+xlab("")+ylab("")

效果展示

示例2 詞頻分析及詞云

# install.packages("wordcloud2")  
library(wordcloud2)  
wordcloud2(demoFreqC, size = 0.7, shape = 'diamond')

效果展示

示例3 查看我國各地空氣質(zhì)量

#install.packages("rvest")  
#install.packages("leafletCN")  
#install.packages("rgeos")  
Sys.setlocale("LC_CTYPE", "eng")  
library(rvest)  
library(leafletCN)  
library(rgeos)  
doc = read_html("http://www.pm25s.com/cn/rank/")  
cities = doc %>% html_nodes(".cityrank a") %>%  
html_text()  
cities = iconv(cities, "UTF-8", "UTF-8")  
AQI = doc %>% html_nodes("span[class^='lv']") %>%  
html_text() %>% .[c(F,F,T)] %>% as.numeric  
dat = data.frame(city = cities, AQI = AQI)  
geojsonMap(dat, "city",  
popup = paste0(dat$city,":",dat$AQI),  
palette = "Reds", legendTitle = "AQI")