Dear all,
I am trying to run the function clmm() on a data table composed as following
:
<http://r.789695.n4.nabble.com/file/n4691592/datatable.png>
I have 187 pea lines assessed on 4 years * 2 locations * 3 blocs * 15
plantes for disease resistance. Disease resistance is assessed with a 0-to-5
scale which made an ordinal variable.
Here is the script a wrote :
>Tab.INRCh=read.csv("INRChamp201013.csv",sep=";")
>Tab.INRCh$Year=as.factor(Tab.INRCh$Year)
>Tab.INRCh$RRI<-factor(Tab.INRCh$RRI,levels=c("5","4","3","2","1","0"),ordered=TRUE)
>Tab.INRCh$RRI_DHW1<-factor(Tab.INRCh$RRI_DHW1,levels=c("5","4","3","2","1","0"),ordered=TRUE)
>Tab.INRCh$RRI_DHW2<-factor(Tab.INRCh$RRI_DHW2,levels=c("5","4","3","2","1","0"),ordered=TRUE)
>Modele=clmm(RRI~Line+RRI_DHW1+RRI_DHW2+(1|Year)+(1|Location)+(1|Not)+(1|Bloc)+(1|Plante),data=Tab.INRCh,na.action=na.omit)
I want to create this modele to run function Anova() to know if disease
resistance is different among lines.
When I run the script, I have an error :
Erreur : all(sapply(gfl, nlevels) > 2) is not TRUE
I can't find anything about it on the internet...
Can someone help me on this?
Thanks
Best regards
Aurore
--
View this message in context:
http://r.789695.n4.nabble.com/Error-clmm-ordinal-tp4691592.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
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.
