I have some data with some errobars that I need to fit to a line. Is there anyway to use scipy.polyfit with the error associated with the data? If not, how can I make a fit routine work with my data? Thanks, Jeremy