1 Leitura dos dados

Para explorar os recursos básicos de visualização, serão usados dados reais sobre a venda do modelo Renault Duster extraídos da WEB. Os dados estão disponíveis online em http://www.leg.ufpr.br/~walmes/data/duster_venda_260314.txt.

Os dados contém informações sobre o preço de venda (R$, valor) e distância percorrida (km, km) além de outras características descritoras do veículo, como ano de fabricação, cor, tipo de câmbio.

#-----------------------------------------------------------------------
# Dados de carros Duster à venda no webmotors em 26/03/2014.

# Importa a tabela de dados da web.
url <- "http://www.leg.ufpr.br/~walmes/data/duster_venda_260314.txt"
dus <- read.table(file = url,
                  header = TRUE,
                  sep = "\t",
                  encoding = "utf-8")
str(dus)

## 'data.frame':    699 obs. of  10 variables:
##  $ modelo: chr  "RENAULT DUSTER 1.6 EXPRESSION 4X2 16V FLEX 4P MANUAL" "RENAULT DUSTER 1.6 4X2 16V FLEX 4P MANUAL" "RENAULT DUSTER 1.6 DYNAMIQUE 4X2 16V FLEX 4P MANUAL" "RENAULT DUSTER 1.6 DYNAMIQUE 4X2 16V FLEX 4P MANUAL" ...
##  $ cor   : chr  "Prata" "Azul" "Prata" "Prata" ...
##  $ km    : int  31442 40800 56000 NA 45000 50000 44000 30000 41000 55000 ...
##  $ ano   : chr  "2011/2012" "2011/2012" "2011/2012" "2011/2012" ...
##  $ valor : num  41990 42500 42900 42990 43800 ...
##  $ cambio: chr  "MANUAL" "MANUAL" "MANUAL" "MANUAL" ...
##  $ poten : num  1.6 1.6 1.6 1.6 1.6 1.6 1.6 1.6 2 1.6 ...
##  $ trac  : chr  "4X2" "4X2" "4X2" "4X2" ...
##  $ cat   : chr  " EXPRESSION " " " " DYNAMIQUE " " DYNAMIQUE " ...
##  $ novo  : chr  "usado" "usado" "usado" "usado" ...

# Cria ano do veículo com extração regex.
dus$ano <- factor(gsub(x = as.character(dus$ano),
                       pattern = "/\\d{4}$",
                       replacement = ""))

# Quantidade de NA em cada coluna.
apply(dus, MARGIN = 2, FUN = function(x) sum(is.na(x)))

## modelo    cor     km    ano  valor cambio  poten   trac    cat   novo 
##      0      0    132      0      0      0      0      0      0      0

#-----------------------------------------------------------------------
# Carrega o pacote `lattice` e `latticeExtra`.

library(lattice)
library(latticeExtra)

ls("package:lattice")

