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é
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.