R語言中g(shù)gplot2繪制雙坐標(biāo)軸圖

第R語言中g(shù)gplot2繪制雙坐標(biāo)軸圖R包ggplot2繪圖精美，可以做出很復(fù)雜的圖形，深受用戶喜愛。它的作者h(yuǎn)adley并不推薦使用ggplot2繪制雙坐標(biāo)軸圖，認(rèn)為這樣會(huì)增加讀圖的難度，但是目前需要雙坐標(biāo)軸圖應(yīng)用的場(chǎng)景還是很多，如下圖形直方圖和曲線分別有不同的坐標(biāo)軸，表示不同的意思，圖片展示內(nèi)容比較豐富

今天我們通過我們的汽車銷售數(shù)據(jù)（公眾號(hào)回復(fù)：汽車銷售，可以獲得數(shù)據(jù)）來演示一下ggplot2繪制雙坐標(biāo)軸圖，主要是通過sec.axis函數(shù)進(jìn)行轉(zhuǎn)換，進(jìn)行第二個(gè)坐標(biāo)軸的轉(zhuǎn)換。我們先來導(dǎo)入數(shù)據(jù)

library(foreign)

library(ggplot2)

library(plyr)

bc-read.spss("E:/r/test/tree_car.sav",

use.value.labels=F,to.data.frame=T)

我們來看下數(shù)據(jù)，car就是汽車售價(jià)，age是年齡，gender是性別，inccat是收入，這里分成4個(gè)等級(jí)，ed是教育程度。

我們用這個(gè)數(shù)據(jù)生成兩個(gè)子數(shù)據(jù)，展示不同教育程度人群購車的區(qū)別，繼續(xù)使用原來的summarySE函數(shù)生成數(shù)據(jù)，想了解這個(gè)函數(shù)的看《R語言繪制帶誤差和可信區(qū)間的折線圖》這篇文章，主要是利用這個(gè)函數(shù)生成不同類型的數(shù)據(jù)分析指標(biāo)，函數(shù)的代碼如下：

summarySE-function(data=NULL,measurevar,groupvars=NULL,na.rm=FALSE,

erval=.95,.drop=TRUE){

library(plyr)

#NewversionoflengthwhichcanhandleNA's:ifna.rm==T,don'tcountthem

length2-function(x,na.rm=FALSE){

if(na.rm)sum(!is.na(x))

elselength(x)

#Thisdoesthesummary.Foreachgroup'sdataframe,returnavectorwith

#N,mean,andsd

datac-ddply(data,groupvars,.drop=.drop,

.fun=function(xx,col){

c(N=length2(xx[[col]],na.rm=na.rm),

mean=mean(xx[[col]],na.rm=na.rm),

sd=sd(xx[[col]],na.rm=na.rm)

measurevar

#Renamethe"mean"column

datac-rename(datac,c("mean"=measurevar))

datac$se-datac$sd/sqrt(datac$N)#Calculatestandarderrorofthemean

#Confidenceintervalmultiplierforstandarderror

#Calculatet-statisticforconfidenceinterval:

#e.g.,ifervalis.95,use.975(above/below),andusedf=N-1

ciMult-qt(erval/2+.5,datac$N-1)

datac$ci-datac$se*ciMult

return(datac)

這個(gè)函數(shù)并不是今天的主要內(nèi)容，主要是看我們生成的數(shù)據(jù)

carss-summarySE(bc,measurevar="car",groupvars=c("gender","ed"))

carss1-summarySE(bc,measurevar="car",groupvars=c("inccat","ed"))

我們來看看生成的兩個(gè)數(shù)據(jù)

Carss

carss1

carss1和carss都有ed這個(gè)數(shù)據(jù)，carss1的數(shù)據(jù)比carss多，carss1有inccat這個(gè)指標(biāo)carss沒有

我們先畫一個(gè)carss數(shù)據(jù)的教育和汽車售價(jià)的直方圖

ggplot()+

geom_col(data=carss,aes(x=ed,y=car),fill="#6794a7")

上圖展示了不同教育程度買車的情況，假設(shè)我們想把carss1這個(gè)數(shù)據(jù)的inccat和ci這個(gè)指標(biāo)映射進(jìn)上圖中（不管臨床意義，只是展示數(shù)據(jù)映射），我們可以看到ci這個(gè)指標(biāo)很小0-5這樣，而汽車這個(gè)指標(biāo)的范圍是0-80，不能通過同一坐標(biāo)軸展示，所以只能通過分軸展示。主要是通過rescale函數(shù)先把ci這個(gè)指標(biāo)在在0-80區(qū)間進(jìn)行標(biāo)準(zhǔn)化轉(zhuǎn)換，然后在投射于右坐標(biāo)軸。

先把inccat轉(zhuǎn)換成因子

carss1$inccat-as.factor(carss1$inccat)

繪圖

ggplot()+

geom_col(data=carss,aes(x=ed,y=car),fill="#6794a7")+

geom_line(data=carss1,aes(x=ed,y=rescale(ci,c(0,80)),colour=inccat,group=inccat),size=1.5)+

geom_point(data=carss1,aes(x=ed,y=rescale(ci,c(0,80))),shape=21,fill="white",size=4)+

scale_y_continuous(breaks=pretty_breaks(4),

sec.axis=sec_axis(~rescale(.,c(0,5)),name="ci",labels=sprintf('%.5f',(0:5))))

圖形左軸表示直方圖中購車情況，右軸表示不同收入人群的ci情況，把圖形進(jìn)一步美化一下

ggplot()+

geom_col(data=carss,aes(x=ed,y=car),fill="#6794a7")+

geom_line(data=carss1,aes(x=ed,y=rescale(ci,c(0,80)),colour=inccat,group=inccat),size=1.5)+

geom_point(data=carss1,aes(x=ed,y=rescale(ci,c(0,80))),shape=21,fill="white",size=4)+

scale_y_continuous(breaks=pretty_breaks(4),

sec.axis=sec_axis(~rescale(.,c(0,5)),name="ci",labels=sprintf('%.5f',(0:5))))+

labs(title="ThisisaTitle!",subtitle="ThisisaSubtitle",caption="ThisisaCaption")+

theme_minimal(base_size=16)%+r