Fit Gnuplot Error Bars
fit f(x) "data_set.dat" using 1:2:3 yerrors via par1, par2, par3
Wild Whale
fit f(x) "data_set.dat" using 1:2:3 yerrors via par1, par2, par3