Last data update: 2014.03.03

R: plot the compositional bias for the sub-taxinomic rank
compbiasPlotR Documentation

plot the compositional bias for the sub-taxinomic rank

Description

.

Usage

compbiasPlot(x, pair=nv(levels(x$samples$group)[1:2] , c('ref','obs')), meta.lev='phylum',  meta.lev.lim=nrow(x$meta.sum[[meta.lev]]), breaks=10, xlim=c(0,.5), type='box', col=NULL, ...)

Arguments

x

a MANTA object.

pair

named vector of the pair of conditions of interest

meta.lev

the taxinomic rank to plot.

meta.lev.lim

the taxinomic rank for diversity measure.

breaks

breaks for the histogram.

xlim

x axis limits.

type

type of plot.

col

colors

...

additional params.

Value

a compositional bias plot

See Also

DGEList, manta

Examples


manta.path <- system.file("extdata","PapaGO-BLAST.results-diatoms.Rdata", package="manta")
load(manta.path)
			
compbiasPlot(x, meta.lev='genus_sp')

Results


R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(manta)
Loading required package: edgeR
Loading required package: limma
> png(filename="/home/ddbj/snapshot/RGM3/R_BC/result/manta/compbiasPlot.Rd_%03d_medium.png", width=480, height=480)
> ### Name: compbiasPlot
> ### Title: plot the compositional bias for the sub-taxinomic rank
> ### Aliases: compbiasPlot
> 
> ### ** Examples
> 
> 
> manta.path <- system.file("extdata","PapaGO-BLAST.results-diatoms.Rdata", package="manta")
> load(manta.path)
> 			
> compbiasPlot(x, meta.lev='genus_sp')
> 
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>