data = Table[{3 + i + RandomReal[{-3, 7}], i + RandomReal[{-2, 5}]}, {i, 1, 20}];
model=LinearModelFit[data,x,x]

Show[ListPlot[data], Plot[model["BestFit"], {x, 0, 30}]]