n_params_junction_tree: Determines the number of free parameters

In nvihrs14/tcherry: Learning the structure of tcherry trees

Description

Determines the number of free parameters in a graphical model based on a junction tree for the graph.

Usage

n_params_junction_tree(cliques, separators, data)

Arguments

cliques	A list containing the cliques of the junction tree.
separators	A list containing the separators of the junction tree.
data	A data frame with observations of the variables.

Value

The number of free parameters.

Author(s)

Katrine Kirkeby, [e-mail]

Maria Knudsen, [e-mail]

Ninna Vihrs, [e-mail]

Examples

set.seed(43)
var1 <- c(sample(c(1, 2), 100, replace = TRUE))
var2 <- var1 + c(sample(c(1, 2), 100, replace = TRUE))
var3 <- var1 + c(sample(c(0, 1), 100, replace = TRUE,
                        prob = c(0.9, 0.1)))
var4 <- c(sample(c(1, 2), 100, replace = TRUE))
var5 <- var2 + var3
var6 <- var1 - var4 + c(sample(c(1, 2), 100, replace = TRUE))
var7 <- c(sample(c(1, 2), 100, replace = TRUE))

data <- data.frame("var1" = as.character(var1),
                   "var2" = as.character(var2),
                   "var3" = as.character(var3),
                   "var4" = as.character(var4),
                   "var5" = as.character(var5),
                   "var6" = as.character(var6),
                   "var7" = as.character(var7))

cliques <- list(c("var1", "var2", "var3"),
                c("var2", "var3", "var5"),
                c("var5", "var6", "var7"),
                c("var1", "var4"),
                c("var2", "var5", "var6"))

separators <- list(c("var1"),
                   c("var2", "var3"),
                   c("var2", "var5"),
                   c("var5", "var6"))

n_params_junction_tree(cliques, separators, data)

nvihrs14/tcherry documentation built on Aug. 1, 2020, 6:25 p.m.