Thank you very much Jay.

2010/3/19 G. Jay Kerns <gke...@ysu.edu>

> Dear Randall,
>
> On Fri, Mar 19, 2010 at 10:24 AM, Randall Wrong <randall.wr...@gmail.com>
> wrote:
> > Dear R users,
> >
> > ?rgamma gives me :
> >
> >               rgamma(n, shape, rate = 1, scale = 1/rate)
> >
> >               rate: an alternative way to specify the scale.
> >
> >               The Gamma distribution with parameters ‘shape’ = a and
> >               ‘scale’ = s has density
> >               f(x)= 1/(s^a Gamma(a)) x^(a-1) e^-(x/s)
> >
> > Should I understand that scale=1/rate ? Is it written somewhere ?
>
> You are kidding, right?  It is written 8 lines above your question, by
> my count.  :-)
>
> Perhaps you meant rate = 1/scale.
>
> >
> > Then rgamma(n, shape=a, scale = s) should be equivalent to rgamma(n,
> > shape=a, rate =1/s).
>
> Yep:
> dgamma(2, shape = 3, scale = 4)
> dgamma(2, shape = 3, rate = 1/4)
>
> >
> > I don't find this very clear.
> >
> > Thanks for your help.
> >
> > Randall
> >
>
> The point is that some books (and software) parameterize by the
> 'scale', and a whole other bunch parameterize by the 'rate'.  The
> reader (and user) always needs to be careful that the version used is
> the one expected. And the help file says that S doesn't have a 'scale'
> parameter at all.
>
> Just be careful, and you should be fine.  And IMHO, given that the PDF
> of the density is shown it is reasonably clear as-is.
>
> Best,
> Jay
>
>
>
>
>
> ***************************************************
> G. Jay Kerns, Ph.D.
> Associate Professor
> Department of Mathematics & Statistics
> Youngstown State University
> Youngstown, OH 44555-0002 USA
> Office: 1035 Cushwa Hall
> Phone: (330) 941-3310 Office (voice mail)
> -3302 Department
> -3170 FAX
> VoIP: gjke...@ekiga.net
> E-mail: gke...@ysu.edu
> http://people.ysu.edu/~gkerns/
>

        [[alternative HTML version deleted]]

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to