qp.search {qp}R Documentation

Matrix of non-rejection rates

Description

Calculates the estimates of the non-rejection rates for every pair of variables.

Usage

qp.search(S, N, q=0, T=500, significance=0.05, binary=TRUE)

Arguments

S sample variance-covariance matrix
N sample size
T number of tests per adjacency
q partial-correlation order
significance significance level of each test
binary flag to switch to the compiled C code

Details

By default binary=TRUE and the compiled and faster C code corresponding to this function will be executed. If binary=FALSE is set, then the R code will be executed.

Value

A matrix with the acceptance test counts
T number of tests per adjacency (copied from the input parameter)

Author(s)

Robert Castelo and Alberto Roverato

References

Castelo, R. and Roverato, A. (2006). A robust procedure for Gaussian graphical model search from microarray data with p larger than n, J. Mach. Learn. Res., accepted

See Also

qp.edge.prob,qp.analyse,qp.hist

Examples

data(jmlr06data)
S <- S.bd5.N20
N <- 20
q <- 6
T <- 100
qp.out <- qp.search(S,N,q,T)

[Package qp version 0.1-1 Index]