diff --git a/R/BanzattoQd6.2.2.R b/R/BanzattoQd6.2.2.R
new file mode 100644
index 0000000000000000000000000000000000000000..a61d3b3d933fead34ab3da2c5ae388c404ac3173
--- /dev/null
+++ b/R/BanzattoQd6.2.2.R
@@ -0,0 +1,52 @@
+#' @name BanzattoQd6.2.2
+#' @title Espaçamento de Plantio na Produção de Cana-de-açúcar
+#' @description O experimento avaliou o efeito do espaçamento de plantio
+#' na produção de variedades de cana-de-açúcar. Para estudar os dois
+#' fatores, espaçamento e variedades, considerou-se o delineamento
+#' de blocos casualizados com arranjo de tratamentos de parcelas
+#' subdivididas com o espaçamento sendo casualizado dentro dos
+#' níveis de variedade.
+#' @format Um \code{data.frame} com 40 observações e 4 variáveis.
+#'
+#' \describe{
+#'
+#' \item{\code{varied}}{Fator de níveis nominais que representa as
+#' variedades de cana-de-açúcar.}
+#'
+#' \item{\code{espac}}{Fator de níveis nominais que representa o
+#' espaçamento entre linhas para o plantio da cana-de-açúcar.}
+#'
+#' \item{\code{bloco}}{Fator de níveis nominais que são os blocos do
+#' experimento.}
+#'
+#' \item{\code{prod}}{Produção de cana-de-açúcar, em ton ha\eqn{^{-1}}.}
+#'
+#' }
+#' @keywords DBC PS
+#' @source Banzatto, D. A., Kronka, S. D. (2013). Experimentação
+#' Agrícola (4th ed.). Jaboticabal, SP: Funep. (Quadro 6.2.2,
+#' pág. 138)
+#' @examples
+#'
+#' library(lattice)
+#'
+#' data(BanzattoQd6.2.2)
+#'
+#' str(BanzattoQd6.2.2)
+#'
+#' ftable(with(BanzattoQd6.2.2,
+#' tapply(prod,
+#' list(varied = varied,
+#' espac = espac,
+#' bloco = bloco),
+#' FUN = mean)))
+#'
+#' xyplot(prod ~ varied, groups = espac, data = BanzattoQd6.2.2,
+#' as.table = TRUE, jitter.x = TRUE, type = c("p", "a"),
+#' auto.key = list(title = "Espaçamento de cultivo",
+#' cex.title = 1, columns = 2,
+#' corner = c(0.85, 0.95)),
+#' xlab = "Variedades de cana-de-açúcar",
+#' ylab = expression("Produção de cana-de-açúcar"~(ton~ha^{-1})))
+#'
+NULL
diff --git a/R/BanzattoQd6.2.5.R b/R/BanzattoQd6.2.5.R
new file mode 100644
index 0000000000000000000000000000000000000000..0200fb9dd01462cbe5087aa7fda55e78862293e0
--- /dev/null
+++ b/R/BanzattoQd6.2.5.R
@@ -0,0 +1,62 @@
+#' @name BanzattoQd6.2.5
+#' @title Fontes e Doses para Adubação Nitrogenada no Milho
+#' @description Experimento que estudou a produção de milho como função
+#' de fontes e doses de adubação nitrogenada. O experimento foi
+#' instalado em delineamento de blocos casualizados com arranjo de
+#' tratamento em parcela subdividida, sendo as doses casualizadas
+#' dentro dos níveis de adubo.
+#' @format Um \code{data.frame} com 48 observações e 4 variáveis.
+#'
+#' \describe{
+#'
+#' \item{\code{adubo}}{Fator de níveis nominais que são os adubos fontes
+#' de nitrogênio.}
+#'
+#' \item{\code{dose}}{Fator de níveis métricos que são as doses
+#' aplicadas, em kg de nitrogênio por hectare.}
+#'
+#' \item{\code{bloco}}{Fator de níveis nominais que representa os blocos
+#' do experimento.}
+#'
+#' \item{\code{prod}}{Produção de milho, em kg ha \eqn{^{-1}}}.
+#'
+#' }
+#' @keywords DBC PS
+#' @source Banzatto, D. A., Kronka, S. D. (2013). Experimentação
+#' Agrícola (4th ed.). Jaboticabal, SP: Funep. (Quadro 6.2.5,
+#' pág. 144)
+#' @examples
+#'
+#' library(lattice)
+#'
+#' data(BanzattoQd6.2.5)
+#'
+#' str(BanzattoQd6.2.5)
+#'
+#' ftable(with(BanzattoQd6.2.5,
+#' tapply(prod,
+#' list(adubo = adubo,
+#' dose = dose,
+#' bloco = bloco),
+#' FUN = mean)))
+#'
+#' xyplot(prod ~ dose, groups = adubo, data = BanzattoQd6.2.5,
+#' as.table = TRUE, jitter.x = TRUE, type = c("p", "a"),
+#' auto.key = list(title = "Adubos fontes de nitrogênio",
+#' cex.title = 1, columns = 2),
+#' xlab = expression(
+#' "Quantidade aplicada de nitrogênio"~(kg~ha^{-1})),
+#' ylab = expression("Produção de milho"~(kg~ha^{-1})))
+#'
+#' xyplot(prod ~ dose | adubo, data = BanzattoQd6.2.5,
+#' as.table = TRUE, jitter.x = TRUE,
+#' xlab = expression(
+#' "Quantidade aplicada de nitrogênio"~(kg~ha^{-1})),
+#' ylab = expression("Produção de milho"~(kg~ha^{-1})),
+#' panel = function(x, y, ...){
+#' panel.xyplot(x, y, ...)
+#' m0 <- lm(y ~ poly(x, degree = 2))
+#' panel.curve(predict(m0, newdata = list(x = x)))
+#' })
+#'
+NULL
diff --git a/data-raw/BanzattoQd6.2.2.txt b/data-raw/BanzattoQd6.2.2.txt
new file mode 100644
index 0000000000000000000000000000000000000000..3757fbc5cfc41b94431688a5dfa269c22994e5ac
--- /dev/null
+++ b/data-raw/BanzattoQd6.2.2.txt
@@ -0,0 +1,41 @@
+varied espac bloco prod
+1 1 1 105
+1 2 1 94.3
+2 1 1 101.8
+2 2 1 86.3
+3 1 1 83.4
+3 2 1 90.7
+4 1 1 72.6
+4 2 1 65.7
+5 1 1 57.3
+5 2 1 51.6
+1 1 2 101.4
+1 2 2 91.7
+2 1 2 96.3
+2 2 2 90.7
+3 1 2 60.7
+3 2 2 58.4
+4 1 2 54.2
+4 2 2 56.3
+5 1 2 50.2
+5 2 2 53.4
+1 1 3 97.2
+1 2 3 93.5
+2 1 3 110
+2 2 3 92.4
+3 1 3 71.3
+3 2 3 65.2
+4 1 3 60.5
+4 2 3 51.3
+5 1 3 61.3
+5 2 3 51.9
+1 1 4 89.6
+1 2 4 81.8
+2 1 4 90.5
+2 2 4 85.8
+3 1 4 62.6
+3 2 4 58.9
+4 1 4 57.4
+4 2 4 52.6
+5 1 4 65.2
+5 2 4 48.7
diff --git a/data-raw/BanzattoQd6.2.5.txt b/data-raw/BanzattoQd6.2.5.txt
new file mode 100644
index 0000000000000000000000000000000000000000..d0d4891dd8049ef20d2ea385375f79975ed8286d
--- /dev/null
+++ b/data-raw/BanzattoQd6.2.5.txt
@@ -0,0 +1,49 @@
+adubo dose bloco prod
+Salitre do Chile 10 1 2747
+Salitre do Chile 40 1 2889
+Salitre do Chile 70 1 3578
+Sulfato de amonio 10 1 3164
+Sulfato de amonio 40 1 3831
+Sulfato de amonio 70 1 4107
+Ureia 10 1 1902
+Ureia 40 1 2547
+Ureia 70 1 3364
+Calnitro 10 1 2978
+Calnitro 40 1 3451
+Calnitro 70 1 3742
+Salitre do Chile 10 2 2702
+Salitre do Chile 40 2 2731
+Salitre do Chile 70 2 3387
+Sulfato de amonio 10 2 2658
+Sulfato de amonio 40 2 3049
+Sulfato de amonio 70 2 3031
+Ureia 10 2 1773
+Ureia 40 2 2642
+Ureia 70 2 2427
+Calnitro 10 2 2769
+Calnitro 40 2 2258
+Calnitro 70 2 2498
+Salitre do Chile 10 3 2671
+Salitre do Chile 40 3 2633
+Salitre do Chile 70 3 3858
+Sulfato de amonio 10 3 3600
+Sulfato de amonio 40 3 4182
+Sulfato de amonio 70 3 3791
+Ureia 10 3 3440
+Ureia 40 3 3347
+Ureia 70 3 4053
+Calnitro 10 3 2640
+Calnitro 40 3 2478
+Calnitro 70 3 3458
+Salitre do Chile 10 4 2547
+Salitre do Chile 40 4 2756
+Salitre do Chile 70 4 4284
+Sulfato de amonio 10 4 2760
+Sulfato de amonio 40 4 3102
+Sulfato de amonio 70 4 3547
+Ureia 10 4 2347
+Ureia 40 4 2924
+Ureia 70 4 3111
+Calnitro 10 4 2542
+Calnitro 40 4 2562
+Calnitro 70 4 3076
diff --git a/data/BanzattoQd6.2.2.rda b/data/BanzattoQd6.2.2.rda
new file mode 100644
index 0000000000000000000000000000000000000000..56638a19e1d48ea1d8c03a30391e8805003a0eab
Binary files /dev/null and b/data/BanzattoQd6.2.2.rda differ
diff --git a/data/BanzattoQd6.2.5.rda b/data/BanzattoQd6.2.5.rda
new file mode 100644
index 0000000000000000000000000000000000000000..40b5802ab10fa9e381fa0588139ff28231b0dcdf
Binary files /dev/null and b/data/BanzattoQd6.2.5.rda differ
diff --git a/man/BanzattoQd6.2.2.Rd b/man/BanzattoQd6.2.2.Rd
new file mode 100644
index 0000000000000000000000000000000000000000..5517b5b6d4869e16fcf7e6abb59fdca2bdd9a8b9
--- /dev/null
+++ b/man/BanzattoQd6.2.2.Rd
@@ -0,0 +1,61 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/BanzattoQd6.2.2.R
+\name{BanzattoQd6.2.2}
+\alias{BanzattoQd6.2.2}
+\title{Espaçamento de Plantio na Produção de Cana-de-açúcar}
+\format{Um \code{data.frame} com 40 observações e 4 variáveis.
+
+\describe{
+
+\item{\code{varied}}{Fator de níveis nominais que representa as
+ variedades de cana-de-açúcar.}
+
+\item{\code{espac}}{Fator de níveis nominais que representa o
+ espaçamento entre linhas para o plantio da cana-de-açúcar.}
+
+\item{\code{bloco}}{Fator de níveis nominais que são os blocos do
+ experimento.}
+
+\item{\code{prod}}{Produção de cana-de-açúcar, em ton ha\eqn{^{-1}}.}
+
+}}
+\source{
+Banzatto, D. A., Kronka, S. D. (2013). Experimentação
+ Agrícola (4th ed.). Jaboticabal, SP: Funep. (Quadro 6.2.2,
+ pág. 138)
+}
+\description{
+O experimento avaliou o efeito do espaçamento de plantio
+ na produção de variedades de cana-de-açúcar. Para estudar os dois
+ fatores, espaçamento e variedades, considerou-se o delineamento
+ de blocos casualizados com arranjo de tratamentos de parcelas
+ subdivididas com o espaçamento sendo casualizado dentro dos
+ níveis de variedade.
+}
+\examples{
+
+library(lattice)
+
+data(BanzattoQd6.2.2)
+
+str(BanzattoQd6.2.2)
+
+ftable(with(BanzattoQd6.2.2,
+ tapply(prod,
+ list(varied = varied,
+ espac = espac,
+ bloco = bloco),
+ FUN = mean)))
+
+xyplot(prod ~ varied, groups = espac, data = BanzattoQd6.2.2,
+ as.table = TRUE, jitter.x = TRUE, type = c("p", "a"),
+ auto.key = list(title = "Espaçamento de cultivo",
+ cex.title = 1, columns = 2,
+ corner = c(0.85, 0.95)),
+ xlab = "Variedades de cana-de-açúcar",
+ ylab = expression("Produção de cana-de-açúcar"~(ton~ha^{-1})))
+
+}
+\keyword{DBC}
+\keyword{PS}
+
diff --git a/man/BanzattoQd6.2.5.Rd b/man/BanzattoQd6.2.5.Rd
new file mode 100644
index 0000000000000000000000000000000000000000..948a6a1ae31b61c448526ec6a8651c8eb5683313
--- /dev/null
+++ b/man/BanzattoQd6.2.5.Rd
@@ -0,0 +1,71 @@
+% Generated by roxygen2: do not edit by hand
+% Please edit documentation in R/BanzattoQd6.2.5.R
+\name{BanzattoQd6.2.5}
+\alias{BanzattoQd6.2.5}
+\title{Fontes e Doses para Adubação Nitrogenada no Milho}
+\format{Um \code{data.frame} com 48 observações e 4 variáveis.
+
+\describe{
+
+\item{\code{adubo}}{Fator de níveis nominais que são os adubos fontes
+ de nitrogênio.}
+
+\item{\code{dose}}{Fator de níveis métricos que são as doses
+ aplicadas, em kg de nitrogênio por hectare.}
+
+\item{\code{bloco}}{Fator de níveis nominais que representa os blocos
+ do experimento.}
+
+\item{\code{prod}}{Produção de milho, em kg ha \eqn{^{-1}}}.
+
+}}
+\source{
+Banzatto, D. A., Kronka, S. D. (2013). Experimentação
+ Agrícola (4th ed.). Jaboticabal, SP: Funep. (Quadro 6.2.5,
+ pág. 144)
+}
+\description{
+Experimento que estudou a produção de milho como função
+ de fontes e doses de adubação nitrogenada. O experimento foi
+ instalado em delineamento de blocos casualizados com arranjo de
+ tratamento em parcela subdividida, sendo as doses casualizadas
+ dentro dos níveis de adubo.
+}
+\examples{
+
+library(lattice)
+
+data(BanzattoQd6.2.5)
+
+str(BanzattoQd6.2.5)
+
+ftable(with(BanzattoQd6.2.5,
+ tapply(prod,
+ list(adubo = adubo,
+ dose = dose,
+ bloco = bloco),
+ FUN = mean)))
+
+xyplot(prod ~ dose, groups = adubo, data = BanzattoQd6.2.5,
+ as.table = TRUE, jitter.x = TRUE, type = c("p", "a"),
+ auto.key = list(title = "Adubos fontes de nitrogênio",
+ cex.title = 1, columns = 2),
+ xlab = expression(
+ "Quantidade aplicada de nitrogênio"~(kg~ha^{-1})),
+ ylab = expression("Produção de milho"~(kg~ha^{-1})))
+
+xyplot(prod ~ dose | adubo, data = BanzattoQd6.2.5,
+ as.table = TRUE, jitter.x = TRUE,
+ xlab = expression(
+ "Quantidade aplicada de nitrogênio"~(kg~ha^{-1})),
+ ylab = expression("Produção de milho"~(kg~ha^{-1})),
+ panel = function(x, y, ...){
+ panel.xyplot(x, y, ...)
+ m0 <- lm(y ~ poly(x, degree = 2))
+ panel.curve(predict(m0, newdata = list(x = x)))
+ })
+
+}
+\keyword{DBC}
+\keyword{PS}
+