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é