agjack: Averaged generalized jackknife estimate of pi0

In pi0: Estimating the proportion of true null hypotheses for FDR

Description

Averaged generalized jackknife estimate of pi0

Usage

agjack.pi0(subtobj,mean.n=c('mean','harmean','geomean') ,
    pointpair=FALSE, trunc=TRUE,tol=1e-5)

Arguments

subtobj	A subt object
mean.n	A character: name of the function to compute the average sample sizes
pointpair	logical: if TRUE, then gjack is called for each pair of rows in subtobj; otherwise, gjack is called for each pair of unique average sample sizes.
trunc	logical, indicating if each gjack estiamte is truncated to [0,1].
tol	A small tolerance number.

Details

When pointpair is FALSE, the rows in subtobj are first grouped by combination of sample sizes and the estimates are averaged for each group. Then gjack is called for each pair of groups.

Value

A numeric scalar of estiamted pi0.

Author(s)

Long Qu

See Also

gjack

Examples

    data(simulatedSubt)
    agjack.pi0(simulatedSubt)  ##

pi0 documentation built on May 2, 2019, 4:47 p.m.