brm (version 1.1.1)

getProbScalarRD: 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

Value

a vector \((P(y=1|x=0),P(y=1|x=1))\)

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\).

Examples

# NOT RUN {
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

# }