qbHyper {qtlbim} | R Documentation |
Retrieve MCMC samples for hyper dataset from R/qtl. Code provided to recreate qbHyper if desired. These samples are used in vignettes and demos.
data(qbHyper)
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)