fdbk_dt_hits_uncert: Calculates stratified hit rates for uncertain obs/fcst

In rfxf/Rfdbk: Handling NetCDF feedback files

Description

Calculates stratified hit rates for uncertain obs/fcst

Usage

fdbk_dt_hits_uncert(DT, thrs, by, cores = 1, incores = 1)

Arguments

DT	data.table with relevant information (at least varno, obs and veri_data)
thrs	list of variable having each a list of lower/upper limit, relative to observation
by	stratify contingency entries by these DT columns
cores	computing cores for the outer loop (splits computation by varnos)
incores	computing cores for the outer loop (splits computation by thresholds)(available cores have to be of number cores x incores)

Value

data.table with columns varno,interval, hits, total and the arguments of 'by'

Author(s)

Felix <[e-mail]>

Examples

#EXAMPLE (CSI for quantile thresholds) 
require(ggplot2)
fnames                = system("ls ~/examplesRfdbk/*/synop/verSYNOP.*",intern=T)
cond                  = list(veri_description="grepl('forecast',veri_description)",
                             veri_forecast_time="veri_forecast_time%in%c(2400,4800,7200,9600,12000,14400,16800)")
columnnames           = c("obs","veri_data","varno","veri_model","veri_forecast_time","statid","veri_initial_date")
DT                    = fdbk_dt_multi_large(fnames,cond,columnnames,20)
thrs                  = list('29'=list('lower'=c(-1/6),'upper'=c(1/6)),
                              '3'=list('lower'=c(-1,-2),'upper'=c(1,2)))
xx                    = fdbk_dt_hits_uncert(DT,thrs,by=c("veri_model","veri_forecast_time","varno"),cores=2,incores=3)
PEC                   = xx[,list(PEC =(hit)/(total) ),by=c("veri_forecast_time","veri_model","interval","varno")]
PEC[,varno:=varno_to_name(varno,T)]
ggplot(PEC,aes(x=veri_forecast_time,y=PEC,group=interaction(veri_model,interval),linetype=veri_model,color=interval))+
geom_line()+
geom_point()+
facet_grid(~varno)+
theme_bw()+
ggtitle("Percent Correct (within interval)")

rfxf/Rfdbk documentation built on May 27, 2019, 7:22 a.m.