varplot: Variance component plot

Description Usage Arguments Author(s) See Also Examples

View source: R/varplot.R

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

1
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

1
2
3
4
5
6
7
# 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 May 2, 2018, 2:46 a.m.