similar to: Optim()

Hi all, I have made an overdispersion test for a data set and get the following result Overdispersion test Obs.Var/Theor.Var Statistic p-value poisson data 16.24267 47444.85 0 after deleting the outliers from the data set I get the following result Overdispersion test Obs.Var/Theor.Var Statistic p-value poisson data 16.27106 0 1 The

Dear R users, I am trying to fully understand the difference between estimating overdispersion with glm.nb() from MASS compared to glm(..., family = quasipoisson). It seems that (i) the coefficient estimates are different and also (ii) the summary() method for glm.nb suggests that overdispersion is taken to be one: "Dispersion parameter for Negative Binomial(0.9695) family taken to be

Hello, I'm a grad student in the Intelligent Transportation Systems lab at Portland State Univ. in Portland, OR, USA. I'm trying to learn the basics of R to run a negative binomial in the near future, and so I ran a test regression on roadway crash data obtained from "Statistical and Econometric Methods for Transportation Data Analysis" by Washington et al (p. 250). I ran the

Hi, there are a lot of messages dealing with overdispersion, but I couldn't find anything about how to test for overdispersion. I applied a GAM with binomial distribution on my presence/absence data, and would like to check for overdispersion. Does anyone know the command? Many thanks, Anna -- View this message in context:

Dear list members, while I appreciate the possibility to deal with overdispersion for count data either by specifying the family argument to be quasipoisson() or negative.binomial(), it estimates just one overdispersion parameter for the entire data set. In my applications I often would like the estimate for overdispersion to depend on the covariates in the same manner as the mean. For example,

dear list! i am running an anlysis on proportion data using binomial (quasibinomial family) error structure. My data comprises of two continuous vars, body size and range size, as well as of feeding guild, nest placement, nest type and foragig strata as factors. I hope to model with these variables the preference of primary forests (#successes) by certain bird species. My code therefore looks

I have a few questions relating to overdispersion in a sex ratio data set that I am working with (note that I already have an analysis with GLMMs for fixed effects, this is just to estimate dispersion). The response variable is binomial because nestlings can only be male or female. I have samples of 1-5 nestlings from each nest (individuals within a nest are not independent, so the response

Hello! I am trying to model data on species abundance (count data) with a poisson error distribution. I have a fixed and a random variables and thus needs a mixed model. I strongly doubt that my model is overdispersed but I don't know how to get the overdispersion parameter in a mixed model. Maybe someone can help me on this point. Secondly, it seems that quasi models cannot be implemented

I fitted a GAM model with Poisson distribution to a data with about 200 observations. I noticed that the plot of the residuals versus fitted values show a trend. Residuals tend to be lower for higher fitted values. Because, I'm dealing with count data, I'm thinking that this might be due to overdispersion. Is there a way to account for overdispersion in any of the packages MGCV or GAM?

Hi, I am performing glm with binomial family and my data show slight overdispersion (HF<1.5). Nevertheless, in order to take into account for this heterogeneity though weak, I use F-test rather than Chi-square (Krackow & Tkadlec, 2001). But surprisingly, outputs of this two tests are exactly similar. What is the reason and how can I scale the output by overdispersion ?? Thank you,

I am looking for the correct commands to do the following things: 1. I have a binomial logistic regression model and i want to test for overdispersion. 2. If I do indeed have overdispersion i need to then run a quasi-binomial model, but I'm not sure of the command. 3. I can get the residuals of the model, but i need to then apply a shapiro wilk test to test them. Does anyone know the command

Dear All, I'm trying to use waldtest to test poolability (parameter stability) between two logistic regressions. Because I need to use robust standard errors (using sandwich), I cannot use anova. anova has no problems running the test, but waldtest does, indipendently of specifying vcov or not. waldtest does not appear to see that my models are nested. H0 in my case is the the vector of

How can I eliminate the overdispersion for binary data apart the use of the quasibinomial? help me Eva Iannario ------------------------------------------------------ Passa a Infostrada. ADSL e Telefono senza limiti e senza canone Telecom http://click.libero.it/infostrada11gen07

Hi all, I was wondering if somebody could give me advice regarding the dispersion parameter in GLMM's. I'm a beginner in R and basically in GLMM's. I've ran a GLMM with a poisson family and got really nice results that conform with theory, as well with results that I've obtained previously with other analysis and that others have obtained in similar studies. But the

Sorry about re-posting this, it never went out to the mailing list when I posted this to r-help forum on Nabble and was pending for a few days, now that I am subscribe to the mailing list I hope that this goes out: I've been a viewer of this forum for a while and it has helped out a lot, but this is my first time posting something. I am running glm models for richness and abundances. For

It seems that MASS suggest to judge on the basis of sum(residuals(mode,type="pearson"))/df.residual(mode). My question: Is there any rule of thumb of the cutpoiont value? The paper "On the Use of Corrections for Overdispersion" suggests overdispersion exists if the deviance is at least twice the number of degrees of freedom. Are there any further hints? Thanks. -- Ronggui

Dear All I recently proposed a simple modification to Wedderburn's 1974 estimate of overdispersion for count and binomial data, which is used in glm for the quasipoisson and quasibinomial families (see the reference below). Although my motivation for the modification arose from considering sparse data, it will be almost identical to Wedderburn's estimate when the data are not sparse.

Dear All I am trying to do a repeated measures analysis using lmer and have a number of issues. I have non-orthogonal, unbalanced data. Count data was obtained over 10 months for three treatments, which were arranged into 6 blocks. Treatment is not nested in Block but crossed, as I originally designed an orthogonal, balanced experiment but subsequently lost a treatment from 2 blocks. My

Hi, I have to analyse the number of provisioning trips to nestlings according to a number of biological and environmental factors. I was thinking of building a mixed-effects model with species and nestid as random effects, using a Poisson distribution, but the data are overdispersed (variance/mean = 5). I then thought of using a mixed-effects model with negative binomial distribution, but I have

I've been poking around with GLMs (on which I am *not* an expert) on behalf of a student, particularly binomial (standard logit link) nested models with overdispersion. I have one possible bug to report (but I'm not confident enough to be *sure* it's a bug); one comment on the general inconsistency that seems to afflict the various functions for dealing with overdispersion in GLMs