[R] Error in ordinal mixed effects

Jeff Newmiller jdnewm|| @end|ng |rom dcn@d@v|@@c@@u@
Tue Jun 9 15:57:26 CEST 2020


Can't tell... example is not reproducible because it is missing "datos'.

On June 9, 2020 4:33:24 AM PDT, SUSANA ALBERICH MESA <SUSANA.ALBERICHMESA using osakidetza.eus> wrote:
>Hi,
>I'm trying to run an ordinal mixed effects model with Mixor command. I
>have 65 cases and repeated visits in 0, 6, 9, 12 and 18 months. My code
>is the following:
>
>cannabis<-c(datos$cannabis0, datos$cannabis6, datos$cannabis9,
>datos$cannabis12, datos$cannabis18)
>time<-c(rep(0, 65), rep(6, 65), rep(9, 65), rep(12, 65), rep(18, 65))
>id<-c(rep(datos$id, 5))
>group<-c(rep(datos$group, 5))
>
>res<-data.frame(cbind(id, group, time, cannabis))
>names(res)<-c("id", "group", "time", "cannabis")
>res<-res[order(res$id),]
>
>cannabismod<-mixor(cannabis~ time + as.factor(group), data=res, id=id,
>na.exclude, which.random.slope=na, link="logit")
>summary(cannabismod)
>
>
>However, I have obtained this error:
>
>Error in xj[i] : invalid subscript type 'closure'
>
>Please, could anyone help me to solve it?
>
>Many thanks,
>Susana
>
>[https://edukiak.osakidetza.net/coronavirus/pie_email7.jpg]
>
>
>
>	[[alternative HTML version deleted]]
>
>______________________________________________
>R-help using r-project.org mailing list -- To UNSUBSCRIBE and more, see
>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.

-- 
Sent from my phone. Please excuse my brevity.



More information about the R-help mailing list