Estimados Hoy en un correo bajo otro tema se utilizaba predict ..., a raíz de esa pregunta realicé una consulta en un buscador y la respuesta me parece interesante, pero no comprendo la parte donde dice object_at_X, Zt, entiendo que en caso de no tener información en algunos puntos pide ejecutar esa parte del código, pero como se produce un error, no se a que hace referencia. Posiblemente al correr el código se interprete eso que no se expresar correctamente. ¿O eso sirve por ejemplo para: hay mediciones en escalas de edades del 1 al 10, pero en el 8 no se pudieron realizar estas mediciones, por lo cuál en realidad tengo 9 (falta la medición del 8), entonces ese object_at_X....Zt permite colocar el valor de 8 y por el modelo lo predice o estima? library(lme4) cbpp$Obs <- factor(seq_len(nrow(cbpp))) (gm1 <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd/Obs), family = binomial, data = cbpp)) str(gm1) predict.lmerBin <- function(object, X, Zt){ if(missing(X)){ X <- object at X } b <- fixef(object) if(missing(Zt)){ Zt <- as.matrix(object at Zt) } z <- unlist(ranef(object)) plogis(X %*% b + t(Zt) %*% z) } all.equal(as.vector(predict.lmerBin(gm1)), fitted(gm1)) [[alternative HTML version deleted]]