Utilizzeremo il pacchetto:<<Statistics`DescriptiveStatistics`
![[Graphics:../Images/index_gr_11.gif]](../Images/index_gr_11.gif)
Eseguiremo ora i vari passi del calcolo della retta di regressione. Calcoliamo:
- le medie di X e di Y;
![[Graphics:../Images/index_gr_12.gif]](../Images/index_gr_12.gif)
- i vettori degli scarti di X e di Y dalla loro media;
![[Graphics:../Images/index_gr_16.gif]](../Images/index_gr_16.gif)
- il vettore delle varianze (quadrato del vettore degli scarti di X);
![[Graphics:../Images/index_gr_20.gif]](../Images/index_gr_20.gif)
- il vettore delle covarianze (vettore dei prodotti degli scarti di X e di Y);
![[Graphics:../Images/index_gr_22.gif]](../Images/index_gr_22.gif)
- la varianza e la covarianza di X e Y (sono le medie dei due vettori precedenti).
![[Graphics:../Images/index_gr_24.gif]](../Images/index_gr_24.gif)
Terminati i calcoli relativi a varianza e covarianza, possiamo utilizzati i risultati cosi' ottenuti per calcolare la pendenza della "retta di covarianza". Tale retta è quella che più si avvicina al grafico della dipendenza da X e Y ottenuto con le misure sperimentali, ed è definita da (covarianza/varianza).
![[Graphics:../Images/index_gr_28.gif]](../Images/index_gr_28.gif)
Possiamo ora disegnare il grafico della retta di regressione utilizzando la formula:
v(x) = pendenza*(x-mediaX)+mediaY
![[Graphics:../Images/index_gr_30.gif]](../Images/index_gr_30.gif)
Confrontiamo il risultato con il grafico della dipendenza da X e Y ottenuto con le misure sperimentali. Definiamo la matrice 2xl Mregr, contenente i punti della retta di regressione in corrispondenza dei valori del vettore X.
![[Graphics:../Images/index_gr_31.gif]](../Images/index_gr_31.gif)
Disegnamo il grafico di tale retta.
![[Graphics:../Images/index_gr_34.gif]](../Images/index_gr_34.gif)
Sovrapponiamo il grafico ottenuto dalle misure empiriche e quello della retta che lo approssima. Occorre usare il pacchetto:
<<Graphics`MultipleListPlot`
![[Graphics:../Images/index_gr_37.gif]](../Images/index_gr_37.gif)
![[Graphics:../Images/index_gr_38.gif]](../Images/index_gr_38.gif)