def linreg(x,y,units='PW/Sv'):
""" Return linear regression model and plot label """
if len(x) > 1:
slope, intercept, r_value, p_value, std_err = stats.linregress(x,y)
y_model = x * slope + intercept
label = '(%5.3f %s)' % (slope, units)
else:
y_model = None
label = ''
return y_model, label
评论列表
文章目录