Ukyeon

This change should be handled with #455 together.

We have identified 769 remaining errors, but many of them would be resolved by upcoming pull requests that address docstring or formatting issues. Therefore, we will handle these errors later...

The resulting parameters are the same after the updates.. est = ss_estimation() # After est2 = ss_estimation_prev() # Before np.array_equal(est.parameters["gamma"], est2.parameters["gamma"]) True np.array_equal(est.aux_param["gamma_intercept"], est2.aux_param["gamma_intercept"]) True np.array_equal(est.aux_param["gamma_r2"], est2.aux_param["gamma_r2"]) True np.array_equal(est.aux_param["gamma_logLL"], est2.aux_param["gamma_logLL"])...