Hi,
I have a GAM quasipoisson that I'd like to run through MuMIn package
- dredge
- gettop.models
- model.avg
However, I'm having no luck with script from an example in MuMIn help file.
In MuMIn help they advise "include only models with smooth OR linear term
(but not both) for each variable". Their example is:
# Example with gam models (based on "example(gam)")
require(mgcv)
dat <- gamSim(1, n = 500, dist="poisson", scale=0.1)
gam1 <- gam(y ~ s(x0) + s(x1) + s(x2) + s(x3) + (x1+x2+x3)^2, family =
poisson, data = dat, method = "REML")
cat(dQuote(getAllTerms(gam1)), "\n")
# include only models with smooth OR linear term (but not both) for each
variable:
dd <- dredge(gam1, subset=(!`s(x1)` | !x1) & (!`s(x2)` | !x2) & (!`s(x3)` |
!x3))
But I'm not sure at all how to apply these instructions to my data
My formula is
GAM<-gam(Species richness (count) ~ Categorical + Continous + Continous + *
s*(Continous ) + Continous : Continous + Continous : Continous,
family=quasipoisson, data =)
Thanks for any advice on script
Karen
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
