similar to: Overdispersed GLM

I would like to try a simple parametric bootstrap, but unfortunately (stupidly?) my models are "overdispersed" gams & glms. I'm hoping for a function that generates overdispersed poisson or negative binomial data with a given mean, scale (& shape parameter). The loose definition I'm using is overdispersed poisson produces integer values with variance=const*mean &

I am wondering if stepAIC in the MASS library may be used for model selection in an overdispersed Poisson situation. What I thought of doing was to get an estimate of the overdispersion parameter phi from fitting a model with all or most of the available predictors (we have a large number of observations so this should not be problematical) and then use stepAIC with scale = phi. Should this

Dear all, I am fitting an nonlinear glm() using optim() by first minimising glm(resp~ var1 + var2, family=binomial, data=data)$deviance where var1= exp(-a1*dist1), and var2= exp(-a2*dist2), where a1 and a2 are parameters and dist1 and dist2 are independent variables. Next, I calculate the value of var1 (and var2) by plugging in the value of al1 (and al2) that minimises deviance, and fit

Dear R users I have been looking for functions that can deal with overdispersed poisson models. According to actuarial literature (England & Verall, Stochastic Claims Reserving in General Insurance , Institute of Actiuaries 2002) this can be handled through the use of quasi likelihoods instead of normal likelihoods. However, we see them frequently in this type of data, and we would like to

Information on package 'aod' Description: Package: aod Version: 1.1-2 Date: 2005-06-08 Title: Analysis of Overdispersed Data Author: Matthieu Lesnoff <matthieu.lesnoff at cirad.fr> and Renaud Lancelot <renaud.lancelot at cirad.fr> Maintainer: Renaud Lancelot <renaud.lancelot at cirad.fr> Depends: R (>=

Information on package 'aod' Description: Package: aod Version: 1.1-2 Date: 2005-06-08 Title: Analysis of Overdispersed Data Author: Matthieu Lesnoff <matthieu.lesnoff at cirad.fr> and Renaud Lancelot <renaud.lancelot at cirad.fr> Maintainer: Renaud Lancelot <renaud.lancelot at cirad.fr> Depends: R (>=

Hi all, First, a quick thank you for R; it's amazing. I am trying to fit models for a count dataset following the overdispersed logisitic regression approach outlined in Baggerly et al. (BMC Bioinformatics, 5:144; Annotated R code is given at the end of the paper) but R is returning an error with the data below. Any help in understanding or overcoming this obstacle is appreciated.

Dear R-users, Earlier this year I posted a message to this list regarding negative binomial mixed models in R. It was suggested that the program I had written should be turned into an R-package. This has now been done, in collaboration with David Fournier and Anders Nielsen. The R-package glmmADMB provides the following GLMM framework: - Negative binomial or Poisson responses. - Zero-inflation

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

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

Dear R users I have been looking for functions that can deal with overdispersed poisson models. Some (one) of the observations are negative. According to actuarial literature (England & Verall, Stochastic Claims Reserving in General Insurance , Institute of Actiuaries 2002) this can be handled through the use of quasi likelihoods instead of normal likelihoods. The presence of negatives is not

Dear all, I am using gamlss (Package gamlss version 4.0-0, R version 2.10.1, Windows XP Service Pack 3 on a HP EliteBook) to relate bird counts to habit variables. However, most models fail because “the global deviance is increasing” and I am not sure what causes this behaviour. The dataset consists of counts of birds (duck) and 5 habit variables measured in the field (n= 182). The dependent

I am creating habitat selection models for caribou and other species with data collected from GPS collars. In my current situation the radio-collars recorded the locations of 30 caribou every 6 hours. I am then comparing resources used at caribou locations to random locations using logistic regression (standard habitat analysis). The data is therefore highly autocorrelated and this causes Type

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

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

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

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:

Hi there, I am using the geepack package and wish to test if my data (family=poisson) is overdispersed so that I may adjust the "scale.value" if necessary. Is there a specific function or method I should be using? Thanks for the advice, Bruce _____________________________ Bruce Catton MSc Candidate Forest Sciences University of British Columbia Vancouver, British

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 R-helpers: I have a question related to fitting overdispersed count data using lmer. Basically, I simulate an overdispsed data set by adding an observation-level normal random shock into exp(....+rnorm()). Then I fit a lmer quasipoisson model. The estimation results are very off (see model output of fit.lmer.over.quasi below). Can someone kindly explain to me what went wrong? Many thanks in