dat <- structure(list(yearID = c(2014L, 2015L, 2016L, 2017L, 2018L, 2019L, 2013L), yield = c(0.328575295095694, 0.0853721681507231,