<html><head></head><body><div class="ydp868c509eyahoo-style-wrap" style="font-family:bookman old style, new york, times, serif;font-size:16px;"><div><div dir="ltr" data-setdir="false"><div><div dir="ltr" data-setdir="false">Boa noite a todos! <br>Estou procurando obter a métrica RMSE e R2 entre os objetos <span><span style="color: rgb(0, 0, 0); font-family: bookman old style, new york, times, serif; font-size: 16px;"><b>pred e </b><span><span style="color: rgb(0, 0, 0); font-family: bookman old style, new york, times, serif; font-size: 16px;"><b>diameter</b> no modelo abaixo, mas está ocorrendo um erro. Parece que reclama de <span style="color: rgb(0, 0, 0); font-family: bookman old style, new york, times, serif; font-size: 16px;">missing, mas ainda não consegui resolver, mesmo depois de muitas pesquisas. <span><span style="color: rgb(0, 0, 0); font-family: bookman old style, new york, times, serif; font-size: 16px;">O interesse é pegar o predito levando em consideração os dois efeitos do modelo. </span></span>Alguém teria alguma dica e sugestão? </span><br></span></span></span></span><br><br><b># CMDR</b><br><br>rm(list=ls(all=TRUE))</div><div>require(lme4)</div><div>data(Penicillin, package="lme4")</div><div>attach(Penicillin)</div><div><br></div><div>library(MCMCglmm)</div><div>model <- MCMCglmm(diameter ~ 1, random = ~ plate + sample, verbose= T, nitt=200000,burn=10000,thin=10, rcov = ~ units, data=Penicillin)</div><div><br></div><div>Penicillin$pred=predict(model, data=Penicillin, marginal = NULL)</div><div><br></div><div><b>Error in object$Sol[, missing.random, FALSE] <- rnorm(nrow(object$Sol) * : </b></div><div><b> incorrect number of subscripts</b></div></div><br></div><div dir="ltr" data-setdir="false">Pedir os valores preditos como abaixo resolve, mas, demora muito tempo em caso de bases de dados grandes.</div><div dir="ltr" data-setdir="false"><br></div><div dir="ltr" data-setdir="false"><span><span style="color: rgb(0, 0, 0); font-family: bookman old style, new york, times, serif; font-size: 16px;">Penicillin$pred=predict(model, data=Penicillin, marginal = NULL, <span><span style="color: rgb(0, 0, 0); font-family: bookman old style, new york, times, serif; font-size: 16px;">interval="prediction"</span></span>)</span></span><br><br></div><div class="ydp868c509esignature"><div style="font-family:old times, serif;font-size:16px;"><div dir="ltr">att,.</div><div>André </div></div></div></div></div></body></html>