similar to: Zero inflated Models- pscl package

I have some questions regarding Zero Inflation Poisson models. I am using count data to analyze abundance trends of salamanders. However, I have surveys which differ in the amount of effort (i.e. the number of people searching and amount of time - I am using a museum database so not all surveys were conducted by me). Therefore I need to account for the effort. If change the count (response

This request is related to the following post from last year: https://stat.ethz.ch/pipermail/r-help/2011-June/279752.html After reading the thread, the idea is still not clear. I have fitted a model using HURDLE from the PSCL package. I am trying to get marginal effects / slopes by multiplying the coefficients by the mean of the marginal effects (I think this is right). To my understanding, this

Dear list, I'm currently trying to develop a model to assess clam yield potential in a lagoon. I'm using the zeroinfl function of the pscl package to fit a Zero-inflated negative binomial model, given the high occurrence of zero counts. I don't understand from the sentence in the pscl guide "Zero-inflated count models are a type of two-component mixture model, with a component

Dear R crew: I think I am in the right mailing list. I have a very simple dataset consisting of two variables: cestode intensity and chick size (defined as CAPI). Intensity is clearly overdispersed, with way too many zeroes. I'm interested in looking at the association between these two variables, i.e. how well does chick size predict tape intensity? I fit a zero inflated negat. binomial

Hello, I'm using GLMM on the pscl package and i'm not getting the AIC on the summary. The code i'm using is (example) : mmall3 <-glmmPQL(allclues ~ cycloc + male, data=dados, family=poisson, random=~1|animal/idfid) and the results: Linear mixed-effects model fit by maximum likelihood Data: dados AIC BIC logLik NA NA NA Random effects: Formula: ~1 | animal

Hi, This may be such a general question that my searches are just failing. I installed the pscl lib, all appears fine, installed it several different ways to be sure, but I am getting: Error: could not find function "zeroinfl" I double checked my spelling of the function and that it had not been evolved out of the package. It is in the same location as the other libraries that are

Dear All, I am using the function vuong from pscl package to compare 2 non nested models NB1 (negative binomial I ) and Zero-inflated model. NB1 <-  glm(, , family = quasipoisson), it is an object of class: "glm" "lm" zinb <- zeroinfl( dist = "negbin") is an object of class: "zeroinfl"   when applying vuong function I get the following: vuong(NB1,

Dear R crew: I am sorry this question has been posted before, but I can't seem to solve this problem yet. I have a simple dataset consisting of two variables: cestode intensity and chick size (defined as CAPI). Intensity is a count and clearly overdispersed, with way too many zeroes. I'm interested in looking at the association between these two variables, i.e. how well does chick

I am interested in utilizing this so called "double hurdle" model in my study. We can write the model in the following way: if (z'a + u > 0 & x'b + e > 0) y = x'b + e, else y = 0 In the model, consumption y is the (left-) censored dependent variable. e and u are the normally distributed error terms. z'a is the participation equation and x'b is the

Hi, I am trying to replicate Lambert (1992)'s simulation with zero-inflated Poisson models. The citation is here: @article{lambert1992zero, Author = {Lambert, D.}, Journal = {Technometrics}, Pages = {1--14}, Publisher = {JSTOR}, Title = {Zero-inflated {P}oisson regression, with an application to defects in manufacturing}, Year = {1992}} Specifically I am trying to recreate Table 2. But my

Dear R-Helpers, Can anyone please help me to interpret warning messages from zeroinfl (package pscl) while fitting a zero inflated negative binomial model? The console reports convergence and the parameters seam reasonable, but these <<Warning messages: 1: algorithm did not converge in: glm.fit(X, Y, family = poisson()) 2: fitted rates numerically 0 occurred in: glm.fit(X, Y, family =

Hi all, I'm trying to run some monte carlo simulation for my roll call data using the ideal() function, which resides in the pscl package. However, I'm receiving an error message that I don't understand. Error in ideal(a, maxiter = 1000, thin = 10, burnin = 50, store.item = TRUE, : NA/NaN/Inf in foreign function call (arg 13) my code is simple the following: > m_a <-

Are you referring to the zeroinfl() function in the countreg package? If so, I think predict(fm_zinb2, type = "zero", newdata = some.new.data) will give you pi for each combination of covariate values that you provide in some.new.data where pi is the probability to observe a zero from the point mass component. As to your second question, I'm not sure that's possible, for any

Hi! I'm running a zero-inflated negative binomial regression on a large (n=54822) set of confidential data. I'm using the code: ZerNegBinRegress<-zeroinfl(Paper~.|., data=OvsP, dist="negbin", EM=TRUE) And keep getting the error: Warning message: glm.fit: fitted probabilities numerically 0 or 1 occurred I've done enough reading about this error to realize that I have

Hello, I have a dataset with 4 variables, each consisting of a vector, all with the same length. I start by interpolating the first three variables using the function "interp", and plot the interpolation successfully using "filled.contour". I then interpolate the first two variables and a fourth using "interp" again, but when I try to overlay the contour lines

Hello, I have count data that illustrate the presence or absence of individuals in my study population. I created a grid cell across the study area and calcuated a count value for each individual per season per year for each grid cell. The count value is the number of time an individual was present in each grid cell. For illustration my data columns look something like this and are repeated for

Good Day All, I have a negative binomial model which I have developed using the MASS library. I now would like to develop some predictions from it. Running the predict.glm (stats library) using type="response" gives me a non-integer value which was rather puzzling. I would like to confirm that this is actually the mean predicted value of the probability mass function as opposed

Dear fellows, I'm trying to extract the AIC statistic from a GLM model with quasipoisson link. The formula I'm referring to is AIC = -2(maximum loglik) + 2df * phi with phi the overdispersion parameter, as reported in: Peng et al., Model choice in time series studies os air pollution and mortality. J R Stat Soc A, 2006; 162: pag 190. Unfortunately, the function logLik

Hello, I have been using mgcv to run GAM hurdle models, analyzing presence/absence data with GAM logistic regressions, and then analyzing the data conditional on presence (e.g. without samples with no zeros) with GAMs with a negative binomial distribution. It occurs to me that using the negative binomial distribution on data with no zeros is not right, as the negative binomial allows zeros.

Dear all.. I have a respon variable y. Predictor variable are x1, x2, x3, x4, x5 (1) What is the syntax to get paramater estimation of ZINB Model by Newton Raphson (not BFGS) (2) What syntax to plot probability of observed & predicted of ZINB Thx. Regards Krist. [[alternative HTML version deleted]]