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é