gnuplot> f(x) = m*x +n
gnuplot> fit f(x) "toreview.dat" u 2:3 via m,n
iter      chisq       delta/lim  lambda   m             n            
   0 3.6936711672e+01   0.00e+00  1.41e+00    1.000000e+00   1.000000e+00
   1 3.6992605848e-01  -9.88e+06  1.41e-01    6.105751e-01   4.583953e-01
   2 7.9414981509e-02  -3.66e+05  1.41e-02    8.423745e-01   2.635659e-02
   3 7.9329087727e-02  -1.08e+02  1.41e-03    8.467400e-01   1.870027e-02
   4 7.9329087724e-02  -3.50e-06  1.41e-04    8.467408e-01   1.869890e-02
iter      chisq       delta/lim  lambda   m             n            

After 4 iterations the fit converged.
final sum of squares of residuals : 0.0793291
rel. change during last iteration : -3.5012e-11

degrees of freedom    (FIT_NDF)                        : 22
rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 0.0600489
variance of residuals (reduced chisquare) = WSSR/ndf   : 0.00360587

Final set of parameters            Asymptotic Standard Error
=======================            ==========================
m               = 0.846741         +/- 0.02881      (3.403%)
n               = 0.0186989        +/- 0.0497       (265.8%)

correlation matrix of the fit parameters:
                m      n      
m               1.000 
n              -0.969  1.000