Hazlo con n = 10000 Est谩s modelizando una variable binaria con otras 2 binarias: "muy chungo"
Un saludo. Isidro Hidalgo Arellano Observatorio Regional de Empleo Consejer铆a de Empleo y Econom铆a http://www.jccm.es > -----Mensaje original----- > De: R-help-es [mailto:r-help-es-boun...@r-project.org] En nombre de > Emilio Torres Manzanera > Enviado el: jueves, 22 de enero de 2015 12:28 > Para: r-help-es@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 > > _______________________________________________ > R-help-es mailing list > R-help-es@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-help-es _______________________________________________ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es
