plotRsq: Plot R square data

In eqtl: Tools for analyzing eQTL experiments: A complementary to Karl Broman's 'qtl' package for genome-wide analysis

Description

Draw histograms of R square value distribution for rsq object.

Usage

plotRsq( rsq, par=c(2,2), ...)

Arguments

rsq	An object of class rsq. See calc.Rsq for function details.
par	A vector of two integers corresponding to the mfrow parameter of the par() function.
...	Passed to the function hist and par when they are called.

Details

Draw histograms of R square value distribution from an object of class rsq. Three histograms are drawn: the first one shows the R square value distribution of single QTLs. The second shows the distribution for QTL interactions. The last one shows all R square values distribution.

Value

none

Author(s)

Hamid A. Khalili

See Also

calc.Rsq,peak.2.array

Examples

data(seed10);

seed10 <- calc.genoprob( cross=seed10, step=2, off.end=0, error.prob=0,
	map.function='kosambi', stepwidth='fixed');
seed10 <- sim.geno( cross=seed10, step=2, off.end=0, error.prob=0,
	map.function='kosambi', stepwidth='fixed');

out.em <- scanone( seed10, pheno.col=1:50, model='normal', method='hk');
out.peak <- define.peak(out.em,'all');
out.rsq <- calc.Rsq(seed10,out.peak);

# plotting R quare data
plotRsq(out.rsq);
plotRsq(out.rsq,par=c(1,3));

eqtl documentation built on May 2, 2019, 5:42 a.m.