tests/testthat/_problems/test_gge-151.R

In gge: Genotype Plus Genotype-by-Environment Biplots

# Extracted from test_gge.R:151

# prequel ----------------------------------------------------------------------
require(gge)
mat1 <- matrix(c(50, 55, 65, 50, 60, 65, 75,
                 67, 71, 76, 80, 82, 89, 95,
                 90, 93, 95, 102, 97, 106, 117,
                 98, 102, 105, 130, 135, 137, 133,
                 120, 129, 134, 138, 151, 153, 155),
               ncol=5, byrow=FALSE)
colnames(mat1) <- c("E1","E2","E3","E4","E5")
rownames(mat1) <- c("G1","G2","G3","G4","G5","G6","G7")
mat2 <- mat1
mat2[1,1] <- NA
bar <- transform(lattice::barley, env=paste0(site,year))

# test -------------------------------------------------------------------------
require(agridat)
dat <- yan.winterwheat
m1 <- gge(dat, yield ~ gen*env, scale=FALSE)