BDgraph (version 2.3)

plotLinks: Plot of posterior distribution for graphs according to number of their links

Description

This function is for ploting the posterior distribution of all possible graphs according to number of their links.

Usage

plotLinks(output, xlim = c(0, output$p * (output$p - 1) / 2))

Arguments

output

a list which is the result of BDMCMC algorithm from the 'bdmcmc', 'bdmcmc.low', or 'bdmcmc.high' functions.

xlim

the x limits (min, max) of the plot.

References

Mohammadi, A. and E. C. Wit (2012). Gaussian graphical model determination based on birth-death MCMC inference, arXiv:1210.5371v4. http://arxiv.org/abs/1210.5371v4

See Also

bdmcmc

Examples

p <- 8 # number of nodes 
  # "truK" is the precision matrix of true graph
  truK <- diag(p)
  for (i in 1:(p-1)) truK[i,i+1] <- truK[i+1,i] <- 0.5
  truK[1,p] <- truK[p,1] <- 0.4
  truK
  
  # generate the data (200 observations) from multivariate normal 
  # distribution with mean zero and percision matrix "truK"
  data <- mvrnorm(200, c(rep(0,p)), solve(truK))  
  
  output <- bdmcmc(data, meanzero = T, iter = 500)
  plotLinks(output)