qpHist: Histograms of non-rejection rates

In rcastelo/qpgraph: Estimation of genetic and molecular regulatory networks from high-throughput genomics data

Histograms of non-rejection rates

Description

Plots the distribution of non-rejection rates.

Usage

qpHist(nrrMatrix, A=NULL,
       titlehist = "all estimated\nnon-rejection rates", freq=TRUE)

Arguments

nrrMatrix	matrix of non-rejection rates.
A	adjacency matrix of an undirected graph whose present and missing edges will be employed to show separately the distribution of non-rejection rates.
titlehist	main title of the histogram(s).
freq	logical; if TRUE, the histograms show frequencies (counts) of occurrence of the different non-rejection rate values; if FALSE, then probability densities are plotted

Details

This function plots histograms using the R-function hist and therefore the way they are displayed follows that of this R-function.

Value

None

Author(s)

R. Castelo and A. Roverato

References

Castelo, R. and Roverato, A. A robust procedure for Gaussian graphical model search from microarray data with p larger than n, J. Mach. Learn. Res., 7:2621-2650, 2006.

See Also

qpNrr qpAvgNrr qpEdgeNrr qpGraphDensity qpClique

Examples

require(mvtnorm)

nVar <- 50  ## number of variables
maxCon <- 5 ## maximum connectivity per variable
nObs <- 30  ## number of observations to simulate

A <- qpRndGraph(p=nVar, d=maxCon)
Sigma <- qpG2Sigma(A, rho=0.5)
X <- rmvnorm(nObs, sigma=as.matrix(Sigma))

nrr.estimates <- qpNrr(X, q=5, verbose=FALSE)

qpHist(nrr.estimates, A)

rcastelo/qpgraph documentation built on Oct. 28, 2024, 5:15 a.m.