similar to: Convergence warnings from zeroinfl (package pscl)

Hello all: I have a question regarding the fitted.values returned from the zeroinfl() function. The values seem to be nearly identical to those fitted.values returned by the ordinary glm(). Why is this, shouldn't they be more "zero-inflated"? I construct a zero-inflated series of counts, called Y, like so: b= as.vector(c(1.5, -2)) g= as.vector(c(-3, 1)) x <- runif(100) # x

I have four files that map populations to appropriate keys: statepop.csv, countypop.csv, zip3.csv, and zip5.csv. I'm using the TIGER/Lines shape files: tl_2009_12_state, tl_2009_12_county tl_2009_12_zcta3, and  tl_2009_12_zcta5. I have carefully followed a number of tutorials and can reproduce those, and can also produce scatter and bubble plots of my data (using long and lat coordinates to

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 <-

Dear R users, I want to apply zero inflated models with continuous and categorical variables and I used pscl package from R and the zeroinf() function. My question are the follow: a) The value of fitted.values is mu or (1-p)*mu? where p is the probability of zero came form a zero point mass b) If mu is zero, how do i know if it is a zero from the zero point mass or from the count process?

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 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

I am attempting to fit a full zero-inflated Poisson model then use backward elimination to arrive at the best-fitting model. When I try to fit the model with zeroinfl() I get this error: Error in while (abs((ll_old - ll_new)/ll_old) > control$reltol) { : missing value where TRUE/FALSE needed Any suggestions? Thanks for your help! Paige Barlow MS Student Virginia Tech Dept Fish

Hello. I am using the zeroinfl package to fit a zero-inflated negative binomial. The explanatory variables are Year and Depth x STemp (interaction). I am in need of guidance for extracting the year effect and the associated precision. Thank you for your time. Cheers, Laura

Hi, I am wondering if there is an easy way to access the hat matrix for zeroinfl and hurdle objects in the pscl library? Thanks, Chris [[alternative HTML version deleted]]

I am fitting a model using zeroinfl() and it runs without errors, returning results that are generally consistent with my hypotheses. One of my variables is percent black (pblack). This variable was highly significant in some of the other count models I ran on the way to my current formulation. It is not significant in this model. As such I decided to try adding pblack^2 to the model to see

Dear list, I'm wondering if anyone can help me calculate the deviance of either a zeroinfl or hurdle model from package pscl? Even if someone could point me to the correct formula for calculating the deviance, I could do the rest on my own. I am trying to calculate a pseudo-R-squared measure based on the R^{2}_{DEV} of [1], so I need to be able to calculate the deviance of the full and null

Hi, I am a little confused at the output from predict() for a zeroinfl object. Here's my confusion: ## From zeroinfl package fm_zinb2 <- zeroinfl(art ~ . | ., data = bioChemists, dist = "negbin") ## The raw zero-inflated overdispersed data > table(bioChemists$art) 0 1 2 3 4 5 6 7 8 9 10 11 12 16 19 275 246 178 84 67 27 17 12 1 2 1 1

I'm trying to fit a zero-inflated poisson model using zeroinfl() from the pscl library. It works fine for most models I try, but when I include either of 2 covariates, I get an error. When I include "PopulationDensity", I get this error: Error in solve.default (as.matrix(fit$hessian)) : system is computationally singular: reciprocal condition number = 1.91306e-34 When I

Hi all, I have been trying to use zeroinfl() with the pscl package with R version 2.1.1. and with the newest versions of the contrib packages compatible with R 2.1.1. I have read the examples, the vignette and all the posts relating to zeroinfl() but I am still confused as to how to structure the model. Here is a small example; the error message is the same for big data sets

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

Hi! I have three models. In the first model, everything is fine. However, in the second and third models, I have NA's for standard errors: The hessians also have NaN's (same for m2 and m3). What should I do about it? It there a way to obtain the hessian without transforming my variables? I will greatly appreciate your help! -- View this message in context:

Hello, I have a data frame with 1425 observations, 539 of which are zeros. I am trying to fit the following ZINB: f3<-formula(Nbr_Abs~ Zone * Year + Source) ZINB2<-zeroinfl(f3, dist="negbin", link= "logit", data=TheData, offset=log(trans.area), trace=TRUE) Zone is a factor with 4 levels, Year a factor with 27 levels, and Source a factor with 3 levels. Nbr_Abs is counts

Dear R Users, When using zeroinfl() function to fit a Zero-Inflated Negative Binomial (ZINB) model to a dataset, the summary() gives an estimate of log(theta) and its standard error, z-value and Pr(>|z|) for the count component. Additionally, it also provided an estimate of Theta, which I believe is the exp(estimate of log(theta)). However, if I would like to have an standard error of Theta

Hi Just a quick problem that I hope is simple to resolve. I'm doing some work with zero inflated poisson models using the pscl package. I can build models using zeroinfl and get outputs fom them with no problem, but when I try to use the predict.zeroinfl function, I get Error: could not find function "predict.zeroinfl". I was using an older version of R, but still had the same