def elastic_evans(tau, alpha, tau_sd, alpha_sd): fit = fitting.odrlin(tau, alpha, tau_sd, alpha_sd) slope, intercept = fit.beta slope_sd, intercept_sd = fit.sd_beta elas = 1/slope elas_sd = slope_sd/(slope*slope) return slope, intercept, elas, elas_sd
def elastic_evans(tau, alpha, tau_sd, alpha_sd): fit = fitting.odrlin(tau, alpha, tau_sd, alpha_sd) slope, intercept = fit.beta slope_sd, intercept_sd = fit.sd_beta elas = 1 / slope elas_sd = slope_sd / (slope * slope) return slope, intercept, elas, elas_sd
def bending_evans(tau, alpha, tau_sd, alpha_sd): fit = fitting.odrlin(log(tau), alpha, tau_sd/tau, alpha_sd) slope, intercept = fit.beta slope_sd, intercept_sd = fit.sd_beta bend = 1/(8*pi*slope) bend_sd = slope_sd/(8*pi*slope*slope) # tau0 = exp(-intercept/slope) # tau0_sd = sqrt(slope*slope*intercept_sd*intercept_sd + slope_sd*slope_sd*intercept*intercept)*tau0/(slope_sd*slope_sd) return slope, intercept, bend, bend_sd
def bending_evans(tau, alpha, tau_sd, alpha_sd): fit = fitting.odrlin(log(tau), alpha, tau_sd / tau, alpha_sd) slope, intercept = fit.beta slope_sd, intercept_sd = fit.sd_beta bend = 1 / (8 * pi * slope) bend_sd = slope_sd / (8 * pi * slope * slope) # tau0 = exp(-intercept/slope) # tau0_sd = sqrt(slope*slope*intercept_sd*intercept_sd + slope_sd*slope_sd*intercept*intercept)*tau0/(slope_sd*slope_sd) return slope, intercept, bend, bend_sd