agjack.pi0: Averaged generalized jackknife estimate of pi0
Description
Averaged generalized jackknife estimate of pi0Usage
agjack.pi0(subtobj,mean.n=c('mean','harmean','geomean') ,
pointpair=FALSE, trunc=TRUE,tol=1e-5)
Arguments
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 sam trunc
logical, indicating if each gjack
estiamte is truncated to [0,1]. tol
A small tolerance number.
Value
- A numeric scalar of estiamted pi0.
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.Examples
Run this codedata(simulatedSubt)
agjack.pi0(simulatedSubt) ##
Run the code above in your browser using DataLab