| arithscan {qtl} | R Documentation |
These binary operators perform arithmetic on genome scan result objects.
scan1+scan2 scan1-scan2
scan1, scan2 |
Genome scan results on the same set of chromosomes and markers. |
This is used to calculate the sum or difference of LOD scores of two
genome scan
results. The basic arithmetic operators were overrided to work on
objects of class scanone, scantwo, and scanqtl.
It currently only works for scanone results.
The same data structure as the input objects, with LOD scores recalcuated.
Karl W Broman, kbroman@jhsph.edu
data(fake.bc)
fake.bc <- calc.genoprob(fake.bc)
# genome scans without covariates
out.nocovar <- scanone(fake.bc)
# genome scans with covariates
ac <- fake.bc$pheno[,c("sex","age")]
ic <- fake.bc$pheno[,"sex"]
out.covar <- scanone(fake.bc, pheno.col=1,
addcovar=ac, intcovar=ic)
# plot the difference of with and without additive covariate
plot(out.covar-out.nocovar)