[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