16.2.4 图形参数
在lattice图形中,lattice函数默认的图形参数包含在一个很大的列表对象中,你可通过trellis.par.get()函数来获取,并用trellis.par.set()函数来修改。show.settings()函数可展示当前的图形参数设置情况。查看当前的默认设置,并将它们存储到一个mysettings列表中:
> show.settings() > mysettings<-trellis.par.get()
查看叠加点的默认设置值:
> mysettings$superpose.symbol $alpha [1] 1 1 1 1 1 1 1 $cex [1] 0.8 0.8 0.8 0.8 0.8 0.8 0.8 $col [1] "#0080ff" "#ff00ff" "darkgreen" "#ff0000" "orange" "#00ff00" "brown" $fill [1] "#CCFFFF" "#FFCCFF""#CCFFCC" "#FFE5CC" "#CCE6FF" "#FFFFCC""#FFCCCC" $font [1] 1 1 1 1 1 1 1 $pch [1] 1 2 3 4 5 6 7 8 9 10
可以看到分组变量的每个水平都使用空心圆圈(pch = 1)。定义了7个水平后,图形符号
将会被循环使用。最后,我们再做如下声明:
> mysettings$superpose.symbol$pch<-c(1:10) > trellis.par.set(mysettings) > show.settings()
此时lattice图形将对分组变量的第一个水平使用符号1(空心圆圈),第二个使用符号2(空心三角形),以此类推。另外,我们对分组变量的10个水平的符号都进行了定义,而不是7个。这种图形设置效果将会一直存在,直到关闭图形设备。你可以按照此方式对其他任意图形参数进行修改。
16.2.5 页面摆放
lattice函数不识别par()设置,先将lattice图形存储到对象中,然后利用plot()函数中的split =或position =选项来进行控制。split选项将页面分割为一个指定行数和列数的矩阵,然后将图形放置到该矩阵中。split选项的格式为:
split=c(placement row,placement column,total number of rows,total number of columns)
例:
> library(lattice) > graph1<-histogram(~height|voice.part,data=singer, + main="Heights of choral singers by voice part") > graph2<-densityplot(~height,data=singer,group=voice.part, + plot.points=FALSE,auto.key=list(columns=4)) > plot(graph1,split=c(1,1,1,2)) > plot(graph2,split=c(1,2,1,2),newpage=FALSE)
它将把第一幅图放置到第二幅图的上面。具体来讲,第一个plot()函数把页面分割成一列
两行的矩阵,并将图形放置到第一列、第一行中(自上往下、从左至右地计数)。第二个plot()
函数做同样的分割,但是把图形放置到第一列、第二行中。因为plot()函数默认启动一个新的页面,所以你需要禁止该操作,因此设定选项newpage = FALSE
使用position =选项可以对大小和摆放方式进行更多的控制:
> plot(graph1,position=c(0,.3,1,1)) > plot(graph2,position=c(0,0,1,.3),newpage=FALSE)
position = c(xmin, ymin, xmax, ymax),该页面的x-y坐标系统是矩形,x轴和y轴的维度范围都是从0到1,原点(0, 0)在图形左下角。在lattice图形中你还可以改变面板的顺序。高级绘图函数的index.cond =选项可以设定条件变量水平的顺序。以voice.part因子为例:
> levels(singer$voice.part) [1] "Bass 2" "Bass 1" "Tenor 2" "Tenor 1" "Alto 2" "Alto 1" "Soprano 2" [8] "Soprano 1"
16.3 ggplot2 包
ggplot2包提供了一个基于全面而连贯的语法的绘图系统。它弥补了R中创建图形缺乏一致
性的缺点,使得用户可以创建有创新性的、新颖的图形类型。ggplot2中最简单的绘图方式是利用qplot()函数,即快速绘图函数。格式为:
qplot(x, y = NULL, ..., data, facets =NULL, margins = FALSE,
geom = "auto", stat = list(NULL), position = list(NULL), xlim= c(NA,NA), ylim = c(NA, NA), log = "", main = NULL,xlab =deparse(substitute(x)), ylab = deparse(substitute(y)), asp = NA)
选项:描 述
alpha :元素重叠的alpha透明度,数值为0(完全透明)到1(完全不透明)间的分数
color、shape、size、fill :把变量的水平与符号颜色、形状或大小联系起来。对于直线图,colo将把线条颜色与变量水平联系起来,对于密度图和箱线图,fill将把填充颜色与变量联系起来。图例将会被自动绘制
data :指定一个数据框
facets :指定条件变量,创建一个栅栏图。表达式如rowvar~ colvar。为创建一个基于单条件变量的栅栏图,可用rowvar ~ .或. ~ colvar
geom :设定定义图形类型的几何形状。geom选项是一个单条目或多条目的字符型向量,
包括"point" 、"smooth" 、"boxplot" 、"line" 、"histogram"、"density"、"bar"和"jitter"
main、sub :字符向量,设定标题和副标题
method、formula :若geom ="smooth",则会默认添加一条平滑拟合曲线和置信区间。当观测数大于1000时,便需要调用更高效的平滑拟合算法。方法包括回归lm、广义可加模型gam、稳健回归rlm。formula参数指定拟合的形式例如,要添加简单的回归曲线,则设定geom = "smooth", method = "lm",formula = y ~ x。将表达式改为y ~poly(x, 2)将生成二次拟合。注意表达式使用的是字母x和y,而不是变量的名称对于method = "gam",一定要记得加载mgcv包。对于method ="rml",则需加载MASS包
x、y :指定摆放在水平轴和竖直轴的变量。对于单变量图形(如直方图),则省略y
xlab、ylab :字符向量,设定横轴和纵轴标签
xlim、ylim :二元素数值型向量,分别指定横轴和纵轴的最小值和最大值
library(ggplot2) mtcars$cylinder<-as.factor(mtcars$cyl) qplot(cylinder,mpg,data=mtcars,geom=c("boxplot","jitter"), fill=cylinder, main="Box plots with superimposed data points", xlab="Number of cylinders", yalb="Miles per gallon")
型添加回归线和置信区间带
传动类型定义了行分面,而气缸数则定义了列分面。数据点的大小代表了汽车马力
的评分:
library(ggplot2) mtcars$cyl<-factor(mtcars$cyl,levels=c(4,6,8), labels=c("4cylinders","6 cylinders","8 cylinders")) mtcars$am<-factor(mtcars$am,levels=c(0,1), labels=c("Automatic","manual")) qplot(wt,mpg,data=mtcars,facets=am~cyl,size=hp)
library(ggplot2) data(singer,package="lattice") qplot(height,data=singer,geom=c("density"), facets=voice.part~.,fill=voice.part)
16.4 交互式图形
16.4.1 与图形交互:鉴别点
可对散点图中的点进行鉴别和标注的函数:identify()。利用该函数,你可用鼠标对散点图中所选择的点标注行数或者行名称,直到你选择了Stop或者右击图形识别工作才会停止。
plot(mtcars$wt,mtcatrs$mpg) identify(mtcars#wt,mtcars$mpg,labels=row.names(mtcars))
光标将从一个点变成一个十字。单击散点图上的点,可以对它们进行标注,直到你从Graphics Device(图形设备)菜单中选择了Stop,或者右击了图形并从右键菜单中选择了Stop。
16.4.2 playwith
playwith包提供了一个GTK+图形用户界面(GUI),使得用户可以编辑R图形并与其交互。
playwith()函数允许用户识别和标注点、查看一个观测所有的变量值、缩放和旋转图形、
添加标注(文本、箭头、线条、矩形、标题和标签)、修改视觉元素(颜色、文本大小等)、应用先前存储的图形风格,以及以多种格式输出图形结果。
library(playwith) library(lattice) playwith( xyplot(mpg~wt|factor(cyl)*factor(am), data=mtcars,subscripts=TRUE, type=c("r","p")) )
playwith窗口。依赖于GTK+ GUI,用户可以使用鼠标编辑图形
与identify()函数不同,playwith()既对R基础图形有效,也对lattice和ggplot2
图形有效。Theme(主题)菜单上的一些选项仅与基础图形契合的很好,一些则与ggplot2图形契合的较好(如标注),还有些对ggplot2图形无效(如识别点)。
16.4.3 latticist
使用latticist包,可通过栅栏图方式探索数据集。该包不仅提供了图形用户界面,也可以通过vcd包来创建新的图形。如果有需要,latticist还能与playwith整合到一起,例:
library(latticist) mtcars$cyl<-factor(mtcars$cyl) mtcars$gear<-factor(mtcars$gear) latticist(mtcars,use.playwith=TRUE)
拥有latticist函数功能的playwith窗口。用户可创建lattice并与vcd图形交互
16.4.4 iplots包的交互图形
playwith和latticist包只能与单幅图形交互,而iplots包提供的交互方式则有所不同。该包提供了交互式马赛克图、柱状图、箱线图、平行坐标图、散点图和直方图,以及颜色刷,并可将它们结合在一起绘制。这意味着你可通过鼠标对观测点进行选择和识别,并且对其中一幅图形的观测点突出显示时,其他被打开的图形将会自动突出显示相同的观测点。另外,还可通过鼠标来收集图形对象(诸如点、条、线)和箱线图的信息。
函数 描 述
ibar() 交互式柱状图
ibox() 交互式箱线图
ihist() 交互式直方图
imap() 交互式地图
imosaic() 交互式马赛克图
ipcp() 交互式平行坐标图
iplot() 交互式散点图
例:
library(iplots) attach(mtcars) cylinders<-factor(cyl) gears<-factor(gear) transmission<-factor(am) ihist(mpg) ibar(gears) iplot(mpg,wt) ibox(mtcars[c("mpg","wt","qsec","disp","hp")]) ipcp(mtcars[c("mpg","wt","qsec","disp","hp")]) imosaic(transmission,cylinders) detach(mtcars)
现在尝试以下步骤。
单击柱状图(gears)窗口的三号齿轮条,直条将会变红色。另外,其他图形窗口中的所
有三齿轮发动机都会被突出显示。
鼠标下移,并在散点图(wt对mpg)窗口选择一个含有点的矩形区域。这些点将会被突出显示,而其他图形窗口中相应的观测点也将会变成红色。
按住Ctrl键不动,将鼠标移动到某幅图形中点、条、箱线图或线上,该对象的详细信息将会在一个弹出窗口中显示出来。
右击任何对象,便可在右键菜单中获得一些选项。例如,你可以右击箱线图(mpg)窗口,将图形转变为一个平行坐标图(PCP)。
拖动鼠标可选择不止一个对象(点、条等),或使用Shift键通过单击选择不邻接的对象。
你可尝试在柱状图(gears)窗口选择三号和五号齿轮条。
16.4.5 rggobi
GGobi有许多吸引眼球的优点,包括:交互式散点图、柱状图、平行坐标图、时间序列图、散点图矩阵和三维旋转的综合使用;窗口刷和点识别;多变量变换方法;复杂的探索平台,如导向动画的和手动的1维及2维图形。令人振奋的是,rggobi软件包为GGobi和R提供了一个无缝接口。