plotCGpvalues: Plot simultaneous p-values for chain graphs

Description

This function plots the simultaneous p-values for SIN model selection for chain graphs.

Usage

plotCGpvalues(blocks, pvals, legend=TRUE, legendpos=NULL)

Arguments

blocks

a list of integer vectors with entries amongst $1,\ldots,p$ where $p$ is the number of variables. Each one of the integer vectors specifies a set of variables that form a block in the chain graph. Furthermore, a partial ordering of the variables is specified by the convention that variables in one block are ordered smaller than variables in a block succeeding in the list blocks.

pvals

a matrix of p-values.

legend

a Boolean.

legendpos

a numerical vector of length 2.

Value

A plot illustrating the simultaneous p-values. If legend==TRUE, then a legend giving information on the variable labelling is drawn. Optionally, the position of the legend can be specified by legendpos.

References

Drton, M. \& Perlman, M.D. (2004) Model Selection for Gaussian Concentration Graphs. Biometrika 91(3): 591-602. Drton, M. \& Perlman, M.D. (2008) A SINful Approach to Gaussian Graphical Model Selection. J. Statist. Plann. Inference 138(4): 1179-1200.

Examples

Run this code

data(fowlbones)
blocks <- list(1:2,3:4,5:6)
pvals <- holm(sinCG(blocks,fowlbones$corr,fowlbones$n, type="AMP"))
## Not run: plotCGpvalues(blocks, pvals)
## Not run: plotCGpvalues(blocks, pvals, legend=FALSE)
## Not run: plotCGpvalues(blocks, pvals, legendpos=c(7,0.5))
## Not run: plotCGpvalues(blocks, pvals, legend=TRUE, legendpos=c(7,0.5))