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]]

