From 0fee36084fe7e24f2e0e3d92056bd4e413091ce0 Mon Sep 17 00:00:00 2001
From: Walmes Zeviani <walmeszeviani@gmail.com>
Date: Sun, 13 Sep 2015 22:57:22 -0300
Subject: [PATCH] Add the defoliation dataset used in the gamma count article.
---
data-raw/defoliation.R | 51 +++++++++++++++++++++++++++++++++++++++++
data/defoliation.RData | Bin 0 -> 441 bytes
2 files changed, 51 insertions(+)
create mode 100644 data-raw/defoliation.R
create mode 100644 data/defoliation.RData
diff --git a/data-raw/defoliation.R b/data-raw/defoliation.R
new file mode 100644
index 0000000..b79f198
--- /dev/null
+++ b/data-raw/defoliation.R
@@ -0,0 +1,51 @@
+##----------------------------------------------------------------------
+## Data generation.
+
+phenol <- c("vegetative", "flower bud", "blossom", "boll", "boll open")
+defoliation <- expand.grid(rept=1:5,
+ defol=seq(0, 1, length.out=5),
+ phenol=factor(phenol, levels=phenol),
+ KEEP.OUT.ATTRS=FALSE)
+
+defoliation$bolls <- c(10, 9, 8, 8, 10, 11, 9, 10, 10, 10, 8, 8, 10, 8,
+ 9, 9, 7, 7, 8, 9, 8, 6, 6, 5, 6, 7, 8, 8, 9, 10,
+ 9, 12, 7, 10, 9, 8, 9, 9, 10, 8, 11, 10, 7, 8, 8,
+ 7, 7, 7, 7, 8, 10, 9, 8, 12, 8, 7, 5, 5, 7, 5, 6,
+ 5, 7, 4, 7, 8, 5, 7, 6, 4, 5, 5, 4, 4, 5, 8, 10,
+ 7, 8, 10, 9, 6, 6, 8, 6, 9, 7, 11, 8, 9,6, 6, 6,
+ 6, 7, 3, 3, 2, 4, 3, 11, 7, 9, 12 , 11, 9, 13, 8,
+ 10, 10, 9, 7, 7, 9, 9, 8, 8, 10, 8, 10, 9, 8, 10,
+ 8, 10)
+
+defoliation <- defoliation[,c(3,2,1,4)]
+str(defoliation)
+save(defoliation, file="../data/defoliation.RData")
+rm(list=ls())
+load(file="../data/defoliation.RData")
+ls()
+
+##----------------------------------------------------------------------
+## Examples.
+
+library(lattice)
+library(latticeExtra)
+
+x11(width=7, height=2.8)
+xyplot(bolls~defol|phenol, data=defoliation,
+ layout=c(NA, 1), type=c("p", "smooth"),
+ xlab="Artificial defoliation level",
+ ylab="Number of bolls produced",
+ xlim=extendrange(c(0:1), f=0.15), jitter.x=TRUE)
+
+## Sample mean and variance in each treatment cell.
+mv <- aggregate(bolls~phenol+defol, data=defoliation,
+ FUN=function(x) c(mean=mean(x), var=var(x)))
+str(mv)
+
+xlim <- ylim <- extendrange(c(mv$bolls), f=0.05)
+
+## Evidence in favor of the underdispersion.
+xyplot(bolls[,"var"]~bolls[,"mean"], data=mv,
+ aspect="iso", xlim=xlim, ylim=ylim,
+ ylab="Sample variance", xlab="Sample mean")+
+ layer(panel.abline(a=0, b=1, lty=2))
diff --git a/data/defoliation.RData b/data/defoliation.RData
new file mode 100644
index 0000000000000000000000000000000000000000..7591edb21dc916b9c442700ecfd896ea234fd3fa
GIT binary patch
literal 441
zcmb2|=3oE==C|Qz{SG+@#BFz7>gLSE&L=13B=xAVSzm#D{sH0ZOwp>5yOyjiFf#ne
zTFcTeSIhM3TY67T;o~@kiyKaUxjlFO+`qRL%D&4xe*NB!<etK31?P4wUaoh2y7IG4
z%Z_cH*)e5K`^}i6J(h0|t@{}BIQZqZ>W9?=>HLr0IZExGzV?&Fv?SG(Rq1OlMP=@t
zx<h>Cx|vfId{wSSsr-9V7?$m)7PNF_)D`dYMN8H^Kb!P!_0kT-=wz+KQzq5QDsnBK
z>GdUN@zbWz)=bm9INgtDh3@NqJj|DSys@2;-=@yx_DA-Ln*Rb9yXVBUu;2ZEZR-19
zf%B&(Us-*TU8S0#u9kV$qs+U^dstlTqskerHyo?D!(Pd!aIQrD!=4lB40~V7oVa}C
zc1a)O8qIag9aVecQ}QK0EshueV;Fxr!QF2E$0N^QAC`Kp%3y!lu0y1r>qG<lOXiNl
zzYg9%I<-;#)6wVsrO%j??#_I5<U#w|1Q|V+txbAcn@+O}YG)}+KaCfEr7yhk%!z<6
z(=)d&oV@vM*`?X7krq8ddu^BGsAbtM*>beJ-Tt-ryyyCBUfG(@SCL&><Il^$005Vx
B*-!ug
literal 0
HcmV?d00001
--
GitLab