getProbScalarRD: Calculate risks from arctanh RD and log OP

In mclements/brm: Binary Regression Model

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

mclements/brm documentation built on June 20, 2020, 9:22 p.m.