From 9921f63e219af537a57fb420eed0903bfa1493ec Mon Sep 17 00:00:00 2001
From: Walmes Zeviani <walmes@ufpr.br>
Date: Fri, 2 Oct 2015 21:22:13 -0300
Subject: [PATCH] Prepares data set from Pimentel page 156.
---
data-raw/peanutYield2.R | 43 +++++++++++++++++++++++++++++++++++++++++
1 file changed, 43 insertions(+)
create mode 100644 data-raw/peanutYield2.R
diff --git a/data-raw/peanutYield2.R b/data-raw/peanutYield2.R
new file mode 100644
index 0000000..5e2b09f
--- /dev/null
+++ b/data-raw/peanutYield2.R
@@ -0,0 +1,43 @@
+##----------------------------------------------------------------------
+## Data generation.
+
+peanutYield2 <- expand.grid(
+ variety=c("40-Roxo", "54-Roxo", "49-Cateto", "53-Tatu"),
+ loc=c("Pindorama 49/50", "Ribeirão Preto 49/50",
+ "Campinas 48/49", "Campinas 42/43"),
+ KEEP.OUT.ATTRS=FALSE)
+
+peanutYield2$meanYield <-
+ c(2100, 2160, 1570, 870, 2650, 2740, 1890, 1570, 2100, 1830, 1890,
+ 1370, 2710, 2610, 2590, 1590)
+
+addmargins(with(peanutYield2,
+ tapply(meanYield, list(variety, loc), FUN=sum)))
+
+peanutYield2 <- peanutYield2[with(peanutYield2,
+ order(loc, variety)),]
+
+str(peanutYield2)
+
+## Put MSE as an attibute to the data.frame.
+mse <- c(52900, 84700, 3970, 106900)
+names(mse) <- levels(peanutYield2$loc)
+attr(peanutYield2, which="MSE") <- mse
+str(peanutYield2)
+
+## save(peanutYield2, file="../data/peanutYield2.RData")
+
+##----------------------------------------------------------------------
+## Examples.
+
+require(lattice)
+
+xyplot(meanYield~variety, data=peanutYield2,
+ groups=loc, type="o",
+ ylab=expression(Yield~(t~ha^{-1})),
+ xlab="Variety")
+
+rm(list=ls())
+load("../data/peanutYield2.RData")
+ls()
+str(peanutYield2)
--
GitLab