##   [1] "as.factorOrShingle"           "as.shingle"                  
##   [3] "axis.default"                 "banking"                     
##   [5] "barchart"                     "barley"                      
##   [7] "bwplot"                       "canonical.theme"             
##   [9] "cloud"                        "col.whitebg"                 
##  [11] "contourplot"                  "current.column"              
##  [13] "current.panel.limits"         "current.row"                 
##  [15] "densityplot"                  "diag.panel.splom"            
##  [17] "do.breaks"                    "dotplot"                     
##  [19] "draw.colorkey"                "draw.key"                    
##  [21] "environmental"                "equal.count"                 
##  [23] "ethanol"                      "histogram"                   
##  [25] "is.shingle"                   "larrows"                     
##  [27] "lattice.getOption"            "lattice.options"             
##  [29] "latticeParseFormula"          "level.colors"                
##  [31] "levelplot"                    "llines"                      
##  [33] "lplot.xy"                     "lpoints"                     
##  [35] "lpolygon"                     "lrect"                       
##  [37] "lsegments"                    "ltext"                       
##  [39] "ltransform3dMatrix"           "ltransform3dto3d"            
##  [41] "make.groups"                  "melanoma"                    
##  [43] "oneway"                       "packet.number"               
##  [45] "packet.panel.default"         "panel.3dscatter"             
##  [47] "panel.3dwire"                 "panel.abline"                
##  [49] "panel.arrows"                 "panel.average"               
##  [51] "panel.axis"                   "panel.barchart"              
##  [53] "panel.brush.splom"            "panel.bwplot"                
##  [55] "panel.cloud"                  "panel.contourplot"           
##  [57] "panel.curve"                  "panel.densityplot"           
##  [59] "panel.dotplot"                "panel.error"                 
##  [61] "panel.fill"                   "panel.grid"                  
##  [63] "panel.histogram"              "panel.identify"              
##  [65] "panel.identify.cloud"         "panel.identify.qqmath"       
##  [67] "panel.levelplot"              "panel.levelplot.raster"      
##  [69] "panel.linejoin"               "panel.lines"                 
##  [71] "panel.link.splom"             "panel.lmline"                
##  [73] "panel.loess"                  "panel.mathdensity"           
##  [75] "panel.number"                 "panel.pairs"                 
##  [77] "panel.parallel"               "panel.points"                
##  [79] "panel.polygon"                "panel.qq"                    
##  [81] "panel.qqmath"                 "panel.qqmathline"            
##  [83] "panel.rect"                   "panel.refline"               
##  [85] "panel.rug"                    "panel.segments"              
##  [87] "panel.smoothScatter"          "panel.spline"                
##  [89] "panel.splom"                  "panel.stripplot"             
##  [91] "panel.superpose"              "panel.superpose.2"           
##  [93] "panel.superpose.plain"        "panel.text"                  
##  [95] "panel.tmd.default"            "panel.tmd.qqmath"            
##  [97] "panel.violin"                 "panel.wireframe"             
##  [99] "panel.xyplot"                 "parallel"                    
## [101] "parallelplot"                 "prepanel.default.bwplot"     
## [103] "prepanel.default.cloud"       "prepanel.default.densityplot"
## [105] "prepanel.default.histogram"   "prepanel.default.levelplot"  
## [107] "prepanel.default.parallel"    "prepanel.default.qq"         
## [109] "prepanel.default.qqmath"      "prepanel.default.splom"      
## [111] "prepanel.default.xyplot"      "prepanel.lmline"             
## [113] "prepanel.loess"               "prepanel.qqmathline"         
## [115] "prepanel.spline"              "prepanel.tmd.default"        
## [117] "prepanel.tmd.qqmath"          "qq"                          
## [119] "qqmath"                       "rfs"                         
## [121] "Rows"                         "shingle"                     
## [123] "show.settings"                "simpleKey"                   
## [125] "simpleTheme"                  "singer"                      
## [127] "splom"                        "standard.theme"              
## [129] "strip.custom"                 "strip.default"               
## [131] "stripplot"                    "tmd"                         
## [133] "trellis.currentLayout"        "trellis.device"              
## [135] "trellis.focus"                "trellis.grobname"            
## [137] "trellis.last.object"          "trellis.panelArgs"           
## [139] "trellis.par.get"              "trellis.par.set"             
## [141] "trellis.switchFocus"          "trellis.unfocus"             
## [143] "trellis.vpname"               "USMortality"                 
## [145] "USRegionalMortality"          "which.packet"                
## [147] "wireframe"                    "xscale.components.default"   
## [149] "xyplot"                       "xyplot.ts"                   
## [151] "yscale.components.default"

ls("package:latticeExtra")

##  [1] "ancestry"                     "as.layer"                    
##  [3] "asTheEconomist"               "axis.grid"                   
##  [5] "biocAccess"                   "combineLimits"               
##  [7] "custom.theme"                 "custom.theme.2"              
##  [9] "dendrogramGrob"               "doubleYScale"                
## [11] "drawLayer"                    "EastAuClimate"               
## [13] "ecdfplot"                     "flattenPanel"                
## [15] "ggplot2like"                  "ggplot2like.opts"            
## [17] "glayer"                       "glayer_"                     
## [19] "gvhd10"                       "horizonplot"                 
## [21] "layer"                        "layer_"                      
## [23] "mapplot"                      "marginal.plot"               
## [25] "mergedTrellisLegendGrob"      "panel.2dsmoother"            
## [27] "panel.3dbars"                 "panel.3dpolygon"             
## [29] "panel.3dtext"                 "panel.ablineq"               
## [31] "panel.ecdfplot"               "panel.ellipse"               
## [33] "panel.horizonplot"            "panel.key"                   
## [35] "panel.levelplot.points"       "panel.lmlineq"               
## [37] "panel.mapplot"                "panel.qqmath.tails"          
## [39] "panel.qqmath.xyarea"          "panel.quantile"              
## [41] "panel.rootogram"              "panel.scaleArrow"            
## [43] "panel.segplot"                "panel.smoother"              
## [45] "panel.tskernel"               "panel.voronoi"               
## [47] "panel.xblocks"                "panel.xyarea"                
## [49] "postdoc"                      "prepanel.ecdfplot"           
## [51] "prepanel.horizonplot"         "prepanel.mapplot"            
## [53] "prepanel.rootogram"           "prepanel.segplot"            
## [55] "resizePanels"                 "rootogram"                   
## [57] "SeatacWeather"                "segplot"                     
## [59] "simpleSmoothTs"               "theEconomist.axis"           
## [61] "theEconomist.opts"            "theEconomist.theme"          
## [63] "tileplot"                     "USAge.df"                    
## [65] "USAge.table"                  "USCancerRates"               
## [67] "useOuterStrips"               "xscale.components.fractions" 
## [69] "xscale.components.log"        "xscale.components.log10.3"   
## [71] "xscale.components.log10ticks" "xscale.components.logpower"  
## [73] "xscale.components.subticks"   "xyplot.list"                 
## [75] "yscale.components.fractions"  "yscale.components.log"       
## [77] "yscale.components.log10.3"    "yscale.components.log10ticks"
## [79] "yscale.components.logpower"   "yscale.components.subticks"

