Hi,
Got my issues sorted.
Error message solved:
I heard from the guy who developed MuMIn and his suggestion worked.
"As for the error you get, it seems you are running an old version of
MuMIn.
Please update the package first."
I did (I was only 1 version behind in both R and in MuMIn) and error
disappeared!
Running quasipoisson GAM in MuMIn:
As for my questions on GAM and " to include only models with smooth OR
linear term (but not both) for each variable", as Gavin said, this refers
only to global models with both smooth and linear terms for the
*same*variables....which I didn't have.
Regards,
Karen
On 4 February 2011 12:31, Karen Moore <kmoore@tcd.ie> wrote:
>
> 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
>
>
>
>
>
--
Karen Moore
PhD Researcher,
FORESTBIO,
Department of Botany,
Trinity College Dublin
Ireland
Ph: 00 353 (0)87 9422 502
http://www.ucc.ie/en/planforbio
[[alternative HTML version deleted]]
