Usage
DearBeggMonotonePvalSelection(y, u, theta0, sigma0, lam = 2, M = 1000, maxiter = 1000, test.stat = function(x){return(min(x))})
Arguments
y
Normally distributed effect sizes.
u
Associated standard errors.
theta0
Initial estimate for $\theta$.
sigma0
Initial estimate for $\sigma$.
lam
Weight of the first entry of $w$ in the likelihood function. Should be the same as used to generate res.
M
Number of runs to compute $p$-value.
maxiter
Maximum number of iterations of differential evolution algorithm. Increase this number to get higher accuracy.
test.stat
A function that takes as argument a vector and returns a number. Defines the test statistic to be used on the estimated
selection function $w$.