[R] Error in ordinal mixed effects
bgunter@4567 @end|ng |rom gm@||@com
Tue Jun 9 16:51:05 CEST 2020
Not sure, and we don't have your data, datos, but this is almost always a
bad thing to do:
res<-data.frame(cbind(id, group, time, cannabis))
Change it to:
res<-data.frame(id, group, time, cannabis)
## and you then won't need to name them either
and see if that fixes things.
Also, res$id is probably a factor: is order(res$id) really meaningful as a
"The trouble with having an open mind is that people keep coming along and
sticking things into it."
-- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )
On Tue, Jun 9, 2020 at 6:01 AM SUSANA ALBERICH MESA
<SUSANA.ALBERICHMESA using osakidetza.eus> wrote:
> 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")
> cannabismod<-mixor(cannabis~ time + as.factor(group), data=res, id=id,
> na.exclude, which.random.slope=na, link="logit")
> However, I have obtained this error:
> Error in xj[i] : invalid subscript type 'closure'
> Please, could anyone help me to solve it?
> Many thanks,
> [[alternative HTML version deleted]]
> R-help using r-project.org mailing list -- To UNSUBSCRIBE and more, see
> PLEASE do read the posting guide
> and provide commented, minimal, self-contained, reproducible code.
[[alternative HTML version deleted]]
More information about the R-help