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]]