diff --git a/.RData b/.RData
index 665926dd757368411f85a2cc293a10cf916bad63..aeec7646785df75aa52fe3c36aaef2217f8d1d13 100644
Binary files a/.RData and b/.RData differ
diff --git a/DESCRIPTION b/DESCRIPTION
index 4092c8312a7ce2fc2c03c52b8ae8d2b62d3a2062..00310e3410a9391f7111fd811236b2743128663e 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -43,6 +43,7 @@ Suggests:
     reshape2,
     MASS,
     qcc,
-    multcomp
+    multcomp,
+    plyr
 VignetteBuilder: knitr
 RoxygenNote: 5.0.1
diff --git a/R/AndradeTb2.1.R b/R/AndradeTb2.1.R
new file mode 100644
index 0000000000000000000000000000000000000000..2cabff3365cb41b849af0ba4b9d2bf5ad5d8c5f0
--- /dev/null
+++ b/R/AndradeTb2.1.R
@@ -0,0 +1,92 @@
+#' @name AndradeTb2.1
+#' @title Experimento de Competição de Híbridos de Milho
+#' @description Um pesquisador instalou um experimento, na região de
+#'     Chapecó, SC, para avaliação do comportamento de híbridos de
+#'     milho.
+#' @format Um \code{data.frame} com 32 observações e 7 variáveis, em que
+#'
+#' \describe{
+#'
+#' \item{\code{hibr}}{Fator de 32 níveis qualitativos que são os
+#'     híbridos de milho.}
+#'
+#' \item{\code{rm}}{Rendimento médio do híbrido de milho, em
+#'     \eqn{kg.ha^{-1}}.}
+#'
+#' \item{\code{ciclo}}{Ciclo de vida do híbrido de milho, em dias.}
+#'
+#' \item{\code{ap}}{Altura da planta do híbrido de milho, em cm.}
+#'
+#' \item{\code{ae}}{Altura da espiga do híbrido de milho, em cm.}
+#'
+#' \item{\code{grao}}{Fator de 3 níveis qualitativos que indica o tipo
+#'     de grão.}
+#' 
+#' \item{\code{resist}}{Fator de 4 níveis qualitativos que indica a
+#'     resistência à ferrugem, sendo: r = resistente; mr = moderadamente
+#'     resistente; ms = moderadamente susceptível; s = susceptível.}
+#' 
+#' }
+#' @keywords AAS
+#' @source Estatística para as ciências agrárias e biológicas: com
+#'     noções de experimentação / Dalton F. Andrade, Paulo J. Ogliari.
+#'     2. ed. rev. e ampl. - Florianópolis: Ed. da UFSC, 2010. (Tabela
+#'     1.1, pág. 37)
+#' @examples
+#'
+#' library(lattice)
+#' library(plyr)
+#'
+#' data(AndradeTb2.1)
+#' str(AndradeTb2.1)
+#'
+#' graf1 <- xyplot(rm ~ hibr,
+#'                 data = AndradeTb2.1,
+#'                 type = c("p", "r"),
+#'                 xlab = "Híbrido de Milho",
+#'                 ylab = "Rendimento Médio (em kg/ha)")
+#'
+#' graf2 <- xyplot(ciclo ~ hibr,
+#'                 data = AndradeTb2.1,
+#'                 type = c("p", "r"),
+#'                 xlab = "Híbrido de Milho",
+#'                 ylab = "Ciclo (em dias)")
+#'
+#' graf3 <- xyplot(ap + ae ~ hibr,
+#'                 data = AndradeTb2.1,
+#'                 type = c("p", "r"),
+#'                 xlab = "Híbrido de Milho",
+#'                 ylab = "Alturas (em cm)",
+#'                 auto.key = list(columns = 2))
+#'
+#' print(graf1, c(0, 0.66, 1, 1), more = TRUE)
+#'
+#' print(graf2, c(0, 0.33, 1, 0.66), more = TRUE)
+#'
+#' print(graf3, c(0, 0, 1, 0.33))
+#'
+#' d <- count(AndradeTb2.1$grao == "dentado")[[2]][2]
+#'
+#' sd <- count(AndradeTb2.1$grao == "semidentado")[[2]][2]
+#'
+#' s <- count(AndradeTb2.1$grao == "semiduro")[[2]][2]
+#'
+#' barplot(c(Dentado = d, Semidentado = sd, Semiduro = s),
+#'         col = "darkturquoise",
+#'         xlab = "Tipo de Grão",
+#'         ylab = "Frequência")                
+#'
+#' r <- count(AndradeTb2.1$resist == "r")[[2]][2]
+#'
+#' mr <- count(AndradeTb2.1$resist == "mr")[[2]][2]
+#'
+#' ms <- count(AndradeTb2.1$resist == "ms")[[2]][2]
+#'
+#' s <- count(AndradeTb2.1$resist == "s")[[2]][2]
+#'
+#' barplot(c(R = r, mr = mr, MS = ms, S = s),
+#'         col = "darkturquoise",
+#'         xlab = "Resistência à Ferrugem",
+#'         ylab = "Frequência")
+#'
+NULL
\ No newline at end of file
diff --git a/data-raw/.~lock.AndradeEg2.1.csv# b/data-raw/.~lock.AndradeEg2.1.csv#
deleted file mode 100644
index 9fb560acd2f55f00a5e1cce3130eb802ad252aed..0000000000000000000000000000000000000000
--- a/data-raw/.~lock.AndradeEg2.1.csv#
+++ /dev/null
@@ -1 +0,0 @@
-,pet-estatistica2,petestatistica2-G31M-ES2C,26.08.2016 16:41,file:///home/pet-estatistica2/.config/libreoffice/4;
\ No newline at end of file
diff --git a/data-raw/AndradeTb2.1.txt b/data-raw/AndradeTb2.1.txt
new file mode 100644
index 0000000000000000000000000000000000000000..2f470468a6cff95a4b224665d22283245445bf0b
--- /dev/null
+++ b/data-raw/AndradeTb2.1.txt
@@ -0,0 +1,33 @@
+hibr	rm	ciclo	ap	ae	grao	resist
+1	6388	65	242	103	dentado	r
+2	6166	65	258	134	semidentado	r
+3	6047	65	240	104	semidentado	s
+4	5889	66	243	108	semidentado	s
+5	5823	69	257	128	dentado	ms
+6	5513	68	241	108	semidentado	s
+7	5202	64	235	108	dentado	r
+8	5172	68	240	103	dentado	s
+9	5166	69	253	123	dentado	ms
+10	4975	70	250	117	semidentado	ms
+11	4778	70	242	114	dentado	mr
+12	4680	66	245	111	semiduro	ms
+13	4660	69	239	110	semiduro	mr
+14	5403	73	264	138	dentado	ms
+15	5117	76	282	149	dentado	mr
+16	5063	72	274	151	dentado	r
+17	4993	71	279	134	semidentado	r
+18	4980	72	274	140	dentado	ms
+19	4770	73	244	140	dentado	r
+20	4685	71	265	139	semiduro	mr
+21	4614	73	248	110	semidentado	r
+22	4552	73	265	128	semidentado	r
+23	3973	74	261	124	semidentado	mr
+24	4550	71	259	129	semiduro	s
+25	5056	64	252	104	semiduro	mr
+26	4500	70	271	109	dentado	ms
+27	4760	68	243	137	semiduro	r
+28	5110	66	252	141	semidentado	ms
+29	4960	70	262	120	dentado	ms
+30	4769	73	260	118	dentado	r
+31	4849	74	250	119	semidentado	s
+32	5230	71	255	138	semiduro	s
diff --git a/data/AndradeTb2.1.rda b/data/AndradeTb2.1.rda
new file mode 100644
index 0000000000000000000000000000000000000000..a32609e694a411da768e08066cc28ff7e1d64e16
Binary files /dev/null and b/data/AndradeTb2.1.rda differ
diff --git a/man/AndradeTb1.1.Rd b/man/AndradeTb1.1.Rd
index c2aae7165e3a88043bcb5f0d2aedfe294aab4a88..c7ab54d7e41f2c8d6e8e8bb703ad908391516341 100644
--- a/man/AndradeTb1.1.Rd
+++ b/man/AndradeTb1.1.Rd
@@ -19,7 +19,10 @@ Foi feito um experimento, na Universidade Federal de
 
 data(AndradeTb1.1)
 