funs <- c(ls("package:lattice"), ls("package:latticeExtra"))

# Funções de alto nível.
c(grep(x = funs, pattern = "^[^.]+plot$", value = TRUE),
  c("histogram", "wireframe", "splom", "qqmath"))

##  [1] "bwplot"       "contourplot"  "densityplot"  "dotplot"      "levelplot"   
##  [6] "parallelplot" "stripplot"    "xyplot"       "ecdfplot"     "horizonplot" 
## [11] "mapplot"      "segplot"      "tileplot"     "histogram"    "wireframe"   
## [16] "splom"        "qqmath"

# Funções de baixo nível.
grep(x = funs, pattern = "^panel\\.", value = TRUE)

##  [1] "panel.3dscatter"        "panel.3dwire"           "panel.abline"          
##  [4] "panel.arrows"           "panel.average"          "panel.axis"            
##  [7] "panel.barchart"         "panel.brush.splom"      "panel.bwplot"          
## [10] "panel.cloud"            "panel.contourplot"      "panel.curve"           
## [13] "panel.densityplot"      "panel.dotplot"          "panel.error"           
## [16] "panel.fill"             "panel.grid"             "panel.histogram"       
## [19] "panel.identify"         "panel.identify.cloud"   "panel.identify.qqmath" 
## [22] "panel.levelplot"        "panel.levelplot.raster" "panel.linejoin"        
## [25] "panel.lines"            "panel.link.splom"       "panel.lmline"          
## [28] "panel.loess"            "panel.mathdensity"      "panel.number"          
## [31] "panel.pairs"            "panel.parallel"         "panel.points"          
## [34] "panel.polygon"          "panel.qq"               "panel.qqmath"          
## [37] "panel.qqmathline"       "panel.rect"             "panel.refline"         
## [40] "panel.rug"              "panel.segments"         "panel.smoothScatter"   
## [43] "panel.spline"           "panel.splom"            "panel.stripplot"       
## [46] "panel.superpose"        "panel.superpose.2"      "panel.superpose.plain" 
## [49] "panel.text"             "panel.tmd.default"      "panel.tmd.qqmath"      
## [52] "panel.violin"           "panel.wireframe"        "panel.xyplot"          
## [55] "panel.2dsmoother"       "panel.3dbars"           "panel.3dpolygon"       
## [58] "panel.3dtext"           "panel.ablineq"          "panel.ecdfplot"        
## [61] "panel.ellipse"          "panel.horizonplot"      "panel.key"             
## [64] "panel.levelplot.points" "panel.lmlineq"          "panel.mapplot"         
## [67] "panel.qqmath.tails"     "panel.qqmath.xyarea"    "panel.quantile"        
## [70] "panel.rootogram"        "panel.scaleArrow"       "panel.segplot"         
## [73] "panel.smoother"         "panel.tskernel"         "panel.voronoi"         
## [76] "panel.xblocks"          "panel.xyarea"

# Lista de itens que podem ser customizados.
str(trellis.par.get())

