metagenomeSeq (version 1.14.0)

plotClassTimeSeries: Plot abundances by class

Description

Plot the abundance of values for each class using a spline approach on the estimated full model.

Usage

plotClassTimeSeries(res, formula, xlab = "Time", ylab = "Abundance", color0 = "black", color1 = "red", include = c("1", "class", "time:class"), ...)

Arguments

res

Output of fitTimeSeries function

formula

Formula for ssanova.

xlab

X-label.

ylab

Y-label.

color0

Color of samples from first group.

color1

Color of samples from second group.

include

Parameters to include in prediction.

...

Extra plotting arguments.

Value

Plot for abundances of each class using a spline approach on estimated null model.

See Also

fitTimeSeries

Examples

data(mouseData)
res = fitTimeSeries(obj=mouseData,feature="Actinobacteria",
   class="status",id="mouseID",time="relativeTime",lvl='class',B=10)
plotClassTimeSeries(res,pch=21,bg=res$data$class,ylim=c(0,8))