## Load the Switzerland data
data(Switzerland)
## Fit the FG DImodel, with factors density and treatment, and with theta = 1
m1 <- DI(y = "yield", density = "density", prop = 4:7, treat = "nitrogen",
FG = c("G","G","L","L"), DImodel = "FG", data = Switzerland)
## Extract the FG terms
extract(m1, what = "FG")
Run the code above in your browser using DataLab