That line looks weird. Looking at scipy.stats.t — SciPy v1.11.4 Manual, it seems you are passing scr as x and scr_Pred as the degrees of freedom of the StudentT.
If you go back to the Normal model and increase the number of trials, do inferences get closer to the true parameters?