这是我要解决的问题。请参见下面
logit = stats.MNLogit(y, x)
logit_fit = logit.fit(method='newton', full_output = True, maxiter = 1000, tol = 1e-8)
print(logit_fit.params)
logit_fit.summary()
outProb0 = logit_fit.predict(x)
logit_fit.params=avg_par
print('after change')
print(logit_fit.params)
outProb = logit_fit.predict(x)
logit_fit.summary()