V interesting point Greg.
But are you not just suggesting left and right truncation?
It strikes me that if the data are Poission then a mixture is likely to be
better - something akin to zero-deflated.
Neg binomail works for greater variance == mix a gamma and poisson, but I'm
unsure what to mix to deflate the variance - maybe bernoulli.
In any event, one possibility is that the data are dependnent (e.g. family
data) so to get a underdispersion you could draw from the poisson and then
redraw from the result based on some condition such as closeness (e.g. as a
proxy for family size).
Hope this helps!
Gerard
Greg Snow
<[EMAIL PROTECTED]
org> To
Sent by: René Holst <[EMAIL PROTECTED]>,
[EMAIL PROTECTED] "R-help@r-project.org"
project.org <R-help@r-project.org>
cc
04/12/2008 17:21 Subject
Re: [R] Simulating underdispersed
counts
Accurate probably depends on what exactly you want to accomplish with the
simulation, but one simple approach would be to generate more data than you
need from a regular Poisson distribution, then throw away the extreme
values on both sides.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
[EMAIL PROTECTED]
801.408.8111
> -----Original Message-----
> From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
> project.org] On Behalf Of René Holst
> Sent: Thursday, December 04, 2008 6:40 AM
> To: R-help@r-project.org
> Subject: [R] Simulating underdispersed counts
>
> Hello,
>
> Anyone who knows a fast and accurate algorithm for generating draws
> from an underdispersed Poisson distribution. Or even better, if there
> is a package containing such an implementation.
>
> Thanks
>
> Rene
>
> ______________________________________________
> 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.
______________________________________________
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.
**********************************************************************************
The information transmitted is intended only for the person or entity to which
it is addressed and may contain confidential and/or privileged material. Any
review, retransmission, dissemination or other use of, or taking of any action
in reliance upon, this information by persons or entities other than the
intended recipient is prohibited. If you received this in error, please contact
the sender and delete the material from any computer. It is the policy of the
Department of Justice, Equality and Law Reform and the Agencies and Offices
using its IT services to disallow the sending of offensive material.
Should you consider that the material contained in this message is offensive
you should contact the sender immediately and also mailminder[at]justice.ie.
Is le haghaidh an duine nó an eintitis ar a bhfuil sí dírithe, agus le haghaidh
an duine nó an eintitis sin amháin, a bheartaítear an fhaisnéis a tarchuireadh
agus féadfaidh sé go bhfuil ábhar faoi rún agus/nó faoi phribhléid inti.
Toirmisctear aon athbhreithniú, atarchur nó leathadh a dhéanamh ar an
bhfaisnéis seo, aon úsáid eile a bhaint aisti nó aon ghníomh a dhéanamh ar a
hiontaoibh, ag daoine nó ag eintitis seachas an faighteoir beartaithe. Má fuair
tú é seo trí dhearmad, téigh i dteagmháil leis an seoltóir, le do thoil, agus
scrios an t-ábhar as aon ríomhaire. Is é beartas na Roinne Dlí agus Cirt,
Comhionannais agus Athchóirithe Dlí, agus na nOifígí agus na nGníomhaireachtaí
a úsáideann seirbhísí TF na Roinne, seoladh ábhair cholúil a dhícheadú.
Más rud é go measann tú gur ábhar colúil atá san ábhar atá sa teachtaireacht
seo is ceart duit dul i dteagmháil leis an seoltóir láithreach agus le
mailminder[ag]justice.ie chomh maith.
***********************************************************************************
______________________________________________
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.
