Dear Ariane,
You must be referring to the Anova() function in the car package. The function
doesn't have an appropriate method for gls objects. The error was produced by
the default method, which is inappropriate for gls objects, as you discovered.
I agree that it would be nice for Anova() to handle gls objects, and it
shouldn't be too hard to do that. I'll consider this a request and put it on my
to-do list.
On the other hand, it should not be difficult to do what you want, by dropping
terms for the model and comparing alternative models by LR tests. It would
likely make most sense to start with the interaction,
mod.1 <- gls(maturity~morph*month,
weights = varIdent(form = ~ 1 |morph*month),
na.action=na.omit, method="ML")
mod.2 <- update(mod.1, . ~ . - morph:month)
anova(mod.1, mod.2)
Then, if the interaction is negligible, proceed to test the main effects:
mod.3 <- update(mod.2, . ~ . - morph)
anova(mod.2, mod.3)
mod.4 <- update(mod.2, . ~ . - month)
anova(mod.2, mod.4)
I hope this helps,
John
On Thu, 28 Jul 2011 18:37:58 +0000
Ariane Charaoui <ariane_...@hotmail.com> wrote:
>
>
>
> Hello,
>
> I have data
> on the maturity of two morphs of fish. I want to test whether their maturity
> is
> evolving differently or not on a temporal scale (month). The maturity
> variable (independent
> variable) is continuous and the morph and month variables (dependant
> variables)
> are categorical. Because the data show variance heterogeneity, I modeled it
> with the function gls:
>
> kg1 =
> gls(maturity~morph*month, weights = varIdent(form = ~ 1 |morph*month) ,
> na.action=na.omit, method="ML")
>
> Next, I
> want to test if the two effects morph and month are significant so I use
> the function anova
>
> anova(kg1)
>
> Denom. DF: 75
>
> numDF F-value p-value
>
> (Intercept) 1
> 174.20833 <.0001
>
> morph 1
> 5.37109 0.0232
>
> month 5 46.41181
> <.0001
>
> morph:month 5 2.85388
> 0.0206
>
>
>
> The problem is that I
> want the results of a non-sequential anova. I tried also the function Anova,
> but I get an error message:
>
>
>
> Error in !factors[,
> term2] : invalid argument type
>
>
>
> Is there a way to fix
> this problem or is it possible to specify a non-sequential anova?
>
>
>
> Thank you
>
>
>
> Ariane
>
>
>
>
>
> [[alternative HTML version deleted]]
>
------------------------------------------------
John Fox
Sen. William McMaster Prof. of Social Statistics
Department of Sociology
McMaster University
Hamilton, Ontario, Canada
http://socserv.mcmaster.ca/jfox/
______________________________________________
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.
