Create an objective function including even distribution and neighbor balance introduced by Piepho 2018.
Usage
objective_function_piepho(
design,
swap,
spatial_cols,
current_score_obj = NULL,
swapped_items = NULL,
pair_mapping = NULL,
row_column = "row",
col_column = "col",
...
)Arguments
- design
A data frame representing the spatial information of the design
- swap
A column name of the items to be swapped
- spatial_cols
Column name(s) of the spatial factors
- current_score_obj
A named list containing the current score
- swapped_items
The items that had just been swapped
- pair_mapping
A named vector of pairs generated from create_pair_mapping
- row_column
Name of column representing the row of the design (default: "row")
- col_column
Name of column representing the column of the design (default: "col")
- ...
Extra parameters passed from speed
Value
A function which returns a named list of numeric values with one required name score representing
the score of the design (lower is better) with a signature function(design_df, swap, spatial_cols, ...).
See signature
details in objective_function_signature.
References
Piepho, H. P., Michel, V., & Williams, E. (2018). Neighbor balance and evenness of distribution of treatment replications in row-column designs. Biometrical journal. Biometrische Zeitschrift, 60(6), 1172–1189. https://doi.org/10.1002/bimj.201800013
Examples
design_df <- initialise_design_df(
items = c(1, 2, 2, 1, 3, 3, 1, 3, 3),
nrows = 3,
ncols = 3
)
pair_mapping <- create_pair_mapping(design_df$treatment)
objective_function_piepho(design_df, "treatment", c("row", "col"), pair_mapping = pair_mapping)
#> $score
#> [1] 16.36667
#>
#> $ed
#> $ed$`2`
#> $ed$`2`$msts
#> $ed$`2`$msts$`2`
#> [1] 1
#>
#>
#> $ed$`2`$min_mst
#> [1] 1
#>
#> $ed$`2`$min_items
#> [1] "2"
#>
#>
#> $ed$`4`
#> $ed$`4`$msts
#> $ed$`4`$msts$`3`
#> [1] 3
#>
#>
#> $ed$`4`$min_mst
#> [1] 3
#>
#> $ed$`4`$min_items
#> [1] "3"
#>
#>
#> $ed$`3`
#> $ed$`3`$msts
#> $ed$`3`$msts$`1`
#> [1] 2
#>
#>
#> $ed$`3`$min_mst
#> [1] 2
#>
#> $ed$`3`$min_items
#> [1] "1"
#>
#>
#>
#> $bal
#> [1] 8
#>
#> $adj
#> [1] 7
#>
#> $nb
#> $nb$nb
#> sorted_pairs
#> 1,1 1,2 1,3 2,2 2,3 3,3
#> 2 1 2 1 2 4
#>
#> $nb$max_nb
#> [1] 4
#>
#> $nb$max_pairs
#> [1] "3,3"
#>
#> $nb$var
#> [1] 1.2
#>
#>
# usage in speed, speed(..., obj_function = objective_function_piepho, pair_mapping = pair_mapping)