convertOutput: Convert output of anal.leaps, anal.mcmc or perm

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

Convert output of anal.leaps, anal.mcmc or perm

Description

Convert the output of anal.leaps, anal.mcmc or perm for use with the function which.correct.

Usage

convertOutput(object, n.mar=rep(11,9))

Arguments

object	An object of the form output from anal.leaps, anal.mcmc, perm; namely, a vector of marker identifiers, as integers, (or a list whose first component is such a vector).
n.mar	Vector indicating the number of markers on each chromosome.

Value

A matrix whose columns are identified QTLs, and rows are chromosome and marker numbers, respectively.

Author(s)

Karl W Broman, [e-mail]

See Also

anal.leaps, anal.mcmc, perm, which.correct

Examples

dat <- simbc(n.ind=250)
res1 <- anal.leaps(dat, "forward", 2.1, 27)
res2 <- anal.leaps(dat, "backward", 2.1, 27)
res <- lapply(list(res1,res2), convertOutput)
out <- which.correct(res)

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