search for: speciesf

...on is about the 'theta' parameter in specification of a NB GAMM. I have fit a GAM with an optimum structure of: SB.gam4<-gam(count~offset(vol_offset)+ s(Depth_m, by=StnF, bs="cs")+StageF*RegionF, family=negbin(1, link=log), data=Zoop_2011[Zoop_2011$SpeciesF=='SB',]) However, this GAM shows heterogeneity in the residuals, and the model fit is not great. I would like to add a random effect to this GAM, for station (StnF), and a spatial correlation structure, for Depth_m. When trying to fit an equivalent GAMM to start with, I used the formula:...

...h"), par.settings=list(strip.background=list(col="white"), fontzize=list(text=16,points=10), lty = 0), space=0, border=NA, col=c("blue","red","green","purple","orange"), stack=TRUE, groups=Zoo_nocyp$SpeciesF, key= list(title="Species", cex.title=1,text=list(c("BB","BC","CH","SB","VS")),space="right", rectangles=list(size=2, border="white",col=c("blue","red","green...