// Données

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

// régression verticale
[aa, bb, sigma] = reglin(X, Y);

// points projetés sur la droite

Yv = aa*X + bb;

// droite de régression verticale
xx1 = (y1 - bb)/aa;
xx2 = (y2 - bb)/aa;

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

xsegs([X ; X], [Y ; Yv], 14) // segments verticaux
h2 = gce();
h2.thickness = 2;

xpoly([xx1, xx2], [y1, y2]) // droite de régression verticale
h3 = gce();
h3.thickness = 2;
h3.foreground = 2;

plot(X, Y, 'ok') // points
h5 = gce();
h5.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')