Ανάλυση χρονο-σειρών ( Time-series analysis )

1. Ανάλυση χρονο-σειρών(Time-series analysis) H. Isliker, 2006

2. Περιεχόμενο • Εισαγωγή:παραδείγματα ΧΣ, γενικός σκοπός και χρησιμότητα, βασικές έννοιες • Απλές μέθοδοι περιγραφής χρονο-σειρών • Μέθοδοι ανάλυσης χρονο-σειρών:αυτο-συσχέτιση (auto-correlation)φασματικήανάλυση (spectral analysis) • Εφαρμογές και ασκήσεις με Mathematica

3. Βιβλιογραφία • C. Chatfield, The Analysis of Time Series,An Introduction, 4th ed., Chapman & Hall,London, 1989 • The Numerical Recipes(free downloadable from the net) • Σημειώσεις:www.astro.auth.gr→ people → H. Isliker → courses

4. Τί είναι οι χρονο-σείρες σήμα (signal) X(t) time t Χρονο-σειρά:=συλλογή από παρατηρήσεις που έγιναν διαδοχικά στο χρόνο Η ανάλυση χρονο-σειρών ανήκει στη στατιστική

5. Παραδείγματα: Οικονομία • π.χ. η εξέλιξη τιμής του πετρελαίου, η αξία μιας μετοχής (εφημερίδες για την οικονομία !) Σοφοκλέους, μετοχή μιας εταιρίας Ευρώ Κύριο ενδιαφέρον: πρόβλεψη !!!

6. Παραδείγματα: Εμπόριο (marketing) • π.χ. πωλήσεις ανά μήνα • σκοπός 1: πρόβλεψη, αλλά και σχεδιασμός • σκοπός 2:σύγκριση, π.χ. δαπάνες για διαφημίσεις και αυξήσεις στις πωλήσειςσύγκριση 2 χρονο-σειρών εξήγηση: Η μια ΧΣ εξηγεί την άλλη διαφημίσεις πωλήσεις ευρώ ευρώ t t

7. Παραδείγματα • Δημογραφία: εξέλιξη πληθυσμού,ενδιαφέρον:πρόβλεψη • Έλεγχος διαδικασιών παραγωγής,ποιότητα ενός προιόντος,π.χ.πάχος μιας βίδας Προδιαγραφή, επιδιωκόμενη τιμή Τιμή μεταβλητής ελέγχου t

8. Παραδείγματα: Ιατρική • Ηλεκτρο-Εγκεφαλογράφημα (ΗΕΓ, EEG) • Ηλεκτρο-Καρδιογράφημα (ΗΚΓ, ΕCG) • Σκοποί:1. κυρίως: αναγνώριση ανωμαλιών 2. κατανόηση και εξήγηση της δυναμικής της καρδιάς, του εγκεφάλου3. πρόβλεψη διαταραχών …

9. Παραδείγματα:Φυσική • Οπου γίνονται παρατηρήσεις ... • Σεισμολογία (κατανόηση και πρόβλεψη) • Μετεωρολογία (πρόβλεψη) • Αστροφυσική • …

10. Παραδείγματα:Ατμοσφαιρική φυσική ταχύτητα ανέμου(διάνυσμα !) θερμοκρασία σκοπός: πρόβλεψη (Παν/μιο Αθηνών, 2004)

11. Παραδείγματα: Αστροφυσική Ηλιακή έκλαμψη (solar flare) Παρατήρηση σε Ράδιο-συχνότητες (~300MHz) Σκοπός: • περιγραφή για σύγκριση με θεωρητικά μοντέλα, • αλλά και κατανόηση των διαδικασιών που παράγουν τις ΧΣ, καθοδήγηση στην δημιουργία μοντέλων • Τα επιτυχή μοντέλα χρησιμοποιούνται μετά και γιαπροβλέψεις

12. Σκοποί της ΑΧΣ • Περιγραφή (description)μιας ΧΣ • Εξήγηση (explanation) (μια ΧΣ μια άλλη) • Πρόγνωση (prediction) • Έλεγχος (control) • Κατανόησητης διαδικασίας

13. Σκοπός 1: περιγραφή • Πρώτα κάνουμε πάντα τη γραφική παράσταση μιας χρονο-σειράς • Απλή περιγραφήπεριοδικότητα (seasonal effects), τάση (trend), αλλαγή τάσης (με το μάτι, και αριθμητικά) • Προχωρημένη περιγραφή μεστοχαστικά μοντέλα (στοχαστικές διαδικασίες), π.χ. ‘η ΧΣ παριστάνει λευκό θόρυβο’

14. Σκοπός 2: Εξήγηση • Έχουμε 2 ή περισσότερες χρονο-σειρές • Μπορεί η μια να εξηγεί την άλλη? • Χ(t) Y(t) = f(X(t)),X(t): input, Y(t): output, π.χ. γραμμικό σύστημα (linear system) X(t) διαφημίσεις πωλήσεις Y(t) ευρώ ευρώ X(t) Y(t) f

