search for: poisson

I would like to test in R what regression fits my data best. My dependent variable is a count, and has a lot of zeros. And I would need some help to determine what model and family to use (poisson or quasipoisson, or zero-inflated poisson regression), and how to test the assumptions. 1) Poisson Regression: as far as I understand, the strong assumption is that dependent variable mean = variance. How do you test this? How close together do they have to be? Are unconditional or conditional mea...

Hi All, I have a sequence of positive integers, which is right skewed. The mean value is 6 and variance is 11. I suspect it can be fitted by poisson distribution. But I'm not familiar with the function to fit distribution. Could you please help me with it? Also, once I fit the poisson distribution, how can I check the good-ness of this fitting? Thank you! [[alternative HTML version deleted]]

Dear Listers, I have a general statistical question. Are hurdle logit-poisson model and posson model nested? Thank you so much?

Craig, Thanks for your follow-up note on using the asypow package. My problem was not only constructing the "constraints" vector but, for my particular situation (Poisson regression, two groups, sample sizes of (1081,3180), I get very different results using asypow package compared to my other (home grown) approaches. library(asypow) pois.mean<-c(0.0065,0.0003) info.pois <- info.poisson.kgroup(pois.mean, group.size=c(1081,3180)) constraints <- matrix(c(2...

Dear All I would like to know whether it is possible with R to generate random numbers from zero-truncated Poisson distribution. Thanks in advance, Galina

Hello, from a search of the archives and functions, I am looking for information on creating random correlated counts from a multivariate Poisson distribution.  I can not seem to find a function that does this. Perhaps, it has not yet  been created. Has anyone created an R package that does this.   thanks,   Jourdan Gold     [[alternative HTML version deleted]]

Hi! I have fitted two glms assuming a poisson distribution which are: fit1 <- glm(Aids ~ Year, data=aids, family=poisson()) fit2 <- glm(Aids ~ Year+I(Year^2), data=aids, family=poisson()) I am trying to work out how to represent the fitted regression curves of fit1 and fit2 on the one graph. I have tried: graphics.off() plot(Aids ~...

Dear members, here is my trouble: My data consists of counts of trapped insects in different attractive traps. I usually use GLMs with a poisson error distribution to find out the differences between my traitments (and to look at other factor effects). But for some dataset where one traitment contains only zeros, GLM with poisson family fail to find any difference between this particular traitment and anyother one (even with traitment that...

I am announcing the release of the exactci package. It calculates exact tests and confidence intervals for binomial and Poisson tests. Here is an example to motivate the package: Suppose you want to see if the observed rates of 2/17877 for group A are significantly different from the observed rates of 10/20000 for group B assuming Poisson counts. The poisson.test function in the stats package gives a significant test res...

I am announcing the release of the exactci package. It calculates exact tests and confidence intervals for binomial and Poisson tests. Here is an example to motivate the package: Suppose you want to see if the observed rates of 2/17877 for group A are significantly different from the observed rates of 10/20000 for group B assuming Poisson counts. The poisson.test function in the stats package gives a significant test res...

I am using gam in mgcv (1.3-22) and trying to use gcv to help with model selection. However, I'm a little confused by the process of assessing GCV scores based on their magnitude (or on relative changes in magnitude). Differences in GCV scores often seem "obvious" with my poisson gams but with negative binomial, the decision seems less clear. My data represent a similar pattern as below - where I see (seemingly) drastic changes in GCV for the poisson with different model structures, but the negative binomial often seems only to change in the second or third decimal place f...

Hello R-help, According to a research article that covers the topic I'm analyzing, in Stata, a Poisson pseudo-maximum-likelihood (PPML) estimation can be obtained with the command poisson depvar_ij ln(indepvar1_ij) ln(indepvar2_ij) ... ln(indepvarN_ij), robust I looked up Stata help for the command, to understand syntax and such: www.stata.com/help.cgi?poisson Which simply says that the co...

...this is to model the mortality rate (or, more accurately, log of the mortality rate, log_m) as (say) a constant plus a proportion of age, plus time, so: r_1 <- lm(formula=log_m ~ age + time) summary(r_1) However, an alternative approach is to use try and estimate the number of deaths from the poisson mean mortality rate, and the number of people, with the poisson mean being defined in terms of age and time (And a constant). Conceptually I can see how this should work, in terms of linking the poisson probabilities together at each age and optimising the coefficients on age and time to maximise t...

This is a basics beginner question. I attempted fitting a a Poisson GLM to data that is non-integer ( I believe Poisson is suitable in this case, because it is modelling counts of infections, but the data collected are all non-negative numbers with 2 decimal places). My question is, since R doesn't return an error with this glm fitting, is it important that th...

Dear all, I have a very basic question: how does the logLik function work for poisson models? Example: I simulate 20 observations from a Poisson distribution with mean 800. y <- rpois(20,800) model <- glm(y ~ 1, family=poisson()) logLik(model) I would like to know what's the exact formula the function logLik uses. I looked at ?extractAIC but I cannot sort it out. Ca...

Hello all, I found a weird result of the GLM function that seems to be a bug. The code: a=c(rep(1,8),rep(2,8)) b=c(rep(0,8),rep(3,8)) cbind(a,b) model=glm(b~a, family=poisson) summary(model) generates a dataset with two groups. One group consists entirely of zeros, the other of 3's (as happened in a dataset I’m analyzing right now). Since they are count data, one should apply a poisson distribution. A GLM with poisson distribution delivers a p value > 0.99, thu...

Hi -I'm running a GAM with 7 explanatory variables with a Poisson error structure. All of the variables are continuous so I'm getting error messages in R. cod.fall.full.gam.model<-gam(Kept.CPUE~s(HOUR)+s(LAT_dec)+s(LONG_dec)+s(meantemp_C)+s(meandepth_fa)+s(change_depth)+s(seds), data=cod.fall.version2,family=poisson) In dpois(y, mu, log = TRUE) ... :...

Hello to all, I'm a biologist trying to tackle a "fish" (Poisson Regression) which is just too big for my modest understanding of stats!!! Here goes... I want to find good literature or proper mathematical procedure to calculate a confidence interval for an inverse prediction of a Poisson regression using R. I'm currently trying to analyse a "dose-r...

Dear list, This is a question about overdispersion and the ML estimates of the parameters returned by the glm.poisson.disp (L. Scrucca) and glm.nb (Venables and Ripley) functions. Both appear to assume a negative binomial distribution for the response variable. Paul and Banerjee (1998) developed C(alpha) tests for "interaction and main effects, in an unbalanced two-way layout of counts involving two fixed fa...

...s student searching in the huge pile of papers on my desk and on the Internet but I can't find out the solution. Would you mind giving me some help? Please. ######################################### I'm trying to use glm with factors: > Pyr.1.glm<-glm(Pyrale~Trait,DataRav,family=poisson) If I have correctly payed attention to my cyber professor explanations I have, for the variable Pyrale which I suppose Poisson-distributed, the following mathematical expression: P(Pyrale=k)=exp(-m).[(m^k)/k!] with log(m)=Intercept+Trait(i) (link function is log for Poisson distribution) The...