• R in action读书笔记(22)第十六章 高级图形进阶(下)


    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提供了一个无缝接口。

  • 相关阅读:
    以友盟+U-Push为例,深度解读消息推送的筛选架构解决方案应用与实践
    [Python图像处理] 三十三.图像各种特效处理及原理万字详解(毛玻璃、浮雕、素描、怀旧、流年、滤镜等)
    走进PEP8——代码规范
    2020全球C++及系统软件技术大会成功落下帷幕
    逆向工程,调试Hello World !程序(更新中)
    520了,用32做个简单的小程序
    细数那些年我用过的前端开发工具
    细数那些年我用过的前端开发工具
    前端几个常用简单的开发手册拿走不谢
    前端几个常用简单的开发手册拿走不谢
  • 原文地址:https://www.cnblogs.com/jpld/p/4496017.html
Copyright © 2020-2023  润新知