On 23-Feb-06 Robin Hankin wrote: > 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) <[EMAIL PROTECTED]> Fax-to-email: +44 (0)870 094 0861 Date: 23-Feb-06 Time: 09:18:50 ------------------------------ XFMail ------------------------------ ______________________________________________ [email protected] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
