diff --git a/data-raw/cornYield.R b/data-raw/cornYield.R index bca2a49ad0772abb18b6acf311a5246b84fd5464..7bb2f73bec96ac9b10cd3d0a8e645db4936706fd 100644 --- a/data-raw/cornYield.R +++ b/data-raw/cornYield.R @@ -1,15 +1,15 @@ ##---------------------------------------------------------------------- ## Data generation. -kornYield <- expand.grid(block=gl(4, 1), N=c(-1,1), P=c(-1,1), +cornYield <- expand.grid(block=gl(4, 1), N=c(-1,1), P=c(-1,1), K=c(-1,1), KEEP.OUT.ATTRS=FALSE) -kornYield$yield <- c(1.32, 2.12, 1.75, 2.35, 1.80, 2.20, 2.95, 2.96, +cornYield$yield <- c(1.32, 2.12, 1.75, 2.35, 1.80, 2.20, 2.95, 2.96, 1.66, 2.66, 1.73, 2.58, 1.72, 3.85, 2.62, 3.00, 2.58, 3.56, 2.86, 2.75, 2.72, 3.20, 2.25, 2.75, 2.26, 2.08, 1.95, 2.70, 2.95, 3.28, 2.40, 3.35) -str(kornYield) +str(cornYield) -save(kornYield, file="../data/kornYield.RData") +save(cornYield, file="../data/cornYield.RData") ##---------------------------------------------------------------------- ## Examples. @@ -17,21 +17,21 @@ save(kornYield, file="../data/kornYield.RData") library(lattice) library(latticeExtra) -data(kornYield) -str(kornYield) +data(cornYield) +str(cornYield) xyplot(yield~N|P, groups=K, - data=kornYield, type=c("p", "a"), + data=cornYield, type=c("p", "a"), ylab=expression(Yield~(ton~ha^{-1})), xlab="Nutrient level") xyplot(yield~N, groups=interaction(P, K), - data=kornYield, type=c("p", "a"), + data=cornYield, type=c("p", "a"), auto.key=list(columns=2), ylab=expression(Yield~(ton~ha^{-1})), xlab="Nutrient level") -m0 <- lm(yield~block+(N+P+K)^3, data=kornYield) +m0 <- lm(yield~block+(N+P+K)^3, data=cornYield) par(mfrow=c(2,2)); plot(m0); layout(1) anova(m0) @@ -61,6 +61,6 @@ levelplot(mu~N+K, data=pred, aspect=1, color=brewer.pal(n=11, name="Spectral"))) rm(list=ls()) -load("../data/kornYield.RData") +load("../data/cornYield.RData") ls() -str(kornYield) +str(cornYield) diff --git a/data-raw/cornYield2.R b/data-raw/cornYield2.R index 0ea865c2854ab20c599f78f753f4b6a4263c1793..91ea8b902abe37477b93f48e3b9bd355da4d0f9c 100644 --- a/data-raw/cornYield2.R +++ b/data-raw/cornYield2.R @@ -1,17 +1,17 @@ ##---------------------------------------------------------------------- ## Data generation. -## kornYield2 <- read.table("clipboard", header=TRUE, sep="\t") -kornYield2 <- data.frame( +## cornYield2 <- read.table("clipboard", header=TRUE, sep="\t") +cornYield2 <- data.frame( N=c(0L, 0L, 45L, 45L, 45L, 45L, 45L, 45L, 90L), P=c(0L, 45L, 0L, 45L, 45L, 45L, 45L, 90L, 45L), K=c(0L, 30L, 30L, 0L, 30L, 30L, 60L, 30L, 30L), limestone=c(0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L), yield=c(4.192, 4.427, 4.146, 5.029, 5.29, 5.325, 5.275, 5.465, 5.39)) -str(kornYield2) +str(cornYield2) -save(kornYield2, file="../data/kornYield2.RData") +save(cornYield2, file="../data/cornYield2.RData") ##---------------------------------------------------------------------- ## Examples. @@ -19,20 +19,20 @@ save(kornYield2, file="../data/kornYield2.RData") library(lattice) library(latticeExtra) -data(kornYield2) -str(kornYield2) +data(cornYield2) +str(cornYield2) ## Axial triple factorial with 2 controls. -ftable(xtabs(~N+P+K, data=kornYield2)) +ftable(xtabs(~N+P+K, data=cornYield2)) xyplot(yield~N+P+K, groups=as.integer(limestone==1 | (N+P+K)==0), - data=kornYield2, type=c("p", "a"), + data=cornYield2, type=c("p", "a"), auto.key=TRUE, ylab=expression(Yield~(ton~ha^{-1})), xlab="Nutrient content") rm(list=ls()) -load("../data/kornYield2.RData") +load("../data/cornYield2.RData") ls() -str(kornYield2) +str(cornYield2) diff --git a/data-raw/wgPigs.R b/data-raw/wgPigs.R index bb92f44dff4a56375ec2ff385ab6ddc43856086d..387451f96f54408b9e36a14a6a0bfc0d53c9c51d 100644 --- a/data-raw/wgPigs.R +++ b/data-raw/wgPigs.R @@ -1,23 +1,23 @@ ##---------------------------------------------------------------------- ## Data generation. -wgpigs <- read.table("http://www.leg.ufpr.br/~walmes/data/pimentel_racoes.txt", +wgPigs <- read.table("http://www.leg.ufpr.br/~walmes/data/pimentel_racoes.txt", header=TRUE, sep="\t") -names(wgpigs) <- c("ft","wg") -str(wgpigs) +names(wgPigs) <- c("ft","wg") +str(wgPigs) -save(wgpigs, file="../data/wgpigs.RData") +save(wgPigs, file="../data/wgPigs.RData") ##---------------------------------------------------------------------- ## Examples. require(lattice) -xyplot(wg~ft, data=wgpigs, +xyplot(wg~ft, data=wgPigs, ylab="Weight gain (kg)", xlab="Feeding type") rm(list=ls()) -load("../data/wgpigs.RData") +load("../data/wgPigs.RData") ls() -str(wgpigs) +str(wgPigs) diff --git a/data-raw/wgPigs2.R b/data-raw/wgPigs2.R index 7c6fcedab412228ad1b48fd4bcf253dd85847462..1e72f6b275fde66f0ed88e4c09acbb84a12fea4e 100644 --- a/data-raw/wgPigs2.R +++ b/data-raw/wgPigs2.R @@ -1,33 +1,33 @@ ##---------------------------------------------------------------------- ## Data generation. -wgpigs2 <- read.table("http://www.leg.ufpr.br/~walmes/data/pimentel_castracao.txt", +wgPigs2 <- read.table("http://www.leg.ufpr.br/~walmes/data/pimentel_castracao.txt", header=TRUE, sep="\t") -names(wgpigs2) <- c("litter", "size", "age", "wg") -wgpigs2 <- transform(wgpigs2, litter=factor(litter), size=factor(size)) +names(wgPigs2) <- c("litter", "size", "age", "wg") +wgPigs2 <- transform(wgPigs2, litter=factor(litter), size=factor(size)) -aggregate(wg~age, data=wgpigs2, FUN=mean) +aggregate(wg~age, data=wgPigs2, FUN=mean) -wgpigs2$age <- factor(wgpigs2$age, - levels=levels(wgpigs2$age)[c(4,3,1,2)], +wgPigs2$age <- factor(wgPigs2$age, + levels=levels(wgPigs2$age)[c(4,3,1,2)], labels=c("control", "7", "21", "56")) -str(wgpigs2) +str(wgPigs2) -save(wgpigs2, file="../data/wgpigs2.RData") +save(wgPigs2, file="../data/wgPigs2.RData") ##---------------------------------------------------------------------- ## Examples. library(lattice) -data(wgpigs2) -str(wgpigs2) +data(wgPigs2) +str(wgPigs2) -xyplot(wg~age, data=wgpigs2, groups=litter, +xyplot(wg~age, data=wgPigs2, groups=litter, ylab="Weight gain (kg)", xlab="Age at castration (days)") -m0 <- lm(wg~litter+size+age, data=wgpigs2) +m0 <- lm(wg~litter+size+age, data=wgPigs2) par(mfrow=c(2,2)); plot(m0); layout(1) anova(m0) @@ -37,10 +37,10 @@ library(multcomp) summary(glht(m0, linfct=mcp(age="Dunnet")), test=adjusted(type="single-step")) -m1 <- glm(wg~litter+size+age, data=wgpigs2, family=Gamma) -m2 <- glm(wg~litter+size+age, data=wgpigs2, +m1 <- glm(wg~litter+size+age, data=wgPigs2, family=Gamma) +m2 <- glm(wg~litter+size+age, data=wgPigs2, family=Gamma(link="log")) -m3 <- glm(wg~litter+size+age, data=wgpigs2, +m3 <- glm(wg~litter+size+age, data=wgPigs2, family=Gamma(link="identity")) rbind(logLik(m0), @@ -57,7 +57,7 @@ summary(glht(m3, linfct=mcp(age="Dunnet")), test=adjusted(type="single-step")) rm(list=ls()) -load("../data/wgpigs2.RData") +load("../data/wgPigs2.RData") ls() -str(wgpigs2) +str(wgPigs2)