Re : Maar hoe vaak komen deze afwijkingen nu voor?

Bericht van: Ben (Lelystad) , 19-08-2010 16:56 

Gezien voor de laatste 30 jaar? Zijn het niet outliers die aan de laatste 2,5 % van de normaalverdeling liggen? Of zijn deze afwijkingen significant, en komen ze meer voor dan volgens het toeval zou moeten kunnen?

Dat kan je heel snel toetsen. Het ging over Moskou, dus daar ga ik even vanuit. Juli 2010 had een maandgemiddelde van 26°C. Dit is een nieuw record, sinds het begin van de reeks in 1881.

We moeten eerst een periode definiëren waar we een normaalverdeling (Gaussische verdeling) op fitten. Ik heb gekozen voor 1951-1980, een periode van 30 jaar en redelijkerwijs voor de 'vermeende' grootste opwarming. De statistische eigenschappen van deze periode waren:

Gem: 18,28°C
Mediaan: 18,3°C
Stdev: 1,66°C
Max: 22,0°C
Min: 15,5°C

Als je een Gaussische verdeling maakt met gemiddelde 18,28 en standaardafwijking 1,66°C, je veronderstelt vervolgens dat de verdeling gelijk blijft, dan zou je 10% van de jaren (3 jaar in 30 jaar) een julitemperatuur van 16,2°C of lager en 10% een temperatuur van 20,4°C of hoger verwachten. Eens?

Nu kijken we naar 1981-2010, wederom een periode van 30 jaar. De waargenomen aantallen per percentiel (links) en verwacht (rechts):

1% of hoger: 29 / 30
2½% of hoger: 29 / 29
5% of hoger: 29 / 29
10% of hoger: 29 / 27

90% of hoger: 6 / 3
95% of hoger: 6 / 2
97½% of hoger: 3 / 1
99% of hoger: 2 / 0

De verdeling is dus scheef. Vooral aan de warme staart is sprake van oververtegenwoordiging. Dit beeld zie je in alle maanden minus één: juni.

Hieronder weergegeven de omgekeerde cumulatieve frequentieverdeling. Dat wil zeggen: je telt alle julimaanden met Tgem >= x, waar bij x = waarde behorende bij een gegeven percentiel uit de gefitte Gaussische verdeling (op de x-as uitgezet). De verwachte verdeling is dus de 'ideale' Gaussische verdeling en de verdeling die je ruwweg verwacht voor elke 30-jarige periode, als het klimaat niet verandert. Het verschil tussen verwachting en waarneming voor het >= 95% percentiel (6 versus 2) is hoogst significant.

Gr. Ben

Bericht laatst bijgewerkt: 19-08-2010 17:01

Moskovitische extremen   ( 1268)
Cees-Rotterdam -- 19-08-2010 10:00
Re : Moskovitische extremen   ( 709)
Rutger (Meppel) -- 19-08-2010 10:09
Re : Moskovitische extremen   ( 519)
Seppie (Buren GLD.) -- 19-08-2010 10:09
Re : Moskovitische extremen   ( 391)
Pascal (Gent) -- 19-08-2010 10:25
Re : Moskovitische extremen   ( 352)
Remko (Gouda) -- 19-08-2010 10:40
Re : Moskovitische extremen   ( 324)
Pascal (Gent) -- 19-08-2010 10:45
  Mooi! 🙂  
Remko (Gouda) -- 19-08-2010 10:59
  Re : 🙂  
