Le istruzioni per R - vedi - con cui è stata tracciata la precedente retta di regressione:

source("http://macosa.dima.unige.it/r.R")
x = c(31,63,87,109); y = c(3.8,6.9,9.6,11.8)
Plane(0,110, 0,12);  POINT(x,y, 1)  # ovvero: POINT(x,y,"black")
# La retta di regressione senza imporre che passi per O
regression1(x,y)
# 0.1033 * x + 0.5341 
g <- function(x) 0.1033 * x + 0.5341; graph1(g,0,110, "blue")

La retta di regressione vincolata a passare per (0,0):

# La retta di regressione passante per O
Plane(0,110,0,12);  POINT(x,y,1)
regression(x,y,0,0)
# 0.1097 * x 
h <- function(x) 0.1097 * x; graph1(h,0,110, "blue")

Ecco come ottenere la rappresentazione delle rette passanti per (0,0) più e meno inclinate passanti per tutti i rettangoli che rappresentano i dati sperimentali:

ex = rep(2,4); ey = rep(0.3,4)
Plane(0,110,0,12);  Point(x,y,1)
pointDif(0,0, x,y, ex,ey)
# 0.1060606 * x 
# 0.1130841 * x