Plot simultaneous p-values for acyclic directed graphs

Share:

Description

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

Usage

1
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

1
2
3
4
5
6
7
data(fowlbones)
p <- dim(fowlbones$corr)[1]
pvals <- holm(sinDAG(list(1:p),fowlbones$corr,fowlbones$n))
## Not run: plotDAGpvalues(pvals)
## Not run: plotDAGpvalues(pvals, legend=FALSE)
## Not run: plotDAGpvalues(pvals, legendpos=c(7,0.5))
## Not run: plotDAGpvalues(pvals, legend=TRUE, legendpos=c(7,0.5))