inst/extdata/Alfalfa.R

In forsearch: Diagnostic Analysis Using Forward Search Procedure for Various Models

Alfalfa <-
structure(list(Variety = structure(c(2L, 2L, 2L, 2L, 2L, 2L, 
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 3L, 3L, 3L, 3L, 3L, 3L, 
3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 
3L, 3L), .Label = c("Cossack", "Ladak", "Ranger"), class = "factor"), 
    Date = structure(c(1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 
    3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 
    2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 
    1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 
    4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 
    3L, 4L), .Label = c("None", "S1", "S20", "O7"), class = "factor"), 
    Block = structure(c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 
    3L, 3L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 1L, 
    1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 
    5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 
    2L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 6L, 6L, 
    6L, 6L), .Label = c("1", "2", "3", "4", "5", "6"), class = "factor"), 
    Yield = c(2.1699999999999999, 1.5800000000000001, 2.29, 2.23, 
    1.8799999999999999, 1.26, 1.6000000000000001, 2.0099999999999998, 
    1.6200000000000001, 1.22, 1.6699999999999999, 1.8200000000000001, 
    2.3399999999999999, 1.5900000000000001, 1.9099999999999999, 
    2.1000000000000001, 1.5800000000000001, 1.25, 1.3899999999999999, 
    1.6599999999999999, 1.6599999999999999, 0.93999999999999995, 
    1.1200000000000001, 1.1000000000000001, 2.3300000000000001, 
    1.3799999999999999, 1.8600000000000001, 2.27, 2.0099999999999998, 
    1.3, 1.7, 1.8100000000000001, 1.7, 1.8500000000000001, 1.8100000000000001, 
    2.0099999999999998, 1.78, 1.0900000000000001, 1.54, 1.3999999999999999, 
    1.4199999999999999, 1.1299999999999999, 1.6699999999999999, 
    1.3100000000000001, 1.3500000000000001, 1.0600000000000001, 
    0.88, 1.0600000000000001, 1.75, 1.52, 1.55, 1.5600000000000001, 
    1.95, 1.47, 1.6100000000000001, 1.72, 2.1299999999999999, 
    1.8, 1.8200000000000001, 1.99, 1.78, 1.3700000000000001, 
    1.5600000000000001, 1.55, 1.3100000000000001, 1.01, 1.23, 
    1.51, 1.3, 1.3100000000000001, 1.1299999999999999, 1.3300000000000001
    )), row.names = 1:72, class = c("nmGroupedData", "groupedData", 
"data.frame"), formula = Yield ~ Date | Block/Variety, formulaList = list(
    Block = ~Block, Variety = ~Variety), order.groups = list(
    Block = TRUE, Variety = TRUE), FUN = structure(function (x) 
max(x, na.rm = TRUE), source = "function (x) max(x, na.rm = TRUE)"), labels = list(
    y = "Yield in 1944 following third date of cutting in 1943"), units = list(
    y = "(T/Acre)"))