Hoi Alwin,
Kijk, dit is leuk. Je regressionele index zit er (natuurlijk) goed uit, ik zie ook dat er een aardige overlap bestaat met de parameters die wij heb genomen (weliswaar anders gecombineerd). Ik zal zeker je index eens programmeren in GrADS, en dan gaan we eens kijken hoe beide indices het doen. Mss. kunnen bepaalde delen van beide indices met elkaar verweven worden, of als blijkt dat bepaalde factoren een overhand hebben tov. anderen het gewicht verhogen/verlagen van bepaalde parameters. Dit ziet er iig. zeer goed uit, top!
Mijn huidige visie over goeie componenten voor zo'n index is wel iets veranderd en mijn regressieformule (eigenlijk geen lineaire, maar een zogeheten logistieke regressie, vanwege die e-macht vorm) zou denk ik beter zijn als ik lapse rates meegenomen had i.p.v. temperatuurverschillen, zoals (T850-T500)/(Z500-Z850) i.p.v. T850 en T500 afzonderlijk.
Het is dus maar een hele simpele benadering en echt bedoeld voor de gein. Bovendien bevat de formule nauwelijks informatie over grootschalige forcering.
Ik heb twee vragen: ik ben benieuwd hoeveel Z1000 voor invloed heeft, immers deze waarde zal toch nooit zo groot zijn? (althans, voor NL)
De waarde van Z1000 varieerde van -212 (20-11-1996 6z) tot +336 (29-1-1997 12z) meter, dus zal een waarde tussen -0.8 en +0.5 kunnen toevoegen aan X. Het is natuurlijk zo dat al deze grootheden met elkaar verbonden zijn, zodat je geen onafhankelijke deelcomponenten hebt.
Vraag 2: heb je mss. die onweersdata die je hebt gebruikt voor je onderzoek nog? Ik zou dan een soortgelijke statistische verificatie kunnen maken voor onze index. Mail me desnoods even.
Ik mail je wel even.
Het script voor vochtconvergentie kende ik nog niet, hartelijk bedankt! Ook daar ga ik eens wat mee plotten en event. kijken of dat in onze index meegenomen kan worden
Ben
Graag gedaan!
Groet,
Alwin