[R] new dgamma rate argument
Martin Maechler
maechler at stat.math.ethz.ch
Fri Jan 11 14:23:16 CET 2002
>>>>> "Jim" == Jim Lindsey <james.lindsey at luc.ac.be> writes:
Jim> Can someone explain to me in what way the new
Jim> (dpqr)gamma parameter can be interpreted as a rate
Jim> (when shape != 1)? The only gamma rate that I am aware
Jim> of is the hazard rate given by dgamma/(1-pgamma), the
Jim> log of which is returned by my hgamma function (event
Jim> library). Jim
NEWS has
o [dpqr]gamma now has third argument `rate' for S-compatibility
(and for compatibility with exponentials). Calls which use
positional matching may need to be altered.
i.e. one point of view (close to mine) could be:
The authors of R (R&R) called that argument of [dpqr]gamma()
`scale' as it should sensibly be called.
OTOH, (at least one of) the original S authors used `rate'
(for 1/scale) in a loose analogy with the exponential and
weibull distribution quite some time before R was born. Now
that there is an increasing drive for S source compatibility
between the different S dialects -- whenever it's ``easy'' --
the compatible parametrization has been allowed as well.
Martin
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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