bose.multi.uniformity {agridat} | R Documentation |

##
Uniformity trials of barley, wheat, lentils

### Description

Uniformity trials of barley, wheat, lentils in India 1930-1932.

### Usage

data("bose.multi.uniformity")

### Format

A data frame with 1170 observations on the following 5 variables.

`year`

year

`crop`

crop

`row`

row ordinate

`col`

column ordinate

`yield`

plot yield, grams

### Details

A field about 1/4 acre was sown in three consecutive years (beginning
in 1929-1930) with barley, wheat, and lentil.

At harvest, borders 3 feet on east and west and 6 feet on north and
south were removed. The field was divided into plots four feet square,
which were harvested separately, measured in grams.

Fertility contours of the field were somewhat similar across years,
with correlation values across years 0.45, 0.48, 0.21.

Field width: 15 plots * 4 feet = 60 feet.

Field length: 26 plots * 4 feet = 104 feet.

Conclusions:

"An experimental field which may be sensibly uniform for one crop or
for one season may not be so for another crop or in a different
season" p. 592.

### Source

Bose, R. D. (1935).
Some soil heterogeneity trials at Pusa and the size and shape of experimental plots.
Ind. J. Agric. Sci., 5, 579-608. Table 1 (p. 585), Table 4 (p. 589), Table 5 (p. 590).
https://archive.org/details/in.ernet.dli.2015.271739

### References

Shaw (1935).
Handbook of Statistics for Use in Plant-Breeding and Agricultural Problems, p. 149-170.
https://krishikosh.egranth.ac.in/handle/1/21153

### Examples

## Not run:
library(agridat)
data(bose.multi.uniformity)
dat <- bose.multi.uniformity
# match sum at bottom of Bose tables 1, 4, 5
# library(dplyr)
# dat
libs(desplot, dplyr)
# Calculate percent of mean yield for each year
dat <- group_by(dat, year)
dat <- mutate(dat, pctyld = (yield-mean(yield))/mean(yield))
dat <- ungroup(dat)
dat <- mutate(dat, year=as.character(year))
# Bose smoothed the data by averaging 2x3 plots together before drawing
# contour maps. Heatmaps of raw data have similar structure to Bose Fig 1.
desplot(dat, pctyld ~ col*row|year,
tick=TRUE, flip=TRUE, aspect=(26)/(15),
main="bose.multi.* - Percent of mean yield")
# contourplot() results need to be mentally flipped upside down
# contourplot(pctyld ~ col*row|year, dat,
# region=TRUE, as.table=TRUE, aspect=26/15)
## End(Not run)

[Package

*agridat* version 1.18

Index]