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