190 likes | 322 Views
Resources on forecast evaluation/verification. Barbara Brown. Joint Numerical Testbed. NCAR, Boulder, CO. July 2011. WMO Verification web page. Developed and maintained by. working group on. forecast verification research. Includes FAQs,. references and links to
E N D
Resourcesonforecast evaluation/verification BarbaraBrown JointNumericalTestbed NCAR,Boulder,CO July2011
WMOVerificationwebpage Developedand maintainedby workinggroupon forecastverification research IncludesFAQs, referencesandlinksto tutorialpresentations Verificationdiscussion group: vx- discuss@rap.ucar.edu http://www.cawcr.gov.au/projects/verification/
EUMETCALOnlineLearningModule Freeon-line tutorialwith exercises Includesmodules on Categorical Continuous Probabilistic Verification methods http://www.eumetcal.org/resources/ukmeteocal/verific ation/www/english/courses/msgcrs/index.htm
Books Jolliffeand Stephenson: Forecast Verification-A Practitioner’s Guidein Atmospheric Sciences PublishedbyWiley (newversionin 2011) Weatherand Climate Prediction Publishedby Cambridge UniverstityPress Includeschapter Warner(2010): Wilks(2011):Statistical MethodsintheAtmospheric Sciences PublishedbyElsevier Includesextensivechapteron verificationmethods Numerical
ForecastEvaluationTools RStatisticsPackage: Includesmanystatisticaltools, includingtheRVerification Library METisfreelyavailableand supportedtothecommunity Mainfocus:Model verification Includestoolsforpointand griddeddata,ensemble forecasts,spatialmethods http://www.dtcenter.org/met/users/
IntroductiontoR TaraJensen NCAR/RAL/JNT
WhatisR? •Astatisticalprogramminglanguage •Inpart,developedfromtheSProgramming LanguagefromBellLabs(JohnChambers) •Createdtoallowrapiddevelopmentof methodsforuseindifferenttypesofdata. •Createnewgraphics.Manydefault parametersarechosen,butusersretain completecontrol.
WhyR? •Rhasbecomethedominantlanguageinthestatistical researchcommunity. •RisOpenSourceandfree. •Runsonalloperatingsystems •Nearly2,400packagescontributed. •Packagesandapplicationsinnearlyeveryfieldof science,businessandeconomics. •SeeRNotes,RJournalandJournalofStatistical Software.www.jstatsoft.org •Morethan100bookswithaccompanyingcode •Verylarge,activeuserbase.
WhynotR? •NCL,IDL,Matlab,SAS,…areallviable alternativestoR.Ifyouareapartofanactive communityofresearchersusinganother language,dolikewise. •Ifwewerebiostatisticianswewouldbeusing SAS.BookTitle:“AnalyzingReceiverOperatingCharacteristicCurveswithSAS” •Considerbuildingverificationfunctionsand utilitiesaspartofcodedevelopment. Verificationneednotbeanexternalprocessto forecasting.
TheRCommunity •Developers –RCoreGroup(17members),only2haveleftsince 1997 –MajorupdateinApril/October(freezedates,beta versions,bugtracking,...) •Mailinglists –Helplist~150messages/day,archived, searchable. •5InternationalConferences,2US,1China
EverythingaboutRisatwww.r-project.org •Sourcecode •Binarycompilations(Windows,MacOS,Linux •Documentation(Maindocuments,plusnumerous contributed.Someinforeignlanguages.) •Newsletter(replacedbyRJournal.) •Mailinglist(Severalsearchengines) •Packagesoneverytopicimaginable •Wikiwithexamples •ReferencelistofbooksusingR.(morethan100) •TaskManager
UseRwithscripts •InLinux-EmacsSpeaksStatistics –Providessyntax-based –Objectnamecompletion –Keystrokeshortcuts –Commandhistory –Alt-xRtoinvokeRwithXemacs. •InWindows,useeditor –AddedGUIfeatures –<control>RsendsalineorhighlightedsectionintoR. –InstallpackagewithGUIs –Savegraphicsbypointandclick. •MacOS –SimilartoWindowswithadvantagesofsystemcalls.
Codingprinciples •Makeverificationcodetransparentandeasy toread •Commentanddocumentliberally •Archiveyourcode •Shareyourcode •Labelandsaveyourdata •Shareyourdata
PackagesinR •Contributedbypeopleworldwide. •Allowscientistsorstatisticianstopushtheir ideas. •ApplyandextendRcapabilitiestomeetthe needsofspecificcommunities. •Accompanymanystatisticaltextbooks •Accompanyappliedarticles(AdrianRaftery, DougNychka,TilmanGneiting,BarbaraCasati, MattBriggs)
Asampleofusefulpackages • • • • • •Rgis,spatstat(GIS) •ncdf(supportfor netcdffiles) •Rcolorbrewer •randomForests verification fields(spatialstats) radiosondes extRemes BMA(BayesianModel Averaging) •BMAensemble •circular •Rsqlite
Packages •Packagesmustbeinstalledtocall. •Packagesmustbecalledtouse. •Basepackagesareinstalledbydefault.
10mostusefulfunctioninR •aggregate-appliesafunctiontogroupsof datasubsetbycategories. •apply-incrediblyefficientinavoidingloops. Appliesfunctionsacrossdimensionsofarrays. •layout-creativelydivideaprintregion. •xyplot(inthelatticepackage)slightlyadvance graphictechniques •%in%returnslogicalshowingwhichelements inAareinB.(e.gA%in%B)
Moretop10 • • • • • table–createcontingencytabelcounts. boot–applybootstrapfunctioncorrectly read.fwf–readfixedwidthformatdata par–controleverythinginagraph system()–allowsyoutocallsystem commandfromR •pairs–themostunderutilizedplot–plotsa matrixof4columnsina4x4plotlayout
RExercises • • • • • Choosegroupsof3-4–findacomputer Logontomachines Bringupatleast2xterms >cp/wrfhelp/MET/R-packages/example/*. >viintro2R.2011comet.R or >kwriteintro2R.2011comet.R&