histTailPP

0th

Percentile

Histogram plot for tail posterior probability

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

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

Aliases
  • histTailPP
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)


Documentation reproduced from package BGmix, version 1.32.0, License: GPL-2

Community examples

Looks like there are no examples yet.