-hist(AndradeTb1.1, prob = TRUE,
+hist(AndradeTb1.1,
+     prob = TRUE,
+     labels = TRUE,
+     col = "darkturquoise",
      xlab = "Duração do Quinto Estádio Ninfal (em dias)",
      ylab = "Densidade",
      main = NULL,
diff --git a/man/AndradeTb2.1.Rd b/man/AndradeTb2.1.Rd
new file mode 100644
index 0000000000000000000000000000000000000000..e9747606f82361c9089934e449500d8a10a7ebd6
--- /dev/null
+++ b/man/AndradeTb2.1.Rd
@@ -0,0 +1,100 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/AndradeTb2.1.R
+\name{AndradeTb2.1}
+\alias{AndradeTb2.1}
+\title{Experimento de Competição de Híbridos de Milho}
+\format{Um \code{data.frame} com 32 observações e 7 variáveis, em que
+
+\describe{
+
+\item{\code{hibr}}{Fator de 32 níveis qualitativos que são os
+    híbridos de milho.}
+
+\item{\code{rm}}{Rendimento médio do híbrido de milho, em
+    \eqn{kg.ha^{-1}}.}
+
+\item{\code{ciclo}}{Ciclo de vida do híbrido de milho, em dias.}
+
+\item{\code{ap}}{Altura da planta do híbrido de milho, em cm.}
+
+\item{\code{ae}}{Altura da espiga do híbrido de milho, em cm.}
+
+\item{\code{grao}}{Fator de 3 níveis qualitativos que indica o tipo
+    de grão.}
+
+\item{\code{resist}}{Fator de 4 níveis qualitativos que indica a
+    resistência à ferrugem, sendo: r = resistente; mr = moderadamente
+    resistente; ms = moderadamente susceptível; s = susceptível.}
+
+}}
+\source{
+Estatística para as ciências agrárias e biológicas: com
+    noções de experimentação / Dalton F. Andrade, Paulo J. Ogliari.
+    2. ed. rev. e ampl. - Florianópolis: Ed. da UFSC, 2010. (Tabela
+    1.1, pág. 37)
+}
+\description{
+Um pesquisador instalou um experimento, na região de
+    Chapecó, SC, para avaliação do comportamento de híbridos de
+    milho.
+}
+\examples{
+
+library(lattice)
+library(plyr)
+
+data(AndradeTb2.1)
+str(AndradeTb2.1)
+
+graf1 <- xyplot(rm ~ hibr,
+                data = AndradeTb2.1,
+                type = c("p", "r"),
+                xlab = "Híbrido de Milho",
+                ylab = "Rendimento Médio (em kg/ha)")
+
+graf2 <- xyplot(ciclo ~ hibr,
+                data = AndradeTb2.1,
+                type = c("p", "r"),
+                xlab = "Híbrido de Milho",
+                ylab = "Ciclo (em dias)")
+
+graf3 <- xyplot(ap + ae ~ hibr,
+                data = AndradeTb2.1,
+                type = c("p", "r"),
+                xlab = "Híbrido de Milho",
+                ylab = "Alturas (em cm)",
+                auto.key = list(columns = 2))
+
+print(graf1, c(0, 0.66, 1, 1), more = TRUE)
+
+print(graf2, c(0, 0.33, 1, 0.66), more = TRUE)
+
+print(graf3, c(0, 0, 1, 0.33))
+
+d <- count(AndradeTb2.1$grao == "dentado")[[2]][2]
+
+sd <- count(AndradeTb2.1$grao == "semidentado")[[2]][2]
+
+s <- count(AndradeTb2.1$grao == "semiduro")[[2]][2]
+
+barplot(c(Dentado = d, Semidentado = sd, Semiduro = s),
+        col = "darkturquoise",
+        xlab = "Tipo de Grão",
+        ylab = "Frequência")                
+
+r <- count(AndradeTb2.1$resist == "r")[[2]][2]
+
+mr <- count(AndradeTb2.1$resist == "mr")[[2]][2]
+
+ms <- count(AndradeTb2.1$resist == "ms")[[2]][2]
+
+s <- count(AndradeTb2.1$resist == "s")[[2]][2]
+
+barplot(c(R = r, mr = mr, MS = ms, S = s),
+        col = "darkturquoise",
+        xlab = "Resistência à Ferrugem",
+        ylab = "Frequência")
+
+}
+\keyword{AAS}
+