## List of 35
##  $ grid.pars        : list()
##  $ fontsize         :List of 2
##   ..$ text  : num 12
##   ..$ points: num 8
##  $ background       :List of 2
##   ..$ alpha: num 1
##   ..$ col  : chr "transparent"
##  $ panel.background :List of 1
##   ..$ col: chr "transparent"
##  $ clip             :List of 2
##   ..$ panel: chr "on"
##   ..$ strip: chr "on"
##  $ add.line         :List of 4
##   ..$ alpha: num 1
##   ..$ col  : chr "#000000"
##   ..$ lty  : num 1
##   ..$ lwd  : num 1
##  $ add.text         :List of 5
##   ..$ alpha     : num 1
##   ..$ cex       : num 1
##   ..$ col       : chr "#000000"
##   ..$ font      : num 1
##   ..$ lineheight: num 1.2
##  $ plot.polygon     :List of 5
##   ..$ alpha : num 1
##   ..$ col   : chr "#00ffff"
##   ..$ border: chr "black"
##   ..$ lty   : num 1
##   ..$ lwd   : num 1
##  $ box.dot          :List of 5
##   ..$ alpha: num 1
##   ..$ col  : chr "#000000"
##   ..$ cex  : num 1
##   ..$ font : num 1
##   ..$ pch  : num 16
##  $ box.rectangle    :List of 5
##   ..$ alpha: num 1
##   ..$ col  : chr "#0080ff"
##   ..$ fill : chr "transparent"
##   ..$ lty  : num 1
##   ..$ lwd  : num 1
##  $ box.umbrella     :List of 4
##   ..$ alpha: num 1
##   ..$ col  : chr "#0080ff"
##   ..$ lty  : num 2
##   ..$ lwd  : num 1
##  $ dot.line         :List of 4
##   ..$ alpha: num 1
##   ..$ col  : chr "#e6e6e6"
##   ..$ lty  : num 1
##   ..$ lwd  : num 1
##  $ dot.symbol       :List of 5
##   ..$ alpha: num 1
##   ..$ cex  : num 0.8
##   ..$ col  : chr "#0080ff"
##   ..$ font : num 1
##   ..$ pch  : num 16
##  $ plot.line        :List of 4
##   ..$ alpha: num 1
##   ..$ col  : chr "#0080ff"
##   ..$ lty  : num 1
##   ..$ lwd  : num 1
##  $ plot.symbol      :List of 6
##   ..$ alpha: num 1
##   ..$ cex  : num 0.8
##   ..$ col  : chr "#0080ff"
##   ..$ font : num 1
##   ..$ pch  : num 1
##   ..$ fill : chr "transparent"
##  $ reference.line   :List of 4
##   ..$ alpha: num 1
##   ..$ col  : chr "#e6e6e6"
##   ..$ lty  : num 1
##   ..$ lwd  : num 1
##  $ strip.background :List of 2
##   ..$ alpha: num 1
##   ..$ col  : chr [1:7] "#ffe5cc" "#ccffcc" "#ccffff" "#cce6ff" ...
##  $ strip.shingle    :List of 2
##   ..$ alpha: num 1
##   ..$ col  : chr [1:7] "#ff7f00" "#00ff00" "#00ffff" "#0080ff" ...
##  $ strip.border     :List of 4
##   ..$ alpha: num 1
##   ..$ col  : chr [1:7] "#000000" "#000000" "#000000" "#000000" ...
##   ..$ lty  : num [1:7] 1 1 1 1 1 1 1
##   ..$ lwd  : num [1:7] 1 1 1 1 1 1 1
##  $ superpose.line   :List of 4
##   ..$ alpha: num 1
##   ..$ col  : chr [1:7] "#0080ff" "#ff00ff" "darkgreen" "#ff0000" ...
##   ..$ lty  : num [1:7] 1 1 1 1 1 1 1
##   ..$ lwd  : num [1:7] 1 1 1 1 1 1 1
##  $ superpose.symbol :List of 6
##   ..$ alpha: num [1:7] 1 1 1 1 1 1 1
##   ..$ cex  : num [1:7] 0.8 0.8 0.8 0.8 0.8 0.8 0.8
##   ..$ col  : chr [1:7] "#0080ff" "#ff00ff" "darkgreen" "#ff0000" ...
##   ..$ fill : chr [1:7] "#CCFFFF" "#FFCCFF" "#CCFFCC" "#FFE5CC" ...
##   ..$ font : num [1:7] 1 1 1 1 1 1 1
##   ..$ pch  : num [1:7] 1 1 1 1 1 1 1
##  $ superpose.polygon:List of 5
##   ..$ alpha : num [1:7] 1 1 1 1 1 1 1
##   ..$ col   : chr [1:7] "#CCFFFF" "#FFCCFF" "#CCFFCC" "#FFE5CC" ...
##   ..$ border: chr [1:7] "black" "black" "black" "black" ...
##   ..$ lty   : num [1:7] 1 1 1 1 1 1 1
##   ..$ lwd   : num [1:7] 1 1 1 1 1 1 1
##  $ regions          :List of 2
##   ..$ alpha: num 1
##   ..$ col  : chr [1:100] "#FF80FFFF" "#FF82FFFF" "#FF85FFFF" "#FF87FFFF" ...
##  $ shade.colors     :List of 2
##   ..$ alpha  : num 1
##   ..$ palette:function (irr, ref, height, saturation = 0.9)  
##  $ axis.line        :List of 4
##   ..$ alpha: num 1
##   ..$ col  : chr "#000000"
##   ..$ lty  : num 1
##   ..$ lwd  : num 1
##  $ axis.text        :List of 5
##   ..$ alpha     : num 1
##   ..$ cex       : num 0.8
##   ..$ col       : chr "#000000"
##   ..$ font      : num 1
##   ..$ lineheight: num 1
##  $ axis.components  :List of 4
##   ..$ left  :List of 3
##   .. ..$ tck : num 1
##   .. ..$ pad1: num 1
##   .. ..$ pad2: num 1
##   ..$ top   :List of 3
##   .. ..$ tck : num 1
##   .. ..$ pad1: num 1
##   .. ..$ pad2: num 1
##   ..$ right :List of 3
##   .. ..$ tck : num 1
##   .. ..$ pad1: num 1
##   .. ..$ pad2: num 1
##   ..$ bottom:List of 3
##   .. ..$ tck : num 1
##   .. ..$ pad1: num 1
##   .. ..$ pad2: num 1
##  $ layout.heights   :List of 19
##   ..$ top.padding      : num 1
##   ..$ main             : num 1
##   ..$ main.key.padding : num 1
##   ..$ key.top          : num 1
##   ..$ xlab.top         : num 1
##   ..$ key.axis.padding : num 1
##   ..$ axis.top         : num 1
##   ..$ strip            : num 1
##   ..$ panel            : num 1
##   ..$ axis.panel       : num 1
##   ..$ between          : num 1
##   ..$ axis.bottom      : num 1
##   ..$ axis.xlab.padding: num 1
##   ..$ xlab             : num 1
##   ..$ xlab.key.padding : num 0
##   ..$ key.bottom       : num 1
##   ..$ key.sub.padding  : num 1
##   ..$ sub              : num 1
##   ..$ bottom.padding   : num 1
##  $ layout.widths    :List of 15
##   ..$ left.padding     : num 1
##   ..$ key.left         : num 1
##   ..$ key.ylab.padding : num 0
##   ..$ ylab             : num 1
##   ..$ ylab.axis.padding: num 1
##   ..$ axis.left        : num 1
##   ..$ axis.panel       : num 1
##   ..$ strip.left       : num 1
##   ..$ panel            : num 1
##   ..$ between          : num 1
##   ..$ axis.right       : num 1
##   ..$ axis.key.padding : num 1
##   ..$ ylab.right       : num 1
##   ..$ key.right        : num 1
##   ..$ right.padding    : num 1
##  $ box.3d           :List of 4
##   ..$ alpha: num 1
##   ..$ col  : chr "#000000"
##   ..$ lty  : num 1
##   ..$ lwd  : num 1
##  $ par.xlab.text    :List of 5
##   ..$ alpha     : num 1
##   ..$ cex       : num 1
##   ..$ col       : chr "#000000"
##   ..$ font      : num 1
##   ..$ lineheight: num 1
##  $ par.ylab.text    :List of 5
##   ..$ alpha     : num 1
##   ..$ cex       : num 1
##   ..$ col       : chr "#000000"
##   ..$ font      : num 1
##   ..$ lineheight: num 1
##  $ par.zlab.text    :List of 5
##   ..$ alpha     : num 1
##   ..$ cex       : num 1
##   ..$ col       : chr "#000000"
##   ..$ font      : num 1
##   ..$ lineheight: num 1
##  $ par.main.text    :List of 5
##   ..$ alpha     : num 1
##   ..$ cex       : num 1.2
##   ..$ col       : chr "#000000"
##   ..$ font      : num 2
##   ..$ lineheight: num 1
##  $ par.sub.text     :List of 5
##   ..$ alpha     : num 1
##   ..$ cex       : num 1
##   ..$ col       : chr "#000000"
##   ..$ font      : num 2
##   ..$ lineheight: num 1

