Para utilizar o link loglog, encontrei esse script abaixo. Então, como faço para implementá-lo no glm() ?
}, loglog = {
linkfun <- function(mu) -log(-log(mu))
linkinv <- function(eta) exp(-exp(-eta))
mu.eta <- function(eta) exp(-exp(-eta)-eta)
valideta <- function(eta) all(eta != 0)
}, stop(sQuote(link), " link not recognised"))
structure(list(linkfun = linkfun, linkinv = linkinv, mu.eta = mu.eta,
valideta = valideta, name = link), class = "link-glm")
}
glm(y~x1+x2+x3,family=binomial(link=make.link("loglog")),data=X)
Em 30/06/2012 23:26, Walmes Zeviani < walmeszeviani@gmail.com > escreveu:
É só passar uma string que identifica os links para o argumento family da glm(), algo como
glm(..., familu=binomial(link="função de ligação"))
você pode usar uma função de ligação própria, por exemplo, se eu quiser usar o modelo gompertz, basta escrever a função na estrutura correta que a glm() é capaz de usar. Veja os resultados da busca com "[R] writing link function to glm":
À disposição.
Walmes.
==========================================================================
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
Departamento de Estatística - Universidade Federal do Paraná
fone: (+55) 41 3361 3573
VoIP: (3361 3600) 1053 1173
e-mail: walmes@ufpr.br
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================