BGmix (version 1.32.0)

histTailPP: Histogram plot for tail posterior probability

Description

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

Usage

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

Arguments

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

References

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

See Also

TailPP, FDRplotTailPP,EstimatePi0

Examples

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)