discrete probability distribution

Hi

I want to sample from a discrete random variable X,  defined on
the  non-negative integers, with

prob(X=0) = (1-p)
prob(X=1) = (1-p)*p
...
prob(X=i)=(1-p)*p^i
...



Before reinventing the wheel,
is there a ready-made R function to give me such a thing?

Use rgeom(n,(1-p))

Note (from ?rgeom) that these functions use the definition

  P(n) = p*(1-p)^n

i.e. the complementary value of p to the one in your description.

Best wishes,
Ted.

--------------------------------------------------------------------
E-Mail: (Ted Harding) <Ted.Harding at nessie.mcc.ac.uk>
Fax-to-email: +44 (0)870 094 0861
Date: 23-Feb-06                                       Time: 09:18:50
------------------------------ XFMail ------------------------------