And we determine the linear relationship from a best fit command in Mathematica.

Input := 

pfit[x_] = Fit[data,{1,x},x]

Output =

0.00998536 - 0.000994001 x