[R-br] Erro predição - Função predict

[Andre Oliveira]

Boa noite a todos! 
Estou procurando obter a métrica RMSE e R2 entre os objetos pred e diameter no modelo abaixo, mas está ocorrendo um erro. Parece que reclama de missing, mas ainda não consegui resolver, mesmo depois de muitas pesquisas. O interesse é pegar o predito levando em consideração os dois efeitos do modelo. Alguém teria alguma dica e sugestão?  


# CMDR

rm(list=ls(all=TRUE))require(lme4)data(Penicillin, package="lme4")attach(Penicillin)
library(MCMCglmm)model <- MCMCglmm(diameter ~ 1, random = ~ plate + sample, verbose= T, nitt=200000,burn=10000,thin=10, rcov = ~ units, data=Penicillin)
Penicillin$pred=predict(model, data=Penicillin, marginal = NULL)
Error in object$Sol[, missing.random, FALSE] <- rnorm(nrow(object$Sol) *  :                                         incorrect number of subscripts
Pedir os valores preditos como abaixo resolve, mas,  demora muito tempo em caso de bases de dados grandes.
Penicillin$pred=predict(model, data=Penicillin, marginal = NULL, interval="prediction")

att,.André 
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20191217/2024001b/attachment.html>