= The ggplot2 Package =
 * ggplot2 is - like lattice based on the grid graphics system (Paul Murrell)
 * graphics and parts of graphics are objects and they are manipulable

== Structure of a ggplot Object ==
begin with an empty object to see the structure:

{{{#!highlight r
> po <- ggplot()
> summary(po)
data: [x]
faceting: facet_null()
}}}
 * what we see are empty place holders
 * when we use str() to explore the structure of the object we see that it is a list with length 9

{{{#!highlight r
> str(po)
List of 9
  List of 9
  $ data       : list()
  ..- attr(*, "class")= chr "waiver"
  $ layers     : list()
  $ scales     :Reference class 'Scales' [package "ggplot2"] with 1 fields
  ..$ scales: NULL
  ..and 21 methods, of which 9 are possibly relevant:
  ..  add, clone, find, get_scales, has_scale, initialize, input, n,
  ..  non_position_scales
  $ mapping    : list()
  $ theme      : list()
  $ coordinates:List of 1
  ..$ limits:List of 2
  .. ..$ x: NULL
  .. ..$ y: NULL
  ..- attr(*, "class")= chr [1:2] "cartesian" "coord"
  $ facet      :List of 1
  ..$ shrink: logi TRUE
  ..- attr(*, "class")= chr [1:2] "null" "facet"
  $ plot_env   :<environment: R_GlobalEnv>
  $ labels     : list()
  - attr(*, "class")= chr [1:2] "gg" "ggplot"
}}}
== Structure of a ggplot Object ==
Now we fill this structure - first the three main steps:

 * the first argument to ggplot is data
 * then specify what graphics shapes you are going to use to view the data (e.g. geom_line() or geom_point()).
 * specify what features (or aesthetics) will be used (e.g. what variables will determine x- and y-locations) with the aes() function
 * if these aesthetics are intented to be used in all layers it is more convenient to specify them in the ggplot object

== Feed the Object ==
 * first we create a little sample data frame\small

{{{#!highlight r
> x1 <- 1:10; y1 <- 1:10; z1 <- 10:1
> l1 <- LETTERS[1:10]
> a <- 10; b <- (0:-9)/10:1
> ex <- data.frame(x=x1,y=y1,z=z1,l=l1,a=a,b=b)
> ex
x  y  z l  a          b
1   1  1 10 A 10  0.0000000
2   2  2  9 B 10 -0.1111111
3   3  3  8 C 10 -0.2500000
4   4  4  7 D 10 -0.4285714
5   5  5  6 E 10 -0.6666667
6   6  6  5 F 10 -1.0000000
7   7  7  4 G 10 -1.5000000
8   8  8  3 H 10 -2.3333333
9   9  9  2 I 10 -4.0000000
10 10 10  1 J 10 -9.0000000
}}}
 * then we create a ggplot object containing the data and some standard aesthetics (here we define the x and the y positions)
 * add one or more geoms, we begin with  geom_point

{{{#!highlight r
> po <- ggplot(ex,aes(x=x1,y=y1))
> summary(po)
> p1 <- po + geom_point()
}}}

[[ggp1.pdf]]

== Layers ==
 * ggplot() creates an object - every "+" adds something to this object (change the object)
 * the default method of ggplot() is print(), which creates the plot
 * it is better to store the object - so you can change it (e.g. you can change the data frame)

== Layers ==
 * so we add another layer, which adds a label to the points (use geom\_text)

{{{#!highlight r
> p2 <- po +
+    geom_point() +
+        geom_text(aes(label=l), hjust=1.1, vjust=-0.2)
> p2
}}}
 * aes(label=l) maps the l variable to the label aesthetic, and hjust and vjust define where our labels are placed

[[ggp2.pdf]]

