x <- seq(0.8, 1.5, by = 0.1)
logloglink(x) # Has NAs
logloglink(x, bvalue = 1.0 + .Machine$double.eps) # Has no NAs
x <- seq(1.01, 10, len = 100)
logloglink(x)
max(abs(logloglink(logloglink(x), inverse = TRUE) - x)) # 0?
Run the code above in your browser using DataLab