diff --git a/R/legTools.R b/R/legTools.R
index e063ffe004fcc55f393db350b56edd7a0f2e775d..9a00dd0aab275a54f8974e29748dc948f52a8152 100644
--- a/R/legTools.R
+++ b/R/legTools.R
@@ -700,6 +700,10 @@ NULL
#' @source Pimentel Gomes, F. (2009). Curso de Estatística Experimental
#' (15th ed.). Piracicaba, São Paulo: FEALQ. (page 126)
#'
+#' Straus, E. (1951). Experimentos de adubação na zona canavieira de
+#' Pernambuco. In: Anais da III Reunião Brasileira de Ciência do
+#' Solo. 1:336-443.
+#'
#' @examples
#'
#' library(lattice)
@@ -1505,3 +1509,65 @@ NULL
#' }
#'
NULL
+
+#' @name sugarcaneYield5
+#'
+#' @title A set of 38 essays on the effect of P in sugar cane yield
+#'
+#' @description Sugar cane yield evaluated in a set of 38
+#' experiments. In each one was a triple factorial of NPK
+#' dosage. These data correspond to the mean values by level of P
+#' and the mean squares of each complete analysis.
+#'
+#' \itemize{
+#'
+#' \item \code{essay} an integer that represent each single experiment.
+#'
+#' \item \code{sugarcane} indicates which harvest the yield corresponts
+#' to. \code{"plant"} represents the first harvest of the crop and
+#' \code{"ratoon"} represents the second harvest.
+#'
+#' \item \code{P} phosphorus levels in kg/ha of \eqn{\textrm{P}_2
+#' \textrm{O}_5}.
+#'
+#' \item \code{yield} sugar cane yield (ton/ha).
+#'
+#' }
+#'
+#' @docType data
+#'
+#' @keywords datasets
+#'
+#' @usage data(sugarcaneYield5)
+#'
+#' @format a \code{data.frame} with 228 records and 4 variables. An
+#' attribute \code{MSE} is a \code{data.frame} (76 by 3) with mean
+#' square error estimate from the analysis of the complete data in
+#' each experiment.
+#'
+#' @source Pimentel Gomes, F. (2009). Curso de Estatística Experimental
+#' (15th ed.). Piracicaba, São Paulo: FEALQ. (page 259)
+#'
+#' Straus, E. (1951). Experimentos de adubação na zona canavieira de
+#' Pernambuco. In: Anais da III Reunião Brasileira de Ciência do
+#' Solo. 1:336-443.
+#'
+#' @examples
+#'
+#' library(lattice)
+#'
+#' data(sugarcaneYield5)
+#' str(sugarcaneYield5)
+#'
+#' xyplot(yield~P|sugarcane, groups=essay,
+#' data=sugarcaneYield5, type="o",
+#' ylab=expression(Yield~(ton~ha^{-1})),
+#' xlab=expression(P[2]*O[5]~(kg~ha^{-1})))
+#'
+#' xyplot(yield~P|essay, groups=sugarcane,
+#' data=sugarcaneYield5, auto.key=TRUE,
+#' type="o", strip=FALSE,
+#' ylab=expression(Yield~(ton~ha^{-1})),
+#' xlab=expression(P[2]*O[5]~(kg~ha^{-1})))
+#'
+NULL
diff --git a/data-raw/sugarcaneYield5.R b/data-raw/sugarcaneYield5.R
new file mode 100644
index 0000000000000000000000000000000000000000..9a5a0724b6dd14c7bcc011997be64993a9239751
--- /dev/null
+++ b/data-raw/sugarcaneYield5.R
@@ -0,0 +1,83 @@
+##----------------------------------------------------------------------
+## Data generation. Pimentel page 259.
+
+sugarcaneYield5 <- expand.grid(essay=1:38,
+ P=c(0, 60, 120),
+ sugarcane=c("plant", "ratoon"),
+ KEEP.OUT.ATTRS=FALSE)
+sugarcaneYield5$yield <-
+ c(602, 409, 290, 610, 841, 157, 502, 870, 883, 420, 473, 531, 132,
+ 204, 231, 436, 328, 321, 436, 275, 349, 240, 280, 284, 318, 212,
+ 253, 315, 373, 564, 245, 449, 552, 402, 323, 453, 300, 300, 706,
+ 561, 470, 675, 974, 452, 613, 972, 929, 566, 563, 987, 864, 548,
+ 323, 716, 672, 495, 481, 573, 459, 401, 436, 469, 355, 558, 505,
+ 602, 462, 746, 348, 564, 630, 578, 507, 527, 497, 388, 748, 627,
+ 556, 724, 1051, 543, 772, 966, 964, 599, 585, 1085, 1067, 539,
+ 525, 714, 718, 575, 515, 650, 553, 434, 459, 488, 377, 715, 589,
+ 614, 502, 807, 350, 598, 652, 575, 666, 474, 488, 395, 308, 402,
+ 106, 427, 463, 105, 167, 435, 647, 241, 238, 232, 17, 93, 133,
+ 463, 79, 421, 614, 454, 462, 157, 577, 352, 265, 180, 289, 359,
+ 350, 645, 161, 446, 288, 538, 301, 332, 122, 185, 319, 508, 110,
+ 353, 516, 282, 263, 517, 681, 347, 308, 414, 269, 126, 241, 599,
+ 272, 446, 596, 636, 521, 189, 541, 419, 268, 400, 516, 500, 369,
+ 684, 298, 577, 330, 582, 437, 353, 200, 180, 355, 574, 168, 390,
+ 563, 411, 371, 506, 635, 380, 304, 539, 422, 191, 401, 654, 348,
+ 423, 600, 734, 566, 236, 579, 490, 283, 537, 515, 497, 383, 687,
+ 359, 514, 309, 625, 493, 331, 227, 206)/10
+str(sugarcaneYield5)
+
+qmr <- expand.grid(essay=1:38, sugarcane=c("plant", "ratoon"))
+qmr$mse <- c(66.4, 65.0, 187.5, 83.4, 185.7, 64.4, 254.6, 88.6, 85.9,
+ 41.0, 65.9, 35.7, 363.5, 112.2, 120, 67.8, 140.8, 125.4,
+ 49.6, 45.5, 108.8, 34.4, 53.2, 35.1, 23.4, 37.2, 130.6,
+ 65.0, 71.6, 60.6, 93.8, 63.1, 131.4, 27.3, 56.6, 123.2,
+ 230.6, 98.6, 22.6, 85, 63.2, 124.4, 170.5, 49, 80.8, 71.1,
+ 103.6, 51.6, 40.2, 21.5, 124.80, 19.4, 50.8, 167.1, 26,
+ 62.2, 117, 129, 183, 82.9, 99.4, 64.1, 13.9, 51.5, 278.2,
+ 102.5, 61.7, 74.1, 120.5, 135.5, 92.8, 53.7, 57.7, 248.2,
+ 185.2, 14.1)
+attr(sugarcaneYield5, "MSE") <- qmr
+
+save(sugarcaneYield5, file="../data/sugarcaneYield5.RData")
+
+##----------------------------------------------------------------------
+
+aggregate(yield~P+sugarcane, data=sugarcaneYield5, FUN=mean)
+
+m0 <- lm(yield~factor(essay)+factor(P),
+ data=subset(sugarcaneYield5, sugarcane=="plant"))
+anova(m0)
+
+with(subset(attr(sugarcaneYield5, "MSE"), sugarcane=="plant"),
+ mean(mse)/9)
+
+m0 <- lm(yield~factor(essay)+factor(P),
+ data=subset(sugarcaneYield5, sugarcane=="ratoon"))
+anova(m0)
+
+with(subset(attr(sugarcaneYield5, "MSE"), sugarcane=="ratoon"),
+ mean(mse)/9)
+
+##----------------------------------------------------------------------
+## Examples.
+
+library(lattice)
+
+data(sugarcaneYield5)
+str(sugarcaneYield5)
+
+xyplot(yield~P|sugarcane, groups=essay,
+ data=sugarcaneYield5, type="o",
+ ylab=expression(Yield~(ton~ha^{-1})),
+ xlab=expression(P[2]*O[5]~(kg~ha^{-1})))
+
+xyplot(yield~P|essay, groups=sugarcane,
+ data=sugarcaneYield5, auto.key=TRUE,
+ type="o", strip=FALSE,
+ ylab=expression(Yield~(ton~ha^{-1})),
+ xlab=expression(P[2]*O[5]~(kg~ha^{-1})))
+
+rm(list=ls())
+load("../data/sugarcaneYield5.RData")
+ls()
+str(sugarcaneYield5)
diff --git a/data/sugarcaneYield5.RData b/data/sugarcaneYield5.RData
new file mode 100644
index 0000000000000000000000000000000000000000..a4fdfccc38ecd5c2f9169b5b8eb383073ea003d4
Binary files /dev/null and b/data/sugarcaneYield5.RData differ
diff --git a/man/sugarcaneYield4.Rd b/man/sugarcaneYield4.Rd
index d86df2ed284bcaed511dc6b781d9aec4ac990e2b..6ee5559b4c71b69df3f8308e484af39597f2e58c 100644
--- a/man/sugarcaneYield4.Rd
+++ b/man/sugarcaneYield4.Rd
@@ -8,6 +8,10 @@
\source{
Pimentel Gomes, F. (2009). Curso de Estatística Experimental
(15th ed.). Piracicaba, São Paulo: FEALQ. (page 126)
+
+Straus, E. (1951). Experimentos de adubação na zona canavieira de
+ Pernambuco. In: Anais da III Reunião Brasileira de Ciência do
+ Solo. 1:336-443.
}
\usage{
data(sugarcaneYield4)
diff --git a/man/sugarcaneYield5.Rd b/man/sugarcaneYield5.Rd
new file mode 100644
index 0000000000000000000000000000000000000000..5c320283756981cfbe389799a5c541310194bad7
--- /dev/null
+++ b/man/sugarcaneYield5.Rd
@@ -0,0 +1,61 @@
+% Generated by roxygen2 (4.1.1): do not edit by hand
+% Please edit documentation in R/legTools.R
+\docType{data}
+\name{sugarcaneYield5}
+\alias{sugarcaneYield5}
+\title{A set of 38 essays on the effect of P in sugar cane yield}
+\format{a \code{data.frame} with 228 records and 4 variables. An
+ attribute \code{MSE} is a \code{data.frame} (76 by 3) with mean
+ square error estimate from the analysis of the complete data in
+ each experiment.}
+\source{
+Pimentel Gomes, F. (2009). Curso de Estatística Experimental
+ (15th ed.). Piracicaba, São Paulo: FEALQ. (page 259)
+
+Straus, E. (1951). Experimentos de adubação na zona canavieira de
+ Pernambuco. In: Anais da III Reunião Brasileira de Ciência do
+ Solo. 1:336-443.
+}
+\usage{
+data(sugarcaneYield5)
+}
+\description{
+Sugar cane yield evaluated in a set of 38
+ experiments. In each one was a triple factorial of NPK
+ dosage. These data correspond to the mean values by level of P
+ and the mean squares of each complete analysis.
+
+\itemize{
+
+\item \code{essay} an integer that represent each single experiment.
+
+\item \code{sugarcane} indicates which harvest the yield corresponts
+ to. \code{"plant"} represents the first harvest of the crop and
+ \code{"ratoon"} represents the second harvest.
+
+\item \code{P} phosphorus levels in kg/ha of \eqn{\textrm{P}_2
+\textrm{O}_5}.
+
+\item \code{yield} sugar cane yield (ton/ha).
+
+}
+}
+\examples{
+library(lattice)
+
+data(sugarcaneYield5)
+str(sugarcaneYield5)
+
+xyplot(yield~P|sugarcane, groups=essay,
+ data=sugarcaneYield5, type="o",
+ ylab=expression(Yield~(ton~ha^{-1})),
+ xlab=expression(P[2]*O[5]~(kg~ha^{-1})))
+
+xyplot(yield~P|essay, groups=sugarcane,
+ data=sugarcaneYield5, auto.key=TRUE,
+ type="o", strip=FALSE,
+ ylab=expression(Yield~(ton~ha^{-1})),
+ xlab=expression(P[2]*O[5]~(kg~ha^{-1})))
+}
+\keyword{datasets}
+