diff --git a/R/PimentelPg382.R b/R/PimentelPg382.R
new file mode 100644
index 0000000000000000000000000000000000000000..8d33f28f5ac77d8904e2842d11c22ff50e583b28
--- /dev/null
+++ b/R/PimentelPg382.R
@@ -0,0 +1,36 @@
+#' @name PimentelPg382
+#' @title Métodos de Enxertia no Pegamento de Mudas
+#' @description Experimento com 3 métodos de enxertia em que haviam 200
+#' estacas para cada método e pegaram 180, 150 e 145 desses
+#' enxertos, respectivamente.
+#' @format Um \code{data.frame} com 3 observações e 3 variáveis, em que
+#'
+#' \describe{
+#'
+#' \item{\code{metod}}{Fator de 3 níveis qualitativos que são os métodos
+#' de enxertia das estacas.}
+#'
+#' \item{\code{morta}}{Quantidade de estacas mortas, ou seja, que não
+#' pegaram com a enxertia, de um total de 200 estacas.}
+#'
+#' \item{\code{viva}}{Quantidade de estacas vivas, ou seja, que pegaram
+#' com a enxertia, de um total de 200 estacas. A soma das vivas com
+#' as mortas é 200, portanto.}
+#'
+#' }
+#' @keywords DIC
+#' @source Pimentel-Gomes, F. (2009). Curso de Estatística Experimental
+#' (15th ed.). Piracicaba, SP: FEALQ. (Página 382)
+#' @examples
+#'
+#' library(lattice)
+#'
+#' data(PimentelPg382)
+#' str(PimentelPg382)
+#'
+#' barchart(morta + viva ~ metod, data = PimentelPg382,
+#' stack = TRUE, auto.key = TRUE,
+#' xlab = "Método de enxertia",
+#' ylab = "Quantidade de estacas")
+#'
+NULL
diff --git a/R/PimentelTb18.2.1.R b/R/PimentelTb18.2.1.R
new file mode 100644
index 0000000000000000000000000000000000000000..f5dda230f4a9b87780e710b7c5d8ce279a779c47
--- /dev/null
+++ b/R/PimentelTb18.2.1.R
@@ -0,0 +1,63 @@
+#' @name PimentelTb18.2.1
+#' @title Ensaio Fatorial com Tratamentos Adicionais de Adubação de
+#' Milho
+#' @description Ensaio de adubação NPK de milho, fatorial de
+#' \eqn{3^{3}}, com confundimento de 2 graus de liberdade da
+#' interação tripla. A cada três blocos de 9 parcelas juntaram-se
+#' mais 5 com tratamentos adicionais combinando calcário e
+#' micronutrientes.
+#' @format Um \code{data.frame} com 42 observações e 8 variáveis, em que
+#'
+#' \describe{
+#'
+#' \item{\code{bloc}}{Fator de 3 níveis qualitativos, usado para
+#' controle local. Os blocos tem 14 parcelas, 9 da porção fatorial e
+#' 5 da porção adicional.}
+#'
+#' \item{\code{trat}}{Fator de 30 níveis qualitativos, que são os
+#' tratamentos aplicados em cada parcela, sendo que cada algarismo
+#' possui um significado diferente conforme sua posição: A posição 1
+#' indica os níveis de nitrogênio, a posição 2 indica os níveis de
+#' fósforo e a posição 3 indica os níveis de potássio. A letra C
+#' indica a adição de calcário e M a adição de micronutrientes.}
+#'
+#' \item{\code{N}}{Fator de níveis codificados que representa a dose de
+#' nitrogênio.}
+#'
+#' \item{\code{P}}{Fator de níveis codificados que representa a dose de
+#' fósforo.}
+#'
+#' \item{\code{K}}{Fator de níveis codificados que representa a dose de
+#' potássio.}
+#'
+#' \item{\code{calc}}{Fator de níveis codificados que indica a presença
+#' (1) ou ausência de calcário.}
+#'
+#' \item{\code{micro}}{Fator de níveis codificados que indica a presença
+#' (1) ou ausência de micronutrientes.}
+#'
+#' \item{\code{prod}}{Produção de milho, em kg ha\eqn{^{-1}}.}
+#'
+#' }
+#' @keywords FAT3 adicional
+#' @source Pimentel-Gomes, F. (2009). Curso de Estatística Experimental
+#' (15th ed.). Piracicaba, SP: FEALQ. (Tabela 18.2.1, pág. 330)
+#' @examples
+#'
+#' library(lattice)
+#'
+#' data(PimentelTb18.2.1)
+#' str(PimentelTb18.2.1)
+#'
+#' xtabs(~trat + bloc, data = PimentelTb18.2.1)
+#'
+#' xyplot(prod ~ N | factor(P), groups = K, data = PimentelTb18.2.1,
+#' layout = c(NA, 1), type = c("p", "a"),
+#' xlab = "Nitrogênio (codificado)",
+#' ylab = expression("Produção de milho"~(ka~ha^{-1})),
+#' auto.key = list(title = "Potássio (codificado)",
+#' cex.title = 1.1, columns = 3),
+#' strip = strip.custom(strip.names = TRUE,
+#' var.name = "Fósforo"))
+#'
+NULL
diff --git a/R/PimentelTb20.2.1.R b/R/PimentelTb20.2.1.R
new file mode 100644
index 0000000000000000000000000000000000000000..6e9428adba8d7f2ad9e871645f356739ac3ec9cf
--- /dev/null
+++ b/R/PimentelTb20.2.1.R
@@ -0,0 +1,36 @@
+#' @name PimentelTb20.2.1
+#' @title Ensaio Fatorial de Adubação de Cana-de-açúcar
+#' @description Ensaio fatorial, de \eqn{3^{2}}, de adubação de
+#' cana-de-açúcar com P e K, em 6 blocos casualizados.
+#' @format Um \code{data.frame} com 9 observações e 3 variáveis, em que
+#'
+#' \describe{
+#'
+#' \item{\code{P}}{Variável que indica os níveis de Fósforo (P) em cada
+#' parcela.}
+#'
+#' \item{\code{K}}{Variável que indica os níveis de Potássio (K) em cada
+#' parcela.}
+#'
+#' \item{\code{totais}}{Produção total nos 6 blocos, em ton
+#' ha\eqn{^{-1}}. Valores individuais não disponíveis.}
+#'
+#' }
+#' @keywords FAT2
+#' @source Pimentel-Gomes, F. (2009). Curso de Estatística Experimental
+#' (15th ed.). Piracicaba, SP: FEALQ. (Tabela 20.2.1, pág. 369)
+#' @examples
+#'
+#' library(lattice)
+#'
+#' data(PimentelTb20.2.1)
+#' str(PimentelTb20.2.1)
+#'
+#' xyplot(totais ~ P, data = PimentelTb20.2.1,
+#' groups = K, type = "o",
+#' auto.key = list(title = "Níveis de Potássio (K)",
+#' cex.title = 1.1, columns = 3),
+#' ylab = "Totais de tratamentos (ton/ha)",
+#' xlab = "Níveis de fósforo (P)")
+#'
+NULL
diff --git a/R/PimentelTb21.5.1.R b/R/PimentelTb21.5.1.R
new file mode 100644
index 0000000000000000000000000000000000000000..58eac7df8f995cb1534cd005cbea1a7255797d1d
--- /dev/null
+++ b/R/PimentelTb21.5.1.R
@@ -0,0 +1,34 @@
+#' @name PimentelTb21.5.1
+#' @title Porcentagem de Plantas Doentes
+#' @description Ensaio inteiramente casualizado para avaliar a
+#' porcentagem de plantas doentes em um experimento de tomateiros
+#' onde foi estudado um fator de 3 níveis. Para análise dos
+#' resultados pode-se considerar o de Kruskal- Wallis.
+#' @format Um \code{data.frame} com 12 observações e 2 variáveis, em que
+#'
+#' \describe{
+#'
+#' \item{\code{trat}}{Fator de 3 níveis qualitativos, que são os
+#' tratamentos.}
+#'
+#' \item{\code{doentes}}{Porcentagem de plantas doentes em um ensaio de
+#' tomateiros.}
+#'
+#' }
+#' @keywords DIC
+#' @source Pimentel-Gomes, F. (2009). Curso de Estatística Experimental
+#' (15th ed.). Piracicaba, SP: FEALQ. (Tabela 21.5.1, pág. 384)
+#' @examples
+#'
+#' library(lattice)
+#'
+#' data(PimentelTb21.5.1)
+#' str(PimentelTb21.5.1)
+#'
+#' xyplot(doentes ~ trat,
+#' data = PimentelTb21.5.1,
+#' type = c("a", "p"),
+#' ylab = "Porcentagem de plantas doentes",
+#' xlab = "Tratamento")
+#'
+NULL
diff --git a/data-raw/PimentelPg382.txt b/data-raw/PimentelPg382.txt
new file mode 100644
index 0000000000000000000000000000000000000000..affe6c620e8308c5f2fe91871db94e8e340ba718
--- /dev/null
+++ b/data-raw/PimentelPg382.txt
@@ -0,0 +1,4 @@
+metod morta viva
+A 20 180
+B 50 150
+C 56 144
diff --git a/data-raw/PimentelTb18.2.1.txt b/data-raw/PimentelTb18.2.1.txt
new file mode 100644
index 0000000000000000000000000000000000000000..edabb753249c6059a8a0be59863ae600af542a52
--- /dev/null
+++ b/data-raw/PimentelTb18.2.1.txt
@@ -0,0 +1,43 @@
+bloc trat N P K calc micro prod
+1 000 0 0 0 0 0 981
+1 012 0 1 2 0 0 1216
+1 021 0 2 1 0 0 1183
+1 101 1 0 1 0 0 1213
+1 110 1 1 0 0 0 1028
+1 122 1 2 2 0 0 1574
+1 202 2 0 2 0 0 1111
+1 211 2 1 1 0 0 1333
+1 220 2 2 0 0 0 1775
+2 001 0 0 1 0 0 1213
+2 010 0 1 0 0 0 1510
+2 022 0 2 2 0 0 901
+2 102 1 0 2 0 0 1404
+2 111 1 1 1 0 0 1602
+2 120 1 2 0 0 0 1214
+2 200 2 0 0 0 0 1221
+2 212 2 1 2 0 0 1694
+2 221 2 2 1 0 0 1593
+3 002 0 0 2 0 0 808
+3 011 0 1 1 0 0 942
+3 020 0 2 0 0 0 998
+3 100 1 0 0 0 0 1327
+3 112 1 1 2 0 0 1019
+3 121 1 2 1 0 0 1401
+3 201 2 0 1 0 0 1120
+3 210 2 1 0 0 0 1444
+3 222 2 2 2 0 0 1343
+1 000 0 0 0 0 0 1290
+1 111 1 1 1 0 0 1314
+1 111+C 1 1 1 1 0 1296
+1 111+M 1 1 1 0 1 1237
+1 111+C+M 1 1 1 1 1 1324
+2 000 0 0 0 0 0 1039
+2 111 1 1 1 0 0 680
+2 111+C 1 1 1 1 0 1324
+2 111+M 1 1 1 0 1 1333
+2 111+C+M 1 1 1 1 1 1079
+3 000 0 0 0 0 0 947
+3 111 1 1 1 0 0 1178
+3 111+C 1 1 1 1 0 947
+3 111+M 1 1 1 0 1 1241
+3 111+C+M 1 1 1 1 1 949
diff --git a/data-raw/PimentelTb20.2.1.txt b/data-raw/PimentelTb20.2.1.txt
new file mode 100644
index 0000000000000000000000000000000000000000..adabf06031b8aa48557edfb1ca7d306e4173c99f
--- /dev/null
+++ b/data-raw/PimentelTb20.2.1.txt
@@ -0,0 +1,10 @@
+P K totais
+0 0 230.5
+0 50 239.6
+0 100 187
+60 0 301.5
+60 50 331.6
+60 100 318.4
+120 0 341.2
+120 50 363.3
+120 100 376.9
diff --git a/data-raw/PimentelTb21.5.1.txt b/data-raw/PimentelTb21.5.1.txt
new file mode 100644
index 0000000000000000000000000000000000000000..adf012b05d210747a984940737b18537f0ee3e11
--- /dev/null
+++ b/data-raw/PimentelTb21.5.1.txt
@@ -0,0 +1,13 @@
+trat doentes
+1 10
+1 15
+1 5
+1 30
+2 27
+2 28
+2 44
+2 35
+3 40
+3 70
+3 55
+3 60
diff --git a/data/PimentelPg382.rda b/data/PimentelPg382.rda
new file mode 100644
index 0000000000000000000000000000000000000000..353f8c60efa1dbf121a63f42a1db54fa6cca8451
Binary files /dev/null and b/data/PimentelPg382.rda differ
diff --git a/data/PimentelTb18.2.1.rda b/data/PimentelTb18.2.1.rda
new file mode 100644
index 0000000000000000000000000000000000000000..7936b4cad3ed6d585b83352d48df893fc9b7826e
Binary files /dev/null and b/data/PimentelTb18.2.1.rda differ
diff --git a/data/PimentelTb20.2.1.rda b/data/PimentelTb20.2.1.rda
new file mode 100644
index 0000000000000000000000000000000000000000..51e73a11734f6b4c37672835e10233d57c3b49eb
Binary files /dev/null and b/data/PimentelTb20.2.1.rda differ
diff --git a/data/PimentelTb21.5.1.rda b/data/PimentelTb21.5.1.rda
new file mode 100644
index 0000000000000000000000000000000000000000..920015be0338ef2deeedd2b5d82ae4be1ebd2568
Binary files /dev/null and b/data/PimentelTb21.5.1.rda differ
diff --git a/man/PimentelPg382.Rd b/man/PimentelPg382.Rd
new file mode 100644
index 0000000000000000000000000000000000000000..bc188bf6f3b1ad01098fe950195bd83cf88c872c
--- /dev/null
+++ b/man/PimentelPg382.Rd
@@ -0,0 +1,44 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/PimentelPg382.R
+\name{PimentelPg382}
+\alias{PimentelPg382}
+\title{Métodos de Enxertia no Pegamento de Mudas}
+\format{Um \code{data.frame} com 3 observações e 3 variáveis, em que
+
+\describe{
+
+\item{\code{metod}}{Fator de 3 níveis qualitativos que são os métodos
+ de enxertia das estacas.}
+
+\item{\code{morta}}{Quantidade de estacas mortas, ou seja, que não
+ pegaram com a enxertia, de um total de 200 estacas.}
+
+\item{\code{viva}}{Quantidade de estacas vivas, ou seja, que pegaram
+ com a enxertia, de um total de 200 estacas. A soma das vivas com
+ as mortas é 200, portanto.}
+
+}}
+\source{
+Pimentel-Gomes, F. (2009). Curso de Estatística Experimental
+ (15th ed.). Piracicaba, SP: FEALQ. (Página 382)
+}
+\description{
+Experimento com 3 métodos de enxertia em que haviam 200
+ estacas para cada método e pegaram 180, 150 e 145 desses
+ enxertos, respectivamente.
+}
+\examples{
+
+library(lattice)
+
+data(PimentelPg382)
+str(PimentelPg382)
+
+barchart(morta + viva ~ metod, data = PimentelPg382,
+ stack = TRUE, auto.key = TRUE,
+ xlab = "Método de enxertia",
+ ylab = "Quantidade de estacas")
+
+}
+\keyword{DIC}
+
diff --git a/man/PimentelTb18.2.1.Rd b/man/PimentelTb18.2.1.Rd
new file mode 100644
index 0000000000000000000000000000000000000000..abce79d9794cec20698dd47a23f31d4b995cb5f2
--- /dev/null
+++ b/man/PimentelTb18.2.1.Rd
@@ -0,0 +1,72 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/PimentelTb18.2.1.R
+\name{PimentelTb18.2.1}
+\alias{PimentelTb18.2.1}
+\title{Ensaio Fatorial com Tratamentos Adicionais de Adubação de
+ Milho}
+\format{Um \code{data.frame} com 42 observações e 8 variáveis, em que
+
+\describe{
+
+\item{\code{bloc}}{Fator de 3 níveis qualitativos, usado para
+ controle local. Os blocos tem 14 parcelas, 9 da porção fatorial e
+ 5 da porção adicional.}
+
+\item{\code{trat}}{Fator de 30 níveis qualitativos, que são os
+ tratamentos aplicados em cada parcela, sendo que cada algarismo
+ possui um significado diferente conforme sua posição: A posição 1
+ indica os níveis de nitrogênio, a posição 2 indica os níveis de
+ fósforo e a posição 3 indica os níveis de potássio. A letra C
+ indica a adição de calcário e M a adição de micronutrientes.}
+
+\item{\code{N}}{Fator de níveis codificados que representa a dose de
+ nitrogênio.}
+
+\item{\code{P}}{Fator de níveis codificados que representa a dose de
+ fósforo.}
+
+\item{\code{K}}{Fator de níveis codificados que representa a dose de
+ potássio.}
+
+\item{\code{calc}}{Fator de níveis codificados que indica a presença
+ (1) ou ausência de calcário.}
+
+\item{\code{micro}}{Fator de níveis codificados que indica a presença
+ (1) ou ausência de micronutrientes.}
+
+\item{\code{prod}}{Produção de milho, em kg ha\eqn{^{-1}}.}
+
+}}
+\source{
+Pimentel-Gomes, F. (2009). Curso de Estatística Experimental
+ (15th ed.). Piracicaba, SP: FEALQ. (Tabela 18.2.1, pág. 330)
+}
+\description{
+Ensaio de adubação NPK de milho, fatorial de
+ \eqn{3^{3}}, com confundimento de 2 graus de liberdade da
+ interação tripla. A cada três blocos de 9 parcelas juntaram-se
+ mais 5 com tratamentos adicionais combinando calcário e
+ micronutrientes.
+}
+\examples{
+
+library(lattice)
+
+data(PimentelTb18.2.1)
+str(PimentelTb18.2.1)
+
+xtabs(~trat + bloc, data = PimentelTb18.2.1)
+
+xyplot(prod ~ N | factor(P), groups = K, data = PimentelTb18.2.1,
+ layout = c(NA, 1), type = c("p", "a"),
+ xlab = "Nitrogênio (codificado)",
+ ylab = expression("Produção de milho"~(ka~ha^{-1})),
+ auto.key = list(title = "Potássio (codificado)",
+ cex.title = 1.1, columns = 3),
+ strip = strip.custom(strip.names = TRUE,
+ var.name = "Fósforo"))
+
+}
+\keyword{FAT3}
+\keyword{adicional}
+
diff --git a/man/PimentelTb20.2.1.Rd b/man/PimentelTb20.2.1.Rd
new file mode 100644
index 0000000000000000000000000000000000000000..89ca4882d70730d8838ad4ebd85410c1223a9cf2
--- /dev/null
+++ b/man/PimentelTb20.2.1.Rd
@@ -0,0 +1,44 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/PimentelTb20.2.1.R
+\name{PimentelTb20.2.1}
+\alias{PimentelTb20.2.1}
+\title{Ensaio Fatorial de Adubação de Cana-de-açúcar}
+\format{Um \code{data.frame} com 9 observações e 3 variáveis, em que
+
+\describe{
+
+\item{\code{P}}{Variável que indica os níveis de Fósforo (P) em cada
+ parcela.}
+
+\item{\code{K}}{Variável que indica os níveis de Potássio (K) em cada
+ parcela.}
+
+\item{\code{totais}}{Produção total nos 6 blocos, em ton
+ ha\eqn{^{-1}}. Valores individuais não disponíveis.}
+
+}}
+\source{
+Pimentel-Gomes, F. (2009). Curso de Estatística Experimental
+ (15th ed.). Piracicaba, SP: FEALQ. (Tabela 20.2.1, pág. 369)
+}
+\description{
+Ensaio fatorial, de \eqn{3^{2}}, de adubação de
+ cana-de-açúcar com P e K, em 6 blocos casualizados.
+}
+\examples{
+
+library(lattice)
+
+data(PimentelTb20.2.1)
+str(PimentelTb20.2.1)
+
+xyplot(totais ~ P, data = PimentelTb20.2.1,
+ groups = K, type = "o",
+ auto.key = list(title = "Níveis de Potássio (K)",
+ cex.title = 1.1, columns = 3),
+ ylab = "Totais de tratamentos (ton/ha)",
+ xlab = "Níveis de fósforo (P)")
+
+}
+\keyword{FAT2}
+
diff --git a/man/PimentelTb21.5.1.Rd b/man/PimentelTb21.5.1.Rd
new file mode 100644
index 0000000000000000000000000000000000000000..9c6c2b07b2d963c9c76d8e3f6ad76015aa5d4d8d
--- /dev/null
+++ b/man/PimentelTb21.5.1.Rd
@@ -0,0 +1,42 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/PimentelTb21.5.1.R
+\name{PimentelTb21.5.1}
+\alias{PimentelTb21.5.1}
+\title{Porcentagem de Plantas Doentes}
+\format{Um \code{data.frame} com 12 observações e 2 variáveis, em que
+
+\describe{
+
+\item{\code{trat}}{Fator de 3 níveis qualitativos, que são os
+ tratamentos.}
+
+\item{\code{doentes}}{Porcentagem de plantas doentes em um ensaio de
+ tomateiros.}
+
+}}
+\source{
+Pimentel-Gomes, F. (2009). Curso de Estatística Experimental
+ (15th ed.). Piracicaba, SP: FEALQ. (Tabela 21.5.1, pág. 384)
+}
+\description{
+Ensaio inteiramente casualizado para avaliar a
+ porcentagem de plantas doentes em um experimento de tomateiros
+ onde foi estudado um fator de 3 níveis. Para análise dos
+ resultados pode-se considerar o de Kruskal- Wallis.
+}
+\examples{
+
+library(lattice)
+
+data(PimentelTb21.5.1)
+str(PimentelTb21.5.1)
+
+xyplot(doentes ~ trat,
+ data = PimentelTb21.5.1,
+ type = c("a", "p"),
+ ylab = "Porcentagem de plantas doentes",
+ xlab = "Tratamento")
+
+}
+\keyword{DIC}
+