getProbRD {brm}

Calculate risks from arctanh RD and log OP (vectorised)

Description

Calculate risks from arctanh RD and log OP (vectorised)

Usage

getProbRD(atanhrd, logop)

Arguments

atanhrd

arctanh of risk difference

logop

log of odds product

Details

The log OP is defined as log OP = log[(P(y=1|x=0)/P(y=0|x=0))*(P(y=1|x=1)/P(y=0|x=1))]. The inverse hyperbolic tangent function arctanh is defined as arctanh(z) = [log(1+z) - log(1-z)] / 2.

Value

a matrix (P(y=1|x=0),P(y=1|x=1)) with two columns

Examples

getProbRD(0,0)

set.seed(0)
logrr = rnorm(10,0,1)
logop = rnorm(10,0,1)
probs = getProbRD(logrr, logop)
colnames(probs) = c("P(y=1|x=0)","P(y=1|x=1)")
probs