sdnet (version 2.4.1)

cnFindBIC-method: Find Network by BIC

Description

This is a model selection routine that finds a network in a set of networks using the BIC criteria.

Usage

cnFindBIC(object, numsamples)

Arguments

object

A list of catNetworkNode objects or catNetworkEvaluate or dagEvaluate

numsamples

The number of samples used for estimating object

Value

A catNetwork object with optimal BIC value.

Details

The function returns the network with maximal BIC value from a list of networks as obtained from one of the search-functions cnSearchOrder, cnSearchSA and cnSearchSAcluster. The formula used for the BIC is log(Likelihood) - 0.5*Complexity*log(numNodes).

See Also

cnFindAIC, cnFind

Examples

# NOT RUN {
  cnet <- cnRandomCatnet(numnodes=12, maxpars=3, numcats=2)
  psamples <- cnSamples(object=cnet, numsamples=10)
  nodeOrder <- sample(1:12)
  nets <- cnSearchOrder(data=psamples, pert=NULL, 
	maxParentSet=2, maxComplexity=36, nodeOrder)
  bicnet <- cnFindBIC(object=nets, numsamples=dim(psamples)[2])
  bicnet
# }