Normal-Block Model with Unknown Clustering

R6 class for a normal-block model with a fixed number of clusters (but unknown clustering).

Super classes

NormalBlockBase -> NormalBlockVarBase -> NormalBlockVarUnknownClusters

Public fields

fixed_tau

whether tau should be fixed at clustering_init during optimization, useful for stability selection

Active bindings

model_par

a list with the matrices of the model parameters: B (covariates), dm1 (species variance), Omega (groups precision matrix))

nb_param

number of parameters in the model

var_par

a list with the matrices of the variational parameters: M (means), S (variances), tau (posterior group probabilities)

entropy

Entropy of the conditional distribution

fitted

Y values predicted by the model, in Y's original units

who_am_I

a method to print what model is being fitted

Methods

Public methods

• NormalBlockVarUnknownClusters$new()

• NormalBlockVarUnknownClusters$clone()

Inherited methods

• NormalBlockBase$best_of_inits()
• NormalBlockBase$candidates_merge()
• NormalBlockBase$candidates_split()
• NormalBlockBase$latent_network()
• NormalBlockBase$optimize()
• NormalBlockBase$plot()
• NormalBlockBase$plot_loglik()
• NormalBlockBase$plot_network()
• NormalBlockBase$predict()
• NormalBlockBase$print()
• NormalBlockBase$update()
• NormalBlockVarBase$merge()
• NormalBlockVarBase$split()
• NormalBlockVarBase$warm_start_from()

---

NormalBlockVarUnknownClusters$new()

Create a new [`NormalBlockVarUnknownClusters`] object.

Usage

NormalBlockVarUnknownClusters$new(
  data,
  q,
  sparsity = 0,
  control = NB_control()
)

Arguments

data

contains the matrix of responses (Y) and the design matrix (X).

q

required number of groups

sparsity

sparsity penalty to add on blocks precision matrix for sparsity

control

structured list for specific parameters

Returns

A new [`NormalBlockVarUnknownClusters`] object

---

NormalBlockVarUnknownClusters$clone()

The objects of this class are cloneable with this method.

Usage

NormalBlockVarUnknownClusters$clone(deep = FALSE)

Arguments

deep

Whether to make a deep clone.

Examples

ex <- generate_normal_block_var_data(n = 50, p = 20, d = 1, q = 3)
data <- NormalBlockData$new(ex$Y, ex$X)
model <- normal_block(data, blocks = 3, control = NB_control(verbose = FALSE))
model$clustering
#>  [1] 1 2 1 3 1 2 3 3 2 2 1 1 1 3 3 1 3 1 2 3