OptimalDesign (version 1.0.1)

Fx_blocks: Matrix of candidate regressors for a block size-two model

Description

Creates the matrix of all candidate regressors of a linear regression model corresponding to the problem of the optimal block size-two design.

Usage

Fx_blocks(n.treats, blocks=NULL, echo=TRUE)

Value

the n times m matrix of all candidate regressors of an auxiliary linear regression model corresponding to the problem of the optimal block size-two design (n is ncol(blocks), m is n.treats-1).

Arguments

n.treats

the number of "treatments" in the block experiment.

blocks

the 2 times n matrix of all permissible blocks (that is, permissible pairings of treatments). If blocks=NULL, blocks is set to combn(n.treats, 2), which means that all treatment pairings are permissible.

echo

Print the call of the function?

Author

Radoslav Harman, Lenka Filova

Details

Creates the matrix Fx of artificial regressors, such that the D- and A-optimal designs for the corresponding artificial LRM are are the same as what is called the D- and A-optimal design in the original block model with blocks of size two.

References

Harman R, Filova, L: Computing efficient exact designs of experiments using integer quadratic programming, Computational Statistics and Data Analysis 71 (2014) 1159-1167.

Sagnol G, Harman R: Computing Exact D-optimal designs by mixed integer second-order cone programming, The Annals of Statistics 43 (2015), 2198-2224.

See Also

Fx_cube, Fx_simplex, Fx_glm, Fx_dose, Fx_survival

Examples

if (FALSE) {
# Compute a D-efficient block size-two design
# with 15 treatments and 10 blocks of size two

Fx <- Fx_blocks(10)
w <- od_KL(Fx, 15, t.max = 5)$w.best
des <- combn(10, 2)[, as.logical(w)]
print(des)


# We can visualize the design as a graph
library(igraph)
grp <- graph_(t(des), from_edgelist(directed = FALSE))
plot(grp, layout=layout_with_graphopt)
}