From 70f5e877392e88e67101f22feb5a410f1e8cf327 Mon Sep 17 00:00:00 2001 From: Walmes Zeviani <walmeszeviani@gmail.com> Date: Thu, 17 Sep 2015 22:36:30 -0300 Subject: [PATCH] Add dataset from a LSD studing 5 sugar cane varieties, Pimentel pg 96. --- R/legTools.R | 71 +++++++++++++++++++++++++++++++++---- data-raw/sugarcaneYield2.R | 43 ++++++++++++++++++++++ data/sugarcaneYield2.RData | Bin 0 -> 343 bytes 3 files changed, 107 insertions(+), 7 deletions(-) create mode 100644 data-raw/sugarcaneYield2.R create mode 100644 data/sugarcaneYield2.RData diff --git a/R/legTools.R b/R/legTools.R index 9dde281..734116e 100644 --- a/R/legTools.R +++ b/R/legTools.R @@ -38,7 +38,7 @@ NULL #' #' @examples #' -#' require(lattice) +#' library(lattice) #' data(wgpigs) #' #' xyplot(wg~ft, data=wgpigs, @@ -75,7 +75,7 @@ NULL #' ed.). Piracicaba, São Paulo: FEALQ. (page 76) #' #' @examples -#' require(lattice) +#' library(lattice) #' data(potatoyield) #' #' plot(yield~variety, data=potatoyield, @@ -115,7 +115,7 @@ NULL #' ed.). Piracicaba, São Paulo: FEALQ. (page 91) #' #' @examples -#' require(lattice) +#' library(lattice) #' data(potatoyield) #' #' xyplot(yield~plow|block, data=plowing, type=c("p", "a"), @@ -237,8 +237,8 @@ NULL #' ed.). Piracicaba, São Paulo: FEALQ. (page 93) #' #' @examples -#' -#' require(lattice) +#' +#' library(lattice) #' data(cassavaYield) #' #' plot(yield~variety, data=cassavaYield, @@ -276,8 +276,8 @@ NULL #' ed.). Piracicaba, São Paulo: FEALQ. (page 92) #' #' @examples -#' -#' require(lattice) +#' +#' library(lattice) #' data(sugarcaneYield) #' #' plot(yield~variety, data=sugarcaneYield, @@ -286,3 +286,60 @@ NULL #' xlab="Variety") #' NULL + +#' @name sugarcaneYield2 +#' +#' @title Sugarcane variety competition experiment +#' +#' @description These data are from an experiment done in a latin square +#' design of size 5. Sugarcane yield (kg/plot) was recorded in each +#' experimental unit. +#' +#' \itemize{ +#' \item \code{row} the rows of the latin square that controls in +#' one dimention. A categorical unordered factor with 5 levels. +#' \item \code{col} the columns of the latin square that controls in +#' one dimention perpendicular to the previus. A categorical +#' unordered factor with 5 levels. +#' \item \code{variety} a categorical unordered factor with 5 +#' levels. +#' \item \code{yield} sugarcane yield (kg/plot). +#' } +#' +#' @docType data +#' +#' @keywords datasets +#' +#' @usage data(sugarcaneYield2) +#' +#' @format a \code{data.frame} with 28 records and 3 variables. +#' +#' @source Frederico, P. (2009). Curso de Estatística Experimental (15th +#' ed.). Piracicaba, São Paulo: FEALQ. (page 96) +#' +#' @examples +#' +#' library(lattice) +#' library(latticeExtra) +#' +#' xyplot(yield~variety|col, groups=row, data=sugarcaneYield2, +#' ylab=expression(Yield~(kg~plot^{-1})), +#' xlab="Variety") +#' +#' ## display.brewer.all() +#' +#' levelplot(yield~row+col, +#' data=sugarcaneYield2, aspect="iso", +#' xlab="Row", ylab="Column", +#' main=expression(Yield~(kg~plot^{-1})), +#' col.regions=colorRampPalette( +#' colors=brewer.pal(n=11, name="Spectral")))+ +#' layer(with(sugarcaneYield2, +#' panel.text(x=row, y=col, +#' label=paste(variety, yield)))) +#' +#' aggregate(yield~row, data=sugarcaneYield2, FUN=mean) +#' aggregate(yield~col, data=sugarcaneYield2, FUN=mean) +#' aggregate(yield~variety, data=sugarcaneYield2, FUN=mean) +#' +NULL diff --git a/data-raw/sugarcaneYield2.R b/data-raw/sugarcaneYield2.R new file mode 100644 index 0000000..418ac36 --- /dev/null +++ b/data-raw/sugarcaneYield2.R @@ -0,0 +1,43 @@ +##---------------------------------------------------------------------- +## Data generation. + +sugarcaneYield2 <- read.table("http://www.leg.ufpr.br/~walmes/data/pimentel_canadeacucar2.txt", + header=TRUE, sep="\t") +names(sugarcaneYield2) <- c("row", "col", "variety", "yield") +sugarcaneYield2 <- transform(sugarcaneYield2, row=factor(row), col=factor(col)) +str(sugarcaneYield2) + +sugarcaneYield2 <- sugarcaneYield2[with(sugarcaneYield2, order(row, col)),] + +save(sugarcaneYield2, file="../data/sugarcaneYield2.RData") + +##---------------------------------------------------------------------- +## Examples. + +library(lattice) +library(latticeExtra) + +xyplot(yield~variety|col, groups=row, data=sugarcaneYield2, + ylab=expression(Yield~(kg~plot^{-1})), + xlab="Variety") + +## display.brewer.all() + +levelplot(yield~row+col, + data=sugarcaneYield2, aspect="iso", + xlab="Row", ylab="Column", + main=expression(Yield~(kg~plot^{-1})), + col.regions=colorRampPalette( + colors=brewer.pal(n=11, name="Spectral")))+ + layer(with(sugarcaneYield2, + panel.text(x=row, y=col, + label=paste(variety, yield)))) + +aggregate(yield~row, data=sugarcaneYield2, FUN=mean) +aggregate(yield~col, data=sugarcaneYield2, FUN=mean) +aggregate(yield~variety, data=sugarcaneYield2, FUN=mean) + +rm(list=ls()) +load("../data/sugarcaneYield2.RData") +ls() +str(sugarcaneYield2) diff --git a/data/sugarcaneYield2.RData b/data/sugarcaneYield2.RData new file mode 100644 index 0000000000000000000000000000000000000000..fb8791f39be4527fb6d00f17e94ffbe7302084dc GIT binary patch literal 343 zcmb2|=3oE==CzYIdL3~PXj9)kdop)d(^P#0j-rGshui}Eb>6HwS9J32Nyj}0)aJ#c z=cYbX{=pDt*Yfr1w5Hyv+m_cCSql|AvgqzRv)|Zyv+>DiOYU+u`#HEYEfG*rSt3+X zv1A43?iq1&GNoTOTE3mVLt5zM(;EIaOmDn4mWl;7uW`Kk;QXhxOWktX)ch7s>w0nJ zb?yd^$ec|#k~YR_K6T&dxARi?-+6tTjyU;h1$PBq|M&m*r)x)!%Riem=g()2Q?jpj zi~U&qan@s*$9q%vYd!t3bZhtZ{!Pb|<{xG3HtRlozsD(0<h|VSp2YYA+M?gXO%}{Q zAX;(wm)JMKGPyO*?;F@(C~x4}%Qo+zsRf()1HL`YcVD$iz9>E{?q<1s%SqMK8l7|J zNN3$Pzu95z{_}sZ=5?2}ZtrKtoF=}PZ)V+hx^Q)7-dCNVJF{ohFwNK-P;i8qfdK#| C@vF4} literal 0 HcmV?d00001 -- GitLab