try glm.fit()