2 Gráficos para uma variável categórica

#-----------------------------------------------------------------------
# Gráfico de barras e setores.

# Tabela de frequência.
x <- xtabs(~cambio, data = dus)
class(x)

## [1] "xtabs" "table"

# Gráfico padrão.
barchart(x,
         horizontal = FALSE,
         xlab = "Tipo de câmbio",
         ylab = "Frequência absoluta")

# A lattice também trabalha com fórmulas.
tb <- as.data.frame(x)

# CUIDADO: a origem desse gráfico não é o zero!
barchart(Freq ~ cambio,
         data = tb,
         xlab = "Tipo de câmbio",
         ylab = "Frequência absoluta",
         col = c("seagreen", "yellowgreen"))

barchart(Freq ~ cambio,
         data = tb,
         origin = 0, # <-- ATENÇÃO.
         xlab = "Tipo de câmbio",
         ylab = "Frequência absoluta",
         col = c("seagreen", "yellowgreen"))

Deve-se ter cuidado com a refência das barras. O argumento origin = permite especificar a origem.

# Para as cores do carro.
x <- xtabs(~cor, data = dus)
levels(dus$cor)

## NULL

tb <- as.data.frame(x)
tb

##          cor Freq
## 1       Azul   14
## 2     Branco  193
## 3      Cinza   30
## 4 Indefinida    7
## 5      Prata  223
## 6      Preto  155
## 7      Verde   58
## 8   Vermelho   18
## 9      Vinho    1

