Baseline logit model  
# weights:  9 (4 variable)
initial  value 1846.767257 
final  value 1821.875901 
converged
Call:
multinom(formula = rsat ~ contact, data = housing, weights = freq)

Coefficients:
       (Intercept) contactLow
Medium  -0.1293244 -0.2572367
Top      0.2585748 -0.2174476

Residual Deviance: 3643.752 
AIC: 3651.752 

 Binary response for comparison purpose  

Call:  glm(formula = oksat ~ factor(contact), family = binomial, data = housing, 
    weights = freq)

Coefficients:
       (Intercept)  factor(contact)Low  
           -0.3720             -0.1053  

Degrees of Freedom: 71 Total (i.e. Null);  70 Residual
Null Deviance:	    2259 
Residual Deviance: 2258 	AIC: 2262

 Ordinal Logistic Regression with wrong order  

 Ordinal Logistic Reg. with right order, empty model  
Call:
polr(formula = factor(nsat) ~ 1, data = housing, weights = freq)

No coefficients

Intercepts:
       0|1        1|2 
-0.6752324  0.4164405 

Residual Deviance: 3648.878 
AIC: 3652.878 

 Ordinal Logist. Reg. with added covariates  
Call:
polr(formula = factor(nsat) ~ infl + type + contact, data = housing, 
    weights = freq)

Coefficients:
    inflLow  inflMedium  typeAtrium typeTerrace   typeTower 
 -1.2888275  -0.7224378   0.2061487  -0.5186378   0.5723756 
 contactLow 
 -0.3602933 

Intercepts:
       0|1        1|2 
-1.5728637 -0.3860317 

Residual Deviance: 3479.149 
AIC: 3495.149 

 Complimentary log log  
Call:
polr(formula = factor(nsat) ~ infl + type + contact, data = housing, 
    weights = freq, method = "cloglog")

Coefficients:
    inflLow  inflMedium  typeAtrium typeTerrace   typeTower 
 -0.9153805  -0.5333315   0.1266692  -0.3352559   0.4071971 
 contactLow 
 -0.2092257 

Intercepts:
       0|1        1|2 
-1.5136135 -0.6620318 

Residual Deviance: 3484.053 
AIC: 3500.053