em_lfdr: EM algorithm to estimate local false discovery rate

In STAREG: An Empirical Bayes Approach for Replicability Analysis Across Two Studies

EM algorithm to estimate local false discovery rate

Description

Estimate the local false discovery rate across two studies and apply a step-up procedure to control the FDR of replicability null.

Usage

em_lfdr(pa_in, pb_in, pi0a_in, pi0b_in)

Arguments

pa_in	A numeric vector of p-values from study 1.
pb_in	A numeric vector of p-values from study 2.
pi0a_in	An initial estimate of the null probability in study 1.
pi0b_in	An initial estimate of the null probability in study 2.

Value

Lfdr	The estimated local false discovery rate for replicability null.
fdr	The adjusted values based on local false discovery rate for FDR control.
xi00	An estimate of the prior probability for joint state (0, 0).
xi01	An estimate of the prior probability for joint state (0, 1).
xi10	An estimate of the prior probability for joint state (1, 0).
xi11	An estimate of the prior probability for joint state (1, 1).
f1	A non-parametric estimate for the non-null probability density function in study 1.
f2	A non-parametric estimate for the non-null probability density function in study 2.

STAREG documentation built on June 8, 2025, 10:56 a.m.