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

classic Classic list List threaded Threaded
1 message Options
Reply | Threaded
Open this post in threaded view
|

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

R-br mailing list
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
[hidden email]
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.