if (FALSE) {
library(agridat)
data(morgan.multi.uniformity)
dat <- morgan.multi.uniformity
libs(desplot)
libs(tidyr)
datw <- pivot_wider(dat, names_from=crop, values_from=yield)
desplot(datw, wheat ~ col*plot,
aspect=945/112.5,
main="morgan.multi.uniformity: wheat yield")
desplot(datw, cornfodder ~ col*plot,
aspect=945/112.5,
main="morgan.multi.uniformity: corn fodder yield")
plot(datw$wheat, datw$cornfodder,
xlab="Wheat yield (bu/a)",
ylab="Corn fodder yield (bu/a)",
main="morgan.multi.uniformity")
}
Run the code above in your browser using DataLab