# Dados
require(MCMCglmm)
data(PlodiaPO)
str(PlodiaPO)
# Divisão dos dados
library(caTools)
divisao = sample.split(PlodiaPO$PO, SplitRatio = 0.75)
conj_treinamento = subset(PlodiaPO,divisao==TRUE)
conj_validacao = subset(PlodiaPO$PO,divisao==FALSE)
# Modelo
require(sommer)
fit=mmer(PO ~ 1, random = ~ FSfamily, data=conj_treinamento)
summary(fit)
# Validação ???
pred=predict(fit, new_data = conj_validacao, classify = "FSfamily")
pred