agridat: Agricultural Datasets

# tulaikow.wheat.uniformity.R

  Note: The Russian word “саж” refers to a unit of
  measurement. Specifically, it represents the sazhen (саже́нь or
  са́жень), which was used in traditional Russian systems of
  measurement. The sazhen itself is approximately 3 meters (7 feet)
  long. Google Translate sometimes converts "sazhen" into "soot",
  "meter" or "fathom". 


libs(desplot,dplyr,kw,lattice,magrittr,readxl,readr,reshape2,tibble)

setwd("c:/drop/rpack/agridat/data-raw/")
datw <- read_excel("tulaikow.wheat.uniformity.xlsx","winter", col_names=FALSE)
dats <- read_excel("tulaikow.wheat.uniformity.xlsx","summer", col_names=FALSE)

datw %<>% as.matrix %>% `colnames<-`(1:ncol(datw)) %>% melt %>% rename(row=Var1,col=Var2,yield=value) %>% mutate(season="winter")
dats %<>% as.matrix %>% `colnames<-`(1:ncol(dats)) %>% melt %>% rename(row=Var1,col=Var2,yield=value) %>% mutate(season="summer")

dat <- dplyr::bind_rows(datw, dats)
head(dat)


tulaikow.wheat.uniformity = dat
kw::agex(tulaikow.wheat.uniformity)

group_by(dat, season) |>
  summarize(mn=mean(yield)) # Roemer has 532, 700
# 1 summer  532.
# 2 winter  700.

libs(desplot)
desplot(dat, yield~col*row, subset=season=="winter",
        aspect=10/24, flip=TRUE, tick=TRUE,
        main="tulaikow.wheat.uniformity (winter)")
desplot(dat, yield~col*row, subset=season=="summer",
        aspect=16/15, flip=TRUE, tick=TRUE,
        main="tulaikow.wheat.uniformity (summer)")