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