Hans (Eindhoven) ( 17m) -- 19-08-2010 10:55
Re : Oooh, gelukkig maar 0,73   ( 263)
Cees-Rotterdam -- 19-08-2010 14:22
Re : Tweede decimaal   ( 163)
Ridwan (Maasland) -- 19-08-2010 15:19
Re : Tweede decimaal   ( 137)
Cees-Rotterdam -- 19-08-2010 15:37
Re : Tweede decimaal   ( 154)
Ben (Lelystad) ( 13m) -- 19-08-2010 18:17
Re : Tweede decimaal   ( 107)
Cees-Rotterdam -- 19-08-2010 21:42
Re : Tweede decimaal   ( 187)
Ben (Lelystad) ( 13m) -- 19-08-2010 22:45
Re : Tweede decimaal   ( 87)
Cees-Rotterdam -- 20-08-2010 09:25
Re : Tweede decimaal   ( 188)
Ben (Lelystad) ( 13m) -- 20-08-2010 11:46
Re : Tweede decimaal   ( 73)
Cees-Rotterdam -- 20-08-2010 13:06
Re : Tweede decimaal   ( 164)
Ben (Lelystad) ( 13m) -- 20-08-2010 13:27
Re : Tweede decimaal   ( 101)
Seppie (Buren GLD.) -- 20-08-2010 00:23
Re : Dank...   ( 85)
Cees-Rotterdam -- 20-08-2010 09:27
Re : Oooh, gelukkig maar 0,73   ( 218)
Seppie (Buren GLD.) -- 19-08-2010 15:01
Re : Oooh, gelukkig maar 0,73   ( 206)
Ben (Lelystad) ( 13m) -- 19-08-2010 15:02
Re : Oooh, gelukkig maar 0,73   ( 187)
Seppie (Buren GLD.) -- 19-08-2010 15:11
Re : Oooh, gelukkig maar 0,73   ( 172)
Ben (Lelystad) ( 13m) -- 19-08-2010 15:20
Wie zei er ook alweer..   ( 215)
Eric, Rotterdam -- 19-08-2010 15:31
Re : Wie zei er ook alweer..   ( 173)
Ben (Lelystad) ( 13m) -- 19-08-2010 15:38
Re : Reken maar eens uit   ( 185)
Cees-Rotterdam -- 19-08-2010 15:31
Re : Reken maar eens uit   ( 188)
Seppie (Buren GLD.) -- 19-08-2010 15:51
Re : Ik geef het op.   ( 134)
Cees-Rotterdam -- 19-08-2010 16:47
Re : Reken maar eens uit   ( 185)
Ben (Lelystad) ( 13m) -- 19-08-2010 16:02
En dat is het andere belangrijke argument.   ( 178)
Remko (Gouda) -- 19-08-2010 16:18
Ho maar wacht...   ( 147)
Bart (O.-L.-V.-Waver) -- 19-08-2010 18:20
Re : Ho maar wacht...   ( 128)
Ben (Lelystad) ( 13m) -- 19-08-2010 18:43
Re : Ho maar wacht...   ( 75)
Remko (Gouda) -- 20-08-2010 09:19
Re : Reken maar eens uit   ( 123)
Dieter Lombaert (Asse) ( 220m) -- 19-08-2010 18:39
Re : Reken maar eens uit   ( 164)
Remko (Gouda) -- 19-08-2010 16:09
Maar hoe vaak komen deze afwijkingen nu voor?   ( 143)
Marcel (Den Haag) -- 19-08-2010 16:29
Re : Maar hoe vaak komen deze afwijkingen nu voor?   ( 127)
Remko (Gouda) -- 19-08-2010 16:46
Re : Maar hoe vaak komen deze afwijkingen nu voor?   ( 140)
Ben (Lelystad) ( 13m) -- 19-08-2010 16:56
Re : Maar hoe vaak komen deze afwijkingen nu voor?   ( 119)
Remko (Gouda) -- 19-08-2010 17:18
Lijkt me helder.   ( 120)
Marcel (Den Haag) -- 19-08-2010 17:30
Re : Gecompenseerd?   ( 150)
Cees-Rotterdam -- 19-08-2010 15:21
Re : Moskovitische extremen   ( 466)
Ben (Lelystad) ( 13m) -- 19-08-2010 10:16
Re : Moskovitische extremen   ( 401)
Jan (Workum, FRL) -- 19-08-2010 10:22
🙂 🙂   ( 225)
Jorrit (Amsterdam / Diever) -- 19-08-2010 10:27
Re : Moskovitische extremen   ( 494)
Rutger (Meppel) -- 19-08-2010 10:17
Re : jaarlijkse gang ca. 65% groter   ( 472)
Ben (Lelystad) ( 13m) -- 19-08-2010 11:02
Re : jaarlijkse gang ca. 65% groter   ( 302)
Lennart (Ottersum) ( 14m) -- 19-08-2010 11:11
Re : jaarlijkse gang ca. 65% groter   ( 295)
Ben (Lelystad) ( 13m) -- 19-08-2010 11:27