From ef1e5d0dd28176d0a50249de45da8ee7a1011ae1 Mon Sep 17 00:00:00 2001 From: Walmes Zeviani <walmeszeviani@gmail.com> Date: Thu, 17 Sep 2015 20:44:24 -0300 Subject: [PATCH] Add dataset of a CRBD studing cassava varieties, Pimentel pg 93. --- R/legTools.R | 41 ++++++++++++++++++++++++++++++++++++++++ data-raw/cassavaYield.R | 38 +++++++++++++++++++++++++++++++++++++ data/cassavaYield.RData | Bin 0 -> 411 bytes man/cassavaYield.Rd | 40 +++++++++++++++++++++++++++++++++++++++ 4 files changed, 119 insertions(+) create mode 100644 data-raw/cassavaYield.R create mode 100644 data/cassavaYield.RData create mode 100644 man/cassavaYield.Rd diff --git a/R/legTools.R b/R/legTools.R index 9cbf85e..1302641 100644 --- a/R/legTools.R +++ b/R/legTools.R @@ -206,3 +206,44 @@ NULL #' layer(panel.abline(a=0, b=1, lty=2)) #' NULL + +#' @name cassavaYield +#' +#' @title Cassava variety competition experiment +#' +#' @description These data are from an experiment done by The Brazilian +#' Agricultural Research Corporation in Cassava & Tropical Fruits +#' (Centro Nacional de Pesquisa da Mandioca e Fruticultura da +#' Embrapa) at Cruz das Almas, Bahia. The experiment was done in a +#' randomized complete block design with 4 blocks. Cassava yield +#' (t/ha) was recorded in each experimental unit. +#' +#' \itemize{ +#' \item \code{block} a categorical unordered factor with 4 levels. +#' \item \code{variety} a categorical unordered factor with 6 +#' levels. +#' \item \code{yield} cassava yield (t/ha). +#' } +#' +#' @docType data +#' +#' @keywords datasets +#' +#' @usage data(cassavaYield) +#' +#' @format a \code{data.frame} with 24 records and 3 variables. +#' +#' @source Frederico, P. (2009). Curso de Estatística Experimental (15th +#' ed.). Piracicaba, São Paulo: FEALQ. (page 93) +#' +#' @examples +#' +#' require(lattice) +#' data(cassavaYield) +#' +#' plot(yield~variety, data=cassavaYield, +#' groups=block, type="o", +#' ylab=expression(Yield~(t~ha^{-1})), +#' xlab="Variety") +#' +NULL diff --git a/data-raw/cassavaYield.R b/data-raw/cassavaYield.R new file mode 100644 index 0000000..2491776 --- /dev/null +++ b/data-raw/cassavaYield.R @@ -0,0 +1,38 @@ +##---------------------------------------------------------------------- +## Data generation. + +cassavaYield <- read.table("http://www.leg.ufpr.br/~walmes/data/pimentel_mandioca.txt", + header=TRUE, sep="\t") +names(cassavaYield) <- c("block", "variety", "yield") +str(cassavaYield) + +cassavaYield <- cassavaYield[with(cassavaYield, order(block, variety)),] + +save(cassavaYield, file="../data/cassavaYield.RData") + +##---------------------------------------------------------------------- +## Examples. + +require(lattice) + +xyplot(yield~variety, data=cassavaYield, + groups=block, type="o", + ylab=expression(Yield~(t~ha^{-1})), + xlab="Variety") + +rm(list=ls()) +load("../data/cassavaYield.RData") +ls() +str(cassavaYield) + +m0 <- lm(yield~block+variety, data=cassavaYield) + +## Accessing the model assumptions. +par(mfrow=c(2,2)); plot(m0); layout(1) + +## ANOVA table. +anova(m0) + +library(doBy) +ls <- LSmatrix(m0, effect="variety") +dput(ls) diff --git a/data/cassavaYield.RData b/data/cassavaYield.RData new file mode 100644 index 0000000000000000000000000000000000000000..07cdb35446a418b6d57821875406864be2fff9a8 GIT binary patch literal 411 zcmb2|=3oE==Cwg*{SG+@#BKZI;Tj@fa7;qz=>bNOJr!#@<ve?0K3HsN&C$8CWoz5) z^$Wx054t~CwZ7%tw#^(Yd48Wt)A!GuJM-Q4@8)~A-McKAJSTIG@w$bxpZxILpmO@r zEaQXU6<gR7t;Ld~1TB3ol@-|be7e@0->{AG<kyH<s&8*BeEZ|YTi1&-d^w7C1k6)A zA}_fv#B(O!KOxbX)8!f@cgR(3o;v?AH}kTTXa7GI>HKq^yd-nd+w6}jZddoQFZ2%i zcFp|mMpc~&(lc*2Ui^|Q_l|?TZQkqE1^c-9dfG3ae>FFq-63=4yJCf7&K<wAR=vL^ zolsfH?AGj>7Q#MJve!KGz1vUkTiXK__LWG^Ec7d!8n3?3>5t(@^^$*Guda9I3&$(} z>-m;0k}r6_s`ongE&0GbZnb`o_pLtEpO~Wk(A4#t?Dd+EaFHhUHwlZ9CjYpoA9}Lt zop*3(RE<Wa(6q@Xg5K~bsk^%7D3%0y%+pr*_Hxxzm)#voS63_W38cDYHYe<2Wi(-z T>vLLkwvO?LkpIf_%nS?w-Z;ae literal 0 HcmV?d00001 diff --git a/man/cassavaYield.Rd b/man/cassavaYield.Rd new file mode 100644 index 0000000..186defb --- /dev/null +++ b/man/cassavaYield.Rd @@ -0,0 +1,40 @@ +% Generated by roxygen2 (4.1.1): do not edit by hand +% Please edit documentation in R/legTools.R +\docType{data} +\name{cassavaYield} +\alias{cassavaYield} +\title{Cassava variety competition experiment} +\format{a \code{data.frame} with 24 records and 3 variables.} +\source{ +Frederico, P. (2009). Curso de Estatística Experimental (15th + ed.). Piracicaba, São Paulo: FEALQ. (page 93) +} +\usage{ +data(cassavaYield) +} +\description{ +These data are from an experiment done by The Brazilian + Agricultural Research Corporation in Cassava & Tropical Fruits + (Centro Nacional de Pesquisa da Mandioca e Fruticultura da + Embrapa) at Cruz das Almas, Bahia. The experiment was done in a + randomized complete block design with 4 blocks. Cassava yield + (t/ha) was recorded in each experimental unit. + +\itemize{ + \item \code{block} a categorical unordered factor with 4 levels. + \item \code{variety} a categorical unordered factor with 6 + levels. + \item \code{yield} cassava yield (t/ha). +} +} +\examples{ +require(lattice) +data(cassavaYield) + +plot(yield~variety, data=cassavaYield, + groups=block, type="o", + ylab=expression(Yield~(t~ha^{-1})), + xlab="Variety") +} +\keyword{datasets} + -- GitLab