plotRsq {eqtl} | R Documentation |
Draw histograms of R square value distribution for rsq
object.
plotRsq( rsq, par=c(2,2), ...)
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. |
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.
none
Hamid A. Khalili
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='em'); 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));