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

Martin Maechler maechler at stat.math.ethz.ch
Tue Feb 2 14:48:27 CET 1999


>>>>> "Yudi" == Yudi Pawitan <yudi at hermes.ucd.ie> 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 stat.math.ethz.ch>	http://stat.ethz.ch/~maechler/
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 http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list