Bookmark and Share

Notice: On March 31, it was announced that Statalist is moving from an email list to a forum. The old list will shut down on April 23, and its replacement, is already up and running.

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: st: Underdispersion in count data from a survery

From   Steve Samuels <>
Subject   Re: st: Underdispersion in count data from a survery
Date   Thu, 7 Oct 2010 19:48:22 -0400


Your analysis  of under-dispersion could also be incorrect. Your
estimated mean and variance have to be  weighted by the survey
weights.  Your predicted values have to come from the weighted Poisson
regression,  and your  test for under-dispersion isn't valid unless
you used -svy: reg- to do it.

After -svy: mean-, the following will give you the estimate of the
sample variance
scalar var_count = e(N)*el(e(V_srs),1,1)
di var_count
On Thu, Oct 7, 2010 at 7:15 PM, Steve Samuels <> wrote:
> --
> You can use the -glm- command in Stata 9 with survey data (link(log)
> family(poisson)) and follow it by -suest- to get valid standard
> errors. See
>  You could also run -svy reg- on the log of (count +1).
> Under-dispersion is a concern for standard errors when your
> inferences are likelihood based. The survey-based commands will base
> standard errors on the sample design.
> The Poisson or regression  models for predicting or estimating  means
> might be useful even if the data are not Poisson distributed. You
> still have to check model diagnostics. See:
>  for some
> ideas.
> Steve
> On Thu, Oct 7, 2010 at 5:46 PM, Laurie Molina <> wrote:
>> Hello all,
>> I have count data from a survey which comes from a complex survey
>> design (stratification, clusters, two stages, and probabilistic
>> design).
>> The mean is 1.89 and the variance is 1.14.
>> Fallowing "Essentials of count data regression" by Cameron and Trivedi
>> (1999)To test for underdispersion i ran the following auxiliary
>> regression:
>> ((yi - muhati)^2 - yi)/muhati= alpha*muhati + ui
>> where yi is the i observation of  original dependent count variable,
>> muhati is the i observation of the fitted values resulting from the
>> poisson coefficients estimated, ui is the error term, and the
>> objective is to test alpha= 0 Vs alpha<0.
>> Whith a p-value of 1, the hypothesis of alpha < 0 cannot be rejected.
>> So my question is... What now? My data is underdispersed so the
>> standard errors estimated when i run the Poisson regression are not
>> correct (even when the estimates of my coefficents are consistent).
>> Is there any solution for this problem in Stata 9, that can be used
>> with survey data?
>> Thank you all very much in advance.
>> Laurie
>> *
>> *   For searches and help try:
>> *
>> *
>> *

*   For searches and help try:

© Copyright 1996–2016 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   Site index