Histogram plot for tail posterior probability

Share:

Description

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

Usage

1
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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
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)