[R-sig-ME] MCMCglmm troubles

Mike Lawrence Mike.Lawrence at dal.ca
Thu Jun 17 14:43:38 CEST 2010

Hi folks,

I'm trying to do a mixed multinomial model with 3 response variables,
but I'm getting the error "please use idh() or us() error structure".
I believe that idh() and us() are used to specify the nature of the
interaction between multiple random variables, but I don't have
multiple random variables, so I'm not sure why I'm getting that error.
Some example data are below, but note that in my real data, the dv's
are indeed non-independent (they are in fact the number of trials
classified into each of three categories).

a = expand.grid(
	sid = 1:20
	, condition = factor(1:2)
a$group = factor(a$sid%%2)
a$sid = factor(a$sid)
a$dv1 = rnorm(nrow(a))
a$dv2 = rnorm(nrow(a))
a$dv3 = rnorm(nrow(a))

	, random = ~ sid
	, family = 'multinomial3'
	, data = a

Mike Lawrence
Graduate Student
Department of Psychology
Dalhousie University

Looking to arrange a meeting? Check my public calendar:

~ Certainty is folly... I think. ~

More information about the R-sig-mixed-models mailing list