plotMixDensity {BGmix}R Documentation

Plot predictive density of data.

Description

Plot predictive density of data superimposed on histograms of observed data. Separate plots for ybar and sums of squares.

Usage

plotMixDensity(res, predres, ybar, ss)

Arguments

res list object output from 'ccParams'
predres list object output from 'ccPred' (need q.trace=T in 'ccPred')
ybar ybar data (see BGmix help for details)
ss ss data (see BGmix help for details)

Details

Note that this function is written for the unpaired differential expression design.

Author(s)

Alex Lewin

Examples

## Note this is a very short MCMC run!
## For good analysis need proper burn-in period.
data(ybar,ss)
outdir <- BGmix(ybar,ss,c(8,8),niter=100,nburn=0,nthin=1,trace.pred=1)
pred <- ccPred(outdir,q.trace=TRUE)
params <- ccParams(outdir)
plotMixDensity(params,pred,ybar,ss)

[Package BGmix version 1.2.0 Index]