R/loglog.R

In BinaryEPPM: Mean and Variance Modeling of Binary Data

loglog <-
function()
{ 
   linkfun <- function(mu) { -log(-log(mu)) }
   linkinv <- function(eta) { pmax(pmin(exp(-exp(-eta)), 
                              1-.Machine$double.eps), .Machine$double.eps) }
   mu.eta <- function(eta) { eta <- pmin(eta, 700)
                             pmax(exp(exp(-eta)), .Machine$double.eps) }
   valideta <- function(eta) TRUE
   link <- paste0("loglog(",")")
   structure (list(linkfun = linkfun, linkinv = linkinv,
                   mu.eta = mu.eta, valideta = valideta, name = link),
              class = "link-glm")
}