varHistPlot {plw}R Documentation

Variance histogram and density

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)


[Package plw version 1.2.0 Index]