| qbExample {qtlbim} | R Documentation |
Retrieve or recreate MCMC samples used in qtlbim.pdf document and demos.
data(qbExample)
This example is used in vignettes and demos.
See vignette qtlbim.pdf or see qtlbim.Rnw in doc folder of package.
qb.genoprob, qb.mcmc, qb.sim.cross
data(qbExample)
summary(qbExample)
## Not run:
## Fake BC Example.
set.seed(1234)
cross <- qb.sim.cross(len = rep(60,3), n.mar = 7, eq.spacing =FALSE,
n.ind = 100, type = "bc", ordinal = c(0.3,0.3,0.2,0.2),
missing.geno = 0.03, missing.pheno = 0.03,
qtl.pos = rbind(qtl.1=c(chr=1,pos=15),qtl.2=c(1,45),qtl.3=c(2,12),qtl.4=c(3,15)),
qtl.main = rbind(main.1=c(qtl=1,add=1.5),main.2=c(2,0),
main3=c(3,-1),main4=c(4,0)),
qtl.epis = rbind(epis1=c(qtl.a=2,qtl.b=3,aa=-2),
epis2=c(2,4,3)),
covariate = c(fix.cov=0.5,ran.cov=0.07),
gbye = rbind(GxE.1=c(qtl=3,add=2)))
summary(cross$qtl)
cross <- qb.genoprob(cross, step=2)
qbExample <- qb.mcmc(cross, pheno.col = 3, rancov = 2, fixcov = 1,
chr.nqtl = rep(3,nchr(cross)), intcov = 1, seed = 1216,
interval = rep(10,3), n.iter = 1000, n.thin = 20)
save("qbExample", file = "qbExample.RData")
## End(Not run)