genmodel runs 4ti2's genmodel program to compute the configuration matrix A corresponding to graphical statistical models given by a simplicial complex and levels on the nodes.
genmodel(varlvls, facets, dir = tempdir(), quiet = TRUE, shell = FALSE, ...)
varlvls | a vector containing the number of levels of each variable |
---|---|
facets | the facets generating the hierarchical model, a list of vectors of variable indices |
dir | Directory to place the files in, without an ending / |
quiet | If FALSE, messages the 4ti2 output |
shell | Messages the shell code used to do the computation |
... | Additional arguments to pass to the function |
The configuration matrix of the model provided
if (has_4ti2()) { varlvls <- rep(2, 2) facets <- list(1, 2) genmodel(varlvls, facets) genmodel(varlvls, facets, quiet = FALSE) varlvls <- rep(3, 3) facets <- list(1:2, 2:3, c(3,1)) genmodel(varlvls, facets) # compare this to algstat's hmat function }#>#> #> #> #> #> #> #> #>#> [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] #> [1,] 1 0 0 0 0 0 0 0 0 1 0 0 0 #> [2,] 0 1 0 0 0 0 0 0 0 0 1 0 0 #> [3,] 0 0 1 0 0 0 0 0 0 0 0 1 0 #> [4,] 0 0 0 1 0 0 0 0 0 0 0 0 1 #> [5,] 0 0 0 0 1 0 0 0 0 0 0 0 0 #> [6,] 0 0 0 0 0 1 0 0 0 0 0 0 0 #> [7,] 0 0 0 0 0 0 1 0 0 0 0 0 0 #> [8,] 0 0 0 0 0 0 0 1 0 0 0 0 0 #> [9,] 0 0 0 0 0 0 0 0 1 0 0 0 0 #> [10,] 1 1 1 0 0 0 0 0 0 0 0 0 0 #> [11,] 0 0 0 1 1 1 0 0 0 0 0 0 0 #> [12,] 0 0 0 0 0 0 1 1 1 0 0 0 0 #> [13,] 0 0 0 0 0 0 0 0 0 1 1 1 0 #> [14,] 0 0 0 0 0 0 0 0 0 0 0 0 1 #> [15,] 0 0 0 0 0 0 0 0 0 0 0 0 0 #> [16,] 0 0 0 0 0 0 0 0 0 0 0 0 0 #> [17,] 0 0 0 0 0 0 0 0 0 0 0 0 0 #> [18,] 0 0 0 0 0 0 0 0 0 0 0 0 0 #> [19,] 1 0 0 1 0 0 1 0 0 0 0 0 0 #> [20,] 0 0 0 0 0 0 0 0 0 1 0 0 1 #> [21,] 0 0 0 0 0 0 0 0 0 0 0 0 0 #> [22,] 0 1 0 0 1 0 0 1 0 0 0 0 0 #> [23,] 0 0 0 0 0 0 0 0 0 0 1 0 0 #> [24,] 0 0 0 0 0 0 0 0 0 0 0 0 0 #> [25,] 0 0 1 0 0 1 0 0 1 0 0 0 0 #> [26,] 0 0 0 0 0 0 0 0 0 0 0 1 0 #> [27,] 0 0 0 0 0 0 0 0 0 0 0 0 0 #> [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25] #> [1,] 0 0 0 0 0 1 0 0 0 0 0 0 #> [2,] 0 0 0 0 0 0 1 0 0 0 0 0 #> [3,] 0 0 0 0 0 0 0 1 0 0 0 0 #> [4,] 0 0 0 0 0 0 0 0 1 0 0 0 #> [5,] 1 0 0 0 0 0 0 0 0 1 0 0 #> [6,] 0 1 0 0 0 0 0 0 0 0 1 0 #> [7,] 0 0 1 0 0 0 0 0 0 0 0 1 #> [8,] 0 0 0 1 0 0 0 0 0 0 0 0 #> [9,] 0 0 0 0 1 0 0 0 0 0 0 0 #> [10,] 0 0 0 0 0 0 0 0 0 0 0 0 #> [11,] 0 0 0 0 0 0 0 0 0 0 0 0 #> [12,] 0 0 0 0 0 0 0 0 0 0 0 0 #> [13,] 0 0 0 0 0 0 0 0 0 0 0 0 #> [14,] 1 1 0 0 0 0 0 0 0 0 0 0 #> [15,] 0 0 1 1 1 0 0 0 0 0 0 0 #> [16,] 0 0 0 0 0 1 1 1 0 0 0 0 #> [17,] 0 0 0 0 0 0 0 0 1 1 1 0 #> [18,] 0 0 0 0 0 0 0 0 0 0 0 1 #> [19,] 0 0 0 0 0 0 0 0 0 0 0 0 #> [20,] 0 0 1 0 0 0 0 0 0 0 0 0 #> [21,] 0 0 0 0 0 1 0 0 1 0 0 1 #> [22,] 0 0 0 0 0 0 0 0 0 0 0 0 #> [23,] 1 0 0 1 0 0 0 0 0 0 0 0 #> [24,] 0 0 0 0 0 0 1 0 0 1 0 0 #> [25,] 0 0 0 0 0 0 0 0 0 0 0 0 #> [26,] 0 1 0 0 1 0 0 0 0 0 0 0 #> [27,] 0 0 0 0 0 0 0 1 0 0 1 0 #> [,26] [,27] #> [1,] 0 0 #> [2,] 0 0 #> [3,] 0 0 #> [4,] 0 0 #> [5,] 0 0 #> [6,] 0 0 #> [7,] 0 0 #> [8,] 1 0 #> [9,] 0 1 #> [10,] 0 0 #> [11,] 0 0 #> [12,] 0 0 #> [13,] 0 0 #> [14,] 0 0 #> [15,] 0 0 #> [16,] 0 0 #> [17,] 0 0 #> [18,] 1 1 #> [19,] 0 0 #> [20,] 0 0 #> [21,] 0 0 #> [22,] 0 0 #> [23,] 0 0 #> [24,] 1 0 #> [25,] 0 0 #> [26,] 0 0 #> [27,] 0 1