diff --git a/data-raw/castorbeansYield.R b/data-raw/castorbeansYield.R
index 1acd1c66e1ecdf0a4818ca49b6474c08425e333f..c4069683f0cbe4777012583e1182435739d39f45 100644
--- a/data-raw/castorbeansYield.R
+++ b/data-raw/castorbeansYield.R
@@ -27,13 +27,16 @@ names(mse) <- levels(castorbeansYield$loc)
 attr(castorbeansYield, which="MSE") <- mse
 str(castorbeansYield)
 
-## save(castorbeansYield, file="../data/castorbeansYield.RData")
+save(castorbeansYield, file="../data/castorbeansYield.RData")
 
 ##----------------------------------------------------------------------
 ## Examples.
 
 require(lattice)
 
+data(castorbeansYield)
+str(castorbeansYield)
+
 xyplot(meanYield~variety, data=castorbeansYield,
        groups=loc, type="o",
        ylab=expression(Yield~(t~ha^{-1})),
diff --git a/data-raw/peanutYield.R b/data-raw/peanutYield.R
index 1834b7b15085339e0e53d979a4caa8a692763e09..d8d3fe5442048db19da024bebb9e68dd97186505 100644
--- a/data-raw/peanutYield.R
+++ b/data-raw/peanutYield.R
@@ -1,5 +1,5 @@
 ##----------------------------------------------------------------------
-## Data generation.
+## Data generation. Pimentel page 149.
 
 peanutYield <- expand.grid(
     variety=c("40-Roxo", "54-Roxo", "49-Cateto", "53-Tatu"),
@@ -28,6 +28,9 @@ str(peanutYield)
 
 require(lattice)
 
+data(peanutYield)
+str(peanutYield)
+
 xyplot(meanYield~variety|year, data=peanutYield,
        groups=loc, type="o",
        ylab=expression(Yield~(t~ha^{-1})),
diff --git a/data-raw/peanutYield2.R b/data-raw/peanutYield2.R
index 5e2b09fe9d93cdc344ac7ff4ec0b897cf0b2add7..86f2650313d3bf5863c786d09620c9c9f99d2600 100644
--- a/data-raw/peanutYield2.R
+++ b/data-raw/peanutYield2.R
@@ -1,5 +1,5 @@
 ##----------------------------------------------------------------------
-## Data generation.
+## Data generation. Pimentel page 156.
 
 peanutYield2 <- expand.grid(
     variety=c("40-Roxo", "54-Roxo", "49-Cateto", "53-Tatu"),
@@ -25,13 +25,16 @@ names(mse) <- levels(peanutYield2$loc)
 attr(peanutYield2, which="MSE") <- mse
 str(peanutYield2)
 
-## save(peanutYield2, file="../data/peanutYield2.RData")
+save(peanutYield2, file="../data/peanutYield2.RData")
 
 ##----------------------------------------------------------------------
 ## Examples.
 
 require(lattice)
 
+data(peanutYield2)
+str(peanutYield2)
+
 xyplot(meanYield~variety, data=peanutYield2,
        groups=loc, type="o",
        ylab=expression(Yield~(t~ha^{-1})),
diff --git a/data-raw/potatoYield2.R b/data-raw/potatoYield2.R
index 5103bbf81864c3c3a8b49ca9ce7e14b8a145be11..1d8b166746df8cebdf9d59c12282ebc76e2584b9 100644
--- a/data-raw/potatoYield2.R
+++ b/data-raw/potatoYield2.R
@@ -1,5 +1,5 @@
 ##----------------------------------------------------------------------
-## Data generation.
+## Data generation. Pimentel page 147.
 
 potatoYield2 <- expand.grid(
     variety=c("Kennebec", "B 25-50 E", "B 1-52", "Huinkul",
@@ -26,13 +26,16 @@ names(mse) <- paste0("loc:", 1:length(mse))
 attr(potatoYield2, which="MSE") <- mse
 str(potatoYield2)
 
-## save(potatoYield2, file="../data/potatoYield2.RData")
+save(potatoYield2, file="../data/potatoYield2.RData")
 
 ##----------------------------------------------------------------------
 ## Examples.
 
 require(lattice)
 
+data(potatoYield2)
+str(potatoYield2)
+
 xyplot(sumYield/4~variety, data=potatoYield2,
        groups=loc, type="o",
        ylab=expression(Yield~(t~ha^{-1})),