Dispense

238 CAPITOLO 4. ANALISI E PREVISIONE DI SERIE STORICHE
Intanto commentiamo i comandi. Il comando par(mfrow=c(h,k)) crea …nestre gra…che
con h righe e k colonne di sotto-disegni. Attende che gli vengano dati h k comandi di plot,
che metterà nella griglia. Attenzione che il comando resta inserito, per così dire, anche nelle
elaborazioni successive, per cui se nella …nestra gra…ca successiva si vuole ad esempio un plot
singolo, lo si deve richiedere esplicitamente con par(mfrow=c(1,1)).
Il comando window(IT4,...) ritaglia una …nestra, quella relativa alle speci…che scritte
nel comando. Si veda ?window per ulteriori dettagli.
La speci…ca type=b (si veda ?plot per tutte le alterntive possibili) disegna il gr…co con
both lines and points (altrimenti con type=l si ottenevano solo lines e con type=p si ottenevano
solo points). Quando si vuole vedere un gr…co al mocroscopio conviene usare
entrambi perché danno entrambi informazioni: le linee danno la visione d’insieme, i punti
permettono di vedere i singoli valori (si noti che a volte due linee successive di un gra…co a
linee sono allineate o quasi e quindi non si percepisce dove sta il punto si suddivisione, se non
lo si disegna esplicitamente).
Fatti questi commenti sul software, torniamo all’analisi della serie. I due anni, 2004 e
2005, mostrano assai poco in comune. Hanno un generico trend crescente, ma questo è solo
il sottoprodotto del trend generale crescente, non è una stagionalità. Forse gli ultimi 3 valori
annuali hanno un pattern simile, ma è poca cosa (volendo la si potrebbe estrapolare). Si
provino altri confronti, per capire di più ad esempio
ts.plot(window(IT4, start=c(2003,1), end=c(2003,12)),type=b)
ts.plot(window(IT4, start=c(2005,1), end=c(2005,12)),type=b)
da cui emerge che agosto è il mese peggiore (così non è stato però nel 2004!), mentre
ad esempio il pattern degli ultimi tre valori si è perso, quindi forse non vale la pena di
evidenziarlo. Si facciano altre prove. E’essenziale conoscere la serie che si studia nei minimi
dettagli.