My outcome variable (y) is 3 categories. Is there anything bad about using
the following code to get a parameter estimate for my bivariate model?

publicfit = glm(y ~ public, data=dataSPSS.vmj, family=binomial)

	[[alternative HTML version deleted]]

