data("HMDA", package = "AER")
fitBinomial <- walsGLM(deny ~ pirat + hirat + lvrat + chist + mhist + phist |
selfemp + afam, data = HMDA,
family = binomialWALS(),
prior = weibull(),
controlInitGLM = controlGLM(restricted = TRUE,
controlGLMfit = list(trace = TRUE)))
Run the code above in your browser using DataLab