varHistPlot {plw} | R Documentation |
Will produce a histogram of observed variance estimators (logged) together with the fitted density.
varHistPlot(model,main="Histogram variance estimators", histCol=8,densityCol=2,drawLegend=TRUE)
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. |
Magnus Astrand
plw, lmw
# ------------------------------------------ # 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)