% TP5ex4 : Régression linéaire gaussienne simple 
%          avec utilisation de linreg

clear;
clf;
x=[28 50 196 55 190 110 60 48 90 35 86 65 32 52 40 70 28 30 105 52 80 60 20 100];
Y=[130 280 800 268 790 500 320 250 378 250 350 300 155 245 200 325 85 78 375 200 270 295 85 495];
[tetac,Itetac,e,s,Is]=linreg(Y,x);
ac=tetac(2);
bc=tetac(1);
sig2c=s*s;
figure(2)
hold on
plot(x,Y,'*')
x=sort(x);
plot(x,ac+bc*x,'r')
hold off
