Helper function: squeezed logit

Usage

flogit(p, sqz = 0.000001)

Arguments

p: a vector of values between 0 and 1 inclusive
sqz: the amount by which to 'squeeze', default is .000001

Value

   a vector of values between -Inf and +Inf

Examples


p <- runif(n = 1000)
summary(p)
#>     Min.  1st Qu.   Median     Mean  3rd Qu.     Max. 
#> 0.001446 0.255711 0.517020 0.511831 0.759035 0.997649 

sqz <- 1 / (10**6)
x <- flogit(p, sqz = sqz)
summary(x)
#>     Min.  1st Qu.   Median     Mean  3rd Qu.     Max. 
#> -6.53706 -1.06839  0.06811  0.07389  1.14739  6.05032 

all(abs(p - fexpit(x, sqz = sqz)) < sqz)
#> [1] TRUE
all(abs(p - fexpit(flogit(p, sqz = sqz), sqz = sqz)) < sqz)
#> [1] TRUE