Yes, that's because I use a small value for `size`. As said in the PR description, maybe implementing analytical ML estimators for Nakagami would be preferred. I am reading about other issues and PRs about scipy.stats to know what the current direction is regarding MLE, especially analytical solutions or not.