margLik {ppRep}

Marginal likelihood of replication effect estimate

Description

This function computes the marginal likelihood of the replication effect estimate tr under the power prior model

f(\code{tr}|\code{to}, \code{so}, \code{sr}, \code{x}, \code{y}) = \int_0^1 \int_{-\infty}^{\infty} \mathrm{N}(\code{tr}; \theta, \code{sr}^2) \times \mathrm{N}(\theta; \mu, \phi) \times \mbox{Beta}(\alpha; \code{x}, \code{y}) ~\mbox{d}\theta~ \mbox{d}\alpha

with \phi = 1/(1/\code{v} + \alpha/\code{so}^2) and \mu = \phi\{(\alpha\times\code{to})/\code{so}^2 + \code{m}/\code{v}\} using numerical integration.

Usage

margLik(tr, to, sr, so, x = 1, y = 1, m = 0, v = Inf, ...)

Arguments

tr

Effect estimate of the replication study.

to

Effect estimate of the original study.

sr

Standard error of the replication effect estimate.

so

Standard error of the replication effect estimate.

x

Number of successes parameter of beta prior for \alpha. Defaults to 1.

y

Number of failures parameter of beta prior for \alpha. Defaults to 1.

m

Mean parameter of initial normal prior for \theta. Defaults to 0.

v

Variance parameter of initial normal prior for \theta. Defaults to Inf (uniform prior).

...

Additional arguments passed to stats::integrate.

Value

Marginal likelihood

Author(s)

Samuel Pawel