语言扩展包q工具箱讲义

1、R数据可视化ggplot2包 第5周2013.2.8DATAGURU专业数据分析法律【】本和幻灯片为炼数网络课程的教学资料，所有资料只能在课程内使用，不得在课程以外范围散播，违者将可能被责任。法律和课程详情炼数培训http:2013.2.82DATAGURU专业数据分析ggthemes包2013.2.83DATAGURU专业数据分析例子n 传统ggplot2风格ggplot(mtcars,aes(x=mpg,y=wt,size=cyl,colour=factor(gear)+ geom_point(alpha=.5)+ #scale_size_area()+ #区域和数值成比例scale_co

2、lour_brewer(palette="Set1")+ggtitle("Motor Trend Car Road Tests")2013.2.84DATAGURU专业数据分析例子n Excel风格ggplot(mtcars,aes(x=mpg,y=wt,size=cyl,colour=factor(gear)+ geom_point()+theme_excel() + scale_fill_excel()+ggtitle("Motor Trend Car Road Tests")2013.2.85DATAGURU专业数据分析例子n

3、Economist风格ggplot(mtcars,aes(x=mpg,y=wt,size=cyl,colour=factor(gear)+ geom_point(alpha=.5)+theme_economist() + scale_colour_economist()+ggtitle("Motor Trend Car Road Tests")2013.2.86DATAGURU专业数据分析例子n 深色风格ggplot(mtcars,aes(x=mpg,y=wt,size=cyl,colour=factor(gear)+ geom_point()+theme_solarize

4、d(light=FALSE) + scale_colour_solarized("red")+ggtitle("Motor Trend Car Road Tests")2013.2.87DATAGURU专业数据分析输出到文件ggsave()qplot(mpg, wt, data = mtcars)ggsave(file = "output.pdf")pdf(file = "output.pdf", width = 6, height = 6)# If inside a script, you will need t

5、o explicitly print() plots qplot(mpg, wt, data = mtcars)qplot(wt, mpg, data = mtcars)dev.off()2013.2.88DATAGURU专业数据分析结果2013.2.89DATAGURU专业数据分析在同一页面上画多幅图n 先产生三幅图(a <- qplot(date, unemploy, data = economics, geom = "line") (b <- qplot(uempmed, unemploy, data = economics) + geom_smooth(

6、se = F)(c <- qplot(uempmed, unemploy, data = economics, geom="path")2013.2.810DATAGURU专业数据分析viewport()Library(grid)# A viewport that takes up the entire plot devicevp1 <- viewport(width = 1, height = 1, x = 0.5, y = 0.5) vp1 <- viewport()# A viewport that takes up half the width a

7、nd half the height, # located in the middle of the plot.vp2 <- viewport(width = 0.5, height = 0.5, x = 0.5, y = 0.5) vp2 <- viewport(width = 0.5, height = 0.5)# A viewport that is 2cm x 3cm located in the centervp3 <- viewport(width = unit(2, "cm"), height = unit(3, "cm"

8、)2013.2.811DATAGURU专业数据分析实验pdf("polishing-subplot-1.pdf", width = 4, height = 4)subvp <- viewport(width = 0.4, height = 0.4, x = 0.75, y = 0.35) bprint(c, vp = subvp)dev.off()2013.2.812DATAGURU专业数据分析改进csmall <- c +theme_gray(9) +labs(x = NULL, y = NULL) + theme(plot.margin = unit(rep

9、(0, 4), "lines")pdf("polishing-subplot-2.pdf", width = 4, height = 4) bprint(csmall, vp = subvp)dev.off()2013.2.813DATAGURU专业数据分析结果2013.2.814DATAGURU专业数据分析指定比例pdf("polishing-layout.pdf", width = 8, height = 6)grid.newpage()pushViewport(viewport(layout = grid.layout(2, 2

10、) vplayout <- function(x, y)viewport(layout.pos.row = x, layout.pos.col = y) print(a, vp = vplayout(1, 1:2)print(b, vp = vplayout(2, 1) print(c, vp = vplayout(2, 2)dev.off()2013.2.815DATAGURU专业数据分析结果2013.2.816DATAGURU专业数据分析画函数图像p <- ggplot(data.frame(x=c(-3,3), aes(x=x)p + stat_function(fun =

11、dnorm)2013.2.817DATAGURU专业数据分析画函数图像myfun <- function(xvar) 1/(1 + exp(-xvar + 10)ggplot(data.frame(x=c(0, 20), aes(x=x) + stat_function(fun=myfun)2013.2.818DATAGURU专业数据分析画微学中常见的梯形# Return dnorm(x) for 0 < x < 2, and NA for all other xdnorm_limit <- function(x) y <- dnorm(x)yx < 0 |

12、 x > 2 <- NAreturn(y)# ggplot() with dummy datap <- ggplot(data.frame(x=c(-3, 3), aes(x=x)p + stat_function(fun=dnorm_limit, geom="area", fill="blue", alpha=0.2)+stat_function(fun=dnorm)2013.2.819DATAGURU专业数据分析结果2013.2.820DATAGURU专业数据分析ggplot2画向量场2013.2.821DATAGURU专业数据分析

13、向量场2013.2.822DATAGURU专业数据分析向量场2013.2.823DATAGURU专业数据分析让图形动起来n 并非ggplot2包的功能library(rgl)plot3d(mtcars$wt, mtcars$disp, mtcars$mpg, type="s", size=0.75, lit=FALSE) play3d(spin3d()2013.2.824DATAGURU专业数据分析在R中实现动态气泡图nn 使用ggplot2包，animation包，ffmpeg软件n 数据：/mbostock//blob/gh-pages/mik

14、e/nations/nations.jsonn 代码：2013.2.825DATAGURU专业数据分析代码#Import the Strange Datasetrawdata=readLines("nations.json")#Get the useful wordsdata=unlist(strsplit(rawdata,""|":"|","|":|,|,"|,"|") data=datawhich(data!="")n=length(data)2013

15、.2.826DATAGURU专业数据分析代码#Set indexnameind=which(data="name") regind=which(data="region")incind=which(data="e")popind=which(data="population")lifind=which(data="lifeExpectancy") endpoint=nameind-1 endpoint=endpoint2:180endpoint180=45999name=datanameind+

16、1region=dataregind+12013.2.827DATAGURU专业数据分析代码#An empty data.frame typeinitdfr=data.frame(name=rep(0,210),region=rep(0,210),yee=rep(0,210),pop=rep(0,210),life=rep(0,210)#An string split function only working for the commasplt=function(x) return(as.numeric(unlist(strsplit(x,",")#Linear inte

17、rpolation function, for those zero dataitpl=function(a)ind=which(a>0) if (ind1>1)a1:(ind1-1)=rep(aind1,ind1-1) n=length(ind)if (indn<length(a)a(indn+1):length(a)=rep(aindn,length(a)-indn) for (i in 1:(n-1)aindi:indi+1=rep(aindi,indi+1-indi+1)+(seq(indi,indi+1,1)-indi)*(aindi+1-aindi)/(indi+

18、1-indi)return(a)2013.2.828DATAGURU专业数据分析代码#Drag data information from wordstbls=NULL for (i in 1:180)dfr=initdfr dfr$name=rep(namei,210) dfr$region=rep(regioni,210)for ( j in (incindi+1):(popindi-1)tmp=spata j)ii=which(dfr$year=tmp1)dfr$eii=tmp22013.2.829DATAGURU专业数据分析代码for ( j in (popindi+1):(lifin

19、di-1)tmp=spata j)ii=which(dfr$year=tmp1)dfr$popii=tmp2for ( j in (lifindi+1):endpointi)tmp=spata j)ii=which(dfr$year=tmp1)dfr$lifeii=tmp2tbls=rbind(tbls,dfr)2013.2.830DATAGURU专业数据分析代码#Two country with only one record, meaninglessind=which(tbls$name="Mayotte") tbls=tbls-ind, ind=which(tbls$

20、name="Tokelau") tbls=tbls-ind,name=namec(-28,-177)region=regionc(-28,-177)2013.2.831DATAGURU专业数据分析代码#Linear interpolationfor (i in 1:178)ind=(210*(i-1)+1):(i*210)tbls$eind=itpl(tbls$eind)tbls$popind=itpl(tbls$popind)tbls$lifeind=itpl(tbls$lifeind)2013.2.832DATAGURU专业数据分析代码require(ggplot2)#

21、Draw function with ggplot2drawit=function(yr,scl=15)ind=which(tbls$year=yr) d.f=data.frame(yr=yr)p=ggplot(aes(x=log(e),y=life,size=pop,colour=as.factor(region),pch=21,data=tblsind,)p+geom_point(show_guide = FALSE)+geom_point(shape = 1,colour = "black",show_guide = FALSE)+ xlim(5.5,11.7)+yl

22、im(10,83)+scale_area(range = c(1, scl)+annotate("text", x=10, y=15, label = yr,size=30,color="grey")#drawit(1800)2013.2.833DATAGURU专业数据分析代码#Automatically repeat the drawing procedurefinaldraw=function(a,b)for (i in 1:10)print(drawit(a) for (i in a:b)print(drawit(i) for (i in 1:10)print(drawit(b)2013.2