Para utilizar o link <em>loglog</em>, encontrei esse script abaixo. Então, como faço para implementá-lo no <em>glm()</em> ?<br /><br />}, loglog = {<br /> linkfun <- function(mu) -log(-log(mu))<br /> linkinv <- function(eta) exp(-exp(-eta))<br /> mu.eta <- function(eta) exp(-exp(-eta)-eta)<br /> valideta <- function(eta) all(eta != 0)<br /> }, stop(sQuote(link), " link not recognised"))<br /> structure(list(linkfun = linkfun, linkinv = linkinv, mu.eta = mu.eta,<br /> valideta = valideta, name = link), class = "link-glm")<br />}<br /><br />glm(y~x1+x2+x3,family=binomial(link=make.link("loglog")),data=X)<br /> <br /><br />
<hr style="border-top: 1px solid #ccc;" />
Em 30/06/2012 23:26, <strong>Walmes Zeviani < walmeszeviani@gmail.com ></strong> escreveu:<br />
<div><span style="font-family: trebuchet ms,sans-serif;">É só passar uma string que identifica os links para o argumento family da glm(), algo como</span></div>
<div><span style="font-family: trebuchet ms,sans-serif;"><br /></span></div>
<div><span style="font-family: trebuchet ms,sans-serif;">glm(..., familu=binomial(link="função de ligação"))</span></div>
<div><span style="font-family: trebuchet ms,sans-serif;"><br /></span></div>
<div><span style="font-family: trebuchet ms,sans-serif;">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 "</span><span style="font-family: 'trebuchet ms',sans-serif;">[R] writing link function to glm</span><span style="font-family:'trebuchet ms',sans-serif">":</span></div>
<div><span style="font-family: trebuchet ms,sans-serif;"><br /></span></div>
<div><a href="http://www.google.com.br/search?sourceid=chrome&client=ubuntu&channel=cs&ie=UTF-8&q=%5BR%5D+writing+link+function+to+glm" target="_blank">http://www.google.com.br/search?sourceid=chrome&client=ubuntu&channel=cs&ie=UTF-8&q=%5BR%5D+writing+link+function+to+glm</a></div>
<div><span style="font-family: trebuchet ms,sans-serif;"><br /></span></div>
<div><span style="font-family: 'trebuchet ms',sans-serif;">À disposição.</span></div>
<div><span style="font-family: 'trebuchet ms',sans-serif;">Walmes.</span></div>
<span style="font-family: trebuchet ms,sans-serif;"><br /> </span><span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br /><span style="font-family:trebuchet ms,sans-serif">Walmes Marques Zeviani</span><br /> <span style="font-family:trebuchet ms,sans-serif">LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)</span><br /><span style="font-family:trebuchet ms,sans-serif">Departamento de Estatística - Universidade Federal do Paraná</span><br /> <span style="font-family:trebuchet ms,sans-serif">fone: (+55) 41 3361 3573</span><br /><span style="font-family:trebuchet ms,sans-serif">VoIP: (3361 3600) 1053 1173</span><br /> <span style="font-family:trebuchet ms,sans-serif">e-mail: <a href="http://mce_host/compose?to=walmes@ufpr.br" target="_blank">walmes@ufpr.br</a></span><br /><span style="font-family:trebuchet ms,sans-serif">twitter: @walmeszeviani</span><br /> <span style=
 "font-family:trebuchet ms,sans-serif">homepage: <a href="http://www.leg.ufpr.br/%7Ewalmes" target="_blank">http://www.leg.ufpr.br/~walmes</a></span><br /><span style="font-family:trebuchet ms,sans-serif">linux user number: 531218</span><br /> <span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br /> <br />