I would like to adopt the hyperbolic arcsin (asinh) for some of my modeling. I have written the following function:
link.asinh <- function(mu)
{
linkfun<-function(mu) asinh(mu)
linkinv<-function(eta) sinh(eta)
mu.eta<-function(eta) sinh(eta)
valideta<-function(eta) TRUE
link <- paste0("lasinh")
structure(list(linkfun=linkfun, linkinv=linkinv, mu.eta=mu.eta,
valideta=valideta, name=link), class="link-glm")
}
But I'm at a loss from that point, how do I use this link in a Gaussian family glm, for example? I have read the ?family and ?make.link, no no avail.