plotMixDensity: Plot predictive density of data.
In BGmix: Bayesian models for differential gene expression

Description Usage Arguments Details Author(s) Examples

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

1	plotMixDensity(res, predres, ybar, ss)

`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)

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

Alex Lewin

## 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)