It seems that an error exists in the code

[fangzhen-li]

In poisson_structure.py, the fit() function of Class PM2, I think it may be more suitable to convert all X_ to X.