So there is no maximum value, just a minimum (greater than 0)? Correct?  In 
that case the sample option will not work (unless you choose some really high 
value and say that you won't go above that), but the rejection sampling would 
still work.  How efficiently it works will depend on how much probability a 
regular poisson would put into the truncated region.

Another possibility is to find the probability of being in the truncated 
region, then generate a uniform between that value and 1, then feed that 
uniform into the qpois function.

From: cassie jones [mailto:cassiejone...@gmail.com]
Sent: Monday, May 16, 2011 7:46 PM
To: Greg Snow
Cc: r-help@r-project.org
Subject: Re: [R] simulation from truncated poisson

It is truncated from left.
On Mon, May 16, 2011 at 6:33 PM, Greg Snow 
<greg.s...@imail.org<mailto:greg.s...@imail.org>> wrote:
Which direction is it truncated?  (only values less than a allowed or only 
greater?).

One simple approach is rejection sampling, just generate from a regular poisson 
distribution, then throw away any values in the truncated region.  Another 
approach if the legal values are those from 0 to a, so that there is a finite 
number of possibilities, then you can use the sample function with replace=TRUE 
and using probabilities from the poisson in the legal range.

-----Original Message-----
From: r-help-boun...@r-project.org<mailto:r-help-boun...@r-project.org> 
[mailto:r-help-boun...@r-project.org<mailto:r-help-boun...@r-project.org>] On 
Behalf Of cassie jones
Sent: Monday, May 16, 2011 5:28 PM
To: r-help@r-project.org<mailto:r-help@r-project.org>
Subject: [R] simulation from truncated poisson

Dear all,

I need to simulate values from a Poisson distribution which is truncated at
certain value 'a'. Can anyone tell me if there is in-built package in R
which can simulate from a truncated Poisson? If not, what should be the
steps to write a function which would do that?

Thanks in advance.


Cassie
       [[alternative HTML version deleted]]

______________________________________________
R-help@r-project.org<mailto: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.


        [[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