diff --git a/R/BanzattoQd5.3.1.R b/R/BanzattoQd5.3.1.R
new file mode 100644
index 0000000000000000000000000000000000000000..48677548268581ec0ae389ea35aecbad5803ef9c
--- /dev/null
+++ b/R/BanzattoQd5.3.1.R
@@ -0,0 +1,54 @@
+#' @name BanzattoQd5.3.1
+#' @title Produção de café em função de adubação NPK
+#' @description Resultados de um experimento em arranjo fatorial triplo
+#' 2 \eqn{\times} 2 \eqn{\times} 2, instalado em delineamento de
+#' blocos casualizados, para estudar o efeito da combinação de
+#' nitrogênio, fósforo e potássio na produção de café coco.
+#' @format Um \code{data.frame} com 48 observações e 5 variáveis.
+#' \describe{
+#'
+#' \item{\code{N}}{Variáveis númerica que indica ausência (0) ou
+#' presença (1) de nitrogênio.}
+#'
+#' \item{\code{P}}{Variáveis númerica que indica ausência (0) ou
+#' presença (1) de fósforo.}
+#'
+#' \item{\code{K}}{Variáveis númerica que indica ausência (0) ou
+#' presença (1) de potássio.}
+#'
+#' \item{\code{bloco}}{Fator de 6 níveis nominais, presente para
+#' controle local.}
+#'
+#' \item{\code{prod}}{Variável resposta contínua que é a produção de
+#' café coco (kg ha\eqn{^{-1}}) obtida em cada uma das parcelas para
+#' cada uma das combinações de N, P e K na adubação.}
+#'
+#' }
+#' @keywords DBC FAT3
+#' @source Banzatto, D. A., Kronka, S. D. (2013). Experimentação
+#' Agrícola (4th ed.). Jaboticabal, SP: Funep. (Quadro 5.3.1,
+#' pág. 113)
+#' @examples
+#'
+#' library(lattice)
+#'
+#' data(BanzattoQd5.3.1)
+#'
+#' addmargins(with(BanzattoQd5.3.1,
+#' tapply(X = prod,
+#' INDEX = list(interaction(N, P, K), bloco),
+#' FUN = sum)))
+#'
+#' xyplot(prod ~ N, data = BanzattoQd5.3.1,
+#' groups = interaction(P, K, sep = ":"), type = c("p", "a"),
+#' auto.key = list(title = "P:K", cex.title = 1,
+#' columns = 1, corner = c(0.05, 0.95)),
+#' xlab = "Nitrogênio (0 = ausente, 1 = presente)",
+#' ylab = expression("Produção de café coco"~(kg~ha^{-1})))
+#'
+#' xyplot(prod ~ interaction(N, P, K, sep = ":"),
+#' data = BanzattoQd5.3.1,
+#' xlab = "N:P:K (0 = ausente, 1 = presente)",
+#' ylab = expression("Produção de café coco"~(kg~ha^{-1})))
+#'
+NULL
diff --git a/R/BanzattoQd5.3.7.R b/R/BanzattoQd5.3.7.R
new file mode 100644
index 0000000000000000000000000000000000000000..8980ad72481a7e34f8a89e975609b3efbe6fa015
--- /dev/null
+++ b/R/BanzattoQd5.3.7.R
@@ -0,0 +1,50 @@
+#' @name BanzattoQd5.3.7
+#' @title Eficiência na utilização de P por cultivares de trigo
+#' @description Experimento em delineamento inteiramente casualizado que
+#' avaliou o efeito a efeiciência na utilização de fósforo por três
+#' cultivares de trigo cultivadas em solo com e sem calagem.
+#' @format Um \code{data.frame} com 48 observações e 5 variáveis.
+#' \describe{
+#'
+#' \item{\code{cult}}{Cultivares de trigo usadas no experimento: BR 20
+#' Guató - tolarante ao alumínio, BR 36 Ianomami - sensível ao
+#' alumínio, BR 40 Tuiúca - moderadamente sensível ao alumínio.}
+#'
+#' \item{\code{calag}}{Quantidade aplicada de calcário feita para
+#' correção da acidez do solo (calagem): 0 e 4,4 ton ha\eqn{^{-1}}.}
+#'
+#' \item{\code{P}}{Quantidade aplicada de fósforo: 0 e 87 mg de P por kg
+#' de solo.}
+#'
+#' \item{\code{rept}}{Identifica a repetição de cada cela experimental.}
+#'
+#' \item{\code{efic}}{Variável resposta que a eficiência na utilização
+#' do fósforo, dado pelo quociente da massa seca da parte aérea
+#' produzida pela quantidade de fósforo absorvida pelas plantas.}
+#'
+#' }
+#' @keywords DIC FAT3
+#' @source Banzatto, D. A., Kronka, S. D. (2013). Experimentação
+#' Agrícola (4th ed.). Jaboticabal, SP: Funep. (Quadro 5.3.7,
+#' pág. 120)
+#' @examples
+#'
+#' library(lattice)
+#'
+#' data(BanzattoQd5.3.7)
+#'
+#' addmargins(with(BanzattoQd5.3.7,
+#' tapply(X = efic,
+#' INDEX = list(interaction(cult, calag, P),
+#' rept),
+#' FUN = sum)))
+#'
+#' xyplot(efic ~ factor(P) | cult, data = BanzattoQd5.3.7,
+#' groups = calag, type = c("p", "a"),
+#' layout = c(NA, 1), jitter.x = TRUE,
+#' auto.key = list(title = expression("Calagem"~(ton~ha^{-1})),
+#' cex.title = 1, columns = 2),
+#' xlab = expression("Fósforo"~(mg~kg^{-1}~de~solo)),
+#' ylab = expression("Eficiência no uso do Fósforo"))
+#'
+NULL
diff --git a/R/BanzattoQd5.5.1.R b/R/BanzattoQd5.5.1.R
new file mode 100644
index 0000000000000000000000000000000000000000..6fce855140dde1a010c9a88c54aad7a9cd6d80ee
--- /dev/null
+++ b/R/BanzattoQd5.5.1.R
@@ -0,0 +1,76 @@
+#' @name BanzattoQd5.5.1
+#' @title Adubação Mineral no Milho com Uso de Confundimento
+#' @description Experimento realizado em Jaboticabal - SP por Vilalta
+#' (1972) que, com um arranjo fatorial triplo, estudou a combinação
+#' de nitrogênio, fósforo e potássio, cada um com 3 níveis (27
+#' combinações ao todo), sobre a produção de grãos de milho. Como o
+#' número de tratamentos é 27, não seria viável, em termos de
+#' controle local, ter blocos desse tamanho. Portanto, usou-se a
+#' técnica de confundimento para distribuir os 27 tratamentos em 3
+#' blocos de tamanho 9. A técnica foi usada para confundir dois
+#' graus de liberdade da interação tripla com o efeito dos blocos. O
+#' grupo de confundimento adotado foi o W.
+#' @format Um \code{data.frame} com 27 observações e 4 variáveis, em que
+#' \describe{
+#'
+#' \item{\code{N}}{Níveis codificados de nitrogênio (N) que representam
+#' 0, 30 e 60 kg ha\eqn{^{-1}}. A fonte do N é o sulfato de amônio
+#' (20\% de N).}
+#'
+#' \item{\code{P}}{Níveis codificados de fósforo (P) que representam 0,
+#' 30 e 60 kg ha\eqn{^{-1}}. A fonte do P é o superfosfato simples
+#' (18\% de P\eqn{_{2}}O\eqn{_{5}}).}
+#'
+#' \item{\code{K}}{Níveis codificados de potássio (K) que representam 0,
+#' 30 e 60 kg ha\eqn{^{-1}}. A fonte do K é o cloreto de potássio
+#' (60\% de K\eqn{_{2}}O).}
+#'
+#' \item{\code{bloc}}{Fator de níveis nominais que identifica os
+#' blocos.}
+#'
+#' \item{\code{prod}}{Variável resposta observada que é a produção de
+#' grãos de milho (kg ha\eqn{^{-1}}).}
+#'
+#' }
+#'
+#' Cada parcela foi constituída de 6 linhas de plantio com 10 metros de
+#' comprimento, num total de 60 m\eqn{^{2}} e, como área útil, foram
+#' consideradas as 4 linhas centrais com 8 metros de comprimento
+#' (por desconsiderar 1 m em cada extremidade), perfazendo 32
+#' m\eqn{^{2}}. O híbrido de milho utilizado foi o HMD-7974.
+#'
+#' @keywords DBC FAT3 confundimento
+#' @source Banzatto, D. A., Kronka, S. D. (2013). Experimentação
+#' Agrícola (4th ed.). Jaboticabal, SP: Funep. (Quadro 5.5.1,
+#' pág. 131)
+#'
+#' Vilalta, O. A. (1972). Avaliação da produção de milho (\emph{Zea
+#' mays} L.) em função da adubação NPK em um latossolo
+#' vermelho-escuro -- fase arenosa. Trabalho de Conclusão de Curso
+#' (Graduação em Agronomia). Faculdade de Ciências Agrárias e
+#' Veterinárias, Universidade Estadual Paulista, Jaboticabal, SP.
+#'
+#' @examples
+#'
+#' library(lattice)
+#'
+#' data(BanzattoQd5.5.1)
+#'
+#' str(BanzattoQd5.5.1)
+#'
+#' a <- with(BanzattoQd5.5.1,
+#' tapply(X = prod, INDEX = list(N, P, K), FUN = sum))
+#' addmargins(apply(a, MARGIN = c(1, 2), FUN = sum))
+#' addmargins(apply(a, MARGIN = c(1, 3), FUN = sum))
+#' addmargins(apply(a, MARGIN = c(2, 3), FUN = sum))
+#'
+#' xyplot(prod ~ factor(N) | factor(P), data = BanzattoQd5.5.1,
+#' groups = K, type = c("p", "a"), layout = c(NA, 1),
+#' auto.key = list(title = "Potássio (K)", cex.title = 1,
+#' columns = 3),
+#' strip = strip.custom(strip.names = TRUE,
+#' var.name = "Fósforo (P)"),
+#' xlab = "Nitrogênio (N)",
+#' ylab = expression("Produção de grãos de milho"~(ka~ha^{-1})))
+#'
+NULL
diff --git a/data-raw/BanzattoQd5.3.1.txt b/data-raw/BanzattoQd5.3.1.txt
new file mode 100644
index 0000000000000000000000000000000000000000..c907246702501692579bb08690a2efc37ef194ef
--- /dev/null
+++ b/data-raw/BanzattoQd5.3.1.txt
@@ -0,0 +1,49 @@
+N P K bloco prod
+0 0 0 1 3.029
+0 0 0 2 3.857
+0 0 0 3 2.448
+0 0 0 4 2.448
+0 0 0 5 3.543
+0 0 0 6 4.314
+0 0 1 1 2.438
+0 0 1 2 3.086
+0 0 1 3 3.771
+0 0 1 4 4.657
+0 0 1 5 1.962
+0 0 1 6 3.21
+0 1 0 1 3.448
+0 1 0 2 3.6
+0 1 0 3 3.895
+0 1 0 4 4.267
+0 1 0 5 3.086
+0 1 0 6 3.657
+0 1 1 1 3.533
+0 1 1 2 5.048
+0 1 1 3 3.467
+0 1 1 4 4.095
+0 1 1 5 1.876
+0 1 1 6 2.895
+1 0 0 1 3.362
+1 0 0 2 3.714
+1 0 0 3 3.429
+1 0 0 4 3.19
+1 0 0 5 2.686
+1 0 0 6 4.038
+1 0 1 1 4.905
+1 0 1 2 6.295
+1 0 1 3 4.924
+1 0 1 4 4.952
+1 0 1 5 5.381
+1 0 1 6 5.543
+1 1 0 1 4.171
+1 1 0 2 3.114
+1 1 0 3 4.124
+1 1 0 4 3.981
+1 1 0 5 3.038
+1 1 0 6 3.59
+1 1 1 1 4.476
+1 1 1 2 4.752
+1 1 1 3 4.848
+1 1 1 4 4.676
+1 1 1 5 6.829
+1 1 1 6 3.771
diff --git a/data-raw/BanzattoQd5.3.7.txt b/data-raw/BanzattoQd5.3.7.txt
new file mode 100644
index 0000000000000000000000000000000000000000..270636115964a2c2db2b31d7e11754b53cbab9f0
--- /dev/null
+++ b/data-raw/BanzattoQd5.3.7.txt
@@ -0,0 +1,49 @@
+cult calag P rept efic
+BR 20 Guató 0 0 1 1255
+BR 20 Guató 0 0 2 1250
+BR 20 Guató 0 0 3 908
+BR 20 Guató 0 0 4 1431
+BR 20 Guató 0 87 1 556
+BR 20 Guató 0 87 2 476
+BR 20 Guató 0 87 3 588
+BR 20 Guató 0 87 4 500
+BR 20 Guató 4.4 0 1 714
+BR 20 Guató 4.4 0 2 770
+BR 20 Guató 4.4 0 3 667
+BR 20 Guató 4.4 0 4 667
+BR 20 Guató 4.4 87 1 417
+BR 20 Guató 4.4 87 2 454
+BR 20 Guató 4.4 87 3 454
+BR 20 Guató 4.4 87 4 385
+BR 36 Ianomami 0 0 1 1428
+BR 36 Ianomami 0 0 2 1444
+BR 36 Ianomami 0 0 3 1667
+BR 36 Ianomami 0 0 4 1428
+BR 36 Ianomami 0 87 1 625
+BR 36 Ianomami 0 87 2 526
+BR 36 Ianomami 0 87 3 667
+BR 36 Ianomami 0 87 4 526
+BR 36 Ianomami 4.4 0 1 769
+BR 36 Ianomami 4.4 0 2 911
+BR 36 Ianomami 4.4 0 3 1000
+BR 36 Ianomami 4.4 0 4 1254
+BR 36 Ianomami 4.4 87 1 370
+BR 36 Ianomami 4.4 87 2 476
+BR 36 Ianomami 4.4 87 3 417
+BR 36 Ianomami 4.4 87 4 357
+BR 40 Tuiúca 0 0 1 1660
+BR 40 Tuiúca 0 0 2 1662
+BR 40 Tuiúca 0 0 3 1667
+BR 40 Tuiúca 0 0 4 1667
+BR 40 Tuiúca 0 87 1 526
+BR 40 Tuiúca 0 87 2 714
+BR 40 Tuiúca 0 87 3 588
+BR 40 Tuiúca 0 87 4 714
+BR 40 Tuiúca 4.4 0 1 625
+BR 40 Tuiúca 4.4 0 2 909
+BR 40 Tuiúca 4.4 0 3 909
+BR 40 Tuiúca 4.4 0 4 667
+BR 40 Tuiúca 4.4 87 1 526
+BR 40 Tuiúca 4.4 87 2 556
+BR 40 Tuiúca 4.4 87 3 400
+BR 40 Tuiúca 4.4 87 4 476
diff --git a/data-raw/BanzattoQd5.5.1.txt b/data-raw/BanzattoQd5.5.1.txt
new file mode 100644
index 0000000000000000000000000000000000000000..f74a3eef04eef548ee4b0e1034bcfa36d7a07c36
--- /dev/null
+++ b/data-raw/BanzattoQd5.5.1.txt
@@ -0,0 +1,28 @@
+N P K prod
+0 0 0 5.719
+0 1 2 5.188
+0 2 1 7.062
+1 0 1 4.812
+1 1 0 5.656
+1 2 2 5.969
+2 0 2 5.938
+2 1 1 5.375
+2 2 0 6.031
+0 0 1 6.406
+0 1 0 6.719
+0 2 2 6.875
+1 0 2 5.812
+1 1 1 6.062
+1 2 0 6.594
+2 0 0 6.031
+2 1 2 6.844
+2 2 1 6.719
+0 0 2 4.781
+0 1 1 6.094
+0 2 0 6.875
+1 0 0 3.5
+1 1 2 7.656
+1 2 1 6.062
+2 0 1 6.812
+2 1 0 5.547
+2 2 2 6.781
diff --git a/data/BanzattoQd5.3.1.rda b/data/BanzattoQd5.3.1.rda
new file mode 100644
index 0000000000000000000000000000000000000000..39228cdea882b634879987bce4652cd80125d7f7
Binary files /dev/null and b/data/BanzattoQd5.3.1.rda differ
diff --git a/data/BanzattoQd5.3.7.rda b/data/BanzattoQd5.3.7.rda
new file mode 100644
index 0000000000000000000000000000000000000000..0e128880c8f7f82513ce9ef38caf8c957905d980
Binary files /dev/null and b/data/BanzattoQd5.3.7.rda differ
diff --git a/data/BanzattoQd5.5.1.rda b/data/BanzattoQd5.5.1.rda
new file mode 100644
index 0000000000000000000000000000000000000000..8c336d4adf19b4fa3fdf9c6b393369cd3dba55a7
Binary files /dev/null and b/data/BanzattoQd5.5.1.rda differ
diff --git a/man/BanzattoQd5.3.1.Rd b/man/BanzattoQd5.3.1.Rd
new file mode 100644
index 0000000000000000000000000000000000000000..7f37fc28447cb68d5b7faa8d5d4a019acbef67e5
--- /dev/null
+++ b/man/BanzattoQd5.3.1.Rd
@@ -0,0 +1,61 @@
+% Generated by roxygen2 (4.1.1): do not edit by hand
+% Please edit documentation in R/BanzattoQd5.3.1.R
+\name{BanzattoQd5.3.1}
+\alias{BanzattoQd5.3.1}
+\title{Produção de café em função de adubação NPK}
+\format{Um \code{data.frame} com 48 observações e 5 variáveis.
+\describe{
+
+\item{\code{N}}{Variáveis númerica que indica ausência (0) ou
+ presença (1) de nitrogênio.}
+
+\item{\code{P}}{Variáveis númerica que indica ausência (0) ou
+ presença (1) de fósforo.}
+
+\item{\code{K}}{Variáveis númerica que indica ausência (0) ou
+ presença (1) de potássio.}
+
+\item{\code{bloco}}{Fator de 6 níveis nominais, presente para
+ controle local.}
+
+\item{\code{prod}}{Variável resposta contínua que é a produção de
+ café coco (kg ha\eqn{^{-1}}) obtida em cada uma das parcelas para
+ cada uma das combinações de N, P e K na adubação.}
+
+}}
+\source{
+Banzatto, D. A., Kronka, S. D. (2013). Experimentação
+ Agrícola (4th ed.). Jaboticabal, SP: Funep. (Quadro 5.3.1,
+ pág. 113)
+}
+\description{
+Resultados de um experimento em arranjo fatorial triplo
+ 2 \eqn{\times} 2 \eqn{\times} 2, instalado em delineamento de
+ blocos casualizados, para estudar o efeito da combinação de
+ nitrogênio, fósforo e potássio na produção de café coco.
+}
+\examples{
+library(lattice)
+
+data(BanzattoQd5.3.1)
+
+addmargins(with(BanzattoQd5.3.1,
+ tapply(X = prod,
+ INDEX = list(interaction(N, P, K), bloco),
+ FUN = sum)))
+
+xyplot(prod ~ N, data = BanzattoQd5.3.1,
+ groups = interaction(P, K, sep = ":"), type = c("p", "a"),
+ auto.key = list(title = "P:K", cex.title = 1,
+ columns = 1, corner = c(0.05, 0.95)),
+ xlab = "Nitrogênio (0 = ausente, 1 = presente)",
+ ylab = expression("Produção de café coco"~(kg~ha^{-1})))
+
+xyplot(prod ~ interaction(N, P, K, sep = ":"),
+ data = BanzattoQd5.3.1,
+ xlab = "N:P:K (0 = ausente, 1 = presente)",
+ ylab = expression("Produção de café coco"~(kg~ha^{-1})))
+}
+\keyword{DBC}
+\keyword{FAT3}
+
diff --git a/man/BanzattoQd5.3.7.Rd b/man/BanzattoQd5.3.7.Rd
new file mode 100644
index 0000000000000000000000000000000000000000..5e7f5b366e21928b9ea271f4622f1b68673d0499
--- /dev/null
+++ b/man/BanzattoQd5.3.7.Rd
@@ -0,0 +1,57 @@
+% Generated by roxygen2 (4.1.1): do not edit by hand
+% Please edit documentation in R/BanzattoQd5.3.7.R
+\name{BanzattoQd5.3.7}
+\alias{BanzattoQd5.3.7}
+\title{Eficiência na utilização de P por cultivares de trigo}
+\format{Um \code{data.frame} com 48 observações e 5 variáveis.
+\describe{
+
+\item{\code{cult}}{Cultivares de trigo usadas no experimento: BR 20
+ Guató - tolarante ao alumínio, BR 36 Ianomami - sensível ao
+ alumínio, BR 40 Tuiúca - moderadamente sensível ao alumínio.}
+
+\item{\code{calag}}{Quantidade aplicada de calcário feita para
+ correção da acidez do solo (calagem): 0 e 4,4 ton ha\eqn{^{-1}}.}
+
+\item{\code{P}}{Quantidade aplicada de fósforo: 0 e 87 mg de P por kg
+ de solo.}
+
+\item{\code{rept}}{Identifica a repetição de cada cela experimental.}
+
+\item{\code{efic}}{Variável resposta que a eficiência na utilização
+ do fósforo, dado pelo quociente da massa seca da parte aérea
+ produzida pela quantidade de fósforo absorvida pelas plantas.}
+
+}}
+\source{
+Banzatto, D. A., Kronka, S. D. (2013). Experimentação
+ Agrícola (4th ed.). Jaboticabal, SP: Funep. (Quadro 5.3.7,
+ pág. 120)
+}
+\description{
+Experimento em delineamento inteiramente casualizado que
+ avaliou o efeito a efeiciência na utilização de fósforo por três
+ cultivares de trigo cultivadas em solo com e sem calagem.
+}
+\examples{
+library(lattice)
+
+data(BanzattoQd5.3.7)
+
+addmargins(with(BanzattoQd5.3.7,
+ tapply(X = efic,
+ INDEX = list(interaction(cult, calag, P),
+ rept),
+ FUN = sum)))
+
+xyplot(efic ~ factor(P) | cult, data = BanzattoQd5.3.7,
+ groups = calag, type = c("p", "a"),
+ layout = c(NA, 1), jitter.x = TRUE,
+ auto.key = list(title = expression("Calagem"~(ton~ha^{-1})),
+ cex.title = 1, columns = 2),
+ xlab = expression("Fósforo"~(mg~kg^{-1}~de~solo)),
+ ylab = expression("Eficiência no uso do Fósforo"))
+}
+\keyword{DIC}
+\keyword{FAT3}
+
diff --git a/man/BanzattoQd5.5.1.Rd b/man/BanzattoQd5.5.1.Rd
new file mode 100644
index 0000000000000000000000000000000000000000..754d17697e3e6435a2ed1573f7d1f763fb2211bd
--- /dev/null
+++ b/man/BanzattoQd5.5.1.Rd
@@ -0,0 +1,82 @@
+% Generated by roxygen2 (4.1.1): do not edit by hand
+% Please edit documentation in R/BanzattoQd5.5.1.R
+\name{BanzattoQd5.5.1}
+\alias{BanzattoQd5.5.1}
+\title{Adubação Mineral no Milho com Uso de Confundimento}
+\format{Um \code{data.frame} com 27 observações e 4 variáveis, em que
+ \describe{
+
+\item{\code{N}}{Níveis codificados de nitrogênio (N) que representam
+ 0, 30 e 60 kg ha\eqn{^{-1}}. A fonte do N é o sulfato de amônio
+ (20\% de N).}
+
+\item{\code{P}}{Níveis codificados de fósforo (P) que representam 0,
+ 30 e 60 kg ha\eqn{^{-1}}. A fonte do P é o superfosfato simples
+ (18\% de P\eqn{_{2}}O\eqn{_{5}}).}
+
+\item{\code{K}}{Níveis codificados de potássio (K) que representam 0,
+ 30 e 60 kg ha\eqn{^{-1}}. A fonte do K é o cloreto de potássio
+ (60\% de K\eqn{_{2}}O).}
+
+\item{\code{bloc}}{Fator de níveis nominais que identifica os
+ blocos.}
+
+\item{\code{prod}}{Variável resposta observada que é a produção de
+ grãos de milho (kg ha\eqn{^{-1}}).}
+
+}
+
+Cada parcela foi constituída de 6 linhas de plantio com 10 metros de
+ comprimento, num total de 60 m\eqn{^{2}} e, como área útil, foram
+ consideradas as 4 linhas centrais com 8 metros de comprimento
+ (por desconsiderar 1 m em cada extremidade), perfazendo 32
+ m\eqn{^{2}}. O híbrido de milho utilizado foi o HMD-7974.}
+\source{
+Banzatto, D. A., Kronka, S. D. (2013). Experimentação
+ Agrícola (4th ed.). Jaboticabal, SP: Funep. (Quadro 5.5.1,
+ pág. 131)
+
+Vilalta, O. A. (1972). Avaliação da produção de milho (\emph{Zea
+ mays} L.) em função da adubação NPK em um latossolo
+ vermelho-escuro -- fase arenosa. Trabalho de Conclusão de Curso
+ (Graduação em Agronomia). Faculdade de Ciências Agrárias e
+ Veterinárias, Universidade Estadual Paulista, Jaboticabal, SP.
+}
+\description{
+Experimento realizado em Jaboticabal - SP por Vilalta
+ (1972) que, com um arranjo fatorial triplo, estudou a combinação
+ de nitrogênio, fósforo e potássio, cada um com 3 níveis (27
+ combinações ao todo), sobre a produção de grãos de milho. Como o
+ número de tratamentos é 27, não seria viável, em termos de
+ controle local, ter blocos desse tamanho. Portanto, usou-se a
+ técnica de confundimento para distribuir os 27 tratamentos em 3
+ blocos de tamanho 9. A técnica foi usada para confundir dois
+ graus de liberdade da interação tripla com o efeito dos blocos. O
+ grupo de confundimento adotado foi o W.
+}
+\examples{
+library(lattice)
+
+data(BanzattoQd5.5.1)
+
+str(BanzattoQd5.5.1)
+
+a <- with(BanzattoQd5.5.1,
+ tapply(X = prod, INDEX = list(N, P, K), FUN = sum))
+addmargins(apply(a, MARGIN = c(1, 2), FUN = sum))
+addmargins(apply(a, MARGIN = c(1, 3), FUN = sum))
+addmargins(apply(a, MARGIN = c(2, 3), FUN = sum))
+
+xyplot(prod ~ factor(N) | factor(P), data = BanzattoQd5.5.1,
+ groups = K, type = c("p", "a"), layout = c(NA, 1),
+ auto.key = list(title = "Potássio (K)", cex.title = 1,
+ columns = 3),
+ strip = strip.custom(strip.names = TRUE,
+ var.name = "Fósforo (P)"),
+ xlab = "Nitrogênio (N)",
+ ylab = expression("Produção de grãos de milho"~(ka~ha^{-1})))
+}
+\keyword{DBC}
+\keyword{FAT3}
+\keyword{confundimento}
+