Statistica.

Schede di studio e quiz

Lezione 27Macroargomento 4 · Statistica

Lezione 27. Correlazione e regressione lineare semplice

In questa lezione

  1. Che cosa misurano covarianza e correlazione?
  2. Perché la correlazione non implica causalità?
  3. Come si trova la retta di regressione?
  4. Come si misura l'associazione tra caratteri qualitativi?
  5. Quali indici misurano la dipendenza tra caratteri qualitativi?
  6. Quali sono le trappole tipiche per i concetti di questa lezione?
  7. Quiz della lezione (10 domande)

1. Che cosa misurano covarianza e correlazione?

IndiceFormulaValori
Covarianza campionariaC(x,y)=1n−1∑i(xi−xˉ)(yi−yˉ)C(x,y)=\dfrac1{n-1}\sum_i(x_i-\bar x)(y_i-\bar y)Qualsiasi; dipende dalle unità di misura
Indice di correlazionerxy=C(x,y)σxσyr_{xy}=\dfrac{C(x,y)}{\sigma_x\sigma_y}Tra −1 e 1

Perché misurano relazioni lineari. Per Cauchy-Schwarz ∣C(x,y)∣≤σxσy|C(x,y)|\le\sigma_x\sigma_y, con uguaglianza se e solo se yi=axi+by_i=ax_i+b. Quindi r=±1r=\pm1 solo per una relazione lineare perfetta (segno di aa).

rxyr_{xy}Lettura
Vicino a 1Forte correlazione positiva: al crescere di xx cresce yy
Vicino a −1Forte correlazione negativa: al crescere di xx cala yy
Vicino a 0Nessuna relazione lineare (ma può esserci una relazione quadratica o di altro tipo)

2. Perché la correlazione non implica causalità?

Correlazione osservataSpiegazione
Agenti di polizia e criminiVariabile nascosta: la dimensione della popolazione
Diffusione del nome «Sarah» e disboscamento in BrasileCoincidenza
Film di Nicolas Cage e annegamenti in piscinaCoincidenza

La correlazione serve a validare una previsione teorica, non a dimostrare una causa.

3. Come si trova la retta di regressione?

Metodo dei minimi quadrati. Si cercano a,ba,b che minimizzano F(a,b)=∑i=1n(yi−axi−b)2.F(a,b)=\sum_{i=1}^n(y_i-ax_i-b)^2.

  1. ∂bF=0\partial_bF=0 dà yˉ=axˉ+b\bar y=a\bar x+b: la retta passa per il punto medio (xˉ,yˉ)(\bar x,\bar y).
  2. ∂aF=0\partial_aF=0 e la sostituzione danno a=∑xiyi−nxˉyˉ∑xi2−nxˉ2a=\dfrac{\sum x_iy_i-n\bar x\bar y}{\sum x_i^2-n\bar x^2}.
  3. Numeratore e denominatore sono (n−1)C(x,y)(n-1)C(x,y) e (n−1)σx2(n-1)\sigma_x^2.
a=rxyσyσx=C(x,y)σx2,b=yˉ−axˉ.a=r_{xy}\frac{\sigma_y}{\sigma_x}=\frac{C(x,y)}{\sigma_x^2},\qquad b=\bar y-a\bar x.

Esempio aggiunto. x=1,2,3,4,5x=1,2,3,4,5; y=2,4,5,4,5y=2,4,5,4,5. xˉ=3\bar x=3, yˉ=4\bar y=4.

QuantitàCalcoloValore
∑(xi−xˉ)(yi−yˉ)\sum(x_i-\bar x)(y_i-\bar y)4+0+0+0+24+0+0+0+26
∑(xi−xˉ)2\sum(x_i-\bar x)^24+1+0+1+44+1+0+1+410
∑(yi−yˉ)2\sum(y_i-\bar y)^24+0+1+0+14+0+1+0+16
Pendenza aa6/106/100,6
Intercetta bb4−0,6⋅34-0{,}6\cdot32,2
rxyr_{xy}6/10⋅66/\sqrt{10\cdot6}≈0,77\approx0{,}77
02460123456xy
Scatter plot e retta di regressione y = 2,2 + 0,6x (r ≈ 0,77)

Con r≈0,4r\approx0{,}4 i punti sono dispersi intorno alla retta; con r≈0,9r\approx0{,}9 sono stretti intorno ad essa.

4. Come si misura l'associazione tra caratteri qualitativi?

Con dati non numerici (per esempio il colore di un prodotto) non esistono correlazione né retta. Si usano le frequenze relative condizionate: fl∣j=#{(xi,yi)=(aj,bl)}∣Yaj∣,Yaj={i:xi=aj}.f_{l|j}=\frac{\#\{(x_i,y_i)=(a_j,b_l)\}}{|Y_{a_j}|},\qquad Y_{a_j}=\{i:x_i=a_j\}.

Idea. Se la distribuzione di yy è la stessa per ogni valore di xx, i caratteri sono indipendenti.

EsperimentoDistribuzioni condizionateConclusione
Primo dado 1–3: secondo lancio con un d6; 4–6: con un d20. yy = sommaMolto diverse per primo lancio 2 o 4Dipendenza
Due dadi a 6 facce, 500 lanciQuella del secondo dado dato «primo = 3» è simile a quella totaleIndipendenza stocastica: fjl≈fajfblf_{jl}\approx f_{a_j}f_{b_l}

5. Quali indici misurano la dipendenza tra caratteri qualitativi?

IndiceChe cosa misuraValori
Distanza χ2\sqrt{\chi^2}Quanto differiscono tra loro le distribuzioni condizionate0 se sono tutte uguali
V di CramérDistanza dall'indipendenza stocasticaTra 0 e 1
Mutual information MIMIDistanza dall'indipendenza, misurata con l'entropia≥0\ge0
MI normalizzata NMINMIVersione normalizzata di MIMITra 0 e 1

Formule.

χ2=∑j<p∑l(fl∣j−fl∣p)2fl∣j+fl∣p\chi^2=\sum_{j<p}\sum_l\frac{(f_{l|j}-f_{l|p})^2}{f_{l|j}+f_{l|p}} V=1min⁡{m,k}−1∑j,l(fjl−fajfbl)2fajfblV=\sqrt{\frac1{\min\{m,k\}-1}\sum_{j,l}\frac{(f_{jl}-f_{a_j}f_{b_l})^2}{f_{a_j}f_{b_l}}} MI=∑j,lfjllog⁡fjlfajfbl,NMI=MIH(x)H(y)MI=\sum_{j,l}f_{jl}\log\frac{f_{jl}}{f_{a_j}f_{b_l}},\qquad NMI=\frac{MI}{\sqrt{H(x)H(y)}}

Qui fjlf_{jl} è la frequenza relativa della coppia (aj,bl)(a_j,b_l); fajf_{a_j}, fblf_{b_l} sono le frequenze marginali. Con pochi dati e molte categorie questi indici sono instabili.

6. Quali sono le trappole tipiche per i concetti di questa lezione?

7. Quiz della lezione (10 domande)

Domande di ripasso create per questa scheda a partire dalla dispensa (non sono il test ufficiale della piattaforma). Scegli una risposta: ricevi subito esito e spiegazione.