cloglog link [was "Re: [R] glm null deviance"]

Martin Maechler maechler at
Tue Feb 2 14:48:27 CET 1999

>>>>> "Yudi" == Yudi Pawitan <yudi at> writes:

      ..... (lots of very helpful notes) ....  [ omitted ]

    Yudi> Question to the devel:
    Yudi> shouldn't the cloglog = log(-log(p))?
    Yudi> rather than   log(-log(1-p)) as implemented?

The current implementation is consistent with  the other 
link functions for the binomial: 
     It is monotonely INcreasing

BTW: Try  
	   plot(function(p) log(-log(1-p)))

     if your R version is >= 0.63.2

Martin Maechler <maechler at>
Seminar fuer Statistik, ETH-Zentrum SOL G1;	Sonneggstr.33
ETH (Federal Inst. Technology)	8092 Zurich	SWITZERLAND
phone: x-41-1-632-3408		fax: ...-1086			<><
r-help mailing list -- Read
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at

More information about the R-help mailing list