loglm1 function

Fit Log-Linear Models by Iterative Proportional Scaling -- Internal function