SIN (version 0.4)

plotDAGpvalues: Plot simultaneous p-values for acyclic directed graphs

Description

This function plots the simultaneous p-values for SIN model selection for acyclic directed graphs (DAGs).

Usage

plotDAGpvalues(pvals, legend=TRUE, legendpos=NULL)

Arguments

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.

See Also

sinDAG

Examples

data(fowlbones)
p <- dim(fowlbones$corr)[1]
pvals <- holm(sinDAG(list(1:p),fowlbones$corr,fowlbones$n))
plotDAGpvalues(pvals)
plotDAGpvalues(pvals, legend=FALSE)
plotDAGpvalues(pvals, legendpos=c(7,0.5))
plotDAGpvalues(pvals, legend=TRUE, legendpos=c(7,0.5))