4. La retta che meglio approsima i dati di partenza.

Utilizzeremo il pacchetto:<<Statistics`DescriptiveStatistics`

[Graphics:../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]
[Graphics:../Images/index_gr_13.gif]
[Graphics:../Images/index_gr_14.gif]
[Graphics:../Images/index_gr_15.gif]

- i vettori degli scarti di X e di Y dalla loro media;

[Graphics:../Images/index_gr_16.gif]
[Graphics:../Images/index_gr_17.gif]
[Graphics:../Images/index_gr_18.gif]
[Graphics:../Images/index_gr_19.gif]

- il vettore delle varianze (quadrato del vettore degli scarti di X);

[Graphics:../Images/index_gr_20.gif]
[Graphics:../Images/index_gr_21.gif]

- il vettore delle covarianze (vettore dei prodotti degli scarti di X e di Y);

[Graphics:../Images/index_gr_22.gif]
[Graphics:../Images/index_gr_23.gif]

- la varianza e la covarianza di X e Y (sono le medie dei due vettori precedenti).

[Graphics:../Images/index_gr_24.gif]
[Graphics:../Images/index_gr_25.gif]
[Graphics:../Images/index_gr_26.gif]
[Graphics:../Images/index_gr_27.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]
[Graphics:../Images/index_gr_29.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]

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]
[Graphics:../Images/index_gr_32.gif]
[Graphics:../Images/index_gr_33.gif]

Disegnamo il grafico di tale retta.

[Graphics:../Images/index_gr_34.gif]

[Graphics:../Images/index_gr_35.gif]

[Graphics:../Images/index_gr_36.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]
[Graphics:../Images/index_gr_38.gif]

[Graphics:../Images/index_gr_39.gif]

[Graphics:../Images/index_gr_40.gif]


Converted by Mathematica      May 26, 2003