plotMixDensity

0th

Percentile

Plot predictive density of data.

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

Keywords
hplot
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.

Aliases
  • plotMixDensity
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)
Documentation reproduced from package BGmix, version 1.32.0, License: GPL-2

Community examples

Looks like there are no examples yet.