// Données

X = 1:4;
Y = [6, 5, 7, 10];

// régression
[a, b, sigma] = reglin(Y, X);

// points projetés sur la droite

X1 = a*Y + b;

// droite de régression
y1 = 4; x1 = a*y1 + b;
y2 = 10.25; x2 = a*y2 + b;

// tracé
clf;
couleurs = [get(sdf(), 'color_map') ; 0.75, 0.75, 0.75];
xset('colormap', couleurs);

xsegs([X ; X1], [Y ; Y], 14) // segments
h1 = gce();
h1.thickness = 2;

xpoly([x1, x2], [y1, y2]) // droite de régression
h2 = gce();
h2.thickness = 2;
h2.foreground = 2;

plot(X, Y, 'ok') // points
h3 = gce();
h3.children.mark_background = 5;

axe = gca();
axe.data_bounds = [0, 4 ; 5, 10.25];
axe.grid = [33, 33];
axe.tight_limits = 'on';
axe.isoview = 'on';
xtitle(' ', 'x', 'y')