== Layers ==
 * imagine you have worked a little time on a plot - and then you detect a mistake in your data, so the ''real'' data frame looks different
 * so you can replace the old, wrong data by the new  data (using \%+\% \footnotesize

{{{#!highlight r
> ## the new data
> ex2 <- data.frame(x1=sample(1:20),
+                   y1=sample(1:10),
+                   l=letters[1:20])
> head(ex2,10)
x1 y1 l
1   3  6 a
2   6  2 b
3  14  1 c
4  19 10 d
5  12  4 e
6  15  8 f
7  20  5 g
8  17  7 h
9  13  3 i
10 16  9 j
}}}
<img alt='sesssion2/ggp3.pdf' src='-1' />

== Layers ==
{{{#!highlight r
> p2 %+% ex2
}}}
<img alt='sesssion2/ggp3.pdf' src='-1' />

== Layers ==
 * by using the line geom you can join the points (we use the new data)

{{{#!highlight r
> pn <- p %+% ex2 ## replace data in p
> pn + geom_line()
}}}
<img alt='sesssion2/ggp4.pdf' src='-1' />

== Layers ==
 * you can also join the points in the order of the data fram by using the path geom instead\footnotesize

{{{#!highlight r
> my.text <- geom_text(aes(label=l),
+                          hjust=1.1,
+                          vjust=-0.2)
> pn + geom_path() + my.text
}}}
<img alt='sesssion2/ggp5.pdf' src='-1' />

== Layers ==
Adding extra lines:

 * there are three geoms: abline, vline, hline
 * abline adds one or more lines with specified slope and intercept to the plot\footnotesize

{{{#!highlight r
> ## one line
> p + geom_abline(intercept=10,slope=-1,
+                          colour=rgb(.5,.5,.9))
> ## two lines
> p + geom_abline(intercept=c(10,9),slope=c(-1,-2),
+                              colour=rgb(.5,.5,.9))
> more lines
> p + geom_abline(intercept=10:1,slope=-(10:1)/10,
}}}
<img alt='sesssion2/ggp6.png' src='-1' />

== Layers ==
 * adding lines referring to the data frame

{{{#!highlight r
> p1 +
+   geom_abline(aes(slope=b,intercept=a,colour=x1)) +
+   scale_x_continuous(limits=c(0,10))
}}}
<img alt='sesssion2/ggp7.pdf' src='-1' />

== Layers ==
 * the same works for the hline and the vline geom which add horizonal and vertical line(s)
 * argument: yintercept, xintercept respectively
 * setting and mapping are possible

{{{#!highlight r
> p1 + geom_hline(yintercept=1:10)
> p1 + geom_hline(yintercept=1:10) +
+     geom_vline(xintercept=1:10)
}}}
<img alt='sesssion2/ggp8.pdf' src='-1' />

== Other Common Layers ==
 * some other layers for 1 continuous variable:
  * geom\_boxplot()
  * geom\_histogram()
  * geom\_density()
 * some other layers for 1 discrete variable:
  * geom\_bar()
 * some other layers for 2 or more continuous variables:
  * geom\_smooth()
  * geom\_density2d()
  * geom\_contour()
  * geom\_quantile()

== Exercises ==
 * use our data frame or load it: \texttt{load("20150310data.rdata")}
 * create a new variable EC1 containing the first 2 letters of the Event.Code column, use the function str\_sub() from the stringr package (type \texttt{?str\_sub}  to get help)

== Exercises ==
{{{#!highlight r
> data$EC1 <- factor(str_sub(data$Event.Code,1,2))
> head(data)
Subject Sex Age_PRETEST Trial Event.Type Code   Time TTime Uncertainty
1       1   f        3.11     7   Response    2 103745  2575           1
2       1   f        3.11    12   Response    2 156493  2737           1
3       1   f        3.11    17   Response    2 214772  6630           1
4       1   f        3.11    22   Response    1 262086  5957           1
5       1   f        3.11    27   Response    2 302589   272           1
6       1   f        3.11    32   Response    1 352703  7197           1
Duration Uncertainty.1 ReqTime ReqDur Stim.Type Pair.Index    Type Event.Code
1     2599             3       0   next       hit          7 Picture   RO26.jpg
2     2800             2       0   next incorrect         12 Picture   RO19.jpg
3     6798             2       0   next       hit         17 Picture   RS23.jpg
4     5999             2       0   next incorrect         22 Picture   OF22.jpg
5      400             2       0   next       hit         27 Picture   AT08.jpg
6     7398             2       0   next       hit         32 Picture   AT30.jpg
testid EC1
1  test2  RO
2  test2  RO
3  test2  RS
4  test2  OF
5  test2  AT
6  test2  AT
}}}
== Exercises ==
Create the five plots and save them into a file. * create a plot using ggplot, map the variable EC1 to x and use geom\_bar() * now to the plot again, but add another aesthetic: fill (colour of the filling); map fill to Stim.Type * add the position argument to geom\_bar(), set it to "fill" * now add \texttt{facet\_wrap(~testid)} to show the same graph per time * make a graph facetted per child showing stacked hit/incorrect bars with time on the x axis

== Exercises ==
 * create a plot using ggplot, map the variable EC1 to x and use geom\_bar()

{{{#!highlight r
> ggplot(data,aes(x=EC1)) +
+     geom_bar()
>
> ggsave("plot1.png")
Saving 16 x 9.13 in image
}}}
<img alt='sesssion2/plot1.png' src='-1' />

== Exercises ==
 * now to the plot again, but add another aesthetic: fill (colour of the filling); map fill to Stim.Type

{{{#!highlight r
> ggplot(data,aes(x=EC1,fill=Stim.Type)) +
+     geom_bar()
>
> ggsave("plot2.png")
Saving 16 x 9.13 in image
}}}
<img alt='sesssion2/plot2.png' src='-1' />

== Exercises ==
 * add the position argument to geom\_bar(), set it to "fill"

{{{#!highlight r
> ggplot(data,aes(x=EC1,fill=Stim.Type)) +
+     geom_bar(position = "fill")
>
> ggsave("plot3.png")
Saving 16 x 9.13 in image
}}}
<img alt='sesssion2/plot3.png' src='-1' />

== Exercises ==
 * now add \texttt{facet\_wrap(~testid)} to show the same graph per time

{{{#!highlight r
> ggplot(data,aes(x=EC1,fill=Stim.Type)) +
+     geom_bar(position = "fill") +
+     facet_wrap(~testid)
>
> ggsave("plot4.png")
Saving 16 x 9.13 in image
}}}
<img alt='sesssion2/plot4.png' src='-1' />

== Exercises ==
 * now add \texttt{facet\_wrap(~testid)} to show the same graph per time

{{{#!highlight r
> ggplot(data,aes(x=EC1,fill=Stim.Type)) +
+     geom_bar(position = "fill") +
+     facet_wrap(~testid,scales = "free")
>
> ggsave("plot4a.png")
Saving 16 x 9.13 in image
}}}
<img alt='sesssion2/plot4a.png' src='-1' />

== Exercises ==
 * make a graph facetted per child showing stacked hit/incorrect bars with time on the x axis

{{{#!highlight r
> ggplot(data,aes(x=testid,fill=Stim.Type)) +
+     geom_bar(position = "fill") +
+     facet_wrap(~ Subject)
> ggsave("plot5.png")
Saving 16 x 9.13 in image
}}}
<img alt='sesssion2/plot5.png' src='-1' />

== Introduction ==
The dplyr package makes each of these steps as fast and easy as possible by:

 * Elucidating the most common data manipulation operations, so that your options are helpfully constrained when thinking about how to tackle a problem.
 * Providing simple functions that correspond to the most common data manipulation verbs, so that you can easily translate your thoughts into code.
 * Using efficient data storage backends, so that you spend as little time waiting for the computer as possible.

== Scales ==
What if we want to change the colours?

 * this leads to another important type of component not yet mentioned
 * if you map a variable to a aesthetic is these done in a default way, in this case some reddish colour is mapped to hit while light blue is mapped to incorrect; in addition a discrete range of colours is automatically used
 * these rules of mapping are called scales
 * different type of scales exists for the axes, colours, shapes etc, some of them exists in discrete and continuous versions, some in just one of them (in general one can say, everytime there can be a legend there is a scale)
 * the name convention: scale\_aesthetic\_specification. for example scale\_x\_discrete for customizing a discrete x axis (e.g. in barplots)

== Changing a Scale ==
 * to change our discrete colour scale for the filling we type \footnotesize

{{{#!highlight r
> ggplot(data,aes(x=EC1,fill=Stim.Type)) +
+     geom_bar(position = "fill") +
+     facet_wrap(~testid,scales = "free") +
+     scale_fill_manual(values=c("forestgreen","firebrick"))
}}}
<img alt='sesssion2/ggp10.png' src='-1' />

== Changing a Scale ==
There are other ways to customize a discrete colour/fill scales

 * scale\_colour\_grey()
 * scale\_colour\_hue()
 * scale\_colour\_brewer()

== Changing a Scale ==
{{{#!highlight r
> ggplot(data,aes(x=EC1,fill=Stim.Type)) +
+     geom_bar(position = "fill") +
+     facet_wrap(~testid,scales = "free") +
+         scale_fill_grey()
> ggplot(data,aes(x=EC1,fill=Stim.Type)) +
+     geom_bar(position = "fill") +
+     facet_wrap(~testid,scales = "free") +
+         scale_fill_hue(h=c(180,360))
> ggplot(data,aes(x=EC1,fill=Stim.Type)) +
+     geom_bar(position = "fill") +
+     facet_wrap(~testid,scales = "free") +
+     scale_fill_brewer(type = "div",palette = 2)
}}}