histTailPP {BGmix}R Documentation

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

Author(s)

Natalia Bochkina

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)


[Package BGmix version 1.2.0 Index]