varHistPlot: Variance histogram and density

In plw: Probe level Locally moderated Weighted t-tests.

Description

Will produce a histogram of observed variance estimators (logged) together with the fitted density.

Usage

varHistPlot(model,main="Histogram variance estimators",
           histCol=8,densityCol=2,drawLegend=TRUE)

Arguments

model	On object obtained from the function plw or lmw.
main	Main title of plot.
histCol	Color for histogram bars.
densityCol	Color for density function.
drawLegend	To draw a legend or not.

Author(s)

Magnus Astrand

See Also

plw, lmw

Examples

# ------------------------------------------
# Example using the result of the analysis of
# the 6 arrays in the AffySpikeU95Subset data set

# Loading the data
data(AffySpikeU95Subset)

# Defining design and contrast matrix
group<-factor(rep(1:2,each=3))
design<-model.matrix(~group-1)
contrast<-matrix(c(1,-1),1,2)

# Analyzing using plw
model1<-plw(AffySpikeU95Subset,design=design,contrast=contrast,
            epsilon=0.01)

## Look at fitted vs observed density for log(s2)
varHistPlot(model1)

plw documentation built on April 28, 2020, 6:38 p.m.