LikelihoodHistory: Calculating and Plot the history of log-Likelihood value

Description

This function is a supplementary function for calculating and ploting the history of log-Likelihood value for each iteration. This function is used for tracking the evolution of the MCMC chain obtained from functions Networks.STD() and Networks.Fast()

Usage

LikelihoodHistory(Trace,pvalue,Status=FALSE,True.Node)

Arguments

Trace

the trace (results of each iteration) generated by the functions Networks.STD() and Networks.Fast().Please set as Trace=total$trace where total is the output of the two functions, Networks.STD() and Networks.Fast()

pvalue

a vector of p-values obtained from large scale statistical hypothesis testing. It should be the same vector which you set for the argument pvalue in functions Networks.STD() and Networks.Fast()

Status

a logical variable indicating whether you have the assumed true value. TRUE means you have the assumed true value for plotting and FALSE, otherwise. The default setting is FALSE

True.Node

a binary(0/1) vector providing the assumed information of each node. The length of vector should be equal to the length of the vector of pvalue

Value

a table indicating the history of the history of log-Likelihood value for each iteration

References

Yize Zhao, Jian Kang, Tianwei Yu (2014) A Bayesian nonparameteric model for selecting gene and gene subnetwork, Annals of Applied Statistics, in press.

Zhou Lan, Jian Kang, Tianwei Yu, Yize Zhao, BANFF: an R package for network identifications via Bayesian nonparametric mixture models, working paper.