Description Usage See Also Examples
Retrieve MCMC samples for hyper dataset from R/qtl. Code provided to recreate qbHyper if desired. These samples are used in vignettes and demos.
1 |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 | data(qbHyper)
summary(qbHyper)
## Not run:
## Here is code to generate qbHyper.
## Get data from R/qtl.
data(hyper)
## Restrict to autosomes.
hyper <- subset(hyper, chr = 1:19)
## Calculate genotype probabilities.
hyper <- qb.genoprob(hyper, step=2)
## Create MCMC samples.
qbHyper <- qb.mcmc(hyper, n.thin = 40, seed = 1616)
## The next line saves qbHyper as an external binary file.
save("qbHyper", file = "qbHyper.RData")
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.