[R-es] Simulación de modelo logit con interacción

-----Mensaje original-----
De: R-help-es [mailto:r-help-es-bounces en r-project.org] En nombre de
Emilio Torres Manzanera
Enviado el: jueves, 22 de enero de 2015 12:28
Para: r-help-es en r-project.org
Asunto: [R-es] Simulación de modelo logit con interacción

Hola,
Deseo simular un modelo logit con interacción, estimar sus coeficientes
y comprobar si son o no parecidos al modelo teórico. Con este ejemplo
obtengo que los coeficientes estimados no se asemejan mucho a los
originales. ¿Se le ocurre a alguien cuál es el motivo de esta
discrepancia? ¿y cómo solucionarlo?
Muchas gracias
Emilio

logisticsimulation <- function(n){
  dat <- data.frame(x1=sample(0:1, n,replace=TRUE),
                    x2=sample(0:1, n,replace=TRUE))
  odds <- exp(-1 - 4 * dat$x1 + 7*dat$x2 - 1 *dat$x1* dat$x2 )
  pr <- odds/(1+odds)
  res <- replicate(100, {
    dat$y <- rbinom(n,1,pr)
    coef(glm(y ~ x1*x2, data = dat, family = binomial()))
  })
  t(res)
}

res <- logisticsimulation(100)
apply(res,2,median)
## (Intercept)          x1          x2       x1:x2
## -1.0986123 -18.4674562  20.4823593  -0.0512933

Deberían salir -1, -4, 7, 1