fitqtl
.
makeqtl(cross, chr, pos, qtl.name, what=c("draws","prob"))
cross
. See
read.cross
for details.chr
) indicating the
positions on the chromosome to be taken. If there is no marker or
pseudomarker at a position, the nearest position is used.fitqtl
.
If unspecified, the names will be of the form "Chr1@10"
for a
QTL on Chromsome 1 at 10 cM.qtl
with the following elements (though only
one of geno
and prob
will be included, according to
whether what
is given as "draws"
or "prob"
):
geno
prob
name
"Chr1@10"
.altname
"Q1"
, "Q2"
,
etc.chr
pos
n.qtl
n.ind
n.gen
cross
object. If both
fields are missing in the input object, the function will report an
error. Before running this function, the user must have first run either
sim.geno
(for what="draws"
) or
calc.genoprob
(for what="prob"
).
fitqtl
, calc.genoprob
,
sim.geno
, dropfromqtl
,
replaceqtl
, addtoqtl
, summary.qtl
,
reorderqtl
data(fake.f2)
# take out several QTLs and make QTL object
qc <- c("1", "6", "13")
qp <- c(25.8, 33.6, 18.63)
fake.f2 <- subset(fake.f2, chr=qc)
fake.f2 <- sim.geno(fake.f2, n.draws=8, step=2, err=0.001)
qtl <- makeqtl(fake.f2, qc, qp, what="draws")
summary(qtl)
Run the code above in your browser using DataLab