varplot: Variance component plot

In maanova: Tools for analyzing Micro Array experiments

Description

This function plots the density curve of each variance component of a result from fitmaanova.

If the input is from fixed model ANOVA, it will plot one curve for error variance component. If the input is from mixed model ANOVA, it will plot multiple curves, one for a random term (including error).

Usage

varplot(anovaobj, xlab, ylab, main)

Arguments

anovaobj	An object of class maanova.
xlab	Figure xlab. Default is "Sigma".
ylab	Figure ylab. Default is "Density".
main	Figure title. Default is "Density plot for sqrt of variance".

Author(s)

Hao Wu

See Also

fitmaanova, density

Examples

# load abf1 data
data(abf1)
## Not run: 
fit.full.mix <- fitmaanova(abf1, formula = ~Strain+Sample, 
   random = ~Sample)
varplot(fit.full.mix)
## End(Not run)

Example output

Attaching package: 'maanova'

The following object is masked from 'package:base':

    norm

Calculating variance components for fixed model...
Fitting mixed effect model...
Finish gene number 100 ...
Finish gene number 200 ...
Finish gene number 300 ...
Finish gene number 400 ...
Finish gene number 500 ...
Warning messages:
1: In any(parsed.formula$random) :
  coercing argument of type 'double' to logical
2: In any(parsed.formula$random) :
  coercing argument of type 'double' to logical

maanova documentation built on Nov. 8, 2020, 8:21 p.m.