Skip to content
Snippets Groups Projects
Commit a946c16a authored by Walmes Marques Zeviani's avatar Walmes Marques Zeviani
Browse files

Add Rd files.

parent a1408a7f
No related branches found
No related tags found
No related merge requests found
% Generated by roxygen2 (4.1.1): do not edit by hand
% Please edit documentation in R/legTools.R
\docType{data}
\name{bib1}
\alias{bib1}
\title{A balanced incomplete block design of type I}
\format{a \code{data.frame} with 56 records and 4 variables.}
\source{
Pimentel Gomes, F. (2009). Curso de Estatística Experimental
(15th ed.). Piracicaba, São Paulo: FEALQ. (page 190)
}
\usage{
data(bib1)
}
\description{
This data is under a balanced complete block design
named type I. There are 7 replications that are groups of 4
blocks of size 2 in a such a way that each treatment occurs once
in each replication and 7 times at all. There are 8 treatments,
28 blocks at all. Treatment occur in pairs once. These treatments
are in fact cells of a \eqn{2^3} factorial design from combining
nitrogen (P), phosphorus (P) and potassium (K) fertilizers at two
levels each.
\itemize{
\item \code{rept} a categorical unordered factor with 7 levels. Each
\code{rept} has 4 bloks of size 2.
\item \code{N} content of nitrogen in the fertilizer (low/high).
\item \code{P} content of phosphorus in the fertilizer (low/high).
\item \code{K} content of potassium in the fertilizer (low/high).
\item \code{block} a categorical unordered factor with 4 levels in
each \code{rept}, so 28 at all.
\item \code{treat} a categorical unordered factor with 8 levels, the
treatments studied.
\item \code{y} some response variable. The book doesn't gave details.
}
}
\examples{
require(lattice)
data(bib1)
str(bib1)
xyplot(y~treat|rept, groups=block, data=bib1, type="b",
ylab="Y", xlab="Treatment")
xyplot(y~treat, data=bib1, jitter.x=TRUE,
ylab="Y", xlab="Treatment")
xyplot(y~N|P+K, groups=rept, data=bib1, type="b",
ylab="Y", xlab="Nitrogen")
}
\keyword{datasets}
% Generated by roxygen2 (4.1.1): do not edit by hand
% Please edit documentation in R/legTools.R
\docType{data}
\name{bib2}
\alias{bib2}
\title{A balanced incomplete block design of type II}
\format{a \code{data.frame} with 42 records and 4 variables.}
\source{
Pimentel Gomes, F. (2009). Curso de Estatística Experimental
(15th ed.). Piracicaba, São Paulo: FEALQ. (page 192)
}
\usage{
data(bib2)
}
\description{
This data is under a balanced complete block design
named type II. There are 3 replications that are groups of 7
blocks of size 2 in a such a way that each treatment occurs twice
in each replication and 6 times at all. There are 7 treatments
and 21 blocks at all. Treatments occur in pairs once.
\itemize{
\item \code{rept} a categorical unordered factor with 3 levels. Each
\code{rept} has 7 bloks of size 2.
\item \code{block} a categorical unordered factor with 8 levels in
each \code{rept}, so 21 at all.
\item \code{treat} a categorical unordered factor with 7 levels, the
treatments studied.
\item \code{y} some response variable. The book doesn't gave details.
}
}
\examples{
require(lattice)
data(bib2)
str(bib2)
xyplot(y~treat|rept, groups=block, data=bib2, type="b",
ylab="Y", xlab="Treatment")
xyplot(y~treat, data=bib2, jitter.x=TRUE,
ylab="Y", xlab="Treatment")
}
\keyword{datasets}
% Generated by roxygen2 (4.1.1): do not edit by hand
% Please edit documentation in R/legTools.R
\docType{data}
\name{bib3}
\alias{bib3}
\title{A balanced incomplete block design of type III}
\format{a \code{data.frame} with 30 records and 3 variables.}
\source{
Pimentel Gomes, F. (2009). Curso de Estatística Experimental
(15th ed.). Piracicaba, São Paulo: FEALQ. (page 185)
}
\usage{
data(bib3)
}
\description{
This data is under a balanced complete block design
named type III. There are 5 treatments and 10 blocks of size 3
plots. Each treatment is repeated 6 times and they occour
together (in pairs) 3 times.
\itemize{
\item \code{block} a categorical unordered factor with 10 levels.
\item \code{treat} a categorical unordered factor with 5 levels, the
treatments studied.
\item \code{y} some response variable. The book doesn't gave details.
}
}
\examples{
require(lattice)
data(bib3)
str(bib3)
xyplot(y~treat|block, data=bib3,
ylab="Y",
xlab="Treatment")
g <- nlevels(bib3$treat)
a <- seq(0, by=(2*pi)/(g), length.out=g)
y <- sin(a)
x <- cos(a)
plot(y~x, asp=1, xlim=c(-1,1), ylim=c(-1,1))
for (b in levels(bib3$block)){
cbn <- combn(x=as.integer(bib3$treat[bib3$block==b]),
m=2)
segments(
x0=x[cbn[1,]], y0=y[cbn[1,]],
x1=x[cbn[2,]], y1=y[cbn[2,]], col=b)
}
}
\keyword{datasets}
% Generated by roxygen2 (4.1.1): do not edit by hand
% Please edit documentation in R/legTools.R
\docType{data}
\name{bib3asin}
\alias{bib3asin}
\title{A balanced incomplete block design of type III}
\format{a \code{data.frame} with 52 records and 3 variables.}
\source{
Pimentel Gomes, F. (2009). Curso de Estatística Experimental
(15th ed.). Piracicaba, São Paulo: FEALQ. (page 198)
Fraga Jr., C. G.; Costa, A. S. (1950). Análise de um experimento para
combate de vira-cabeça do tomateiro. Bragantia, 10:305-316.
}
\usage{
data(bib3asin)
}
\description{
This data is under a balanced complete block design
named type III. There are 13 treatments and 13 blocks of size 4
plots. Each treatment is repeated 6 times and they occour
together (in pairs) only once.
\itemize{
\item \code{block} a categorical unordered factor with 13 levels.
\item \code{treat} a categorical unordered factor with 13 levels, the
treatments studied.
\item \code{z} correpond a transformation of the original recorded
variable, the observed percentual of healthy plants, \eqn{p}. So,
\eqn{z = \arcsin{\sqrt{p/100}}}, in radians is applied to
stabilize the variance to comply with the assumption of constant
variance.
}
}
\examples{
require(lattice)
data(bib3asin)
str(bib3asin)
xyplot(z~treat|block, data=bib3asin,
ylab="Arc sin of heathy plants fraction",
xlab="Treatment")
## Why not consider a beta distribution for p?
bib3asin$p <- sin(bib3asin$z*pi/180)^2
xyplot(p~treat|block, data=bib3asin,
ylab="Fraction of healthy plants",
xlab="Treatment")
g <- nlevels(bib3asin$treat)
a <- seq(0, by=(2*pi)/(g), length.out=g)
y <- sin(a)
x <- cos(a)
plot(y~x, asp=1, xlim=c(-1,1), ylim=c(-1,1))
for (b in levels(bib3asin$block)){
cbn <- combn(x=as.integer(bib3asin$treat[bib3asin$block==b]),
m=2)
segments(
x0=x[cbn[1,]], y0=y[cbn[1,]],
x1=x[cbn[2,]], y1=y[cbn[2,]], col=b)
}
}
\keyword{datasets}
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment