similar to: Estimating overdispersion when using glm for count and binomial data

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.

I am attempting to run a glm with a binomial model to analyze proportion data. I have been following Crawley's book closely and am wondering if there is an accepted standard for how much is too much overdispersion? (e.g. change in AIC has an accepted standard of 2). In the example, he fits several models, binomial and quasibinomial and then accepts the quasibinomial. The output for residual

http://biomet.oxfordjournals.org/cgi/reprint/86/3/677 biometrika1999 http://biomet.oxfordjournals.org/cgi/reprint/94/4/1006 biometrika2000 Hi All: I just want to try some luck. I am currenly working on my project,one part of my project is to reanalysis the kenward cattle data by using the method in Mohsen's paper,but I found I really can get the same or close output as he did,so,any

If you mean this paper by Fahrmeir: http://biomet.oxfordjournals.org/cgi/content/abstract/81/2/317 I would recommend BayesX: http://www.stat.uni-muenchen.de/~bayesx/. BayesX interfaces with R and estimates discrete (and continuous) time survival data with penalized regression methods. If you are looking for a bona fide Bayesian survival analysis method and do not wish to spend a lot of time

R Users, The d/p/q/r functions for the bridge distribution are now available in bridgedist. When a random intercept follows the bridge distribution, as detailed in Wang and Louis (2003) <doi:10.1093/biomet/90.4.765 <http://dx.doi.org/10.1093/biomet/90.4.765>>, a marginalized random-intercept logistic regression will still be a logistic regression with marginal coefficients that are

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

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 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, The share of concurring votes (i.e. yes-yes and no-no) in total votes between a pair of voters is a function of their ideological distance (index continuous on [1,2]). I show by other means that the votes typically are highly positively correlated (with an average c=0.6). This is because voters sit together and discuss the issue before taking a vote, but also because they share common

Dear R users, I work with a descrete variable (sclae 0 - 27) which is highly skwed to the right (many zeros and small numbers). I measure this variable on a control and intervention cohort 5 times a year. When I analyze analyze this varoable at each time point separately and use GLM with family quasi-Poisson (descrete outcome and two binary variables, gender and cohort, are predictors), I observe

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

Hello everybody, I have count data and with these data, I would like to build a mixed model by using the function glmer(). In a first time, I calculated the c-hat of a simple model with glm() to verify overdispersion and I found a c-hat = 18. I also verified overdispersion in the mixed model by checking the residuals of random effects via the function glmmPQL and I found a c-hat = 15. Thus,

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

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

Dear friends, As we all know, the usual model selection criteria(e.g.deviance,AIC...) in GLMs isn't very good for selecting the best model when overdispersion exist, so we need to adjust the corresponding statistic,see(Fitzmaurice,G.M. (1997) Model selection with overdispersed

Hello, I am analysing counts data with a mixed model using lmer procedure. I therefore use the quasipoisson distribution but I'm not sure if this is sufficient to account for overdispersion. Actually the results are not very different to what I get when specifying a poisson distribution although my data are clearly overdispersed. this my model: >model <- lmer(NB ~ T + volume +

Dear all, I am using function glmer from package lme4 to fit a generalized linear mixed effect model. My model is as follows: model1 <- glmer(fruitset ~ Dist*wire + (1|Site), data, binomial) summary(model1) Generalized linear mixed model fit by the Laplace approximation Formula: fruitset ~ Dist * wire + (1 | Site) Data: data AIC BIC logLik deviance 68.23 70.65 -29.11 58.23 Random

Dear R users, Overdispersion is often a problem in binomial data. I attempt to model a binary response (sex-ratio) with three categorical explanatory variables, using GLM, which could assume the form: y<-cbind(sexf, sample-sexf) model<-glm(y ~ age+month+year, binomial) summary(model) Output: (Dispersion parameter for binomial family taken to be 1) Null deviance: 8956.7 on 582

Hey Buddy, Hope you have been doing well since last contact. If you have the answer to the following question, please let me know. If you have chance to travel up north. let me know. best, -Sean ---------- Forwarded message ---------- From: Sean Zhang <seanecon@gmail.com> Date: Sat, Apr 11, 2009 at 12:12 PM Subject: question related to fitting overdispersion count data using lmer

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,