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