ANOVA a una via

L'ANOVA a una via serve a testare le differenze tra le medie di gruppo

L'analisi della varianza a una via è un metodo statistico per testare l'ipotesi nulla (H₀) in base alla quale le medie di tre o più popolazioni sarebbero uguali, contro l'ipotesi alternativa (H_a) secondo cui almeno una media sarebbe diversa. Servendoci della notazione formale per le ipotesi statistiche, per k medie scriviamo:

$ H_0:\mu_1=\mu_2=\cdots=\mu_k $

$ H_a:\mathrm{non\mathrm{\ }tutte\ le\ medie\ sono\ uguali} $

dove $\mu_i$ è la media dell'i-esimo livello del fattore.

D'accordo, penserai, ma quando mai mi potrebbe capitare di dover determinare se le medie di più popolazioni sono le stesse o meno? Un caso comune è se sospetti che una determinata variabile di processo indipendente possa influire in maniera determinante sui risultati di tale processo. Per esempio, potresti avere dei dubbi su come diversi lotti produttivi, operatori o partite di materiali grezzi influiscano sui risultati (quindi sulle misurazioni di qualità) di un processo produttivo.

Per verificare i tuoi sospetti, puoi eseguire il processo utilizzando tre o più variazioni (o livelli) di questa variabile indipendente (o fattore) e ricavarne dei campioni per l'osservazione. Se dal confronto delle medie di ogni gruppo di osservazione tramite ANOVA emergono delle differenze, allora (presupponendo che tutto sia stato fatto nella maniera corretta!) hai dimostrato di essere nel giusto e che il fattore che hai analizzato influisce sui risultati ottenuti.

Esempio di ANOVA a una via

Vediamo l'analisi della varianza a una via applicata a un esempio pratico: immagina di lavorare per un'azienda produttrice di un gel adesivo che viene commerciato in barattolini. La viscosità del gel è fondamentale: se eccessiva rischia di rendere difficile l'applicazione, se insufficiente ne diminuisce l'adesività. Di recente, hai ricevuto lamentele da alcuni clienti, i quali sostengono che la viscosità del gel è più varia che in passato. Il tuo capo ti ha chiesto di indagare sulle possibili cause.

Decidi quindi che un primo passo può essere esaminare la viscosità media degli ultimi cinque lotti di produzione. Se dovessi riscontrare differenze tra un lotto e l'altro, vorrebbe dire che c'è davvero un problema. Inoltre, avresti delle basi per ipotizzare quali fattori causino le incongruenze tra un lotto e l'altro.

Per misurare la viscosità, ti servi di uno strumento dotato di un girante immerso nel barattolo di adesivo, che fornisce il valore della resistenza alla torsione. Il test viene effettuato su cinque barattoli presi a caso da ciascuno dei cinque lotti più recenti, per un totale di 25 barattoli. Dopo aver ottenuto la misura della resistenza alla torsione di ciascun barattolo, inserisci i dati in un diagramma.

Grazie al diagramma, ti accorgi che i valori della torsione dei barattoli del lotto numero 3 tendono a essere inferiori rispetto a quelli dei campioni prelevati dagli altri lotti. Facendo una media delle diverse misurazioni, ti accorgi che la resistenza alla torsione media del lotto numero 3 è pari a 26,77: molto più bassa degli altri quattro lotti, che presentano valori medi attorno a 30.

Tabella 1: Media della resistenza alla torsione dall'analisi dei cinque lotti di adesivo

La tabella ANOVA

I risultati dell'analisi della varianza in genere vengono presentati servendosi di una tabella. Una tabella ANOVA include:

• Origine: le origini della varianza, inclusi il fattore in esame (nel nostro caso il lotto), errori e totali.
• DF: gradi di libertà per ogni origine della varianza.
• Devianza (SS): somma dei quadrati per ogni origine della varianza, insieme al totale di tutte le origini.
• Media quadratica: la somma dei quadrati divisa per i relativi gradi di libertà.
• Statistica F: la media quadratica del fattore (lotto) divisa per la media quadratica dell'errore.
• Prob > F: il p-value.

Tabella 2: tabella ANOVA con i risultati ottenuti dalla misurazione della resistenza alla torsione

La derivazione dei componenti della tabella è illustrata di seguito. Un elemento chiave della tabella su cui vogliamo concentrarci è il p-value. Il p-value viene usato per verificare la validità dell'ipotesi nulla secondo cui tutte le medie si equivarrebbero. Nel nostro esempio, il p-value (Prob > F) è pari a 0,0012 e, per quanto piccolo, implica che le medie non sono tutte uguali. Dai nostri campioni emerge quindi che c'è differenza tra la resistenza alla torsione media di uno o più dei lotti in esame.

