Hello All,

I am using the MCMCglmm package to fit a multivariate ordered probit model
for 5 ordinal variables. Parameter estimation was successful, so was the
prediction on the linear scale. However, when I tried prediction on the
response scale I got the following error:

> eq5d.mcmc2bfix.p.pred <- predict(eq5d.mcmc2bfix.p, type="response")
Error in M[, which(rm.v), ] <- 0 : incorrect number of subscripts

When I looked at the predict.MCMCglmm() function from the package, I
noticed that the M matrix was specified as a 3-dimentional array - hence
the incorrect number of substripts error.

 M <- model.matrix(~v.terms[nz] - 1)
                  M[, which(rm.v), ] <- 0

I am using the version:

> packageVersion("MCMCglmm")
[1] '2.17'


Note:
> eq5d.mcmc2bfix.p.pred <- predict(eq5d.mcmc2bfix.p, type="terms")
works with no error.

Regards,
Shamil.

	[[alternative HTML version deleted]]