# Vetor que cores para usar com cada cor de veiculo.
cols <- c("blue", "white", "gray50", "Yellow", "gray90", "black",
          "green4", "red", "red4")

cbind(levels(dus$cor), cols)

##       cols    
##  [1,] "blue"  
##  [2,] "white" 
##  [3,] "gray50"
##  [4,] "Yellow"
##  [5,] "gray90"
##  [6,] "black" 
##  [7,] "green4"
##  [8,] "red"   
##  [9,] "red4"

# Com anotação sobre as barras.
barchart(Freq ~ reorder(cor, -Freq),
         data = tb,
         origin = 0,
         col = cols,
         xlab = "Cores",
         ylab = "Frequência absoluta") +
    layer({
        panel.text(x = x, y = y, pos = 3, label = y)
    })

3 Gráfico para duas variáveis categóricas

#-----------------------------------------------------------------------
# Gráficos de barras emplilhadas (stacked) e lado a lado.

x <- xtabs(~cambio + ano, data = dus)
x

##             ano
## cambio       2011 2012 2013 2014
##   AUTOMÁTICO   51   40  100    4
##   MANUAL       88  204  180   32

cols <- c("#660d32", "#bc1a5e")

# Barras empilhadas.
barchart(x,
         xlab = "Ano",
         stack = TRUE,
         auto.key = TRUE,
         ylab = "Frequência absoluta")

# Barras lado a lado.
barchart(x,
         xlab = "Ano",
         stack = FALSE,
         horizontal = FALSE,
         auto.key = list(space = "right"),
         ylab = "Frequência absoluta")

tb <- as.data.frame(x)
tb

##       cambio  ano Freq
## 1 AUTOMÁTICO 2011   51
## 2     MANUAL 2011   88
## 3 AUTOMÁTICO 2012   40
## 4     MANUAL 2012  204
## 5 AUTOMÁTICO 2013  100
## 6     MANUAL 2013  180
## 7 AUTOMÁTICO 2014    4
## 8     MANUAL 2014   32

barchart(Freq ~ ano,
         groups = cambio,
         auto.key = list(space = "right", title = "Ano"),
         data = tb,
         origin = 0,
         xlab = "Câmbio",
         ylab = "Frequência absoluta")

cols <- c("#07B1FF", "#0CE8BB", "#00FF40", "#BEE80C")
barchart(Freq ~ cambio,
         groups = ano,
         auto.key = list(space = "right", title = "Ano"),
         data = tb,
         origin = 0,
         xlab = "Ano",
         ylab = "Frequência absoluta",
         par.settings = list(superpose.polygon = list(col = cols)))

#-----------------------------------------------------------------------
# Anotações nas barras.

cols <- c("#07B1FF", "#0CE8BB", "#00FF40", "#BEE80C")
barchart(Freq ~ cambio,
         groups = ano,
         auto.key = list(corner = c(0, 1), title = "Ano"),
         data = tb,
         origin = 0,
         xlab = "Ano",
         ylab = "Frequência absoluta",
         par.settings = list(superpose.polygon = list(col = cols))) +
    layer({
        # print(ls()) # Descomente para ver os objetos do escopo.
        g <- scale(as.integer(groups)/(nlevels(groups) + box.ratio),
                   scale = FALSE)
        print(g)
        panel.text(x = as.integer(x) + g, y = y, pos = 3, label = y)
    }) +
    layer({
        panel.abline(h = seq(0, 200, by = 25),
                     lty = 2,
                     col = "gray50")
    }, under = TRUE)

##             [,1]
## [1,] -0.25000000
## [2,] -0.25000000
## [3,] -0.08333333
## [4,] -0.08333333
## [5,]  0.08333333
## [6,]  0.08333333
## [7,]  0.25000000
## [8,]  0.25000000
## attr(,"scaled:center")
## [1] 0.4166667

# Mais dimensões.
x <- xtabs(~cambio + ano + cor, data = dus)

tb <- as.data.frame(x)
head(tb)