In base ai risultati ottenuti, decidi di sottoporre il lotto n° 3 a ulteriori test. Nel report, potresti scrivere: in seguito all'analisi ANOVA della resistenza alla torsione di cinque barattoli di prodotto per ognuno dei cinque lotti più recenti, supportata dall'osservazione, è stata rilevata una differenza di resistenza alla torsione media tra i diversi lotti (p = 0,0012). Un diagramma dei dati così ottenuti ha mostrato che il lotto n° 3 aveva una resistenza alla torsione media inferiore (26,77) rispetto agli altri quattro lotti, pertanto è stato trattenuto per essere sottoposto a ulteriori valutazioni.

Calcolo dell'ANOVA a una via

Ora, vediamo l'esempio della resistenza alla torsione più nel dettaglio: come ricorderai, avevamo cinque lotti di materiale, da cui avevamo prelevato cinque barattoli scelti a caso per ciascuno. In questo caso, quindi, abbiamo un piano fattoriale a un fattore. Il fattore, cioè il lotto, ha cinque livelli, che vengono replicati (testati) cinque volte. I risultati dei test sono elencati di seguito.

Tabella 3: Misure della resistenza alla torsione per lotto

Per meglio analizzare i calcoli risultanti nella tabella ANOVA di cui sopra (Tabella 2), dobbiamo prima stabilire le seguenti definizioni:

$n_i$ = numero di osservazioni del trattamento $i$ (nel nostro esempio il lotto $i$)

$N$ = numero di osservazioni totale

$Y_{ij}$ = j-esima osservazione dal livello i-esimo

$\overline{Y}_i$ = media campione per l'i-esimo trattamento

$\overline{\overline{Y}}$ = media delle osservazioni (media generale)

Somma dei quadrati

Tenendo a mente queste definizioni, possiamo affrontare la colonna della somma dei quadrati nella tabella ANOVA. La somma dei quadrati ci dà modo di quantificare la variabilità di un data set focalizzandoci sulla differenza tra ogni punto di dati e la media di tutti i punti di dati in quel data set. La formula riportata di seguito divide la variabilità generale in due parti: la variabilità dovuta al modello, o ai livelli del fattore, e la variabilità dovuta a un errore casuale.  

$$ \sum_{i=1}^{a}\sum_{j=1}^{n_i}(Y_{ij}-\overline{\overline{Y}})^2\;=\;\sum_{i=1}^{a}n_i(\overline{Y}_i-\overline{\overline{Y}})^2+\sum_{i=1}^{a}\sum_{j=1}^{n_i}(Y_{ij}-\overline{Y}_i)^2 $$

$$ SS(Totale)\;     =     \;SS(Fattore)\;     +     \;SS(Errore) $$

Se vista così l'equazione può sembrare complicata, analizzarne i singoli elementi uno per uno può aiutare a capirla più facilmente. La Tabella 4 qui di seguito elenca i diversi componenti della formula, riportandoli ai quadrati che vanno a costituire la somma dei quadrati. La prima colonna di dati ($Y_{ij}$) contiene le misurazioni della resistenza alla torsione raccolte nella Tabella 3 di cui sopra.

Tabella 4: Calcolo della somma dei quadrati

Gradi di libertà (DF)

A ogni somma dei quadrati corrisponde un valore chiamato grado di libertà (DF). Il grado di libertà indica il numero delle variabili indipendenti usate per calcolare ogni somma dei quadrati. Nei piani fattoriali a un fattore con un fattore a k livelli (nel nostro esempio i cinque lotti) e un totale di N osservazioni (cinque barattoli per lotto per un totale di 25) i gradi di libertà corrispondono a:

Tabella 5: Calcolo dei gradi di libertà

Media quadratica e statistica F

Per ottenere la media quadratica (MS), ogni somma dei quadrati viene divisa per i gradi di libertà corrispondenti. Se l'ipotesi nulla è vera (quindi le medie sono uguali), MS (Fattore) e MS (Errore) sono entrambe stime di errore nella varianza, per cui i valori risultanti dovrebbero essere molto simili e il rapporto tra i due, detto statistica F, dovrebbe essere vicino a 1. Quando invece l'ipotesi nulla è falsa, la MS (Fattore) è maggiore della MS (Errore) e il rapporto tra le due è maggiore di 1. Nel nostro esempio di test sugli adesivi la statistica F risultava pari a 6,90, chiara dimostrazione che l'ipotesi nulla secondo cui le medie sarebbero uguali era da confutare.

Tabella 6: Calcolo di medie quadratiche e statistica F