Learn R Programming

agridat (version 1.27)

baena.bean.uniformity: Uniformity trial of beans in Columbia

Description

Uniformity trial of beans in Columbia.

Usage

data("baena.bean.uniformity")

Arguments

Format

A data frame with 1296 observations on the following 3 variables.

row

row ordinate

col

column ordinate

yield

yield per plot, in grams

Details

The experiment was carried out during the period between Sep 21 and Dec 16, 1976. The sowing was carried out in plot Q2 of the experimental field at the International Center for Tropical Agriculture, CIAT-Palmira, Columbia. Latitude of 3 deg 31 min N and a longitude of 76 deg 19 min W.

The document says (translated): "For the drawing of the fertility contour map, a standard deviation of unit plot yields (S.D. = 40 g) was taken as the grouping criterion for plots with homogeneous yields. The seven classes of homogeneous plots range between 80 and 320 g; plots with yields above or below the limits were considered in the first or last class. Figure 1 shows the contour map."

Note: The source document has a heatmap on page 22. I cannot match the heatmap to the data values given on page 20.

Field width: 36 columns x 1m = 36 m

Field length: 36 rows x 1m = 36 m

Data provenance: OCR performed by Claude. Hand-corrected by K.Wright.

References

None

Examples

Run this code
if (FALSE) {
library(agridat)
data(baena.bean.uniformity)
dat <- baena.bean.uniformity

# Source has mean 183.80, variance 1573.34. We get close by omitting the
# one outlier. (Checked the source data which shows the yield is 450 g.)
filter(dat, yield < 400) |> pull(yield) |> mean() # 183.8672
filter(dat, yield < 400) |> pull(yield) |> var() # 1576.738

libs(desplot)
desplot(dat, yield ~ col*row,
        flip=TRUE, tick=TRUE, aspect=1,
        at=c(0,seq(80,320,40), 500),
        main="baena.bean.uniformity")

}

Run the code above in your browser using DataLab