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