getProbScalarRD {brm}

Calculate risks from arctanh RD and log OP

Description

Calculate risks from arctanh RD and log OP

Usage

getProbScalarRD(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 vector (P(y=1|x=0),P(y=1|x=1))

Examples

getProbScalarRD(0,0)

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