Introduzione
Come noto i trasmettitori di Radio Vaticana situati a S.M. di
Galeria, presso Cesano (provincia di Roma), sono da tempo sotto accusa per i
casi di luecemia (particolarmente leucemia infantile) che si sono verificati
nella zona intorno alla stazione radio. Ne è scaturito un processo ancora in
corso.
Uno studio epidemiologico dei casi di leucemia occorsi dal 1987 nella zona, in un raggio di 10 km dall'emittente, è stato condotto dalla Azienda di Sanità Pubblica (ASP) del Lazio. I risultati, secondo gli autori, non sono conclusivi a causa del piccolo numero di casi e dalla carente misura dell'esposizione dei residenti nella zona (circa 60'000 abitanti). La stampa è invece giunta a conclusioni opposte, ritenendo evidentemente di essere più capace degli stessi autori dello studio.
Uno studio epidemiologico dei casi di leucemia occorsi dal 1987 nella zona, in un raggio di 10 km dall'emittente, è stato condotto dalla Azienda di Sanità Pubblica (ASP) del Lazio. I risultati, secondo gli autori, non sono conclusivi a causa del piccolo numero di casi e dalla carente misura dell'esposizione dei residenti nella zona (circa 60'000 abitanti). La stampa è invece giunta a conclusioni opposte, ritenendo evidentemente di essere più capace degli stessi autori dello studio.
Vista del centro trasmittente con alcune delle antenne indicate
![]() |
![]() |
| Una delle antenne rotanti per onde corte (C) | Le 4 torri direttive del trasmettitore a 1530 kHz (B) |
In questa pagina faremo
passo per passo un'analisi elementare dei dati di Cesano per mostrare (speriamo)
che è alla portata di chiunque mastica un po' di matematica. Cercheremo di
spiegare alcuni concetti di statistica che, benché elementari, non sembrano
molto diffusi. Vogliamo anche mettere in evidenza quali sono le difficoltà
oggettive nell'interpretare i risultati e come si potrebbe raggiungere una
maggiore sicurezza nelle conclusioni. Lo scopo di questa pagina è essenzialmente
didattico, per cercare di avvicinare anche chi non è del mestiere alle
problematiche con cui ci si scontra in questo genere di studi, studi che
sembrano appassionare l'opinione pubblica, se si giudica dai dibattiti a volte
feroci. Per questo motivo la terminologia è necessariamente semplificata e non
tecnica.
I dati
L'analisi dei casi di leucemia
nella zona di Cesano in un raggio R=10 km intorno all'emittente di Radio
Vaticana è stata effettuata dall'ASP Lazio [1-2]. Gli stessi dati sono stati
rianalizzati da una commissione dell'Istituto Superiore di Sanità [3]. I dati
raccolti sono il numero di decessi per gli adulti e il numero di casi di
malattia per i bambini (0-14 anni). Il numero di casi attesi di morte o
incidenza della malattia è stato calcolato prendendo le popolazioni residenti
nel cerchio e moltiplicandole per i relativi tassi del Comune di Roma nell'anno
in esame. L'analisi copre 12 anni (1987-1998) per gli adulti e 13 anni
(1987-1999) per i bambini.
Bibliografia
(1) "Mortalità per leucemia nella popolazione adulta ed incidenza di leucemia infantile in prossimità della stazione Radio Vaticana di Roma" P. Michelozzi et al., Epidemiologia e Prevenzione no. 6, p. 249 (2001)(2) "Adult and Chilhdood Leukemia near a High-Power Radio Station in Rome, Italy" P. Michelozzi et al., Am. J. Epidemiol. Vol. 155, No. 12, p. 1096 (2002)
(3) "Esposizione a campi di radiofrequenza e leucemia infantile: stato attuale delle conoscenze scientifiche in rapporto alle problematiche dell'area di Cesano" D. Greco et al., Rapporto ISTISAN 01/25.
Per una discussione fra gli autori vedi anche Tempo Medico 02/2002
|
La cosiddetta
"ipotesi nulla" H0 è che nell'area di Cesano la mortalità
"naturale" per leucemia sia la stessa che nel Comune di Roma.
|
Per
evidenziare una possibile dipendenza dalla distanza dal centro (ossia dai
trasmettitori radio) il cerchio è stato suddiviso in 5 anelli concentrici di
ampiezza radiale 2 km. I risultati (da [2]) sono riportati nelle due tabelle che
seguono.
Decessi da leucemia negli adulti nella zona di Cesano nel periodo 1987-1998
| 0-2 km | 2-4 km | 4-6 km | 6-8 km | 8-10 km | Totale | ||
| Donne |
Oss. | 0 | 3 | 5 | 6 | 5 | 19 |
| Att. | 0,4 | 2,3 | 4,7 | 4,1 | 3,6 | 15,1 | |
| Uomini | Oss. | 2 | 6 | 7 | 5 | 1 | 21 |
| Att. | 0,7 | 3,7 | 7,2 | 6,2 | 4,4 | 22,2 | |
| Totale | Oss. | 2 | 9 | 12 | 11 | 6 | 40 |
| Att. | 1,1 | 6,0 | 11,9 | 10,4 | 8,00 | 37,4 |
Casi di leucemia infantile nella zona di Cesano nel periodo 1987-1999
| 0-2 km | 2-4 km | 4-6 km | 6-8 km | 8-10 km | Totale | ||
| Bambini | Oss. | 1 | 2 | 5 | 0 | 0 | 8 |
| Att. | 0,2 | 0,9 | 2,6 | 1,7 | 1,4 | 6,8 |
E' da notare che, essendoci sostanziali differenze fra la
leucemia che colpisce adulti e bambini, i casi sono stati analizzati
separatamente.
Per rispondere a questa domanda divideremo i campioni di dati in due categorie:
Un po' di statistica
La domanda a cui si deve rispondere è questa: è plausibile che i casi di leucemia osservati non abbiano nulla a che vedere con la presenza dell'emittente radio oppure, viceversa, questa è responsabile di un aumento dei casi nelle sue vicinanze? La risposta a questa domanda si deve basare necessariamente sul concetto di probabilità dato che il contrarre una malattia come il cancro è il risultato di un processo aleatorio ovvero probabilistico.Per rispondere a questa domanda divideremo i campioni di dati in due categorie:
- casi vicini all'antenna
- casi lontani dall'antenna
- chiameremo vicini i primi tre anelli
- chiameremo lontani gli ultimi due
Tabella per gli adulti vicini e lontani
| Vicini (0-6 km) | Lontani (6-10 km) | Totale | ||
| Adulti | Oss. | 23 | 17 | 40 |
| Att. | 19 | 18,4 | 37,4 |
Il numero di casi attesi è ovviamente calcolato come se l'antenna non avesse effetto alcuno. Si tratta
quindi dell'incidenza "naturale" della malattia. Naturalmente non ha senso parlare di casi frazionari
se non in maniera probabilistica.
Nel lancio della moneta si ammette che testa e croce abbiano la stessa probabilità di verificarsi, entrambe pari al 50 % (0,50). La nostra moneta non è però perfettamente simmetrica: infatti ci attendiamo alla lunga di avere 19 teste su 37,4 (oppure 190 su 374):
Costruiamo una tabella per il lancio della moneta
|
Ma avere 40 casi di malattia suddivisi in due categorie (23 vicini e 17 lontani) è esattamente lo stesso
che avere 40 lanci di una moneta in cui si siano verificate 23 teste e 17 croci. Troppi
casi "vicini" è esattamente come troppe "teste". Ma che si intende per "troppo"? |
|
Nel lancio della moneta si ammette che testa e croce abbiano la stessa probabilità di verificarsi, entrambe pari al 50 % (0,50). La nostra moneta non è però perfettamente simmetrica: infatti ci attendiamo alla lunga di avere 19 teste su 37,4 (oppure 190 su 374):
- P(testa) = 19/37,4 = 0,51
- P(croce) = 18,4/37,4 = 0,49
Costruiamo una tabella per il lancio della moneta
Tabella per il lancio della moneta
| No. teste | No. croci | Totale | ||
| Risultati | Osservati | 23 | 17 | 40 |
| Probabilità | 0,51 | 0,49 | 1,00 |
e domandiamoci se l'osservazione di 23 teste e 17 croci sia compatibile con le rispettive
probabilità (che corrispondono all'ipotesi H0) oppure se le teste (o le croci)
non siano "troppe".
Per fare questo calcoleremo la probabilità di osservare esattamente 23 teste su 40 lanci partendo dalla cosiddetta distribuzione binomiale con probabilità p=0,51 trovando
Questo risultato significa che, ripetendo l'esperimento (che consiste in 40 lanci della moneta), la probabilità di ottenere 23 o più teste è del 25%. Una probabilità del 25% è considerata molto elevata per cui l'osservazione di 23/17 è ritenuta del tutto accettabile. Ne consegue che il numero di decessi da leucemia fra gli adulti è perfettamente in linea con l'ipotesi che l'antenna non abbia affetto alcuno.
Per fare questo calcoleremo la probabilità di osservare esattamente 23 teste su 40 lanci partendo dalla cosiddetta distribuzione binomiale con probabilità p=0,51 trovando
P(23 |H0) = 0,090.
A questa P sommiamo la probabilità di avere 24 teste, poi 25, poi 26 ecc., fino a 40. Otteniamo così la cosidetta probabilità cumulativa ossia la probabilità di osservare almeno 23 teste, che vale|
|
P(n≥23|H0) = 0,25 = 25 % |
Questo risultato significa che, ripetendo l'esperimento (che consiste in 40 lanci della moneta), la probabilità di ottenere 23 o più teste è del 25%. Una probabilità del 25% è considerata molto elevata per cui l'osservazione di 23/17 è ritenuta del tutto accettabile. Ne consegue che il numero di decessi da leucemia fra gli adulti è perfettamente in linea con l'ipotesi che l'antenna non abbia affetto alcuno.
| La convenzione utilizzata comunemente in biostatistica è che si considerano significative probabilità minori del 5%. Diciamo in tal caso che l'osservazione falsifica l'ipotesi H0. |
Perché il 5% ?
La scelta del 5% per discriminare fra effetto significativo e non è arbitraria e fu fissata da Fisher, un famoso statistico, che si domandava come si potesse valutare l'efficacia della concimazione di un terreno coltivato. Niente impedisce di utilizzare altre convenzioni e sicuramente questo avviene spesso. In genere però si va nella direzione di diminuire la probabilità. Più basso è questo valore e più si abbassa il rischio di rigettare un'ipotesi vera.| In effetti
P=5% è il rischio di rigettare l'ipotesi H0 nel caso che
questa sia vera. Se H0 fosse l'ipotesi che l'imputato
in un processo è innocente P=5% rappresenterebbe il
rischio di condannare un innocente (di non riconoscerlo come tale). In tal
caso forse il 5% ci sembrerebbe un rischio eccessivo e vorremmo avere una
maggiore sicurezza per emettere la sentenza Ma se invece del 5%
scegliessimo l'1% molti veri colpevoli resterebbero impuniti. Il 5% è in
definitiva un compromesso e in questa scelta la scienza non ha molta
voce in capitolo. |
Per fare un altro esempio: in Fisica per rigettare o "falsificare" un'ipotesi o una teoria si chiede che la probabilità di errore sia meno dello 0,1%. E inoltre si richiede sempre la ripetibilità degli esperimenti. Non si da un premio Nobel se un esperimento non è stato confermato da altri esperimenti indipendenti. Se il limite fosse elevato al 5% saremmo ogni anno sommersi da "sensazionali" scoperte del tipo "Einstein sbugiardato" che piacciono tanto ai giornali. |
Un principio di precauzione?
Nel caso di Cesano si potrebbe obbiettare che è meglio rischiare di sbagliarsi accettando una P di errore elevata dato che vi è in gioco la salute delle persone. Tuttavia un ragionamento di questo tipo, che ricorda il Principio di Precauzione, anche se fatto a fin di bene, è pericoloso.In questo contesto specifico siamo prima di tutto di fronte a un processo i cui ci sono delle responsabilità penali in gioco che non si possono attribuire con leggerezza o per "precauzione". Ciò corrisponderebbe a dire che per precauzione è meglio un innocente in galera che un colpevole a piede libero. In questo caso la scelta del valore appropriato di P è una responsabilità che spetta al Tribunale e non alla scienza e pertanto non ci interessa qui. Se invece dimentichiamo gli aspetti processuali e cerchiamo di ragionare scientificamente, ignorando risvolti emotivi e magari anche ideologici, dobbiamo innanzitutto considerare di essere di fronte a uno studio che va ad aggiungersi a una mole di studi precedenti. Da questo punto di vista scartare (per sbaglio) un'ipotesi vera per sostituirla con un'altra falsa non è un servizio reso a nessuno: contribuisce solo a ritardare, mettendoci fuori strada, la corretta comprensione di un fenomeno che riguarda la salute delle persone: è importante capire, non trovare un capro espiatorio. Infine non si deve dimenticare che la scienza è democratica. Questo studio va ad aggiungersi ad una mole di studi già esistenti e non abbiamo nessuna ragione a priori per preferire le sue conclusioni e scartare le altre che non ci piacciono. Tutto va visto nel suo insieme. Se la scienza non dà certezze meno ancora le può dare una singola misura.
I bambini
I casi di leucemia infantile
hanno ricevuto molta più attenzione di quelli degli adulti e ciò è del tutto
comprensibile. I genitori sono alla ricerca di una causa "esterna" responsabile
di un fatto doloroso come la malattia o la morte di un bambino che è difficile
accettare come evento "naturale".
Nel caso dei bambini si sono potuti esaminare i casi di incidenza della malattia e non solo i casi di morte (vedi Tabella precedente).
Ripetiamo l'esercizio precedente dividendo il campione osservato in "teste" e "croci". Abbiamo la tabella seguente
Questa volta lo sbilanciamento dei casi osservati è notevole, mentre le probabilità di
osservare testa o croce sono cambiate di poco. La probabilità di
osservare 8 teste e 0 croci vale
Nel caso dei bambini si sono potuti esaminare i casi di incidenza della malattia e non solo i casi di morte (vedi Tabella precedente).
Ripetiamo l'esercizio precedente dividendo il campione osservato in "teste" e "croci". Abbiamo la tabella seguente
Tabella per i bambini
| No. teste | No. croci | Totale | ||
| Risultati | Osservati | 8 | 0 | 8 |
| Probabilità | 0,54 | 0,46 | 1,00 |
P(8 |H0) = 0,00723
Dato che abbiamo osservato il massimo possibile di teste (8 in 8 lanci)
questa coincide con la probabilità della nostra ipotesi:
|
Tuttavia la limitata statistica (solo 8 casi) non permette una conclusione
robusta: se infatti uno solo dei casi
vicini fosse stato male assegnato avremmo
|
Cosa concludere?
Da quel che si è detto dovrebbe essere chiaro che per gli adulti non vi è nessun motivo
per ritenere che vi sia un aumento di morti di leucemia legati alla vicinanza dell'antenna.
La statistica è già sufficientemente elevata per rendere la conclusione robusta.
Diverso è il caso per i bambini, in cui l'ipotesi nulla non è compatibile con l'osservazione. Tuttavia, dato che la statistica è estremamente limitata, questa conclusione dovrebbe essere corroborata da ulteriori studi. Sarebbe quindi fondamentale che il Tribunale o la ASP prolungassero l'analisi fino ad includere i dati del 2004 o del 2005, una cosa abbastanza semplice visto che il metodo di studio è ormai stato collaudato. Sembra in effetti che sarà richiesta una perizia aggiornata. Cosa si guadagnerebbe?
Se consideriamo i bambini: includendo i dati del 2005 lo studio coprirebbe un arco di tempo di 19 anni. Il numero di casi osservati crescerebbe presumibilmente del rapporto 19/13 = 1,5.
Diverso è il caso per i bambini, in cui l'ipotesi nulla non è compatibile con l'osservazione. Tuttavia, dato che la statistica è estremamente limitata, questa conclusione dovrebbe essere corroborata da ulteriori studi. Sarebbe quindi fondamentale che il Tribunale o la ASP prolungassero l'analisi fino ad includere i dati del 2004 o del 2005, una cosa abbastanza semplice visto che il metodo di studio è ormai stato collaudato. Sembra in effetti che sarà richiesta una perizia aggiornata. Cosa si guadagnerebbe?
Se consideriamo i bambini: includendo i dati del 2005 lo studio coprirebbe un arco di tempo di 19 anni. Il numero di casi osservati crescerebbe presumibilmente del rapporto 19/13 = 1,5.
| Includendo i dati fino al 2005 il numero
totale di casi di leucemia osservati su tutta l'area potrebbe diventare 12.
Per avere una P del 5% o inferiore dovremmo quindi osservare un numero di casi vicini almeno uguale a 10. Questo sarebbe un risultato significativo ma non ancora robusto. Ci vorrebbero 11 o 12 casi su 12 per essere praticamente certi che l'ipotesi H0 non sta in piedi. |
Quanto costa?
Per una nuova perizia epidemiologica che includa i dati recenti l'Istituto Tumori di Milano ha chiesto la bella cifra di 200.000 Euro.Potremmo concludere che la colpa è delle antenne?
La nostra analisi semplicistica può mettere in evidenza che esiste un effetto
che varia con la distanza dall'antenna ma questo non equivale a dire che vi sia
un rapporto causa-effetto fra radiazione E.M. e leucemia. E' però abbastanza ovvio che
una correlazione con la distanza sarebbe motivo sufficiente per ulteriori
investigazioni.
Un'analisi più rigorosa dovrebbe basarsi sulla dose di radiazione ricevuta dalla popolazione (la cosiddetta dosimetria). Di che si tratta?
In base al principio di causa ed effetto, i casi di malattia devono crescere con la quantità di energia E.M. a cui i malati sono esposti (che chiameremo "dose" mutuando un termine usato per la radiazione ionizzante). Senza questa assunzione non è ovviamente possibile fare della scienza sperimentale. Non vogliamo definire una forma della dipendenza malattia-dose se non limitandoci a dire che dev'essere monotona e crescente (più campo più dose più malati). Il problema è allora quello di misurare la dose per gli abitanti dell'area di Cesano e in particolare per i casi esaminati.
Tralasciando la possibilità di utilizzare dei misuratori di campo portatili su un campione di persone (sarebbe la tecnica corretta ma evidentemente dispendiosa e non sarebbe nemmeno semplice trovare dei volontari) una semplificazione ragionevole consiste nel misurare il campo nell'arco delle 24h e dei 365 gg nelle abitazioni che si trovano distribuite entro i cinque anelli esaminati. I casi di incidenza o di decesso andrebbero raggruppati quindi non più secondo la distanza dalle antenne ma in fasce di dose di radiazione ricevuta mediata sull'arco della giornata. E' chiaro comunque che questa tecnica è, a rigore, applicabile se le persone passano tutto il loro tempo nello stesso posto, un'ipotesi irrealistica eccetto forse per i bambini molto piccoli (e in ogni caso, anche per essi sarebbe valida su tempi limitati). Negli studi fatti (inclusa la nostra analisi didattica) si è quindi ricorsi, in manacanza di meglio, a una semplificazione ancora più cruda: dato che, a parità di potenza irraggiata dall'antenna, la dose di radiazione ricevuta decresce col crescere della distanza R, si è utilizzato la distanza R come surrogato della dose. Se ci fosse una correlazione perfetta sarebbe indifferente raggruppare la popolazione per dose o per distanza. Tuttavia una correlazione dose-distanza che persista sull'arco delle 24h e sui 365 gg dell'anno e sull'arco di un decennio è pochissimo plausibile in quanto
E' ovvio dal grafico che la correlazione fra densità di potenza e distanza è completamente inesistente. In siti fra loro vicini ci sono variazioni di oltre 100 volte! La ragione è dovuta all'orografia del terreno, ed agli effetti di schermatura e di riflessione degli edifici e al fatto che (come detto prima) l'emissione non è omnidirezionale. Addirittura, la dipendenza 1/R2 non è verificata neanche misurando lungo una stessa direzione.
Un'analisi più rigorosa dovrebbe basarsi sulla dose di radiazione ricevuta dalla popolazione (la cosiddetta dosimetria). Di che si tratta?
In base al principio di causa ed effetto, i casi di malattia devono crescere con la quantità di energia E.M. a cui i malati sono esposti (che chiameremo "dose" mutuando un termine usato per la radiazione ionizzante). Senza questa assunzione non è ovviamente possibile fare della scienza sperimentale. Non vogliamo definire una forma della dipendenza malattia-dose se non limitandoci a dire che dev'essere monotona e crescente (più campo più dose più malati). Il problema è allora quello di misurare la dose per gli abitanti dell'area di Cesano e in particolare per i casi esaminati.
Tralasciando la possibilità di utilizzare dei misuratori di campo portatili su un campione di persone (sarebbe la tecnica corretta ma evidentemente dispendiosa e non sarebbe nemmeno semplice trovare dei volontari) una semplificazione ragionevole consiste nel misurare il campo nell'arco delle 24h e dei 365 gg nelle abitazioni che si trovano distribuite entro i cinque anelli esaminati. I casi di incidenza o di decesso andrebbero raggruppati quindi non più secondo la distanza dalle antenne ma in fasce di dose di radiazione ricevuta mediata sull'arco della giornata. E' chiaro comunque che questa tecnica è, a rigore, applicabile se le persone passano tutto il loro tempo nello stesso posto, un'ipotesi irrealistica eccetto forse per i bambini molto piccoli (e in ogni caso, anche per essi sarebbe valida su tempi limitati). Negli studi fatti (inclusa la nostra analisi didattica) si è quindi ricorsi, in manacanza di meglio, a una semplificazione ancora più cruda: dato che, a parità di potenza irraggiata dall'antenna, la dose di radiazione ricevuta decresce col crescere della distanza R, si è utilizzato la distanza R come surrogato della dose. Se ci fosse una correlazione perfetta sarebbe indifferente raggruppare la popolazione per dose o per distanza. Tuttavia una correlazione dose-distanza che persista sull'arco delle 24h e sui 365 gg dell'anno e sull'arco di un decennio è pochissimo plausibile in quanto
- le trasmissioni sono direttive e quindi l'energia irraggiata dipende considerevolmente dalla direzione
- la direzione cambia nel corso delle 24h
- vi sono schermature dovute all'orografia del terreno, alle riflessioni e agli edifici
Densità di potenza a 1530 kHz in funzione della distanza dall'antenna trasmittente. La correlazione è inesistente
E' ovvio dal grafico che la correlazione fra densità di potenza e distanza è completamente inesistente. In siti fra loro vicini ci sono variazioni di oltre 100 volte! La ragione è dovuta all'orografia del terreno, ed agli effetti di schermatura e di riflessione degli edifici e al fatto che (come detto prima) l'emissione non è omnidirezionale. Addirittura, la dipendenza 1/R2 non è verificata neanche misurando lungo una stessa direzione.
Qualcuno potrebbe proporre di
sopprimere il dato a 1,85 km perché "anomalo". Se però l'anomalia è solo nel
fatto che questo dato è molto diverso dagli altri per cui non ci piace,
l'operazione si chiama "truccare" i dati. Per scartare un punto dev'esserci un
motivo serio a priori (per esempio sapere che è stato mal misurato) ma non si
deve essere influenzati dal risultato a posteriori.
Purtroppo le misure fatte si fermano a 3 km. Sarebbe fondamentale arrivare fino ad almeno a 10 km in modo da coprire tutto il cerchio. Una misura che si può fare in maniera molto semplice dato che basta limitarla al trasmettitore a 1530 kHz. Vi è anche da dire che, non sapendo che una trasmittente dell'impianto dominava, gli autori dello studio assumevano come centro del cerchio di 10 km un "centro dell'impianto" teorico. Questo peggiorerebbe ancora la correlazione se esistesse. In assenza di correlazione è chiaro che il raggruppamento in anelli di distanza crescente dal "centro" teorico o dall'antenna a 4 torri è del tutto arbitrario. Questo sta a significare che un'eventuale evidenza di un effetto dipendente dalla distanza non avrebbe connessione con la stazione radio.
| Aggiornata | |||


