which.correct: Determine which QTLs were correctly identified in an analysis

In kbroman/qtlsim: Software for Testing QTL Analysis via Simulations

Determine which QTLs were correctly identified in an analysis

Description

Given the results of a model selection procedure and the true underlying locations of QTLs, determine which QTLs were correctly identified and how many were extraneous.

Usage

which.correct(results, truth=rbind(c(1,1,2,2,3,4,5),c(4,8,4,8,6,4,1)),
              within=1)

Arguments

results	A list (or a list of lists), with each component in the inner lists being either numeric(0) or a matrix whose columns are identified QTLs, and rows are chromosome and marker numbers, respectively.
truth	Matrix whose columns are the true QTLs, with rows giving the chromosome and marker numbers, respectively.
within	Distance (in number of markers) allowable between true and inferred QTLs for the inferred QTL to be considered correct.

Value

A list of matrices, with the matrices of size n.sim x (n.qtl+2). The first set of columns in each matrix indicate which QTLs were correctly identified. The last two columns give the number of extraneous loci on the same chromosome as a QTL and the number of unlinked extraneous loci.

Author(s)

Karl W Broman, [e-mail]

See Also

anal.multi, anal.mcmc

Examples

## Not run: dat <- simbc()
## Not run: output1 <- anal.multi(dat)
## Not run: output2 <- which.correct(output1)

kbroman/qtlsim documentation built on May 17, 2023, 11:53 p.m.