15. Σκοπός 3: Πρόγνωση • Από μια ΧΣ μπορούμε να προβλέψουμε το μέλλον, και με πόση ακρίβεια ?

16. Σκοπός 5: Κατανόηση • Η στατιστική δίνει φορμαλιστικές περιγραφές των ΧΣ καθεαυτών • Στη φυσική μπορεί να έχουμε μια ΧΣ από ένα κατά τα άλλα άγνωστο σύστημα.Τί μπορούμε να καταλάβουμε από την ΧΣ για το ίδιο το σύστημα, το οποίο έχει παράγει την ΧΣ? • Πχ. ‘το σύστημα είναι περιοδικό, με περίοδο ...’,ή ‘το σύστημα είναι εντελώς στοχαστικό’ή ‘το σύστημα είναι χαοτικό, με fractal διάσταση ...’(→Μη-γραμμική ανάλυση χρονο-σειρών:fractal διαστάσεις, εκθέτες Lyapounov)

17. Βασικές Έννοιες • Συνεχής (continuous)ΧΣX(t) παρατηρείται συνεχώς (π.χ. kαταγράφηση σε ‘άπειρο’ χαρτί – βλ. σεισμολογία) • διακριτή (discrete) ΧΣ(π.χ. αν είναι σε ηλεκτρονική μορφή)sampling time t, ti = i t X(t) t X(ti), i=1,2,3,… t t1 t2 ti

18. Στον Η/Υ βολεύουν οι διακριτές ΧΣ(πρόβλημα αποθήκευσης !) Sampling (read off, digitize): διαβάζουμε και κρατάμε από συνεχή ΧΣ τιμές μόνο σε σημεία με σταθερή χρονική απόσταση Δt (sampling time/interval) ή μετράμε εξ’αρχής μόνο σε διακριτές χρονικές στιγμές X(t) t X(ti) Αποφεύγουμε μη-σταθερό Δt, οι περισσότερες μέθοδοι γίνονται πιο δύσκολες ! t t t1 t2 t3

19. Σε άλλες περιπτώσεις, το X(ti) είναι άθροισμα ή ολοκλήρωμα για όλο το t • π.χ. X(ti) βροχή ανά m2ανά μέρα,άθροισμα όλων των βροχοπτώσεωνγια κάποιο διάστημα • π.χ. X(ti) ράδιο-ακτινοβολία από τον ήλιο • όπου f(t) η συνεχής ακτινοβολία, καιti-ti-1=t

20. Ειδικά στην στατιστική θεωρία για την ανάλυση χρονο-σειρών:συνήθως, οι διαδοχικέςπαρατηρήσεις δεν είναι ανεξάρτητες, άρα πρέπει να λάβουμευπ’όψιν μας τη σειρά των παρατηρήσεων • Ακριβώς αυτή η εξάρτηση επιτρέπει τηνπρόγνωση του μέλλοντος με βάση το παρελθόν • Ορισμός:Ντετερμινιστική ΧΣ: επιτρέπει πρόγνωση με ακρίβειαΣτοχαστική ΧΣ: επιτρέπει προβλέψεις μόνο εν μέρει, ‘με πιθανότητα p θα συμβεί Α, ...’

21. Δυο βασικές προσεγγίσεις για την ΑΧΣ • Time-domain:οι μέθοδοι είναι συναρτήσεις του χρόνου,π.χ. αυτο-συσχέτιση (auto-correlation) • Frequency domain:κάνουμε μετασχηματισμό Fourier και δουλεύουμε στο χώρο των συχνοτήτων,π.χ. φασματική πυκνότητα (spectral density)

22. Απλές περιγραφικές μέθοδοι • γραφική παράσταση !!!!! • μέσος όρος μ • διασπορά2 • τρέχοντας μέσος όρος • στασιμότητα (stationarity) • ανάλυση: περιοδικότητα, τάσεις, θόρυβος • φιλτράρισμα (filtering)

23. Παράδειγμα X(t) t • Κάπως περιοδικό • Έχει θόρυβο Ιδιότητες?

24. Μέσος όρος μ Εστώ η ΧΣ Ορισμός του μέσου όρου δηλ. η μέση τιμή όλων των τιμών της ΧΣ, ο ορισμός είναι όπως συνήθως, και η σειρά των X(ti) δεν παίζει ρόλο ! Στο παράδειγμα:  =-0.77

25. Τί σημαίνει =-0.77 ? Plot it ! Μέσος όρος Συμπίπτει με την εκτίμησή μας …

26. Διασπορά 2 Εστώ ΧΣ Ορισμός της διασποράς Μέσος όρος των αποκλίσεων από τη μέση τιμή στο τετράγωνο X(ti) X(ti)-  t