##       cambio  ano  cor Freq
## 1 AUTOMÁTICO 2011 Azul    4
## 2     MANUAL 2011 Azul    3
## 3 AUTOMÁTICO 2012 Azul    0
## 4     MANUAL 2012 Azul    6
## 5 AUTOMÁTICO 2013 Azul    0
## 6     MANUAL 2013 Azul    1

barchart(Freq ~ ano | cor,
         groups = cambio,
         auto.key = list(columns = 2, title = "Câmbio", cex.title = 1.1),
         data = tb,
         origin = 0,
         xlab = "Ano",
         ylab = "Frequência absoluta",
         par.settings = list(superpose.polygon = list(col = cols)))

barchart(Freq ~ ano | cor,
         groups = cambio,
         auto.key = list(columns = 2, title = "Câmbio", cex.title = 1.1),
         data = tb,
         origin = 0,
         xlab = "Ano",
         ylab = "Frequência absoluta",
         scales = list(y = list(relation = "free")),
         par.settings = list(superpose.polygon = list(col = cols)))

4 Gráficos para uma variável contínua

#-----------------------------------------------------------------------
# Histograma.

# Gráfico básico.
# hist(dus$valor)

histogram(~valor,
          data = dus,
          xlab = "Preço de venda (R$)",
          ylab = "Frequência relativa",
          col = "orange") +
    layer(panel.rug(x, col = 1))

# Se breaks é um escalar então entende-se que é uma *sugestão* para o
# número de clases.
histogram(~valor,
          data = dus,
          breaks = 15,
          type = c("percent", "count", "density")[2],
          xlab = "Preço de venda (R$)",
          ylab = "Frequência absoluta",
          col = "orange") +
    layer(panel.rug(x, col = 1))

# Se breaks é um vetor então entende-se que são os limites para
# classificação dos valores.
histogram(~valor,
          data = dus,
          breaks = seq(35000, 75000, by = 2500),
          type = c("percent", "count", "density")[2],
          xlab = "Preço de venda (R$)",
          ylab = "Frequência absoluta",
          col = "orange") +
    layer(panel.rug(x, col = 1))

# Gráfico onde a altura é a densidade e não a frequência.
histogram(~valor,
          data = dus,
          type = c("percent", "count", "density")[3],
          xlab = "Preço de venda (R$)",
          ylab = "Frequência absoluta",
          col = "#ba6dff") +
    layer(panel.rug(x, col = 1))

#-----------------------------------------------------------------------
# Gráficos de densidade.

densityplot(~valor,
            data = dus,
            kernel = "rectangular",
            bw = 2000,
            xlab = "Preço de venda (R$)",
            ylab = "Frequência absoluta",
            # plot.points = c("jitter", "rug")[2],
            plot.points = FALSE,
            col = "#ba6dff") +
    layer({
        panel.rug(x, col = 1)
    })

densityplot(~valor,
            groups = cambio,
            auto.key = list(corner = c(0.95, 0.95)),
            data = dus,
            xlab = "Preço de venda (R$)",
            ylab = "Frequência absoluta",
            plot.points = "jitter")

densityplot(~valor | ano,
            groups = cambio,
            auto.key = list(corner = c(0.05, 0.95)),
            data = dus,
            xlab = "Preço de venda (R$)",
            ylab = "Frequência absoluta",
            plot.points = "jitter")

#-----------------------------------------------------------------------
# Gráfico de frequência acumulada empírica.

# Gráfico onde a altura é a densidade e não a frequência.
ecdfplot(~valor,
         data = dus,
         xlab = "Preço de venda (R$)",
         ylab = "Frequência acumulada",
         col = "#ba6dff") +
    layer({
        panel.rug(x, col = 1)
    })

ecdfplot(~valor,
         groups = cambio,
         auto.key = list(corner = c(0.05, 0.8)),
         data = dus,
         xlab = "Preço de venda (R$)",
         ylab = "Frequência acumulada")

ecdfplot(~valor | ano,
         groups = cambio,
         auto.key = list(corner = c(0.05, 0.85)),
         data = dus,
         xlab = "Preço de venda (R$)",
         ylab = "Frequência acumulada")

5 Gráfico para duas variáveis contínuas

#-----------------------------------------------------------------------
# Diagrama de dispersão.

# Diagrama de dispersão básico.
xyplot(valor ~ km, data = dus)

# Adicionar uma linha de tendência suave.
xyplot(valor ~ km,
       data = dus,
       type = c("p", "smooth"))

xyplot(valor ~ km, data = dus) +
    layer({
        panel.smoother(x, y)
    })

