Interessante Beiträge :  Wetterzentrale Forum Wetterzentrale Forum

| Wetterzentrale-Startseite | Wetterzentrale-Chat | Wetterchronik | wetter-foto.de | Das alte WZ-Forum-Archiv (05/1999-04/2007) |

| Nutzungsbedingungen/Forumsregeln/Datenschutz | Die neuesten Postings aus allen Foren | Wettersymbole einbinden | Präfixe in Betreffs verwenden |

KLIMA: Wie prüft man auf Trends
geschrieben von: org: FokkerPlanck am 07.07.03 (IP-Adresse bekannt)
Datum: 14. Juli 2003 14:11

Im folgenden stelle ich eine kleine Anleitung zum Überprüfen von
Zeitreihen auf lineare Trends dar.
Ich habe also eine Zeitreihe mit einer Anzahl n von Messwerten xk (k=1 bis k=n), die immer in gleichen Zeitabständen aufgezeichnet wurden (z.B. Monatsmittelwerte der Temperatur aller Juli-Monate der Jahre 1905 bis 1966, oder den Luftdruck um 13.00 MEZ vom 28.2. bis 21.10.1871,...)
Überprüft werden soll nun, ob diese Zeitreihe {x} einen linearen Trend zeigt.
Wir setzen voraus, daß die Messwerte von allen systematischen Fehlern bereinigt sind (Instrumentenwechsel, Standortverlagerung,...).
Dann bildet man zunächst die üblichen Maßzahlen zur Stichprobenbeschreibung.
1. den Mittelwert m = (x1+x2+...+xn)/n
2. die Standardabweichung s, dazu bilden wir zunächst
a) die Quadrate aller Messwerte x1*x1, x2*x2, ... xn*xn und
b) summieren diese Quadrate SQ = x1^2 + x2^2 +... + xn^2.
c) Die sogenannte Varianz (oder Streuung) V = s^2 ergibt sich dann zu:
V = 1/(n-1) * [SQ - n*m^2] und damit
d) die Standardabweichung s = wurzel (V).
Diese beiden Werte m und s charakterisieren die Zeitreihe schon ziemlich gut.
Nicht nur der Mittelwert ist wichtig, sondern auch die Standardabweichung. Sie wird um so größer, je weiter die Messwerte vom Mittelwert abweichen. Zwei Zeitreihen können zwar denselben Mittelwert haben, aber bei der Zeitreihe mit der größeren Standardabweichung werden extremere Werte eher mal tatsächlich erreicht.
Aber was sagen diese Maßzahlen über die Zukunft der Zeitreihe?
Kann ich behaupten, daß alle zukünftigen Messwerte in einem bestimmten Band um den berechneten Mittelwert liegen?
Und wie breit wäre dieses Band?
Die Antwort lautet schlicht:
Über die Zukunft läßt sich keine Aussage treffen, solange man nicht weiß,
wie die Messwerte entstehen. Alleine aus dem Verlauf der Zeitreihe läßt sich ohne das Hinzuziehen zusätzlicher Annahmen oder Modelle nichts schließen!
Wenn das nicht möglich ist, was läßt sich sonst von der Zeitreihe
lernen?
Naja, nicht viel, wenn man nicht ein paar Annahmen macht. Eine häufige Annahme ist, daß die vorliegende Zeitreihe ein Ausschnitt oder eine Stichprobe aus einer Normalverteilung (allgemeiner einer Grundgesamtheit GG) ist. Diese Normalverteilung ist unbekannt, weil sich klimatische Messwerte auf Dauer verändern (Klimaänderung) und sich das Wetter sich nicht nochmal ablaufen lassen kann. Man kann die Zeitreihe aber überprüfen, ob sie eine Stichprobe
einer Normalverteilung sein könnte.
Wenn das der Fall ist, so wird der Mittelwert interessant. Der stimmt dann nämlich in gewissen Grenzen mit dem Mittelwert M dieser unbekannten Normalverteilung überein. Es ist dann: |m-M| < t * s/wurzel(n), d.h. die Differenz (der Fehler) zwischen Stichproben-Mittelwert m und Mittelwert der GG wird begrenzt durch der Term auf der rechten Seite. Der setzt sich aus der Standardabweichung s, der Anzahl der Messwerte n und einem zusätzlichen Faktor t zusammen. Wenn die Streuung der Messwerte also klein ist oder die Anzahl der Messwerte n groß, so kann man den Mittelwert der GG immer besser mit dem Stichproben-Mittelwert identifizieren.
Wenn zwei Zeitreihen mit jeweils n Werten vorliegen und jeweils ihre Maßzahlen
m1, m2 und s1, s2 berechnet sind, so kann ich auch entscheiden, ob sich die
beiden Zeitreihen signifikant unterscheiden oder ob sie nur zwei Stichproben aus derselben GG sind.
Das geht so: man bildet t1 = |m1-m2|*wurzel(n)/wurzel(s1^2+s2^2), f = 2n-2
und vergleicht dieses t1 mit dem oben schon erwähnten t. Bei t1Stichproben signifikant aus derselben GG.
Zu diesem t (aus dem t-Test):
das t wird größer je sicherer ich mit meiner Aussage sein will, daß
die vorliegende Stichprobe zu einer GG behört bzw. beide Stichproben
zur selben GG. Diese Sicherheit nennt man Signifikanz und wird üblicherweise
in Prozent angegeben. 100% Signifikanz heißt, das meine Aussage todsicher ist.
Die kann ich aber bei nur endlich vielen Messwerte nie erreichen. Typisch
für eine befriedigende Signifikanz sind 95% oder 90%. 95% Signifikanz heißt,
daß nur in 5% aller Fälle von Stichproben meine anhand der Überprüfung
gewonnene Aussage falsch ist. Wenn ich also sage, daß zwei
Stichproben-Mittelwerte mit 95%iger Signikanz gleich sind, so können
sie dennoch aus verschiedenen GGs stammen, aber eben nur selten.
Das t ist dummerweise noch von der Werteanzahl n abhängig
(über das f, den Freiheitsgrad). Für n > 30 ändert sich t aber nicht
mehr stark. Bei 95%-Signifikanz ist t=1.96, bei 90% 1.645, bei 99% 2.576.
3. Wie entscheide ich nun, ob eine Zeitreihe einen (linearen) Trend
aufweist?
Das geht so, man berechnet schrittweise den Korrelationskoeffizienten r
(m Mittelwert, s Standardabweichung der Zeitreihe)
a) bilde die Produkte aus Zeitpunkt * Messwert zu diesem Zeitpunkt: 1*x1, 2*x2,...n*xn
b) summiere die Produkte auf: ST = x1+2*x2+3*x3+...n*xn
c) berechne den Dividenden A = ST - m*n*(n+1)/2
d) berechne den Divisor B = s*(n-1)*wurzel(n*(n+1)/12))
e) der Korrelationskoeffizient r ist dann: r = A/B.
Wenn der signifikant positiv ist, so liegt ein positiver
Trend vor und umgekehrt. Wie entscheidet man, ob er das ist?
Wenn das Intervall zwischen r - t*(1-r*r)/wurzel(n-1)
und r + t*(1-r*r)/wurzel(n-1) nicht die Null
enthält, so liegt ein Trend mit der durch t gegegenen Signifikanz vor.
Das Intervall gibt analog zur Signifikanz des Mittelwertes an, wie
"echt" die Korrelation (der Trend) ist. Wenn das Intervall um r die Null
einschließt, so kann man nicht signifikant ausschließen, daß keine
Korrelation (Trend) vorliegt.
Wer meinen Ausführungen bis hierhin gefolgt ist, der kann sich schon mal
ein paar Beispiele zurechtlegen, die er überprüft haben will.
Ich bitte darum!
Fokker



