[R] Error clmm(){ordinal}
Rune Haubo
rune.haubo at gmail.com
Wed Jun 4 10:21:37 CEST 2014
Aurore,
I don't know if car::Anova is able/should be able to produce anova
tables for clmm objects; I usually use drop1() (and sometimes add1) to
test terms in CLMMS:
> library(ordinal)
> fm1 <- clmm(rating ~ temp + contact + (1|judge), data=wine)
> drop1(fm1, test="Chi")
Single term deletions
Model:
rating ~ temp + contact + (1 | judge)
Df AIC LRT Pr(>Chi)
<none> 177.13
temp 1 209.59 34.464 4.343e-09 ***
contact 1 189.48 14.347 0.000152 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Here you get the more accurate likelihood ratio tests instead of the
faster and more convenient Wald tests that (if I recall correctly)
car::Anova can provide.
HTH,
Rune
On 4 June 2014 09:18, adesgroux <aurore.desgroux at rennes.inra.fr> wrote:
> Dear Rune,
>
> Thanks a lot for your comment!
> I tried with "Location" as a fixed-effect and clmm() ran! It take quiet a
> long time but I have approximatively 67000 observations in my datatable so I
> guess it's quiet normal!
>
> I want to run function Anova{car} on Modele but it gave my another error :
>>Anova(Modele)
> /Error in eval(expr, envir, enclos) : object 'RRI' not found/
>
> any idea?
> Thanks again
>
> Aurore
>
>
>
> --
> View this message in context: http://r.789695.n4.nabble.com/Error-clmm-ordinal-tp4691592p4691693.html
> Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________
> R-help at 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.
More information about the R-help
mailing list