library("glmx")
## data
data("MexicanLabor", package = "glmx")
## visualizations
plot(I(laborforce/total) ~ interaction(income, age), data = MexicanLabor)
plot(I(laborforce/total) ~ interaction(schooling, locality), data = MexicanLabor)
Call: glm(formula = cbind(laborforce, total - laborforce) ~ ., family = binomial,
data = MexicanLabor)
Coefficients:
(Intercept) localityurban age> 24 incomehigh
-2.58575 0.38931 0.29212 -0.01521
schoolingfurther
0.54814
Degrees of Freedom: 15 Total (i.e. Null); 11 Residual
Null Deviance: 73.89
Residual Deviance: 15.96 AIC: 95.61