zoukankan      html  css  js  c++  java
  • ggplot2(2) 从qplot开始入门

    2.1 简介

    qplot的意思是快速作图(quick plot)。

    qplot是一种快捷方式,如果您已习惯于使用基础plot(),则可以使用它。它可以使用一致的调用模式快速创建许多不同类型的图。

    qplot(x, y, ..., data, facets = NULL, margins = FALSE, geom = "auto", xlim = c(NA, NA), ylim = c(NA, NA), log = "", main = NULL, xlab = NULL, ylab = NULL, asp = NA, stat = NULL, position = NULL)

    • x,y:图中对象的x坐标和y坐标;
    • data:可选,用于指定数据框,若进行了指定,那么函数会首先在该数据框内查找变量名,如果没有指定,将创建一个,从当前环境中提取向量;
    • facets:用于图形分面;
    • margins:逻辑值或字符向量。margins是附加的面;
    • geom:指定绘图类型的字符向量。如果指定x和y,默认为“点”;如果指定x,默认为“柱状图”;
    • xlim,ylim:x和y坐标边界;
    • log:哪些变量要进行对数转换(“x”、“y”或“xy”);
    • main:标题;
    • xlab,ylab:x和y坐标轴标签;
    • asp:y/x的长宽比。

    2.2 数据集

    > head(diamonds)
    # A tibble: 6 x 10
      carat cut       color clarity depth table price     x     y     z
      <dbl> <ord>     <ord> <ord>   <dbl> <dbl> <int> <dbl> <dbl> <dbl>
    1 0.23  Ideal     E     SI2      61.5    55   326  3.95  3.98  2.43
    2 0.21  Premium   E     SI1      59.8    61   326  3.89  3.84  2.31
    3 0.23  Good      E     VS1      56.9    65   327  4.05  4.07  2.31
    4 0.290 Premium   I     VS2      62.4    58   334  4.2   4.23  2.63
    5 0.31  Good      J     SI2      63.3    58   335  4.34  4.35  2.75
    6 0.24  Very Good J     VVS2     62.8    57   336  3.94  3.96  2.48

    A dataset containing the prices and other attributes of almost 54,000 diamonds. The variables are as follows:

    • price: price in US dollars ($326–$18,823)
    • carat: weight of the diamond (0.2–5.01)
    • cut: quality of the cut (Fair, Good, Very Good, Premium, Ideal)(均匀、良好、非常好、优质、理想)
    • color: diamond colour, from D (best) to J (worst)
    • clarity: a measurement of how clear the diamond is (I1 (worst), SI2, SI1, VS2, VS1, VVS2, VVS1, IF (best))
    • x: length in mm (0–10.74)
    • y: width in mm (0–58.9)
    • z: depth in mm (0–31.8)
    • depth: total depth percentage = z / mean(x, y) = 2 * z / (x + y) (43–79)
    • table: width of top of diamond relative to widest point (43–95)
    > str(diamonds)
    Classes ‘tbl_df’, ‘tbl’ and 'data.frame':	53940 obs. of  10 variables:
     $ carat  : num  0.23 0.21 0.23 0.29 0.31 0.24 0.24 0.26 0.22 0.23 ...
     $ cut    : Ord.factor w/ 5 levels "Fair"<"Good"<..: 5 4 2 4 2 3 3 3 1 3 ...
     $ color  : Ord.factor w/ 7 levels "D"<"E"<"F"<"G"<..: 2 2 2 6 7 7 6 5 2 5 ...
     $ clarity: Ord.factor w/ 8 levels "I1"<"SI2"<"SI1"<..: 2 3 5 4 2 6 7 3 4 5 ...
     $ depth  : num  61.5 59.8 56.9 62.4 63.3 62.8 62.3 61.9 65.1 59.4 ...
     $ table  : num  55 61 65 58 58 57 57 55 61 61 ...
     $ price  : int  326 326 327 334 335 336 336 337 337 338 ...
     $ x      : num  3.95 3.89 4.05 4.2 4.34 3.94 3.95 4.07 3.87 4 ...
     $ y      : num  3.98 3.84 4.07 4.23 4.35 3.96 3.98 4.11 3.78 4.05 ...
     $ z      : num  2.43 2.31 2.31 2.63 2.75 2.48 2.47 2.53 2.49 2.39 ...

    可以看到cut、color、clarity为有序分类变量。

    抽样:(方便观察,加快运算速度)

    set.seed(1410)  # 让样本可重复
    
    dsmall <- diamonds[sample(nrow(diamonds), 100), ]

    2.3 基本用法

    绘制散点图。

    qplot(carat, price, data = diamonds)

    对数变换(对坐标系进行变换,并没有改变变量的值)。

    qplot(carat, price, data = diamonds, log = 'xy')

    2.4 颜色、大小、形状和其他图形属性

    自动向重量和价格的散点图中添加颜色和切工的信息。

    qplot(carat, price, data = dsmall, colour = color)
    
    qplot(carat, price, data = dsmall, shape = cut)


    不透明度选择,在数据量很大时,太多的点重合,影响观察,使用半透明的颜色可以有效减轻图形元素重叠的现象,观察数据集聚趋势,alpha值可以设定不透明度,取值从0(完全透明)到1(完全不透明),通常用分数来表示,例如1/10或1/20,其分母表示经过多少次重叠之后颜色将变得不透明。

    qplot(carat, price, data = diamonds, log = 'xy', alpha = I(1/10))

    qplot(carat, price, data = diamonds, colour = color, alpha = I(1/5))

    2.5 几何对象

    二维对象:

    • geom = "point":绘制散点图;
    • geom = "smooth":拟合一条平滑曲线,并将标准误展示在图中,如果不需要展现标准误可以使用se = FALSE;
    • geom = "boxplot":绘制箱线图,用以概括一系列点的分布情况;
    • geom = "jitter":绘制扰动点图;
    • geom = "path"或geom = "line":可以在数据点之间绘制连线,这类图形的作用是探索时间和其他变量之间的关系,但连线同样可以用其他的方式将数据点连接起来,线条图只能创建从左到右的连线,而路径图则可以使任意的方向。

    一维对象:

    • geom = "histogram":绘制直方图;
    • geom = "freqpoly":绘制频率多边形;
    • geom = "density":绘制密度曲线;
    • geom = "bar:对离散变量绘制条形图。

    2.5.1 向图中添加平滑曲线

    使用c()将多个对象传递给geom,几何对象会按照指定的顺序进行堆叠。

    qplot(carat, price, data = dsmall, geom = c("point", "smooth"))

    利用method参数可以选择不同的平滑器。

    • method = "loess"使用的是局部回归方法,在当n较小是是默认选项,其平滑程度由span控制,取值范围为0(很不平滑)到1(很平滑)。Loess对于大数据并不合适,因为其内存消耗是$O(n^{2})$。
    • 使用method = "gam", formula = y ~ s(x)可以调用mgcv包拟合一个广义可加模型。对于大数据应使用y ~ s(x, bs = "cs")。
    library(mgcv)
    qplot(carat, price, data = dsmall, geom = c("point", "smooth"), method = "gam", formula = y ~ s(x))

    • method = "lm"拟合的是线性模型,默认得到一条直线,但可以通过formula = y ~ poly(x, 2)来拟合一个二次多项式。
    qplot(carat, price, data = dsmall, geom = c("point", "smooth"), method = "lm") 

    • method = "rlm"与"lm"类似,但采用了一种更为稳健的拟合方法,使得结果对异常值不太敏感,这一方法是MASS包的一部分,因此使用时需要先加载MASS包。

    2.5.2 箱线图和扰动点图

    qplot(color, price/carat, data = diamonds, geom = "jitter", alpha = I(1/30))
    qplot(color, price/carat, data = diamonds, geom = "boxplot")

    每种方法都有它的优势和不足,箱线图只用了5个数字进行描述,更具概括性;扰动点图绘制了全部点,更具全面性。

    2.5.3 直方图和密度曲线图

    直方图和密度曲线图可以展现单个变量的分布,相对于箱线图而言,它们提供了更多的关于单个变量分布的信息,但它们不太容易在不同组之间进行比较。

    qplot(carat, data = diamonds, geom = "histogram")
    qplot(carat, data = diamonds, geom = "density")

    对于密度曲线而言,adjust参数控制了曲线的平滑程度(adjust取值越大,曲线越平滑);对于直方图,binwidth参数通过设置组距来调节平滑度,或者也可以使用breaks对切分位置进行显示的指定。

    在直方图中,应当尝试多种组距:组距较大时图形能够反应总体特性;组距较小时,则能显示出更过细节。

    要在不同的分组之间进行比较,只需在加上一个图形映射。

    qplot(carat, data = diamonds, geom = "density", colour = color)
    qplot(carat, data = diamonds, geom = "histogram", fill = color)

    密度曲线图似乎更吸引人,因为很容易阅读,而且适于在不同的曲线之间进行比较。然而要真正理解密度曲线则比较困难,而且密度曲线有一些隐含的假设,例如曲线应该是无界、连续和平滑的,这些假设不一定适用于真实的数据。

    2.5.4 条形图

    用于绘制离散变量,与直方图类似。条形图会计算每一个水平下观测的数量。使用weight可以进行加权。

    qplot(color, data = diamonds, geom = "bar")
    qplot(color, data = diamonds, geom = "bar", weight = carat) + scale_y_continuous("carat")
    

    第一幅图展现了分组的计数,第二幅图展现了没中颜色钻石的总重量。

    2.5.5 时间序列中的线条图和路径图

    线条图和路径图用于可视化时间序列数据。线条图将点从左到右进行连接,而路径图则按照点在数据集中的顺序进行连接。线条图的x轴一般是时间,它展现了单个变量随时间的变化情况;路径图则展现了两个变量随时间联动的情况,时间反映在点的连接顺序上。

    数据集:economics

    > head(economics)
    # A tibble: 6 x 6
      date         pce    pop psavert uempmed unemploy
      <date>     <dbl>  <dbl>   <dbl>   <dbl>    <dbl>
    1 1967-07-01  507. 198712    12.6     4.5     2944
    2 1967-08-01  510. 198911    12.6     4.7     2945
    3 1967-09-01  516. 199113    11.9     4.6     2958
    4 1967-10-01  512. 199311    12.9     4.9     3143
    5 1967-11-01  517. 199498    12.8     4.7     3066
    6 1967-12-01  525. 199657    11.8     4.8     3018

    该数据集来自美国经济时间序列数据。

    • date:数据收集月份
    • psavert:个人储蓄率
    • pce:个人消费支出
    • unemploy:失业人数(以千计)
    • uempmed:失业持续时间的中位数,以周为单位。
    • pop:总人口,以千为单位

    绘制时间序列线条图。(失业率)

    qplot(date, unemploy/pop, data = economics, geom = "line")

    绘制路径图(失业率&失业持续时间)。将时间date映射到colour属性上更容易看出时间的行进方向。

    qplot(unemploy/pop, uempmed, data = economics, geom = c("point", "path"), colour = date)

    2.6 分面

    分面将数据分隔成若干子集,然后创建一个图形的矩阵,将一个子集绘制到矩阵的窗格中。

    通过row_var ~ col_var的表达式进行指定,如果指向指定一行或者一列,可以使用 . 作为占位符。

    qplot(carat, data = diamonds, facets = color ~ ., geom = "histogram", binwidth = 0.1, 
        xlim = c(0, 3))
    
    qplot(carat, ..density.., data = diamonds, facets = color ~ cut, geom = "histogram", 
        binwidth = 0.1, xlim = c(0, 3))

    ..density..告诉ggplot2将密度而不是频数映射到y轴。

    2.7 与plot函数的区别

    • qplot不是泛型函数,当将不同类型的R对象传入qplot时,它并不会自动匹配默认的函数调用。ggplot()是一个泛型函数;
    • gglpot2中的图形属性名称如colour、shape和size等比基础绘图系统中的名称如col、pch、cex等更直观,容易记忆;
    • 在基础绘图系统中,可以通过points()、lines()和text()函数来向已有的图形中添加更多的元素,而在ggplot2中,你需要在当前的图形中加入额外的图层。

    总结

  • 相关阅读:
    智能实验室-结构化存储浏览器(SSExplorer) 1.5.0.150
    智能实验室-杀马(Defendio) 3.1.0.681
    智能实验室-结构化存储浏览器(SSExplorer) 1.6.0.160
    IT餐馆—第八回 三十
    使用Silverlight Toolkit 绘制图表区域图和冒泡图
    IT餐馆—第十二回 软培
    IT餐馆—第四回 离职
    IT餐馆—第一回 前言
    IT餐馆—第十回 潜伏
    IT餐馆—第十三回 重构
  • 原文地址:https://www.cnblogs.com/dingdangsunny/p/12340894.html
Copyright © 2011-2022 走看看