plotMixDensity {BGmix} | R Documentation |
Plot predictive density of data superimposed on histograms of observed data. Separate plots for ybar and sums of squares.
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)