summary fn zero-inflated MCMC models.

Dear all / Jarrod

Have just updated to the latest version of MCMCglmm but am now  
unable to produce model info using the summary function for zero  
inflated models. It doesn't work on hurdle models fitted on the  
previous version either (but these do on an older version I have on  
my laptop).

Error is

Error in rep(rep(1:length(object$Residual$nrt), object$Residual$nrt),  :
  invalid 'times' argument

Unless something was wrong with my model specification previously  
and something in the update is now showing that?

Best wishses
John Hodsoll

prior and model definitions..

#### Define prior  #####
zal2.prior <-  list(B= list (mu = matrix(c(rep(0,20)),20),V = diag(20)*(20)),
                    R = list(V = diag(2), n = 0.002, fix = 2),
                    G = list(G1 = list(V = diag(2), n = 0.002)))


zal2.prior$B$mu[9] <- 1
diag(zal2.prior$B$V)[9]<-1e-9

zal2.prior$B$mu[10] <- 1
diag(zal2.prior$B$V)[10]<-1e-9

#Level 2 models with complementary log log link.

system.time(
  cf.za.1 <- MCMCglmm(totflct ~ trait-1 + trait:(expcon.r*period.x +  
log.beds + nshift + cday),
                      data = rswAll.df, family = "zapoisson",
                      random = ~us(trait):wardn,
                      rcov = ~idh(trait):units,
                      prior = zal2.prior,
                 #     nitt = 300000, burnin = 50000, thin = 250,
                      verbose = TRUE, pr = TRUE, pl = FALSE)
)

summary(cf.za.1)
save(cf.za.1, file="cf.za.1.rda")

	[[alternative HTML version deleted]]