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]