MCMCglmm predict fails on MCMCglmm output generated 11/15/15
Hi,
I had to change how the predict function works in order to accommodate multivariate responses coming from different distributions. If you use:
model$Residual$family<-rep("gaussian", nrow(model$X))
it should run I think (if your response is gaussian).
Cheers,
Jarrod
On 11/09/2016 04:00, Kerry Samerotte wrote:
Hi mixed modelers,
I could really use some help with the MCMCglmm predict function. It works
for me on MCMCglmm output that I generated recently, but fails on MCMCglmm
output that I generated in November of 2015. In the past, the predict
function worked just fine on the 11/15 data, but now it fails with the
following error:
Error in if (!grepl("hu|zi|za|multinomial", object$Residual$family[i])) { :
argument is of length zero
Here is the model I run in MCMCglmm:
model = MCMCglmm(phen ~ condition, random = ~us(condition):strain,
rcov=~idh(condition):units, data=data, prior = prior_PE, thin = 15, nitt =
30000, burnin = 6000, verbose = FALSE, pr = TRUE)
Then I simply run:
predictions <- predict(model, marginal = NULL)
And it fails with the above error. This also fails with the same error:
predictions <- predict(model)
Help is appreciated!
Thank you
Kerry
[[alternative HTML version deleted]]
_______________________________________________ R-sig-mixed-models at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
The University of Edinburgh is a charitable body, registered in Scotland, with registration number SC005336.