load 'uh5yTm'
M = uh5yTm;

years = M (:,1);
Tm = M (:,2);

bb = robustfit (years,Tm);

b = bb (1);
m = bb (2);

line = years * m + b;

save ('/Users/sarahstuart/Desktop/MK/validated/temperature/overall/bb','bb');


load 'uh5yTm2'
M = uh5yTm2;

years2 = M (:,1);
Tm2 = M (:,2);

bb2 = robustfit (years2,Tm2);

b = bb (1);
m = bb (2);

line2 = years * m + b;

save ('/Users/sarahstuart/Desktop/MK/validated/temperature/overall/bb2','bb2');