BMNJT: Junction tree algorithm for pairwise binary Markov networks

In BMN: The pseudo-likelihood method for pairwise binary markov networks

Description

Calculates the expections and second moments for a pairwise binary Markov network using the junction tree algorithm.

Usage

BMNJT(thetaMat, adjMat = NULL, var = NULL, onlyActive = FALSE, timeout = 60)

Arguments

thetaMat	Theta matrix of the model (assumed to be symmetric).
adjMat	Adjacency matrix of the underlying graph. If not given, inferred from thetaMat. If given, values in thetaMat that do not correspond to an edge in the graph are being ignored.
var	if given, the expectation and second moments of variable var only will be calculated. Otherwise, for all variables.
onlyActive	If TRUE, only values in the second moment matrix corresponding to edges in the graph will be calculated (much faster). Otherwise, all entries will be calculated.
timeout	Number of seconds until the functions terminates automatically.

Details

Uses the junction tree algorithm to calculate the matrix of second moments of the underlying pairwise binary Markov model. This is useful for inference on this type of models as the matrix of second moments is the derivative of the partition function in this class of models. The function provides the option to calculate the second moments only for one variable with all the others (variables are numbered from 1 to p; see var). Also, the juntion tree algorithm is by a factor p more efficient if only second moments corresponding to edges in the underlying graph (non-zero values in the thetaMat or adjMat) are being calculated. This can be done using the switch onlyActive.

Value

The return value is a list with elements:

Expectation	Expectation of the variables.
SecondMomentMatrix	Second moments of the variables (only present if Var=NULL).
SecondMomentVector	Vector of second moments of variable Var with the others if Var != NULL.

Author(s)

Holger Hoefling

See Also

BMNPseudo, BMNExact

Examples

library(BMN)
Theta = matrix(numeric(25), ncol=5)
Theta[1,1]=0.5; Theta[2,2]=0.5; Theta[3,3]=0; Theta[4,4]= -0.5; Theta[5,5]= 0.5;
Theta[1,2]=Theta[2,1]=1; Theta[1,4]=Theta[4,1]=1; Theta[2,3]=Theta[3,2]= -1;

BMNJT(Theta)

BMN documentation built on May 29, 2017, 10:55 a.m.