# Estimation of rate parameter in exponential distribution
T <- rexp(n = 1000, rate = 3)
lambda <- maxlogL(x = T, dist = "dexp", start = 5,
link = list(over = "rate", fun = "NegInv_link"))
summary(lambda)
# Link function name
fun <- NegInv_link()$name
print(fun)
# Link function
g <- NegInv_link()$g
curve(g(x), from = 0.1, to = 1)
# Inverse link function
ginv <- NegInv_link()$g_inv
curve(ginv(x), from = 0.1, to = 1)
Run the code above in your browser using DataLab