## Load the sim1 data
data(sim1)
## Fit the FULL model
m1 <- DI(y = "response", block = "block", prop = 3:6,
DImodel = "FULL", data = sim1)
## Compare with the reference model
DI_compare(m1, test = "F")
Run the code above in your browser using DataLab