search for: obspr

...sh 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 = obsPred$ar), mean) Is this correct? Thanks in advance > version _ platform i386-pc-mingw32...