I was looking at different link functions for binomial glms recently for
the same reason as you (more zeros than ones). I did a bit of reading up
on the various link functions and IIRC, you can use the cloglog link on
your data, just turn your 0's into 1's and vice versa. This was stated
in the
can you please say specifically how to modify the make.link function
--
View this message in context:
http://r.789695.n4.nabble.com/glm-binomial-loglog-NOT-cloglog-link-tp849370p2135909.html
Sent from the R help mailing list archive at Nabble.com.
__
Dear Bill,
Perhaps the cloglog function in the VGAM package might be useful for
you.
HTH,
Jorge
On Fri, Jan 23, 2009 at 11:32 AM, William Simpson
william.a.simp...@gmail.com wrote:
I would like to do an R glm() with
family = binomial(link=loglog)
Right now, the cloglog link exists, which
I wish to write a new link function for a GLM. R's glm routine does
not supply the loglog link. I modified the make.link function adding
the code:
}, loglog = {
linkfun - function(mu) -log(-log(mu))
linkinv - function(eta) exp(-exp(-eta))
mu.eta - function(eta)
I would like to do an R glm() with
family = binomial(link=loglog)
Right now, the cloglog link exists, which is nice when the data have a
heavy tail to the left. I have the opposite case and the loglog link
is what I need. Can someone suggest how to add the loglog link onto
glm()? It would be
William Simpson william.a.simpson at gmail.com writes:
I would like to do an R glm() with
family = binomial(link=loglog)
Right now, the cloglog link exists, which is nice when the data have a
heavy tail to the left. I have the opposite case and the loglog link
is what I need. Can someone
6 matches
Mail list logo