%get cereals data 
Z=csvread('Cereals no alpha.csv');
[rows,cols]=size(Z);
disp([rows,cols])
sugar=Z(:,7);
fiber=Z(:,5);
rating=Z(:,13);
[b0,b1,rsq,StdErrEst,StdRes] = linfit(sugar,rating);
fit=b0+b1*sugar;
plot(sugar,StdRes,'*')