pi0 (version 1.3-250)

histf1: Histogram estimator of p-value density evaluated at 1

Description

Histogram estimator of p-value density evaluated at 1. See references.

Usage

histf1(p,max.bins=20,bin.method=c("max","nmse","bootstrap",
      "Sturges","Scott","FD"),discrete=FALSE,seq.perm=FALSE,
      nboots=200,plotit=FALSE,perm.n,perm.h,...)

Arguments

p

Vector of p-values

max.bins

maximum number of bins

bin.method

binning method

discrete

Whether p-values are discrete

seq.perm

Whether p-values come from sequential permutation tests

nboots

bootstrap sample size

plotit

Whether to plot the histogram

perm.n

n for sequential permutation tests

perm.h

h for sequential permutation tests

...

Other arguments passed to hist

Value

• Numeric scalar

References

Nettleton, Hwang, Caldo, Wise. 2006. Estimating the number of true null hypotheses from a histogram of $p$ values. Journal of Agricultural, Biological, and Environmental Statistics. 11. 337-356. Bancroft and Nettleton. 2009. Estimation of False Discovery Rate Using Permutation P-values with Different Discrete Null Distributions. Iowa State University Department of Statistics Preprint Series, #2009-05. Bancroft and Nettleton. 2009. Computationally Efficient Estimation of False Discovery Rate Using Sequential Permutation P-values. Iowa State University Department of Statistics Preprint Series, #2009-04.

See Also

lastbin, qvalue

Examples

set.seed(9992722)
histf1(runif(5e5)^1.5) ##  [1] 0.6762