gnuplot> f(x) = m*x +n
gnuplot> fit f(x) "toreview.dat" u 2:3 via m,n

After 4 iterations the fit converged.
final sum of squares of residuals : 0.0531383
rel. change during last iteration : -3.75201e-12

degrees of freedom    (FIT_NDF)                        : 32
rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 0.0407501
variance of residuals (reduced chisquare) = WSSR/ndf   : 0.00166057

Final set of parameters            Asymptotic Standard Error
=======================            ==========================

m               = 1.00834          +/- 0.01949      (1.933%)
n               = -0.0311848       +/- 0.02683      (86.03%)


correlation matrix of the fit parameters:

               m      n      
m               1.000 
n              -0.965  1.000