GBNlikelihood: Loglikelihhood

In andreamrau/GBNcausal: Causal estimation for Gaussian Bayesian networks

Description

Compute the loglikelihood between a GBN and data.

Usage

GBNlikelihood(GBN, data)

Arguments

GBN	GBN - an object of type GBN. colnames and rownames of the Weightmatrix must be the same as the colnames and rownamesof elements from data.
data	data - Data is a format obtained by the function dataFormat.

Value

It returns the loglikelihood of GBN and data.

Examples

# Data creation

seed = 1990

p <- 10
m<-rep(0,10)
sigma<-rep(0.1,10) 

W <- 1*upper.tri(matrix(0,p,p))

data <- dataCreate(nbData = 2*p, p = 10,KO = list(1,9), nbKO = c(p,p), W = W , m = m,sigma = sigma, seed = seed)$data

# Initial Value

W1=1*upper.tri(matrix(0,p,p)) 
m1=rep(0,p)
s1=rep(10e-4,p)
colnames(W1)=names(m1)=names(s1)=rownames(W1)=paste("N",1:p,sep="")

GBnetwork = new("GBNetwork",WeightMatrix=W1,resMean=m1,resSigma=s1)
GBNlikelihood(GBnetwork,data)

andreamrau/GBNcausal documentation built on May 12, 2019, 3:34 a.m.