def nu(t, s=0.025, nu0=0.005): return pd.Series(utl.sig(0.5 * s * t + utl.logit(nu0)), index=t) def plotLDDecayTheoretical(sim):
def H(t, s=0.025, nu0=0.005): return 2 * utl.sig(0.5 * s * t + utl.logit(nu0)) * ( 1 - utl.sig(0.5 * s * t + utl.logit(nu0))) def nu(t, s=0.025, nu0=0.005): return pd.Series(utl.sig(0.5 * s * t + utl.logit(nu0)), index=t)