Homepage: https://kwstat.github.io/desplot

Repository: https://github.com/kwstat/desplot

desplot makes it easy to plot experimental designs of field trials in agriculture.

Key features

Flexible options to customize appearance of graphic.
Stable, well-tested using lattice graphics.
Beta version using ggplot2 graphics. Note ggplot2 is about 4 times slower.

Installation

# Install the released version from CRAN:
install.packages("desplot")
# Install the development version from GitHub:
install.packages("devtools")
devtools::install_github("kwstat/desplot")

Example 1

This data is is from a split-plot experiment with 6 replicates. The replicates are shown by the colored regions and outlined by the thick lines. The text codes and the thin lines define the whole-plots. The nitrogen sub-plot treatments are shown by the text colors.

require(agridat)
require(desplot)
desplot(yates.oats, block ~ col+row,
        col=nitro, text=gen, cex=1, aspect=511/176,
        out1=block, out2=gen,
        out2.gpar=list(col = "gray50", lwd = 1, lty = 1))

The default graphics are based on lattice. It is also possible to create graphics based on ggplot2 by adding 'gg=TRUE' to the function. This functionality is in development and the legend for the ggplot version is not as polished as the lattice version.

require(agridat)
require(desplot)
desplot(yates.oats, block ~ col+row,
        col=nitro, text=gen, cex=1, aspect=511/176,
        out1=block, out2=gen,
        out2.gpar=list(col = "gray50", lwd = 1, lty = 1), gg=TRUE)
        # note, out2.gpar is ignored

Example 2

Another very useful technique is to color the cells according to a continuous response variable (such as plot yield).

require(agridat)
require(desplot)
desplot(yates.oats, yield ~ col*row,
        col=gen, num=nitro, cex=1, out1=block, aspect=511/176)

Example 3

Multiple locations can be shown in a single figure. Data quality for each cell can be indicated with a slash drawn across questionable values and an X drawn across bad values.

require(agridat)
require(desplot)
bmet <- agridat::besag.met
# Simulate a data quality flag
bmet$flag <- sample(c("G","Q","B"), nrow(bmet), replace=TRUE, prob=c(0.8,0.1,0.1))
desplot(bmet, yield ~ col*row|county,
        main="besag.met",
        num=gen, out1=rep, dq=flag,
        aspect=1.0)

Logo note

The hexagon logo for this package shows a simple experimental design layout on top of a tan (soil-color) background.