Hi,

I am trying to run a MCMCglmm model using multinomial family but getting an error message.

My prior is:

prior1<-list(
G=list((G1=list(V=0.01, nu = 0.1)),
    (G2=list(V=diag(2)* 0.01, nu = 0.1))),
R=list(V= diag(2)*0.01,nu=0.1))
prior1

my model is following:

m1<-MCMCglmm(cbind(M18, A18) ~ sex + treatment + group,
random = ~ batch + idh(group):mitotype,
rcov = ~ idh(group):units,
family = "multinomial2",
prior  = prior1,
verbose = FALSE,
nitt=150000,thin=100,burnin=50000,
data = SBR)

when I try to run the above model I get the following error message

Error in if (nJ < 1) { : missing value where TRUE/FALSE needed

Can any one please help me to understand what is wrong with my model setting?

Many Thanks

Sumayia

-- 
The University of Stirling has been ranked in the top 12 of UK universities for graduate employment*.
94% of our 2012 graduates were in work and/or further study within six months of graduation.
*The Telegraph
The University of Stirling is a charity registered in Scotland, number SC 011159.


	[[alternative HTML version deleted]]

