histTailPP: Histogram plot for tail posterior probability

Description Usage Arguments Author(s) References See Also Examples

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)


Search within the BGmix package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.