# Usar cores diferentes para identificar o tipo de câmbio.
xyplot(valor ~ km,
       groups = cambio,
       auto.key = TRUE,
       data = dus)

# Adicionar uma linha de tendência suave.
xyplot(valor ~ km,
       groups = cambio,
       auto.key = TRUE,
       data = dus) +
    glayer({
        panel.smoother(x = x, y = y, col = col.line, span = 0.9)
    })

n <- head(names(sort(table(dus$cor), decreasing = TRUE)), n = 4)

xyplot(valor ~ km | cor,
       groups = cambio,
       auto.key = TRUE,
       data = subset(dus, cor %in% n)) +
    glayer({
        panel.smoother(x = x, y = y, col = col.line, span = 0.9)
    })

6 Gráfico para variável contínua e categórica

#-----------------------------------------------------------------------

# Preço em função dos anos.
bwplot(valor ~ ano, data = dus, pch = "|")

# Edita níveis do fator.
levels(dus$cat)

## NULL

levels(dus$cat) <- trimws(levels(dus$cat))
dus2 <- droplevels(subset(dus, cat != ""))
dus2 <- transform(dus2, valor = valor/1000, km = km/1000)

names(trellis.par.get())

##  [1] "grid.pars"         "fontsize"          "background"        "panel.background" 
##  [5] "clip"              "add.line"          "add.text"          "plot.polygon"     
##  [9] "box.dot"           "box.rectangle"     "box.umbrella"      "dot.line"         
## [13] "dot.symbol"        "plot.line"         "plot.symbol"       "reference.line"   
## [17] "strip.background"  "strip.shingle"     "strip.border"      "superpose.line"   
## [21] "superpose.symbol"  "superpose.polygon" "regions"           "shade.colors"     
## [25] "axis.line"         "axis.text"         "axis.components"   "layout.heights"   
## [29] "layout.widths"     "box.3d"            "par.xlab.text"     "par.ylab.text"    
## [33] "par.zlab.text"     "par.main.text"     "par.sub.text"

str(trellis.par.get()$box.dot)

## List of 5
##  $ alpha: num 1
##  $ col  : chr "#000000"
##  $ cex  : num 1
##  $ font : num 1
##  $ pch  : num 16

str(trellis.par.get()$box.rectangle)

## List of 5
##  $ alpha: num 1
##  $ col  : chr "#0080ff"
##  $ fill : chr "transparent"
##  $ lty  : num 1
##  $ lwd  : num 1

str(trellis.par.get()$box.umbrella)

## List of 4
##  $ alpha: num 1
##  $ col  : chr "#0080ff"
##  $ lty  : num 2
##  $ lwd  : num 1

# Customização do gráfico.
bw_ps <- list(plot.symbol = list(col = 1),
              box.dot = list(pch = "|"),
              box.umbrella = list(col = 1, lty = 1),
              box.rectangle = list(col = 1, fill = "seagreen"))

bwplot(valor ~ cat,
       data = dus2,
       xlab = "Modelo",
       ylab = "Preço de venda (R$)",
       par.settings = bw_ps)

bwplot(valor ~ cat | ano,
       data = dus2,
       xlab = "Modelo",
       ylab = "Preço de venda (R$)",
       as.table = TRUE,
       par.settings = bw_ps)

# Indicação do valor da média.
bwplot(valor ~ cat,
       data = dus2,
       xlab = "Modelo",
       ylab = "Preço de venda (R$)",
       par.settings = bw_ps) +
    layer({
        # panel.xyplot(x = x, y = y, jitter.x = TRUE)
        m <- aggregate(y ~ x, FUN = mean)
        panel.points(x = m$x, y = m$y, col = "red", pch = 19)
    })

# Indicação do valor da média.
bwplot(valor ~ cat | cambio,
       data = dus2,
       xlab = "Modelo",
       ylab = "Preço de venda (R$)",
       par.settings = bw_ps) +
    layer({
        m <- aggregate(y ~ x, FUN = mean)
        panel.points(x = m$x, y = m$y, col = "red", pch = 19)
        fiv <- by(y, x, FUN = fivenum, simplify = FALSE)
        fiv <- do.call(cbind, fiv)
        aiq <- fiv[4, ] - fiv[2, ]
        l <- fiv[2, ] - 1.5 * aiq
        u <- fiv[4, ] + 1.5 * aiq
        i <- seq_along(u)
        panel.segments(x0 = i - 0.5, x1 = i + 0.5, y0 = l, y1 = l,
                       col = "gray50", lty = 3)
        panel.segments(x0 = i - 0.5, x1 = i + 0.5, y0 = u, y1 = u,
                       col = "gray50", lty = 3)
    })