library("Zelig")
library("stargazer")
library("car")
data("Prestige")
reg1 <- zelig(prestige ~ income + education, data = Prestige, model = "normal")
stargazer(reg1, type ="html", style = "demography")
prestige
income 0.001***
(0.0002)
education 4.137***
(0.349)
Constant -6.848*
(3.219)
N 102
Log Likelihood -353.863
AIC 713.725
p < .05; p < .01; p < .001

Based on the results, we can see that education is a stronger facter than income when predicting prestige.

```