[R] GAM and Poisson distribution

Ben Bolker bolker at ufl.edu
Sat Nov 15 10:49:31 CET 2008




sallyr wrote:
> 
> 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) ... : non-integer x = 5.325517
> 
> Thanks Sally
> 

If your response variable is an intensity (catch per unit effort) you
probably need to
include an offset for the total effort, e.g. something like
offset=log(effort)


-- 
View this message in context: http://www.nabble.com/GAM-and-Poisson-distribution-tp20508604p20513810.html
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list