Reproductie KNMI homogenisatie

Bericht van: Ben Lankamp (Lelystad) , 06-10-2019 23:42 

In het voorjaar verscheen de publicatie 'Het raadsel van de verdwenen hittegolven', door Rob de Vos, Marcel Crok e.a. In navolging daarvan heb ik contact gehad met zowel de auteurs als het KNMI, om allereerst te kijken of de homogenisatie binnen acceptabele marges kon worden geproduceerd.

Ik heb de methode nagebouwd en uitgevoerd in het statistisch pakket R, net als het KNMI. Het technisch rapport beschrijft de methode grondig, maar een paar cruciale details hebben de auteurs van de hittegolven-publicatie via e-mail van het KNMI gekregen (staan ook in de publicatie). Daarmee kan ik vrijwel exact alles nabootsen.

Wat cijfers voor de periode 1901-1950:

Het gemiddelde verschil in maximumtemperatuur met de officiële homogenisatie is 0,04 graden. Een T-toets op de verschillen tussen de KNMI homogenisatie en mijn reproductie levert p=0.5690 op.

De figuren tonen de correcties per percentiel van mijn reproductie en de KNMI homogenisatie (figuur 16 uit Brandsma, 2016).

Gr. Ben

 



 

Spreidingsdiagram reproductie versies KNMI-reeks:

Bericht laatst bijgewerkt: 07-10-2019 00:07

Reproductie KNMI homogenisatie   ( 845)
Ben Lankamp (Lelystad) -- 06-10-2019 23:42
Re: Reproductie KNMI homogenisatie   ( 307)
Bart (Abcoude) -- 07-10-2019 09:39
Re: Reproductie KNMI homogenisatie   ( 250)
Ben Lankamp (Lelystad) -- 07-10-2019 11:11
Bestand met metingen Pagode - Stevenson - Operationeel   ( 236)
Bart (Abcoude) -- 07-10-2019 12:28
Verschillen Pagode - Stevensonhut - Operationeel nader onderzocht   ( 139)
Bart (Abcoude) -- 10-10-2019 17:33
Maximum correctie Tx in augustus, in plaats van juni   ( 148)
Paul (Odijk) -- 13-10-2019 10:09
Re: Maximum correctie Tx in augustus, in plaats van juni   ( 114)
Ben Lankamp (Lelystad) -- 13-10-2019 13:22
Re: Maximum correctie Tx in augustus, in plaats van juni   ( 133)
Paul (Odijk) -- 14-10-2019 11:31
Re: Maximum correctie Tx in augustus, in plaats van juni   ( 120)
Ben Lankamp (Lelystad) -- 14-10-2019 13:25
Gemiddelde gemeten maximumtemperatuur versus tabel   ( 125)
Ben Lankamp (Lelystad) -- 14-10-2019 13:45
Bonus: percentage dagen met een bepaalde correctie   ( 159)
Ben Lankamp (Lelystad) -- 14-10-2019 14:07
Re: Bonus: percentage dagen met een bepaalde correctie   ( 110)
Michiel (Bunnik) -- 16-10-2019 10:50
Je hebt gelijk, mijn grafiek was iets mis mee   ( 130)
Ben Lankamp (Lelystad) -- 16-10-2019 14:11
Brandsma ( 2016) gaf dat aan in zijn Meteorologica-artikel   ( 104)
sebastiaan (bussum) -- 16-10-2019 17:50
Het was nog niet af (?)   ( 88)
Bart (Abcoude) -- 18-10-2019 09:42
Het is nooit 'af'   ( 82)
Ben Lankamp (Lelystad) -- 18-10-2019 10:22
Grote complimenten Ben! Ik zie veel positieve reacties   ( 224)
Jan (Workum, FRL) -- 07-10-2019 20:08
Extra: waarom ik niet op 100% kom   ( 251)
Ben Lankamp (Lelystad) -- 08-10-2019 00:14
Percentielen, hoe werkt dat?   ( 198)
Bart (Abcoude) -- 08-10-2019 08:45
Interessante tweets   ( 175)
sebastiaan (bussum) -- 12-10-2019 11:00
visie op homogenisatie   ( 157)
Bart (Abcoude) -- 12-10-2019 18:47
Re: visie op homogenisatie   ( 163)
Ben Lankamp (Lelystad) -- 12-10-2019 20:48