27. Στάνταρτ απόκλιση • Η διασπορά έχει μαθηματικάπλεονεκτήματα, δηλ. στην στατιστική θεωρία • Πιο διαισθητική είναι η στάνταρτ απόκλιση (standard deviation)μέση απόκλιση από την μέση τιμή

28. =-0.77 Στο παράδειγμα: 2 = 54.10, = 7.36, +   -  Μεταξύ -και + βρίσκονται τα περισσότερα σημεία της ΧΣ, Το διάστημα αυτό μας δίνει τη διακύμανση των τιμών της ΧΣ

29. Παράδειγμα 2 = 4.24, 2 = 59.50,  = 7.71 X(t) +  + t και §δεν δίνουν καλή περιγραφή της πραγματικής μέσης τιμής εδώ Λόγος: υπάρχει μια τάση (trend), δηλ. σαν να αλλάζει το πραγματικό στο χρόνο:  = (t)

30. Με διαίσθηση: (t) = a + t¢b, συγκεκριμένα (t) = t¢10/512 δηλ. κατάσκευάσαμε ένα μοντέλο για την τάση (trend), αλλά κάπως αυθαίρετα. Πώς μπορούμε όμωςνα το κάνουμε πιο συστηματικά ?

31. Τρέχοντας μέσος όρος (running mean, moving average) • Παράμετρος:K, μήκος του παραθύρου ti-K ti+K ti

32. Τρέχων μέσος όρος, Κ=40 ) Ο (τρέχοντας) μέσος όρος αλλάζει στο χρόνο

33. Τρέχοντας μέσος όρος, Κ=40 ) Ο (τρέχοντας) μέσος όρος αλλάζει στο χρόνο

34. Στασιμότητα (stationarity) • διαισθητικός ορισμός:μια ΧΣ είναι στάσιμη αν δεν υπάρχει συστηματική αλλαγή του μέσου όρου και της διασποράς στο χρόνο • π.χ. τάση )μη-στασιμότητα • η στασιμότητα είναι προυπόθεση για τα περισσότερα εργαλεία της ΑΧΣ (π.χ. αυτο-συσχέτιση, φασματική ανάλυση) • )χρειάζονται εργαλεία μετατροπής μη-στάσιμων σε στάσιμες ΧΣ !

35. Μέθοδος της ανάλυσης (decomposition) • βασικός σκοπός:να διαχωρίσουμε και να απομονώσουμε τα διάφορα χαρακτηριστικά μιας χρονοσειράς:κυρίως τα μη-στάσιμα από τα στάσιμα χαρακτηριστικά • Εργαλεία:εξομάλυνση (smoothing), προσαρμογή (fitting), αφαίρεση, παίρνουμε τη διαφορά …

36. συνυπάρχουν 3 • χαρακτηριστικα: • τάση, • περιοδικότητα, • θόρυβος Σκοπός:απομόνωσητου καθένα + + = περιοδικότητα τάση θόρυβος έτσι μπορούμε να χαρακτηρίσουμε πιο καθαρά το κάθε χαρακτηριστικό

37. Ανάλυση Ι: τάση μέθοδοι για ταυτοποίηση της τάσης: • γραμμικό φιλτράρισμα (smoothing, low-pass filtering) • μέθοδος διαφορών (differencing) • προσαρμογή (fitting)μιας συνάρτησης

38. Ανάλυση Ι.1: Γραμμικό φιλτράρισμα • Από την ΧΖ X(ti) σχήματίζουμε μια καινούργια ΧΣ Υ(t),όπου(για να έχουμε weighted (ανισοβαρή) μέσο όρο) • ενδιαφέρον έχουν 2 ΧΣ:1. Υ(ti), περιέχει την τάση2. περιέχει το ‘υπόλοιπο’ (residual)

39. γραμμικό φιλτράρισμα (i):π.χ. s=q=K, ar=1/(2K+1))τρέχοντας μέσος όρος πρόβλημα: περιθώρια (άκρες) (α) είτε με πρόσθεση μηδενικών ή της μέσης τιμής στις άκρες, (β) είτε

40. τρέχοντας μέσος όρος, Κ=40 με πρόσθεση μηδενικών στις άκρες (zero padding): χωρίς πρόσθεση μηδενικών στις άκρες, με προσαρμογή του συντελεστή Ai

41. X(ti) και Y(ti) Y(ti) τάση R(ti)=X(ti)-Y(ti) περιοδικότητα και θόρυβος

42. ΑΧΣ με Mathematica • άσκηση 1α: α) δημιουργία της ΧΣi = 1,2,3,…,512, ti = i t, t = 1u(ti): θόρυβος με ομοιόμορφη κατανομή στο διάστημα [-2,2]β) γραφική παράσταση ui = Random[ Real, {-2.,2.} ]

43. Άσκηση 1β: α) υπολόγισε τηνμέση τιμή μ, τη διασπορά σ2 και τη στάνταρτ απόκλιση σ της Χ(ti) β) γραφική παράσταση της Χ(ti), μ, μ+σ, μ-σ στην ίδια εικόνα