Displaying 1 result from an estimated 1 matches for "mesocosm".
2013 Mar 15
0
Poisson and negbin gamm in mgcv - overdispersion and theta
Dear R users,
I am trying to use "gamm" from package "mgcv" to model results from a mesocosm experiment. My model is of type
M1 <- gamm(Resp ~ s(Day, k=8) + s(Day, by=C, k=8) + Flow + offset(LogVol),
data=MyResp,
correlation = corAR1(form= ~ Day|Mesocosm),
family=poisson(link=log))
where the response variable is counts, offset by the log of sample...