unknown = np.random.normal(loc=1.1, scale=2.0, size=100)
Loc, Scale = ss.norm.fit_loc_scale(unknown)
unknown_cdf = lambda x: ss.norm.cdf(x, loc=Loc, scale=Scale)
plt.plot(np.linspace(-10, 10), unknown_cdf(np.linspace(-10, 10)), '-')
Is there a method to do arithmetic with SciPy's random variables?