ESEMPI DI ANALISI STATISTICA APPLICATA: LA ‘CERCA’ DI UNICORNI: date di nascita e di morte, correlate alle fasi lunari e ad altro; del dott. Piero Pistoia; intermezzi pitture di Gabriella Scarciglia

CURRICULUM DI PIERO PISTOIA

piero-pistoia-curriculumok (#)

luna1

   DATA DI NASCITA E FASI

LUNARI

INTERMEZZI: LE TRE “LUNE” DI GABRIELLA SCARCIGLIA

STATLgraf

PREMESSA GENERALE A QUESTO POST A TAGLIO STATISTICO

Fasi guidate e TROUBLE del ricercatore che ‘tenta di guardare’ all’interno di un oggetto complesso, per ‘sbrogliare’ la matassa di un migliaio di dati relativi alle date di nascita e di morte in Val di Cecina, da classificare secondo vari criteri, onde trovare correlazioni, per es., con le fasi lunari, i sette giorni della settimana, i dodici mesi dell’anno, maschi e femmine…., attraverso ipotesi di lavoro successive.

LUNA_PREMESSA0001

Seguendo i percorsi tracciati nel link precedente, consideriamo, al termine del processo di ricerca, due soli gruppi di dati, il primo dei quali (CAMPIONE n.1) rappresenta le date di nascita  di persone decedute fra maggio 1976 e marzo 1992 (822 casi), considerando solo quelle comprese fra il 1880 ed il 1925 (risultate 741), venendo così a trascurare 81 casi. Questi dati sono stati ripresi dal prezioso Registro dei Defunti della Parrocchia di S. Giovanni Battista di Pomarance (Pi), messo a disposizione dell’autore dall’allora proposto Don Piero Burlacchini, che si ringrazia per la gentile preziosa collaborazione.

Il secondo Gruppo (CAMPIONE n. 2)  è costituito da 444 date di nascita (11980-1991) riprese dal Registro delle Vaccinazioni dell’USL della Val di Cecina. Un ringraziamento particolare va alla Signora Fulvia Gronchi Borghetti, ostetrica,  per avere ricopiato direttamente a mano dal detto Registro le date di nascita al tempo attuali (444 casi).

Ringraziamo altresì l’amico prof. Vasco Pineschi  per aver riportato in ‘bella copia’ i grafici del risultato.

I diversi gruppi di date scelti in successione per questa ricerca  vengono via via sottoposti a vari programmi di lettura e classificazione dati, in Qbasic (LUNFRE), scritti da Piero Pistoia, simili, per es., al contenuto del link  seguente:

NASCO2

Vi facciamo vedere, per esempio, anche il risultato ottenuto da uno di questi programmi, da sottoporre ad una analisi statistica successiva in cerca di correlazioni.

 

BREVE COMMENTO AL PROGRAMMA DATI1.BAS

Di tratta di una routine in QB per immettere contemporaneamente in tre file i dati relativi a nascita, morte , sesso, età e per richiamarli.

Due di questi file sono di tipo random ed uno è sequenziale. I tre file vengono aperti all’inizio e chiusi alla fine della fase di immissione o della fase richiamo dati. Per i file random è necessario specificare anche il nome ed il tipo delle variabile immesse.

DT0.DT → file random che contiene il giorno-mese-anno di morte, l’età, il sesso, il giorno-mese-anno di nascita. E un file BINARIO

ET.DT → file random che contiene l’anno di morte e l’età. E’ un file BINARIO.

TE.DT → file sequenziale che contiene l’anno di morte e l’età. E’ un file ASCII.

In tal modo di due file ET.DT e TE.DT contengono gli stessi dati, solo che uno è scritto in BINARIO e l’altro in ASCII.

E possibile aggiungere a questi due file anche il SESSO per ogni età.

ALGORITMO PER IL CALCOLO DELLA ETA’

(Subroutine 100 del programma accessibile con il link al termine del commento)

L’età viene calcolata prendendo i GIORNI GIULIANI (Subroutine 1000, vedere programma nel linK sotto), relativi ad ogni data di nascita e di morte, sottraendoli per ottenere i GIORNI DI DURATA DELLA VITA.

Conoscendo l’anno di nascita ed i giorni di vita, l’algoritmo procede secondo le fasi seguenti:

1 – Aumento di 1 l’anno di nascita e controllo se è BISESTILE  (vedere NOTA) con una ulteriore subroutine.

2 – Sottraggo in successione dai giorni di vita, 365 (o 366 se l’anno in oggetto risulta bisestile) e per ogni sottrazione sommo 1 in un contatore.

3 – Ad ogni sottrazione controllo se il numero dei giorni di vita restanti è inferiore a 365 (o 366), giorni dell’anno su cui stiamo lavorando. Finché tale numero risulta maggiore, si torna al punto 1, continuando a sommare 1 all’anno di nascita, passando quindi al punto 2 e 3.

4 – Quando si esce (numero di giorni di vita restanti (< di 365 o 366), nel contatore c’è un numero senza decimali (intero) equivalente agli ANNI DI VITA o ETA’.

5 – I giorni di vita che restano rappresentano la frazione di anno, che posso calcolare dividendo questo numero di giorni di vita per 365 (o 3666, se l’ultimo anno del conto è bisestile).

6 – Sommo infine alla frazione di anno, il numero di anni di vita, ottenendo l’ETA’ cercata, in anni e frazioni di anno.

NOTA: anni bisestili

 Anno Tropico -> è l’intervallo di tempo che trascorre fra due successivi passaggi del sole, considerato che si muova a velocità angolare costante, all’equinozio di primavera. L’anno tropico, che è poi l’anno del senso comune, è più breve dell’anno siderale (perché l’equinozio (gamma) si muove incontro al sole ed è più breve del tempo che il sole impiega a percorrere un arco di 50”. 26
360° -> 1 296 000″
Anno Siderale -> è il periodo di rivoluzione della terra attorno al sole ed è pari a 365.256360 giorni solari medi.
L’Anno Tropico = Anno Siderale * (1 – 50”.27 / 1 296 000″) =365.24220 giorni solari medi
L’Anno Tropico presenta alcuni inconvenienti per le frazioni di giorno di cui bisogna tener conto.
Si è definito così l’anno civile che inizia sempre col principio esatto di un giorno medio ed è composto di un numero intero di giorni medi. Si può così ottenere una buona approssimazione, per l’anno civile all’anno tropico, considerando il civile pari a 365 e un 1/4 giorni solari medi.
Il problema del calendario si risolve facendo uguale a 365 giorni la durata dell’anno civile ed inserendo ogni 4 anni un anno bisestile di 366 giorno.
Ma rispetto all’anno tropico ogni 400 anni si ha:
365 1/4 * 400 -> 365*4000 + 100 giorni medi , anziché:
365.24222 * 400 -> 365 * 400 + 96.88 giorni medi.
Quindi ogni 400 anni si devono sopprimere 3 giorni, se non si vuole che la data del 21 Marzo anticipi.superi l’equinozio di primavera.
Così nel 1582 si dovettero sopprimere 10 giorni (il giorno dopo il 4 ottobre 1582 fu così il 15 invece del 5).
D’allora in poi gli anni indicati da un numero terminante con 2 zero (che prima erano tutti bisestili) fossero bisestili solo se tale numero fosse anche divisibile per 400, effettuando così la soppressione di 3 giorni ogni 400 anni.

Il precedente commento è relativo al programma in basic di Piero Pistoia, accessibile con il link:

NASCITA_MORTE_ETA_DATI1

 

Da continuare…..

STATL1

OLYMPUS DIGITAL CAMERA
STATL12


OLYMPUS DIGITAL CAMERA
STATL13

PRIMA LUNA

marmel0004

STATL2

STATL21

STATL22i

STATL23

SECONDA LUNA

OLYMPUS DIGITAL CAMERA

STATL31i

STATL32

STATL33i

STATL34

TERZA LUNA
goz2

STATL35

NOTA PREMESSA