data-raw/lu.stability.r

In kwstat/agridat: Agricultural Datasets

rm(list=ls(all=TRUE)) 
yield = scan() 
3679    1875    3418 4071 2384   1574 
3175    2016    3829 3650 2071   1502 
3597    2820    3573 3907 1254   1998 
2704    1939    3951 4055 646   2147 
3799    1913    3407 3667 1763   1618 
3246    1522    2793 3910 1672   1829 
3279    1543    2969 3932 3134   1067 
3625    1660    3424 4080 1177   1102 
2497    2007    2792 3526 1361   979 
2750    2015    2130 3035 824   978 
3091    2482    2080 2056 1125   1132 
2267    1481    2496 3062 620   1543 
2125    2038    3209 4007 1612   1266 
2604    2109    2764 3855 2208   1327 
2687    2434    1975 3808 1615   1873 
2452    1800    2350 3302 1948   1489 
2844    1674    4041 3645 1803   2259 
2624    2099    3534 5367 2524   1490 
2424    1476    3846 4061 1237   1429 
3118    1796    4230 3642 1994   1534

env = gl(6,1,length(yield),labels=c('E1','E2','E3','E4','E5','E6') )
gen = gl(5,24, labels=c('G1','G2','G3','G4','G5') )
block = gl(4, 6, length(yield), labels=c('B1','B2','B3','B4') )

dat <- data.frame(yield=yield,gen=gen,env=env,block=block)

lu.stability <- dat