gamma0_evidence: Calculate marginal probability of observed case-control status y under model gamma = 0
Description
Marginal probability calculated exactly by integration.
Usage
gamma0_evidence(y, tau0_shape = c(1, 1))
Arguments
y
Logical vector of case (TRUE) control (FALSE) status.
tau0_shape
Beta shape hyper-priors for prior on rate of case labels
Value
Log marginal likelihood.