# \donttest{
library(metan)
# Genotype data (5 genotypes and 3 replicates)
gen_data <-
g_simula(ngen = 5,
nrep = 3,
seed = 1)
gen_data
inspect(gen_data, plot = TRUE)
aov(V1 ~ GEN + REP, data = gen_data) %>% anova()
# Genotype-environment data
# 5 genotypes, 3 environments, 4 replicates and 2 traits
df <-
ge_simula(ngen = 5,
nenv = 3,
nrep = 4,
nvars = 2,
seed = 1)
ge_plot(df, ENV, GEN, V1)
aov(V1 ~ ENV*GEN + ENV/REP, data = df) %>% anova()
# Change genotype effect (trait 1 with fewer differences among genotypes)
# Define different intercepts for the two traits
df2 <-
ge_simula(ngen = 10,
nenv = 3,
nrep = 4,
nvars = 2,
gen_eff = c(1, 50),
intercept = c(80, 1500),
seed = 1)
ge_plot(df2, ENV, GEN, V2)
# }
Run the code above in your browser using DataLab