histTailPP: Histogram plot for tail posterior probability
In BGmix: Bayesian models for differential gene expression

Description Usage Arguments Author(s) References See Also Examples

Plots a histogram of tail posterior probability with its density under the null hypothesis

1	histTailPP(tpp.res, bw=0.05, xlim=c(0,1),nc=10)

`tpp.res`	output of TailPP
`bw`	bandwidth for kernel estimate of the null density
`xlim`	limits on the x axis
`nc`	number of bins of the histogram

Natalia Bochkina

Bochkina N., Richardson S. (2007) Tail posterior probability for inference in pairwise and multiclass gene expression data. Biometrics.

TailPP, FDRplotTailPP,EstimatePi0

data(ybar, ss)
 nreps <- c(8,8)

## Note this is a very short MCMC run!
## For good analysis need proper burn-in period.
 outdir <- BGmix(ybar, ss, nreps, jstar=-1, nburn=0, niter=100, nthin=1)

 params <- ccParams(outdir)  
 res <-  ccTrace(outdir)
  
 tpp.res <- TailPP(res, nreps, params, plots  = FALSE)
 histTailPP(tpp.res, bw=0.04, xlim=c(0,1), nc=10)