search for: rekt1

...I don't know if anyone has answer my query. Here is the original message: I'm applying a gam model (package mgcv) to predict relative abundances of a fish species. The covariates are year, month, vessel and statistical rectangle. The model looks like this: g1 <- gam(log(cpue) ~ s(rekt1) + s(year) + s(mon) + s(reg1), data = dataTest) Once the model is fitted to the data I want to get the mean model estimates by year. I do the following: obsPred <- data.frame(year = dataTest$year, pred = predict(g1, type = "response")) gamFit <- tapply(obsPred$pred, list(year =...