Plot predictive density of data.

Share:

Description

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

Usage

1
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

1
2
3
4
5
6
7
## 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)