plotMixDensity: Plot predictive density of data.

Description Usage Arguments Details Author(s) Examples

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)


Search within the BGmix package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.