Thema Klicks geschrieben von Datum/Zeit
  Berichte zu Klima / Witterung (oT) 26445 Admin 09.10.08 08:20
  Hat El Niño oder La Niña einen Einfluss auf Europa? (lang) / Org: Lobo 20165 Admin 24.01.08 12:50
  Die Witterung in Deutschland im Dezember und Jahr 2007 (Karten & Beschreibungen) Org: Jörg, Berlin-Malchow (53 m) 6860 Admin 24.01.08 12:27
  Alle Temperaturabweichungen in Deutschland ab 1901 (Monate, Jahreszeiten, Jahre), Org: Badischer Reisbauer, 30.06.2007 8409 Admin 20.09.07 11:24
  KLIMA: Wahrhaft warme Zeiten - einige Diagramme dazu - Org: Astronom am 17.04.2007 10245 Admin 29.06.07 10:19
  KLIMA: Deutschland seit 2000 75 % "zu warme" Monate - Org.: Badischer Reisbauer 17306 Admin 29.06.07 10:09
  VERIFIKATION: Die wilden 70er und unser aktuelles Wetter 9960 derLüneburgerHeide 06.03.07 12:10
  Das "richtige" Klimamittel: 30 Jahre? Oder 10, 40? 6704 Astronom 06.03.07 11:53
  27.11.2006 - Zunahme der Großwetterlage Südwest 6284 Simon (Karlsruhe) 20.12.06 10:14
  KLIMA: Der extreme Winter 1709 10936 Astronom 20.12.06 10:11
  Milde 1. Januardekade = Winter ade? Studie dazu 6052 Astronom 02.01.05 17:52
  KLIMA: Zusammenhang Januar - letzte Dezemberdekade 6249 Astronom 20.12.06 09:56
  KLIMA: Dichtung und Wahrheit 6076 Astronom 20.12.06 09:43
  KLIMA: Größten Temp-Abweichungen in Deutschland 6379 Sebastian, Bielefeld/Uni, 124m 20.12.06 09:05
  Die größten negativen Temp-Abweichungen in D. 7312 Sebastian, Bielefeld/Uni, 124m 20.12.06 09:08
  130 Jahre Hitze, Kälte und Regen auf einen Blick 7971 Astronom 21.11.06 12:00
  Die homogene 300jährige Berliner Klimareihe(lang!) 5913 Berlinjörg 17.08.06 12:12
  Das Rätsel der zweiten Augustdekade 6362 Astronom 10.08.05 10:51
  Sommer zu kalt, zu warm? Hundstage, ja, nein? 6692 Astronom 24.07.05 22:06
  Märchenstunde! Heute: Die Siebenschläferregel 6129 Astronom 09.07.05 22:38
  Tägl. Extremtemperaturen letzte 200 Jahre (Europa) 10035 Astronom 12.04.05 10:09
  KLIMA: Rätselhafte Korrelationen (lang!) 7313 Astronom 26.10.04 15:49
  Der Gruselsommer 1980 - ein Rückblick (lang!) 17021 Astronom am 10.07.2004 13.07.04 10:20
  Ein Langzeit-Test der Siebenschläfer-Regel 5953 Astronom am 02.07.2004 06.07.04 07:12
  Das Märchen von den Singularitäten 5761 Astronom am 25.04.2004 29.04.04 21:34
  KLIMA: Analyse Temperaturverteilung Deutschland 5660 FokkerPlanck, am 08.03.04 28.03.04 22:54
  KLIMA: Nochmal zur Erläuterung: 4241 FokkerPlanck am 09.03.04 28.03.04 23:05
  KLIMA: Jetzt wird's spannend: Ursachenforschung 4212 FokkerPlanck, am 09.03.04 28.03.04 23:08
  Vulkanausbrüche und Klimaschwankungen 4933 Astronom am 22.03.04 28.03.04 22:57
  Vortrag in Plön: Das sagte M.Latif zum Klimawandel 4932 Wetterfuchs, am 06.02.04 11.02.04 21:54
  2003 -Tornado-Rekord in Deutschland? Eine Analyse. 5101 Matthias, Bonn am, 22.12.03 22.12.03 22:27
  Gletscher-Vergleich Schweiz und Norwegen (Bilder) 9995 Christoph (AN-Elpersdorf 501m-Mittelf 22.09.03 22:27
  KLIMA: Wie prüft man auf Trends 5336 org: FokkerPlanck am 07.07.03 14.07.03 14:11
  Jahrhundertchronik 1900-1999 7956 org: W.v.J, 30.04.2003 30.04.03 13:53
  KLIMA: Winterbegriffsdefinitionen??? 4959 org:Peter von Zug 35km südl.von Züric 24.10.02 19:03
  KLIMA: Winteraus- und Rückblick (lang!!) 9065 Christian (Braunschweig) 15.09.02 23:10
  Analyse der Regentage in Karlsruhe von 1888 - 2007 (Astronom) 4381 Admin 23.04.09 10:38
  Schweizermesser Ursache für zu hohe Niederschläge 1869-1887 geschrieben von: Badischer Reisbauer 7475 Admin 30.11.09 13:45
  Das Problem langer Klimareihen - neue Wege zur Wahrheit?(lang!) Org.: Jörg, Berlin-Malchow; Mai 2009 6718 Admin 27.11.09 18:26
  Potsdamer Sonnenschein-Diagramm seit 1893... geschrieben von: Jörg, Berlin-Malchow (September 2008) 6780 Admin 30.11.09 13:47
  Und noch was: Bedeckungsgrad / heitere+trübe Tage geschrieben von: Jörg, Berlin-Malchow 5922 Admin 30.11.09 13:48
  Klimawandel anschaulich: 2 lange Temperaturreihen aus D geschrieben von: Jörg, Berlin-Malchow (Februar 2008) 7776 Admin 30.11.09 13:50
  500 Jahre Eiswinter im Ostseeraum: Gibt es Hoffnung für Winterfans? geschrieben von: Jörg, Berlin-Malchow (Februar 2008) 6300 Admin 30.11.09 13:52
  Kältesummen ab 1750 von Basel, Berlin, Prag usw.! geschrieben von: Jörg, Berlin-Malchow (November 2008) 7886 Admin 30.11.09 13:54
  Abweichungen des Deutschen Flächenmittels vom Mittel 1961-1990 und 1981-2010; geschrieben von: Christoph (Halle/Saale) (März 2015) 3246 Admin 28.03.15 21:09


Du hast nicht die erforderliche Berechtigung, um in diesem Forum zu schreiben.