A Hierarchical Bayes Model for Combining Precipitation Measurements from Different Sources
|
|
- Justin Lindsey
- 5 years ago
- Views:
Transcription
1 A Herarchcal Bayes Model for Combnng Precpaon Measuremens from Dfferen Sources Tamre P. Cardoso and Peer Guorp NRCSE T e c h n c a l R e p o r S e r e s NRCSE-TRS No. 8 March, 7
2 A Herarchcal Bayes Model for Combnng Precpaon Measuremens from Dfferen Sources Tamre P. Cardoso and Peer Guorp ABSTRACT Surface ran rae s an mporan clmac varable and many enes are neresed n obanng accurae ran rae esmaes. Ran rae, however, canno be measured drecly by currenly avalable nsrumenaon. A herarchcal Bayes model s used as he framework for esmang ran rae parameers hrough me, condonal on observaons from mulple nsrumens such as ran gauges, ground radars, and dsromeers. The herarchcal model ncorporaes relaonshps beween physcal ranfall processes and colleced daa. A key feaure of hs model s he evoluon of drop-sze dsrbuons (DSD) as a hdden process. An unobserved DSD s modeled as wo ndependen componen processes: ) an AR() me-varyng mean wh GARCH errors for he oal number of drops evolvng hrough me, and ) a me-varyng lognormal dsrbuon for he sze of drops. From he modeled DSDs, precpaon parameers of neres, ncludng ran rae, are calculaed along wh assocaed uncerany. Ths model formulaon devaes from he common noon of ran gauges as ground ruh ; raher, nformaon from he varous precpaon measuremens s ncorporaed no he parameer esmaes and he esmae of he hdden process. The model s mplemened usng Markov chan Mone Carlo mehods. INTRODUCTION Surface ranfall s an mporan envronmenal varable ha s ncorporaed across many areas of sudy, ncludng meeorology, clmaology, agrculure, land use, and hydrology, for example. The varous felds of sudy requre esmaon of precpaon a a range of emporal and spaal scales. Hydrologss and land use planners may be neresed n shor-erm ranfall n relavely small regonal areas for sudes nvolvng flood and flash flood forecasng. Researchers n clmaology or agrculure may be neresed n clmac sudes ha focus on weekly, monhly, or annual oals over large spaal exens. Alhough of consderable neres, surface precpaon raes and amouns are dffcul o esmae. Unlke oher amospherc varables, ranfall can dsplay exreme heerogeney n space and me. Boh he occurrence of precpaon and he rae a whch falls may be hghly varable, even whn a sngle ran even. As such, due o he mporance of and dffcully wh esmaon, here has been much research n he area of precpaon measuremen and esmaon. Precpaon s measured usng many dfferen nsrumens. Some nsrumenaon measures ranfall drecly, whle ohers make ndrec measuremens of quanes ha can be relaed o ranfall. All of he nsrumenaon, however, make ndrec measuremens of he usual
3 quany of neres, namely, ran rae. The mos common nsrumens nclude ran gauges and ground-based scannng radar. Dsromeer n su measuremens and saelle deployed nsrumenaon are also used, alhough less commonly. Each of he nsrumens has nheren srenghs and weaknesses, whch affec he observed measuremens and subsequen esmaon of surface ranfall. There exs many algorhms, each wh s own se of assumpons, for calculang ranfall relaed parameers from he nsrumenal observaons. Many mes he observaons obaned from one nsrumen are used n he calculaon of esmaes based on he observaons from anoher nsrumen. Ran gauge daa are ofen consdered o be ground ruh and are used o adjus he esmaes based on oher nsrumens. Ths s done despe he known fac ha observaons from ran gauges are nherenly prone o errors and ha hey are no more rue han he observaons from oher nsrumens. Surface precpaon can be defned by a populaon of fallng drops. The dsrbuon of he sze dameers of he drops characerzes he populaon behavor; hus, he drop-sze dsrbuon (DSD) forms a basc descrpor n he modelng of ran mcrophyscs. The DSD s used o compue a varey of derved properes of ranfall hrough mahemacal relaonshps; common compued properes nclude waer conen, ran rae, and reflecvy. DSD daa, however, are no colleced on a roune bass due o he expense of dsromeer nsrumens and he lmed spaal coverage of an ndvdual nsrumen. Ran gauge measuremens and ground radar mages are much more readly avalable. The leraure s full of modelng approaches for he esmaon of surface ranfall, ncludng emprcal sascal models (Sern and Coe 984; Skaugen e al. 996, for example) and a varey of sochasc models (Bell 987; Rodrguez-Iurbe and Eagleson 987; Cox e al. 988; Rodrquez-Iurbe e al. 988; Smh 993, for example). More recenly, herarchcal Bayesan models have been used n a number of envronmenal applcaons (Royle e al. 998, Wkle 998, Berlner, Berlner e al. a, Wkle e al., Hrafnkelsson 3). In hs arcle we presen a generalzed model for surface ranfall ha accommodaes daa from mulple nsrumens and produces esmaes of precpaon parameers and assocaed unceranes. Usng a herarchcal Bayesan approach, he evoluon of DSDs are modeled as an unobserved process, whle also ncorporang nformaon from daa gahered hrough commonly deployed nsrumens. Ths formulaon does no depend on havng expensve DSD daa and depars from
4 he noon of gauge daa as beng ground ruh. Raher, an esmae of he drop-sze dsrbuon a each pon n me, wh aendan uncerany measures, can be used o calbrae all he nsrumens smulaneously. HIERARCHICAL MODEL Surface ranfall parameers hrough me are esmaed usng a flexble fve-sage herarchcal Bayesan formulaon ha s based on a model for generang ranfall drop-sze dsrbuons. Esmaes of he unobserved process, denoed by N(D), are he prmary quany of neres. The frs sage of he herarchy specfes a measuremen error model for he observaonal daa, denoed by and G, boh of whch are observaons of funcons of N(D), wh error. The second sage of he herarchy allows for a me seres formulaon of he unobserved DSD process. Tme seres parameers and emporal dynamc erms of he DSD evoluon process are modeled n he hrd sage. The herarchcal Bayes formulaon s compleed n sages four and fve by specfyng prors on model parameers. Emprcal Analyss of DSD Daa Exploraory daa analyss was conduced on a se of one-mnue drop sze specra colleced usng a Joss-Waldvogel dsromeer (JWD) from a se n Eureka, Calforna from January hrough March 999. Ran raes were calculaed for each of he one-mnue specra o denfy ranfall evens. Seven evens rangng from o 4 hours n duraon were chosen for evaluaon. Barplos of he bnned drops and me seres of he oal number of drops/mnue were evaluaed for each even. Two general observaons were made:. The dsrbuons of drop-szes were generally unmodal and rgh-skewed. Based on vsual nspecon and fve number summares, he dsrbuons of he shapes of he onemnue DSD daa were smlar boh whn and across ranfall evens. Fgure shows barplos for some one-mnue DSDs.. The oal number of drops from one-mnue o he nex was hghly varable boh whn and across evens. Log-ransformed me seres daa suggesed an underlyng auoregressve process. Fgure shows some me-seres plos for oal number of drops. Unobserved Process A DSD, denong he number of drops per mm dameer bn nerval and per defned as 3 m of ar s well ( ) N D = ( ) ( ) n D AV D D () 3
5 where ( ) n D s he observed number of drops (raw couns) n each dameer bn nerval samplng nerval ; A s he sample surface area n D for m ; s he lengh of sample collecon n seconds; V ( D ) s he ermnal fall speed n sll ar n m/s; and, D s he drop dameer nerval n mm. Now we defne he number of drops over dameers D as ( ) ( ) n D = n P D () where n s he oal number of raw couns over all drop dameers for he me nerval and s assumed ndependen of P( D ), a probably densy funcon defnng he shape of he DSD over he nerval ; P( D ) gves he probably of fndng a drop whn a dameer D o D + dd. Subsung equaon no equaon yelds ( ) N D = n P( D ) ( ) AV D D (3) Nex, defne N as he oal number of drops per he DSD formulaed n erms of oal couns as 3 m of ar. Subsuon no equaon 3 gves ( ) N D = N ( ) P D D (4) For a ranfall even of duraon T, he unobserved saes N and N ( D ) are esmaed hrough models for P( D ) ; =,,T. To be conssen wh commonly colleced dsromeer daa, each represens a one-mnue nerval. Model for Toal Number of Drops, Le N N denoe he oal number of drops a me mnues, =,,,T. To accoun for he seral correlaon and burs varably observed n many of he me seres of oal number of drops n conguous DSD specra, we nroduce a model composed of a me-varyng mean process negraed wh a GARCH(,) condonal varance process. Le TN = ln N = +. N 4
6 The s are an ndependen sequence of correlaed errors from an unknown dsrbuon. The values, condonal on all prevous values are assumed o be normally dsrbued:,, N,# ( ). Under a GARCH(,) model = a + a # + b # where he parameers a, a, and b are subjec o he consrans a >, a, b, a <, and a + b (Alexander 998). By subsuon, ( ) ( ) N,a + a + b, N,a + a + b #, #,,, N (,a + a # + b # ) Gven he condonal dsrbuon of,,,,, he jon probably for s: #$,,, % & = # $,, % # &,, % $ 3 & #$ % & # % = (' # $ +,,,, % &) #$ % & $ = & #,. % = (. ' exp - ( * ) ) ( = $ /. + ). * exp, - ( ) & / Assume ha condonal on ndependen and are dsrbued as N (, ) Se N and he parameers N { a, a, b, # } N TN # $ = TN,,TN T N,, NT ;% ' & = ' ' TN N T = =, he ln N s are. Then he jon dsrbuon of ln N s:,,% &T ( ) + - exp,* TN )% & % * N.- & =. Ths reduces ( ) - / /. T, = # $( P o he consan erm % ' exp & $ TN # $ N (' # /# - ( ( exp,* TN - )% $ & % * N.- & ( ) # and yelds: ) ' * 3 +' # ( ) / - - 5
7 Prors on GARCH parameers Snce he GARCH parameers are consraned o be greaer han zero, each coeffcen s assgned a lognormal pror: a a b : Lognormal, : Lognormal, : Lognormal, ( # ) ( # ) ( # ) The hyperparameers,,,, and, # # # are aken as fxed, wh values ha may vary by applcaon. Mean process for number of drops The mean for N N N N TN, = a + where varance N, s modeled as a frs-order auoregressve process. Le N s an ndependen random normal process wh mean and. The parameers are { a N,#} =. Condonal on a N and, N N a N N #$ % &' N a (,( N N ) ); ( ) IG ( ),* ) The parameer a N s esmaed assumng a N N ( N, N ). The hyperparameers,#,, and # are aken as fxed, wh values ha may vary by applcaon. N N DSD shape, P( D ) The DSD shape s modeled as a lognormal dsrbuon wh me varyng parameers and. Condonal on and, he D s are assumed ndependen yeldng he jon dsrbuon for D as ( ln d & ) T $ % ( D, ) ' ', - = exp *& + = # d '. '/ where N (, ); # IG ( $ #,% # ). The hyperparameers fxed, wh values ha may vary by applcaon. # $ % are aken as,,, and 6
8 Daa Componens The herarchcal model ncorporaes wo common daa sources, namely, ran gauge and ground radar observaons. Snce gven a DSD he nsananeous (rue) ran rae and a derved (rue) radar reflecvy can be calculaed, we develop a measuremen error model for boh he ran gauge and radar observaons. Gauge Observaons, G Ran gauges measure he amoun of ranfall accumulang over me a a fxed locaon n space. If he amoun of ranfall measured a a pon s measured whou error, he gauge measuremen, G, can be equaed o he approprae me negral of ran rae. The predomnan error n gauge measuremens s a sysemac bas nduced by wnds, whch ofen resuls n an underesmae of he surface ranfall. Thus, defne a measuremen error model for gauges as a funcon of he rue ran rae and wnd speed (akng oher sysemac bases as neglgble). Defne a gauge observaon as $ % & & G = ' c w + R ds ( + ) * ( ) s G & # & where R s he derved nsananeous ran rae a me ; ( ) c w are gauge ype-specfc coeffcens ( c ) prmarly based on wnd speed, w; and errors modeled as ( N, G) ; G IG( G, G) :. G are ndependen measuremen Gven a DSD, he nsananeous ran rae can be derved based on meeorologcal prncples as a funcon of he hrd power of drop dameer (Baan 973). The ran rae, R, for a gven, s esmaed as 3 = R ( ) ( ) 6 # R c D V D N D dd where D s he drop dameer n mm; V(D) s a deermnsc funcon for he ermnal velocy for drops wh dameer D; per N ( D ) s he DSD gvng he number of drops per mm dameer D 3 m a me ; and c R s an addonal consan o accoun for uns. 7
9 Gven he model for ( ) N ( D ) yelds: where P( D ). N D (equaon 4), summng over all D 3 R 6 N s he oal drop concenraon a me and and subsung N P( D) for R = c ' D V ( D) # % N f ( D) $ & dd (5) f ( D ) s he pdf for D defned hrough Assumng an ndependen Gaussan measuremen error model, and condonal on ran rae as derved from a DSD and parameer { G G} case G (, N m G G) : where measuremens, of he nsananeous ran rae defned above. Ground Radar Observaons, =, he elemens of G are ndependen and n each m G s he negral over he me nerval beween gauge The equvalen radar reflecvy observaons from ground radar,, are a funcon of he average reurned power and range of he radar scan (Baan 973). If here are no measuremen errors, wll conform o he reflecvy derved from a DSD. To accommodae varous facors ha can lead o errors when esmang, we consruc a smple measuremen error model for equvalen radar reflecvy as a funcon of he derved reflecvy and random errors; D s he meeorologcally derved reflecvy a me and errors modeled as ( N, ); IG(, ) :. The reamen of derved reflecvy, = + where D are ndependen measuremen D, s smlar o ha for ran rae. Gven a DSD, D can be derved as he sxh momen of he DSD (Baan 973). Derved reflecvy a me s esmaed as 6 D = ( ) c D N D dd where D s he drop dameer n mm; N ( D ) s he curren drop specra gvng he number of drops per mm dameer D per 3 m a me ; and c s a consan relaed o sample volume and 8
10 3 equals one for a one- m sample. As for gauge observaons, summng over all subsung N P( D) for ( ) N D yelds: D and 6 D = ( ) & c D $ N f D # % dd where varables are as prevously defned. In pracce, D s esmaed as a sum over bnned drop szes. Assumng hs ndependen Gaussan measuremen error model, and condonal on D derved from a DSD and parameer { } each case N ( D, ). Model Summary =, he elemens of are ndependen and n We use he componens presened above o consruc he herarchcal model formulaon. A he hghes level we assume, based on a seres of condonal ndependence assumpons, ha gven a DSD ( N( D ) ), he gauge (G) and radar () observaons are condonally ndependen wh he followng facorzaon: ( ),, = ( ), ( ) # N D G $ # N D $ G N D # $ The componens G N ( D) and N ( D) ( ) ( ) ( ) = # N D $ G N D N D # $ # $ # $ # $ represen lkelhoods based on he gauge and radar daa, respecvely. The frs componen, N ( D) # $ represens he pror probably for he DSD process. The condonal probables for he full model, ncludng he parameers and hyperparameers, are summarzed n Table. 9
11 Table. Herarchcal Model Summary. Varables Observaonal Daa Hdden Process Temporal Dynamcs for mean of TN Model Parameers b G N ( D ) # $ G % ( ) N D # $ % TN N, # $ N % $ D D #% # $ N % Condonal Probables for Model Componens a c G,, N, D, H G # H G # H N # H N # H D # $ # % = $ % $ % $ % $ H $ % D % where H s a collecon of hyperparameers and { } ; G = G = $ G G # % & G ' { } ; = = $ # % & ' $ N = { a, a, b,% } = ' a ; & ( ' a ; & ( ' # # ; & ( # '% ( ) * ) * ) * ) * # {, } ; ; = an $ = & a N N % ' & $ N % ' ( ) ( ) # {, } ; D = = ' $ ( ' % ; & ( ) * ) * Hyperparameers & H = { G, ' G,,',,,, ',',', N,', #,'#, $,( $,%, %} N a The bracke noaon # $ s used as a shor hand for denong condonal dsrbuons. b N P D. c N ( D ) s general noaon for ( ) TN = ln N.
12 Parameer Esmaon We used a Markov Chan Mone Carlo (MCMC) mplemenaon o sample from he poseror dsrbuon of he herarchcal model. The approprae Markov chan s consruced usng a sngle-componen Meropols-Hasngs algorhm based on he herarchcal srucure of he model, and on he condonal probables assocaed wh he varous model componens. Sarng values, hyperparameer values, he lengh of burn-n, and he number of eraons mus be deermned for each applcaon. Sarng values are random, based on he pror dsrbuons of he parameers; he dsrbuons can be alered by he choce of he fxed hyperparameers. The requred number of eraons for a gven applcaon depends on he lengh of he burn-n, he auocorrelaon of he parameers, and he number of samples needed for parameer characerzaon. Convergence s monored usng vsual analyss of ploed parameers n wo ways: ) by performng mulple runs usng dfferen random sars; and ) runnng mulple chans usng he same sarng values. APPLICATION TO EUREKA RAINFALL Daa The Eureka se consss of hree sources of daa: dsromeer, ran gauge, and ground radar daa. The dsromeer daa were colleced usng a JWD nsrumen. The dsromeer was locaed near he Eureka radar se a 4 48'N/4 9'W. Each daa fle conaned one hours worh of raw bn couns colleced a one-mnue nervals. The raw bn couns were convered o one-mnue DSD specra couns 3 N ( m mm ) by adjusng raw couns for he nsrumen collecon area, D me, ermnal fall speed, and drop dameer. A se of processed -mnue dsromeer daa was also avalable. The -mnue daa are drop specra ha have been averaged over conguous one-mnue perods. The observed -mnue drop specra were used o calculae ran raes and reflecvy values o use n place of gauge and radar observaons for model verfcaon. Hourly gauge accumulaon daa (TD-34) were obaned from he Naonal Clmac Daa Cener (NCDC) on-lne servce. The daa ncluded measuremens from 4 gaugng saons n he vcny of Eureka, Calforna over he perod January 999 hrough 3 March 999, of whch nne saons were n he vcny of he radar and dsromeer. Daa were colleced usng Fscher-Porer precpaon gauges wh auomaed readous. Precpaon daa from he NCDC daabase are qualy checked and eded, as necessary, by an auomaed and a manual ed
13 (Hammer and Seurer 998). To concde wh he me doman only mplemenaon of he model, gauge daa were lmed o measuremens from he gaugng saon locaed closes o he dsromeer a 4 49'N/4 'W. Level II WSR-88D radar scans were ordered from NCDC for he Eureka radar se, whch s locaed a 4 9'54''N/4 7'3''W, for perods concdng wh ranfall evens denfed by gauge and dsromeer daa. There was a radar scan approxmaely every sx mnues for several elevaon ncremens. The sphercal coordnae Level II daa were nerpolaed o a km x km Caresan coordnae sysem of equvalen reflecvy (db) values usng sandard programs for NEXRAD radar. The nerpolaed daa were only avalable for one of he denfed ranfall evens; he avalably of radar daa defned, and lmed, he me perod over whch he model was mplemened. Approprae reflecvy values needed o be exraced from he km x km grdded values o correspond wh he fxed spaal locaon of he dsromeer and closes ran gauge. Reflecvy n he me doman was exraced from he spaal daa by frs denfyng he closes pxel locaed vercally above he dsromeer and nearby ran gauge. Gven he uncerany n he vercal flow of ranfall above he surface of he earh o a pon on he ground, an equvalen reflecvy was calculaed as he aeral average of frs and second order neghbors usng he denfed pxel and he egh neghborng pxels. Issues of beam blockage along he pah of he radar o he locaon of he dsromeer were consdered o be neglgble due o he gradually ncreasng naure of he coasal erran of Norhern Calforna. Therefore, he reflecvy values were exraced from he lowes radar beam (l elevaon of.5 o ). We mplemened he herarchcal Bayes model usng daa from a nne hour (55 mnues) ranfall even on 6 February 999. For he course of he model run here were nne hourly gauge readngs and 9 ground radar observaons. Specfcaon of Hyperparameers, Inal Values, and Consrans on Parameers We used exploraory analyss of dsromeer daa from fve February ranfall evens o come up wh nformave means and varances for he fxed hyperparameers. The consan hyperparameer values are summarzed n Table wh model resuls.
14 In addon o he hyperparameer values, we specfed several = sarng values as follows: = = =. We based nal condons for he remanng,., and N 3. parameers on random draws from he approprae pror dsrbuons. We nroduced cuoffs for wo ypes of consrans; namely, consrans due o ranfall dynamcs/nsrumenaon and consrans due o sascal consderaons. Consrans based on ranfall dynamcs ncluded upper bounds for hree parameers:. We lmed he oal number of drops n a one-mnue nerval o 6., abou.5% above he maxmum value observed from he dsromeer daa. Ths corresponds o a oal of exp 6 43 raw coun drops when summed over he drop dameer bns. We ( ) appled hs cuoff o he parameers TN and N.. We lmed he upper bound for drop dameer sze such ha P( D 5).. Ths lmaon was mposed snce he ran even of neres represened lgh ranfall. The dsromeer daa had no raw couns beyond a mean drop dameer of 3.5 mm. Proposals for parameer updang are drawn from unform dsrbuons cenered a he curren parameer values. Ths mplemenaon made possble o draw negave-valued proposals, dependng on he curren value of he parameer. Consrans based on sascal consderaons hus ncluded lower bounds on varance parameers ( ) model consrans were also placed on he hree GARCH coeffcens. Esmang m G and D > and >. Requred A each me perod for whch here were gauge or radar observaons, he daa componens of he model requred compuaon of he wo quanes m G and D, he rue equvalen gauge value and he rue derved reflecvy, respecvely. Boh of hese values are deermnsc quanes based on he curren sae of he parameers ha defne he hdden DSD process. We calculae he rue equvalen gauge value as he sum over he prevous 6 mnues of each of he nsananeous ran raes n mm/mn derved from he sae of he DSD before each pon n me. The nsananeous ran raes defned by equaon 5 are approxmaed by summng over he drop dameer bns ha would ypcally be observed usng a JWD nsrumen as follows: { TN } $ P( D ) 3.6 exp # R D V D D &( ') 3 = * & ' ( )% 6 = % D 3
15 where ( ) P D s based on he cdf of a lognormal dsrbuon wh parameers ermnal fall speed, V ( D ) drop dameer for bn. Then, was calculaed as and. The D (Smh 993), where D s he mean = 6+ m G = True derved reflecvy s calculaed as an average, over he prevous mnues, of he sxh momen of he DSD. As wh he equvalen gauge calculaon, we approxmae derved reflecvy by summng over he drop dameer bns usng 6 R Then D, n uns of { TN }# P( D ) D D %' &( exp 6 = ) % & $ = $ D = + D = 6 3 mm m, are convered o db as log ( D ). Tunng, Thnnng, and Convergence Tunng he chan was opmzed by makng many small runs (5 o 5 eraons n lengh) and monorng he accepance raes for each of he parameers. Adjusmens o he accepance raes were made by eher ncreasng or decreasng he upper and lower bounds on he unformly dsrbued parameer proposals. Afer much expermenaon, bounds on he unform proposals were adjused o yeld accepance raes rangng from abou.4 on he low end o abou.85 on he hgh end. Mos parameers had accepance raes n he range of Whle hese values exceed he commonly recommended values of. o.5 for updang muldmensonal componens (Robers 996), we found ha mxng and convergence performed beer a hese hgher levels. We used hnnng o oban roughly ndependen samples from he jon poseror dsrbuon. We used gbbs (Rafery and Lews 995) o ge esmaes of he deal spacng beween eraons and he number of burn-n eraons, as well as he number of eraons for a desred precson. In mos cases a hnnng value of 5 was more han suffcen o oban 4
16 approxmaely ndependen samples, alhough on occason, a parameer hnnng value would h 3 or so. We monored convergence by vsual nspecon of ploed oupu. The number of burnn eraons recommended by gbbs was generally low and always < 5. Convergence was also assessed by lookng a he pah of several runs usng boh he same and dfferen sarng values. Resuls Model Verfcaon There s no way o specfcally valdae he model n erms of comparng oupus o rue values. Pror o runnng he fully mplemened model, however, we conduced model verfcaon by assessng he ably of he model o realscally reproduce argeed sgnals and negraed parameers assocaed wh surface ranfall. We replace NEXRAD radar observaons wh derved reflecvy calculaed from -mnue dsromeer daa. Ran gauge daa are omed for hs assessmen. To maxmze he npus we use a -mnue me sep, hus nsurng a daa value a each me sep. The verfcaon run ncludes 49 derved reflecvy observaons compued from he one se of dsromeer daa. The MCMC run was ended afer 8, eraons followng, burn-n eraons. Based on graphcal oupu, all parameers converged whn he, burnn perod. Resuls for a ypcal model verfcaon run showng me ndependen parameer esmaes and me varyng parameer esmaes (usng he 8 MCMC samples generaed afer applyng a hnnng value of 35) are shown n Tables 3, respecvely. Furher, usng he poseror parameer esmaes for he hdden process we can compare calculaed negraed ranfall parameers wh values obaned from he -mnue dsromeer daa. Fgures 3 and 4 show he comparsons for ran rae and reflecvy. Noe ha alhough we compare he model oupu wh quanes derved from dsromeer daa, we do no consder he dsromeer daa as ground ruh. Dsromeer daa are subjec o measuremen errors. The poseror means of he hree GARCH coeffcens are all lower han hose specfed for he pror dsrbuons. The varably for he a and b coeffcens s slghly less han ha for he prors. The mean of he AR() coeffcen s slghly hgher han ha specfed n he 5
17 pror. The varance parameer for he AR() errors,, s consderably hgher han ha specfed n he pror, as s he measuremen error varance for radar ( ). The mean values for TN and N are abou he same and poseror esmaes of TN are close o expeced couns alhough he model ends o underesmae he couns. In general, he DSD dsrbuon parameers and are more varable han esmaes based on he dsromeer daa. On average, he poseror esmaes of dsromeer daa. The model ends o slghly underesmae are smlar o he values esmaed from he The poseror esmaes for boh negraed parameers (ran rae and reflecvy) are more varable han values calculaed from he dsromeer daa realzaon. Ths s expeced snce hese esmaes nclude he dsromeer daa realzaon wh some uncerany. Asde from he exra varably, he model esmaes are a good mach wh he negraed quanes calculaed from dsromeer observaons. The model does a good job of capurng he overall means and he sgnals presen n he me seres. Full Model Run The full model run ncludes 9 hourly ran gauge observaons and 9 effecve reflecvy observaons from ground radar. The MCMC run was ended afer 45, eraons followng a burn-n perod of 35, eraons. Afer applyng a hnnng value of 35, here were 7 MCMC samples for poseror esmaes. Convergence for he full model run s much slower han ha for he verfcaon run, akng up o abou 6, eraons for a few parameers. Wh he excepon of he AR() varance erm, eraons. The., all parameers appear o converge by 6, varance hardly moves unl 6, eraons, a whch pon he value sars bouncng around and hen seles down a abou.5, followed by more nermen bouncng. Ths paern could be due o slow mxng or denfably ssues. Resuls for a ypcal full model run showng me ndependen parameer esmaes and me varyng parameer esmaes are shown n Tables 3, respecvely, along sde he verfcaon run resuls. Tme seres resuls for negraed parameers calculaed from he esmaed parameers for he hdden DSD process are shown n Fgures 5 8. Excep for a, he poseror esmaes for he GARCH parameers are smlar o hose from he verfcaon run. The generally larger sandard errors reflec he greaer uncerany n 6
18 he one-mnue model, whch does no have daa a every me sep. The mean error erms,, are no sgnfcanly dfferen from zero based on 95% confdence nervals for he means. The poseror dsrbuons are relavely symmerc and approxmaely normal, whch s n lne wh he GARCH model for. The whn me poseror means for TN and N are essenally he same. The poseror dsrbuons for TN are nearly symmercal and appear o be approxmaely normally dsrbued. The represenave dsrbuons all have means close o he mean for he pror dsrbuon. The poseror dsrbuons are farly symmerc; values for are varable, rangng from -.8 o. and do conan he observed esmaes based on he one realzaon of dsromeer daa. The represenave dsrbuons all have means close o.8, he mean of he pror dsrbuon. The poseror dsrbuons are rgh-skewed; values for are varable, rangng from.5 o.3 and do conan he observed esmaes calculaed from he dsromeer daa. The n-seres varably for TN,, and s greaer han ha for he observed dsromeer daa. Excep for, he me seres values esmaed from he dsromeer daa are nearly cenered n he poseror model resuls. As wh he verfcaon run, he full model underesmaes he DSD shape parameer,. From he poseror parameer esmaes for he hdden DSD process, we can esmae he negraed values for ran rae and derved reflecvy. The me seres of calculaed ran raes show hgher varably han ha produced from he one realzaon of dsromeer daa. The large varably makes dffcul o compare he sgnals, hough here appears o be some amoun of agreemen n he gross sgnal feaures for ran rae (Fgure 5). The me seres plos of derved reflecvy also show hgher varably han ha produced from one realzaon of dsromeer daa, hough no o he same exen as for ran rae (Fgure 6). The model oupu conssenly overesmaes dsromeer derved values. A mean normalzed plo for reflecvy (Fgure 7) removes he bas and clearly shows some smlary n he sgnals compared wh he dsromeer based calculaons. Whle he cause of he bas s unknown, s possble ha he observed bas ndcaes poor calbraon of he radar. 7
19 An alernae vsualzaon of he calculaed ran rae parameer based on model oupu dsplays hsograms for ran raes from pons along he me seres (Fgure 8). Correspondng values calculaed from he dsromeer daa are also shown. The ndvdual hsograms conan he dsromeer-based ran raes, alhough he dsromeer-based values are n he lower ends of he dsrbuons for mos me pons. Table. MCMC Summary Sascs for Tme Independen Parameers. MCMC Poseror Mean (SE) Model Parameer Pror Mean Pror SD Verfcaon Run Full Run GARCH Coeffcens a (.) (.) a (.8) (.) b (.99) (.3) AR() Process a (.8) (.) N (.3) Measuremen Errors.478 (.33) G (.44) (.3) (.865) 8
20 Table 3. MCMC Summary Sascs for Tme Varyng Parameers for hree me pons. For he verfcaon run, mes are = 9, 4, 39. For he full run, = 9, 35, 39. Verfcaon Run Full Model Run Model Parameer Observed Value Poseror Mean (SE) Observed Value Poseror Mean (SE)n Toal # drops TN (.6) (.56) TN TN (.9) 6.9 (.8) (.55) 5.4 (.5) Mean # drops N N N 3 6. (.6) 6.34 (.8) 6.3 (.6) 5. (.5) 5.5 (.49) 5.5 (.48) Errors for TN 3 -. (.8) -.5 (.8) -. (.) -.7 (.4) -.6 (.).3 (.7) Error Varance 3.5 (.5).7.36).8 (.4).3 (.77).39 (.84).39 (.75) Lognormal Scale Parameer Pror mean: -.5 SD:.4 3 Lognormal Shape Parameer Pror mean:.8 SD: a a a Esmaed from -mnue dsromeer daa. b Esmaed from -mnue dsromeer daa (.9) -.46 (.9) -.49 (.9).55 (.8).76 (.3).66 (.3) -.56 b b (.5) -.54 (.5) -.53 (.5).8 (.3).9 (.3).9 (.4) 9
21 DISCUSSION The model resuls for he Eureka ranfall are promsng. Usng herarchcal Bayes mehods we are able o negrae nformaon from mulple daa sources for purposes of esmang ran rae and relaed parameers. Whle we are unable o valdae he poseror esmaes of TN,, and, he man parameers for he unobserved DSD process, we can compare her values wh he one avalable dsromeer realzaon. For TN, many of he MCMC realzaons, hough no all, produce me seres sgnals smlar o ha observed for he dsromeer realzaon wh model oupu ha s always more varable. Ths hgher varably s due, n par, o he hgher uncerany ha s nroduced when observaons from wo dfferen daa sources are consdered. The DSD scale and shape parameers capure he gross feaures of he drop specra, bu he me seres for dsromeer esmaes. and are many mes more varable han hose produced from A comparson of auocorrelaons for hese parameers based on dsromeer daa versus model oupu suggess ha here s more srucure n he me evoluon of he DSD shape han wha s accouned for n he curren model. Ulmaely he esmaed DSD parameers are used o produce esmaes of he negraed quanes ran rae and derved reflecvy. The verfcaon runs do a good job of capurng he feaures of hese quanes; he full model runs capure he gross feaures, bu agan produce hgher varable me seres oupu and overesmaes of derved reflecvy. The bas seen n derved reflecvy could be due o a number of facors; mos noably, dfferences may be relaed o samplng errors assocaed wh dsromeer daa. Derved reflecvy s more adversely affeced by small dfferences n drop-sze dsrbuons snce funcon of D s a funcon of 6 D whle R s a 3 D. Eher overesmaon of larger drops by he model, or under recordng of larger drops by he dsromeer could lead o hese observed dfferences. Verfcaon runs produced beer conssency beween model derved reflecvy and dsromeer derved values; use of - mnue dsromeer daa helps compensae for dsromeer samplng errors ha arse due o he small sample volume colleced n a one-mnue nerval. There s no measurable ruh when comes o deermnng ran raes a he earh s surface. For each nsrumen deployed n he feld, here are clams and belefs as o how well hey perform her respecve asks, n lgh of uncerany. Insead of usng ran gauge daa as ground ruh, he poseror model dsrbuons capure mulple sources of uncerany ha
22 reflec uncerany assocaed wh one or more sources such as: he nably o drecly measure quanes of neres; measuremen and samplng errors assocaed wh he mulple nsrumens; sparse gauge daa; and, ncorrec or nadequae assumpons and specfcaons relaed o equaons used o represen he underlyng physcal processes. As a proof of concep, he herarchcal Bayes model s able o generae reasonable esmaes for mos of he modeled parameers when appled o he Eureka daa se. Usng exsng model runs as baselnes we can furher characerze he observed unceranes by makng sepwse modfcaons o model sages or by ncorporang addonal daa. In he shor erm, a emporal dynamc sage should be developed for he evoluon of DSD hrough me. Also, we should run he model for oher ranfall evens. The parcular Eureka even conssed of relavely lgh ranfall. Deermnng how well he model performs for heaver ranfall evens wll provde some dea as o he porably/usably of hs model under a varey of ranfall scenaros. In he longer erm, some deas for expandng and usng hs modelng approach nclude:. If applcable, ncorporae wnd correcons for gauge observaons n fuure model runs. Gauge cachmens are nherenly affeced by wnds; a wnd effec would be one componen of he overall uncerany observed n he model oupu.. Incorporae a vercal correcon model componen for radar reflecvy observaons, based on an underlyng process model, o compensae for he fac ha ground-based radar provdes nformaon abou precpaon above he earh s surface. Include daa from vercal ponng radar o augmen he vercal correcon model componen. 3. Expand he me only model o a space-me model. Such a model could provde a ool ha allows for he creaon of spaal ranfall maps wh uncerany esmaes, whle also provdng he opporuny for more expanded analyss of he componen unceranes. 4. Use model oupu o esmae beer -R relaonshps dynamcally. 5. Use model oupu o assess he value added by usng oher ypes of precpaon measuremens (or by removng some nsrumens). ACKNOWLEDGEMENTS We would lke o hank Dr. Sandra Yuer, who provded he dsromeer and NEXRAD radar daa ses, for her help wh daa pre-processng and her nvaluable nsghs relaed o nsrumenaon, meeorologcal properes of precpaon, and relaonshps beween he wo. REFERENCES Alexander C. Volaly and correlaon: measuremen, models, and applcaons. In Alexander C, edor, Rsk Managemen and Analyss. Vol. : Measurng and Modellng Fnancal Rsk. John Wley Sons Ld, 998.
23 Baan LJ. Radar observaon of he amosphere. Unversy Chcago Press, Chcago, revsed edon, 973. Bell TL. A space-me sochasc model of ranfall for saelle remoe-sensng sudes. Journal of Geophyscal Research, 9(D8): , 987. Berlner LM. Herarchcal Bayesan modelng n he envronmenal scences. Allgemenes Sasshe Archv, 84:4-53,. Berlner LM, Levne RA, and Shea DJ. Bayesan clmae change assessmen. Journal of Clmae, 3():385-38, a. Cox DR, S. RF, and IshamV. A smple spaal-emporal model of ranfall. Proc. R. Soc. Lond., 45:37-38, 988. Hammer G and Seurer P. Daa documenaon for hourly precpaon daa TD-34. Techncal repor, Naonal Clmac Daa Cener, February, 998. Hrafnkelsson B. Herarchcal modelng of coun daa wh applcaon o nuclear fall-ou. Journal of Envronmen and Ecologcal Sascs, :79-, 3. Rafery A and Lews S. Implemenng MCMC. In Glks W, Rchardson S, and Spegelhaler D, edors, Markov Chan Mone Carlo n Pracce, pages 5-3. Chapman Hall, London, 996. Robers GO. Markov chan conceps relaed o samplng algorhms. In Glks W, Rchardson S, and Spegelhaler D, edors, Markov Chan Mone Carlo n Pracce, pages Chapman Hall, London, 996. Rodrguez-Iurbe I, Cox DR, F.R.S., and Isham V. A pon process model for ranfall: furher developmens. Proc. R. Soc. Lond., 47(83-98), 988. Rodrguez-Iurbe I and Eagleson PS. Mahemacal models of ransorm evens n space and me. Waer Resources Research, 3():8-9, 987. Royle JA, Berlner LM, Wkle CK, and Mllff R. A herarchcal spaal model for consrucng wnd felds from scaeromeer daa n he Labrador Sea. In Case Sudes n Bayesan Sascs IV, pages Sprnger-Verlag, 998. Skaugen T, Creun JD, and Goschald L. Reconsrucon and frequency esmaes of exreme daly areal precpaon. Journal of Geophyscal Research, (D): , 996. Smh JA. Marked pon process models of randrop-sze dsrbuons. Journal of Appled Meeorology, 3:84-96, 993. Sern R and Coe R. A model fng analyss of daly ranfall daa. Journal of he Royal Sascal Socey A, 47(Par ):-34, 984. Wkle CK, Mllff RF, Nychka D, and Berlner M. Spaoemporal herarchcal Bayesan modelng: ropcal ocean surface wnds. Journal of he Amercan Sascal Assocaon, 96(454):38-397,. Wkle CK. Herarchcal Bayesan space-me models. Envronmenal and Ecologcal Sascs, 5:7-54, 998.
24 Fgure. Barplos for nne consecuve mnues of dsromeer DSD daa from a Eureka ranfall even; n s he oal number of drops n each one-mnue nerval. Fgure. Tme seres of he one-mnue oal number of drops (lef) and log ransformed oal number of drops (rgh) for each of fve Eureka ranfall evens. 3
25 Sample: Sample: 5 log(rr) (mm/hr) log(rr) (mm/hr) me (mnues) me (mnues) Sample: 3 Sample: 45 log(rr) (mm/hr) log(rr) (mm/hr) me (mnues) me (mnues) Sample: 5 Sample: 7 log(rr) (mm/hr) log(rr) (mm/hr) me (mnues) me (mnues) Fgure 3. Tmes seres plos of poseror log ran rae esmaes calculaed for some represenave samples from he valdaon run (lgh lnes) wh ran rae esmaes from he dsromeer daa (heavy lnes). Sample: Sample: 5 Reflecvey (db) 3 Reflecvey (db) me (mnues) me (mnues) Sample: 3 Sample: 45 Reflecvey (db) 5 3 Reflecvey (db) me (mnues) me (mnues) Sample: 5 Sample: 7 Reflecvey (db) 3 Reflecvey (db) me (mnues) me (mnues) Fgure 4. Tmes seres plos of poseror reflecvy esmaes calculaed for some represenave samples from he valdaon run (lgh lnes) wh reflecvy esmaes from he dsromeer daa (heavy lnes). 4
26 Sample: 4 log(rr) (mm/hr) me (mnues) Sample: 365 log(rr) (mm/hr) me (mnues) Sample: 4 log(rr) (mm/hr) me (mnues) Fgure 5. Tme seres plos of full model run poseror log ran rae esmaes calculaed for seleced samples (lgh lnes) wh logged ran rae values based on dsromeer daa (heavy lnes). Sample: 4 Reflecvy (db) me (mnues) Sample: 365 Reflecvy (db) me (mnues) Sample: 4 Reflecvy (db) me (mnues) Fgure 6. Tme seres plos of full model run poseror reflecvy esmaes calculaed for seleced samples (lgh lnes) wh reflecvy values based on dsromeer daa (heavy lnes). 5
27 Sample: 4 Reflecvy (db) me (mnues) Sample: 365 Reflecvy (db) me (mnues) Sample: 4 Reflecvy (db) me (mnues) Fgure 7. Tme seres plos of full model run poseror reflecvy (normalzed o mean ) esmaes calculaed for seleced samples (lgh lnes) wh reflecvy values based on dsromeer daa (heavy lnes). Fgure 8. Frequency hsograms of he full model run poseror dsrbuon of log ran rae esmaes for seleced me pons whn he 55 mnue modeled even. All hsograms are scaled wh equvalen y-axes. The whe crcles correspond o he values of ran rae calculaed from one realzaon of dsromeer daa. 6
Correlation of default
efaul Correlaon Correlaon of defaul If Oblgor A s cred qualy deeroraes, how well does he cred qualy of Oblgor B correlae o Oblgor A? Some emprcal observaons are efaul correlaons are general low hough hey
More informationOnline Technical Appendix: Estimation Details. Following Netzer, Lattin and Srinivasan (2005), the model parameters to be estimated
Onlne Techncal Appendx: Esmaon Deals Followng Nezer, an and Srnvasan 005, he model parameers o be esmaed can be dvded no hree pars: he fxed effecs governng he evaluaon, ncdence, and laen erence componens
More informationChain-linking and seasonal adjustment of the quarterly national accounts
Sascs Denmark Naonal Accouns 6 July 00 Chan-lnkng and seasonal adjusmen of he uarerly naonal accouns The mehod of chan-lnkng he uarerly naonal accouns was changed wh he revsed complaon of daa hrd uarer
More informationNormal Random Variable and its discriminant functions
Normal Random Varable and s dscrmnan funcons Oulne Normal Random Varable Properes Dscrmnan funcons Why Normal Random Varables? Analycally racable Works well when observaon comes form a corruped sngle prooype
More informationSection 6 Short Sales, Yield Curves, Duration, Immunization, Etc.
More Tuoral a www.lledumbdocor.com age 1 of 9 Secon 6 Shor Sales, Yeld Curves, Duraon, Immunzaon, Ec. Shor Sales: Suppose you beleve ha Company X s sock s overprced. You would ceranly no buy any of Company
More informationA valuation model of credit-rating linked coupon bond based on a structural model
Compuaonal Fnance and s Applcaons II 247 A valuaon model of cred-rang lnked coupon bond based on a srucural model K. Yahag & K. Myazak The Unversy of Elecro-Communcaons, Japan Absrac A cred-lnked coupon
More informationLab 10 OLS Regressions II
Lab 10 OLS Regressons II Ths lab wll cover how o perform a smple OLS regresson usng dfferen funconal forms. LAB 10 QUICK VIEW Non-lnear relaonshps beween varables nclude: o Log-Ln: o Ln-Log: o Log-Log:
More informationSOCIETY OF ACTUARIES FINANCIAL MATHEMATICS. EXAM FM SAMPLE SOLUTIONS Interest Theory
SOCIETY OF ACTUARIES EXAM FM FINANCIAL MATHEMATICS EXAM FM SAMPLE SOLUTIONS Ineres Theory Ths page ndcaes changes made o Sudy Noe FM-09-05. January 4, 04: Quesons and soluons 58 60 were added. June, 04
More informationInterest Rate Derivatives: More Advanced Models. Chapter 24. The Two-Factor Hull-White Model (Equation 24.1, page 571) Analytic Results
Ineres Rae Dervaves: More Advanced s Chaper 4 4. The Two-Facor Hull-Whe (Equaon 4., page 57) [ θ() ] σ 4. dx = u ax d dz du = bud σdz where x = f () r and he correlaon beween dz and dz s ρ The shor rae
More informationIFX-Cbonds Russian Corporate Bond Index Methodology
Approved a he meeng of he Commee represenng ZAO Inerfax and OOO Cbonds.ru on ovember 1 2005 wh amendmens complan wh Agreemen # 545 as of ecember 17 2008. IFX-Cbonds Russan Corporae Bond Index Mehodology
More informationDifferences in the Price-Earning-Return Relationship between Internet and Traditional Firms
Dfferences n he Prce-Earnng-Reurn Relaonshp beween Inerne and Tradonal Frms Jaehan Koh Ph.D. Program College of Busness Admnsraon Unversy of Texas-Pan Amercan jhkoh@upa.edu Bn Wang Asssan Professor Compuer
More informationThe Financial System. Instructor: Prof. Menzie Chinn UW Madison
Economcs 435 The Fnancal Sysem (2/13/13) Insrucor: Prof. Menze Chnn UW Madson Sprng 2013 Fuure Value and Presen Value If he presen value s $100 and he neres rae s 5%, hen he fuure value one year from now
More informationNetwork Security Risk Assessment Based on Node Correlation
Journal of Physcs: Conference Seres PAPER OPE ACCESS ewor Secury Rs Assessmen Based on ode Correlaon To ce hs arcle: Zengguang Wang e al 2018 J. Phys.: Conf. Ser. 1069 012073 Vew he arcle onlne for updaes
More informationFITTING EXPONENTIAL MODELS TO DATA Supplement to Unit 9C MATH Q(t) = Q 0 (1 + r) t. Q(t) = Q 0 a t,
FITTING EXPONENTIAL MODELS TO DATA Supplemen o Un 9C MATH 01 In he handou we wll learn how o fnd an exponenal model for daa ha s gven and use o make predcons. We wll also revew how o calculae he SSE and
More informationOnline appendices from Counterparty Risk and Credit Value Adjustment a continuing challenge for global financial markets by Jon Gregory
Onlne appendces fro Counerpary sk and Cred alue Adusen a connung challenge for global fnancal arkes by Jon Gregory APPNDX A: Dervng he sandard CA forula We wsh o fnd an expresson for he rsky value of a
More informationImproving Forecasting Accuracy in the Case of Intermittent Demand Forecasting
(IJACSA) Inernaonal Journal of Advanced Compuer Scence and Applcaons, Vol. 5, No. 5, 04 Improvng Forecasng Accuracy n he Case of Inermen Demand Forecasng Dasuke Takeyasu The Open Unversy of Japan, Chba
More informationDynamic Relationship and Volatility Spillover Between the Stock Market and the Foreign Exchange market in Pakistan: Evidence from VAR-EGARCH Modelling
Dynamc Relaonshp and Volaly pllover Beween he ock Marke and he Foregn xchange marke n Paksan: vdence from VAR-GARCH Modellng Dr. Abdul Qayyum Dr. Muhammad Arshad Khan Inroducon A volale sock and exchange
More informationA Framework for Large Scale Use of Scanner Data in the Dutch CPI
A Framework for Large Scale Use of Scanner Daa n he Duch CPI Jan de Haan Sascs Neherlands and Delf Unversy of Technology Oawa Group, 2-22 May 215 The basc dea Ideally, o make he producon process as effcen
More informationDYNAMIC ECONOMETRIC MODELS Vol. 8 Nicolaus Copernicus University Toruń 2008
DYNAMIC ECONOMETRIC MODELS Vol. 8 Ncolaus Coperncus Unversy Toruń 2008 Por Fszeder Ncolaus Coperncus Unversy n Toruń Julusz Preś Szczecn Unversy of Technology Prcng of Weaher Opons for Berln Quoed on he
More informationBaoding, Hebei, China. *Corresponding author
2016 3 rd Inernaonal Conference on Economcs and Managemen (ICEM 2016) ISBN: 978-1-60595-368-7 Research on he Applcably of Fama-French Three-Facor Model of Elecrc Power Indusry n Chnese Sock Marke Yeld
More informationMind the class weight bias: weighted maximum mean discrepancy for unsupervised domain adaptation. Hongliang Yan 2017/06/21
nd he class wegh bas: weghed maxmum mean dscrepancy for unsupervsed doman adapaon Honglang Yan 207/06/2 Doman Adapaon Problem: Tranng and es ses are relaed bu under dfferen dsrbuons. Tranng (Source) DA
More informationBank of Japan. Research and Statistics Department. March, Outline of the Corporate Goods Price Index (CGPI, 2010 base)
Bank of Japan Research and Sascs Deparmen Oulne of he Corporae Goods Prce Index (CGPI, 2010 base) March, 2015 1. Purpose and Applcaon The Corporae Goods Prce Index (CGPI) measures he prce developmens of
More informationFloating rate securities
Caps and Swaps Floang rae secures Coupon paymens are rese perodcally accordng o some reference rae. reference rae + ndex spread e.g. -monh LIBOR + 00 bass pons (posve ndex spread 5-year Treasury yeld 90
More informationDeriving Reservoir Operating Rules via Fuzzy Regression and ANFIS
Dervng Reservor Operang Rules va Fuzzy Regresson and ANFIS S. J. Mousav K. Ponnambalam and F. Karray Deparmen of Cvl Engneerng Deparmen of Sysems Desgn Engneerng Unversy of Scence and Technology Unversy
More informationEstimating intrinsic currency values
Esmang nrnsc currency values Forex marke praconers consanly alk abou he srenghenng or weakenng of ndvdual currences. In hs arcle, Jan Chen and Paul Dous presen a new mehodology o quanfy hese saemens n
More informationEstimation of Optimal Tax Level on Pesticides Use and its
64 Bulgaran Journal of Agrculural Scence, 8 (No 5 0, 64-650 Agrculural Academy Esmaon of Opmal Ta Level on Pescdes Use and s Impac on Agrculure N. Ivanova,. Soyanova and P. Mshev Unversy of Naonal and
More informationAgricultural and Rural Finance Markets in Transition
Agrculural and Rural Fnance Markes n Transon Proceedngs of Regonal Research Commee NC-04 S. Lous, Mssour Ocober 4-5, 007 Dr. Mchael A. Gunderson, Edor January 008 Food and Resource Economcs Unversy of
More informationA Novel Application of the Copula Function to Correlation Analysis of Hushen300 Stock Index Futures and HS300 Stock Index
A Novel Applcaon of he Copula Funcon o Correlaon Analyss of Hushen3 Sock Index Fuures and HS3 Sock Index Fang WU *, 2, Yu WEI. School of Economcs and Managemen, Souhwes Jaoong Unversy, Chengdu 63, Chna
More informationAlbania. A: Identification. B: CPI Coverage. Title of the CPI: Consumer Price Index. Organisation responsible: Institute of Statistics
Albana A: Idenfcaon Tle of he CPI: Consumer Prce Index Organsaon responsble: Insue of Sascs Perodcy: Monhly Prce reference perod: December year 1 = 100 Index reference perod: December 2007 = 100 Weghs
More informationThe UAE UNiversity, The American University of Kurdistan
MPRA Munch Personal RePEc Archve A MS-Excel Module o Transform an Inegraed Varable no Cumulave Paral Sums for Negave and Posve Componens wh and whou Deermnsc Trend Pars. Abdulnasser Haem-J and Alan Musafa
More informationAmerican basket and spread options. with a simple binomial tree
Amercan baske and spread opons wh a smple bnomal ree Svelana orovkova Vre Unverse Amserdam Jon work wh Ferry Permana acheler congress, Torono, June 22-26, 2010 1 Movaon Commody, currency baskes conss of
More informationUC San Diego Recent Work
UC San Dego Recen Work Tle On More Robus Esmaon of Skewness and Kuross: Smulaon and Applcaon o he S&P500 Index Permalnk hps://escholarshp.org/uc/em/7b5v07p Auhors Km, Tae-Hwan Whe, Halber Publcaon Dae
More informationThe Effects of Nature on Learning in Games
The Effecs of Naure on Learnng n Games C.-Y. Cynha Ln Lawell 1 Absrac Ths paper develops an agen-based model o nvesgae he effecs of Naure on learnng n games. In parcular, I exend one commonly used learnng
More informationPricing and Valuation of Forward and Futures
Prcng and Valuaon of orward and uures. Cash-and-carry arbrage he prce of he forward conrac s relaed o he spo prce of he underlyng asse, he rsk-free rae, he dae of expraon, and any expeced cash dsrbuons
More informationFugit (options) The terminology of fugit refers to the risk neutral expected time to exercise an
Fug (opons) INTRODUCTION The ermnology of fug refers o he rsk neural expeced me o exercse an Amercan opon. Invened by Mark Garman whle professor a Berkeley n he conex of a bnomal ree for Amercan opon hs
More informationTerms and conditions for the MXN Peso / US Dollar Futures Contract (Physically Delivered)
The Englsh verson of he Terms and Condons for Fuures Conracs s publshed for nformaon purposes only and does no consue legal advce. However, n case of any Inerpreaon conroversy, he Spansh verson shall preval.
More informationGlobal regional sources of risk in equity markets: evidence from factor models with time-varying conditional skewness
Global regonal sources of rsk n equy markes: evdence from facor models wh me-varyng condonal skewness Aamr R. Hashm a, Anhony S. Tay b, * a Deparmen of Economcs, Naonal Unversy of Sngapore, AS2, Ars Lnk,
More informationNoise and Expected Return in Chinese A-share Stock Market. By Chong QIAN Chien-Ting LIN
Nose and Expeced Reurn n Chnese A-share Sock Marke By Chong QIAN Chen-Tng LIN 1 } Capal Asse Prcng Model (CAPM) by Sharpe (1964), Lnner (1965) and Mossn (1966) E ( R, ) R f, + [ E( Rm, ) R f, = β ] + ε
More informationTax Dispute Resolution and Taxpayer Screening
DISCUSSION PAPER March 2016 No. 73 Tax Dspue Resoluon and Taxpayer Screenng Hdek SATO* Faculy of Economcs, Kyushu Sangyo Unversy ----- *E-Mal: hsao@p.kyusan-u.ac.jp Tax Dspue Resoluon and Taxpayer Screenng
More informationAccuracy of the intelligent dynamic models of relational fuzzy cognitive maps
Compuer Applcaons n Elecrcal Engneerng Accuracy of he nellgen dynamc models of relaonal fuzzy cognve maps Aleksander Jasrebow, Grzegorz Słoń Kelce Unversy of Technology 25-314 Kelce, Al. Tysącleca P. P.
More informationOptimal Combination of Trading Rules Using Neural Networks
Vol. 2, No. Inernaonal Busness Research Opmal Combnaon of Tradng Rules Usng Neural Neworks Subraa Kumar Mra Professor, Insue of Managemen Technology 35 Km Mlesone, Kaol Road Nagpur 44 502, Inda Tel: 9-72-280-5000
More informationGaining From Your Own Default
Ganng From Your Own Defaul Jon Gregory jon@ofranng.com Jon Gregory (jon@ofranng.com), Quan ongress US, 14 h July 2010 page 1 Regulaon s Easy () Wha don lke as a regulaor? Dfferen nsuons valung asses dfferenly
More informationCan Multivariate GARCH Models Really Improve Value-at-Risk Forecasts?
2s Inernaonal Congress on Modellng and Smulaon, Gold Coas, Ausrala, 29 ov o 4 Dec 205 www.mssanz.org.au/modsm205 Can Mulvarae GARCH Models Really Improve Value-a-Rsk Forecass? C.S. Sa a and F. Chan a a
More informationEconomics of taxation
Economcs of axaon Lecure 3: Opmal axaon heores Salane (2003) Opmal axes The opmal ax sysem mnmzes he excess burden wh a gven amoun whch he governmen wans o rase hrough axaon. Opmal axes maxmze socal welfare,
More informationCointegration between Fama-French Factors
1 Conegraon beween Fama-French Facors Absrac Conegraon has many applcaons n fnance and oher felds of scence researchng me seres and her nerdependences. The analyss s a useful mehod o analyse non-conegraon
More informationTime-Varying Correlations Between Credit Risks and Determinant Factors
me-varyng Correlaons Beween Cred Rsks and Deermnan Facors Frs & Correspondng Auhor: Ju-Jane Chang Asssan Professor n he Deparmen of Fnancal Engneerng and Acuaral Mahemacs, Soochow Unversy, awan 56, Sec.
More informationA Hybrid Method to Improve Forecasting Accuracy Utilizing Genetic Algorithm An Application to the Data of Operating equipment and supplies
A Hyrd Mehod o Improve Forecasng Accuracy Ulzng Genec Algorhm An Applcaon o he Daa of Operang equpmen and supples Asam Shara Tax Corporaon Arkne, Shzuoka Cy, Japan, e-mal: a-shara@arkne.nfo Dasuke Takeyasu
More informationCentre for Computational Finance and Economic Agents WP Working Paper Series. Amadeo Alentorn Sheri Markose
Cenre for Compuaonal Fnance and Economc Agens WP002-06 Workng Paper Seres Amadeo Alenorn Sher Markose Removng maury effecs of mpled rsk neural denses and relaed sascs February 2006 www.essex.ac.uk/ccfea
More informationMethodology of the CBOE S&P 500 PutWrite Index (PUT SM ) (with supplemental information regarding the CBOE S&P 500 PutWrite T-W Index (PWT SM ))
ehodology of he CBOE S&P 500 PuWre Index (PUT S ) (wh supplemenal nformaon regardng he CBOE S&P 500 PuWre T-W Index (PWT S )) The CBOE S&P 500 PuWre Index (cker symbol PUT ) racks he value of a passve
More informationEmpirical Study on the Relationship between ICT Application and China Agriculture Economic Growth
Emprcal Sudy on he Relaonshp beween ICT Applcaon and Chna Agrculure Economc Growh Pengju He, Shhong Lu, Huoguo Zheng, and Yunpeng Cu Key Laboraory of Dgal Agrculural Early-warnng Technology Mnsry of Agrculure,
More informationMORNING SESSION. Date: Wednesday, May 4, 2016 Time: 8:30 a.m. 11:45 a.m. INSTRUCTIONS TO CANDIDATES
SOCIETY OF ACTUARIES Exam QFICORE MORNING SESSION Dae: Wednesday, May 4, 016 Tme: 8:30 a.m. 11:45 a.m. INSTRUCTIONS TO CANDIDATES General Insrucons 1. Ths examnaon has a oal of 100 pons. I consss of a
More informationMichał Kolupa, Zbigniew Śleszyński SOME REMARKS ON COINCIDENCE OF AN ECONOMETRIC MODEL
M I S C E L L A N E A Mchał Kolupa, bgnew Śleszyńsk SOME EMAKS ON COINCIDENCE OF AN ECONOMETIC MODEL Absrac In hs paper concep of concdence of varable and mehods for checkng concdence of model and varables
More informationRecall from last time. The Plan for Today. INTEREST RATES JUNE 22 nd, J u n e 2 2, Different Types of Credit Instruments
Reall from las me INTEREST RATES JUNE 22 nd, 2009 Lauren Heller Eon 423, Fnanal Markes Smple Loan rnpal and an neres paymen s pad a maury Fxed-aymen Loan Equal monhly paymens for a fxed number of years
More informationUsing Fuzzy-Delphi Technique to Determine the Concession Period in BOT Projects
Usng Fuzzy-Delph Technque o Deermne he Concesson Perod n BOT Projecs Khanzad Mosafa Iran Unversy of Scence and Technology School of cvl engneerng Tehran, Iran. P.O. Box: 6765-63 khanzad@us.ac.r Nasrzadeh
More informationSTOCK PRICES TEHNICAL ANALYSIS
STOCK PRICES TEHNICAL ANALYSIS Josp Arnerć, Elza Jurun, Snježana Pvac Unversy of Spl, Faculy of Economcs Mace hrvaske 3 2 Spl, Croaa jarnerc@efs.hr, elza@efs.hr, spvac@efs.hr Absrac Ths paper esablshes
More informationAre Taxes Capitalized in Bond Prices? Evidence from the Market for Government of Canada Bonds* Stuart Landon **
PRELIINARY DRAFT Are Taxes Capalzed n Bond Prces? Evdence from he arke for Governmen of Canada Bonds* Suar Landon ** Deparmen of Economcs Unversy of Albera Edmonon, Albera Canada T6G 2H4 14 ay 2008 Absrac
More informationAssessing Long-Term Fiscal Dynamics: Evidence from Greece and Belgium
Inernaonal Revew of Busness Research Papers Vol. 7. No. 6. November 2011. Pp. 33-45 Assessng Long-Term Fscal Dynamcs: Evdence from Greece and Belgum JEL Codes: Ε62 and Η50 1. Inroducon Evangela Kasma 1,2
More informationA Hybrid Method for Forecasting with an Introduction of a Day of the Week Index to the Daily Shipping Data of Sanitary Materials
Journal of Communcaon and Compuer (05) 0-07 do: 0.765/548-7709/05.0.00 D DAVID PUBLISHING A Hyrd Mehod for Forecasng wh an Inroducon of a Day of he Week Inde o he Daly Shppng Daa of Sanary Maerals Dasuke
More informationESSAYS ON MONETARY POLICY AND INTERNATIONAL TRADE. A Dissertation HUI-CHU CHIANG
ESSAYS ON MONETARY POLICY AND INTERNATIONAL TRADE A Dsseraon by HUI-CHU CHIANG Submed o he Offce of Graduae Sudes of Texas A&M Unversy n paral fulfllmen of he requremens for he degree of DOCTOR OF PHILOSOPHY
More informationCorrelation Smile, Volatility Skew and Systematic Risk Sensitivity of Tranches
Correlaon Smle Volaly Skew and Sysemac Rsk Sensvy of ranches Alfred Hamerle Andreas Igl and lan Plank Unversy of Regensburg ay 0 Absac he classcal way of eang he correlaon smle phenomenon wh cred ndex
More informationMixtures of Normal Distributions: Application to Bursa Malaysia Stock Market Indices
World Appled Scences Journal 6 (6): 78-790, 0 ISSN 88-495 IDOSI Publcaons, 0 Mxures of Normal Dsrbuons: Applcaon o Bursa Malaysa Sock Marke Indces Zey An Kamaruzzaman, Zad Isa and Mohd ahr Ismal School
More informationPermanent Income and Consumption
roceedngs of 30h Inernaonal onference Mahemacal Mehods n Economcs ermanen Income and onsumpon Václava ánková 1 Absrac. A heory of consumer spendng whch saes ha people wll spend money a a level conssen
More informationOnline Data, Fixed Effects and the Construction of High-Frequency Price Indexes
Onlne Daa, Fxed Effecs and he Consrucon of Hgh-Frequency Prce Indexes Jan de Haan* and Rens Hendrks** * ascs eherlands / Delf Unversy of Technology ** ascs eherlands EMG Worksho 23 Ams of he aer Exlan
More informationPrediction of Oil Demand Based on Time Series Decomposition Method Nan MA * and Yong LIU
2017 2nd Inernaonal Conference on Sofware, Mulmeda and Communcaon Engneerng (SMCE 2017) ISBN: 978-1-60595-458-5 Predcon of Ol Demand Based on Tme Seres Decomposon Mehod Nan MA * and Yong LIU College of
More informationOnline appendices from The xva Challenge by Jon Gregory. APPENDIX 14A: Deriving the standard CVA formula.
Onlne appendces fro he xa Challenge by Jon Gregory APPNDX 4A: Dervng he sandard CA forla We wsh o fnd an expresson for he rsky vale of a need se of dervaves posons wh a ax ary dae Denoe he rsk-free vale
More informationSUMMARY INTRODUCTION. Figure 1: An illustration of the integration of well log data and seismic data in a survey area. Seismic cube. Well-log.
Perophyscal propery esmaon from sesmc daa usng recurren neural neworks Moaz Alfarraj, and Ghassan AlRegb Cener for Energy and Geo Processng CeGP, Georga Insue of Technology SUMMARY Reservor characerzaon
More informationQuarterly Accounting Earnings Forecasting: A Grey Group Model Approach
Quarerly Accounng Earnngs Forecasng: A Grey Group Model Approach Zheng-Ln Chen Deparmen of Accounng Zhongnan Unversy of Economcs and Law # Souh Nanhu Road, Wuhan Cy, 430073 Hube People's Republc of Chna
More informationReturn Calculation Methodology
Reurn Calculaon Mehodology Conens 1. Inroducon... 1 2. Local Reurns... 2 2.1. Examle... 2 3. Reurn n GBP... 3 3.1. Examle... 3 4. Hedged o GBP reurn... 4 4.1. Examle... 4 5. Cororae Acon Facors... 5 5.1.
More informationThe Keynesian micro-foundations of the business cycle: some implications of globalisation
The Keynesan mcro-foundaons of he busness cycle: some mplcaons of globalsaon Paul Ormerod, Volerra Consulng Ld., London e-mal pormerod@volerra.co.uk el: 44 0208 878 6333 I am graeful o Rod Gbson n parcular
More informationHoldings-based Fund Performance Measures: Estimation and Inference 1
1 Holdngs-based Fund Performance Measures: Esmaon and Inference 1 Wayne E. Ferson Unversy of Souhern Calforna and NBER Junbo L. Wang Lousana Sae Unversy Aprl 14, 2018 Absrac Ths paper nroduces a panel
More information(Version ) I. SOURCE, DEFINITIONS AND APPLICABILITY
Threh meeng Repor Page 1 1. Sources Draf afforesaon and reforesaon baselne and monorng mehodology AR-AM00XX Afforesaon or reforesaon of degraded or abandoned agrculural lands (Verson 01.0.0) I. SOURE,
More informationConditional Skewness of Aggregate Market Returns: Evidence from Developed and Emerging Markets
Global Economy and Fnance Journal Vol. 7. No.. March 04. Pp. 96 Condonal Skewness of Aggregae Marke Reurns: Evdence from Developed and Emergng Markes Anchada Charoenrook and Hazem Daouk Ths paper examnes
More informationNumerical Evaluation of European Option on a Non Dividend Paying Stock
Inernaonal Journal of Compuaonal cence and Mahemacs. IN 0974-389 olume Number 3 (00) pp. 6--66 Inernaonal Research Publcaon House hp://www.rphouse.com Numercal Evaluaon of European Opon on a Non Dvdend
More informationOption-Implied Currency Risk Premia
Opon-Impled Currency Rsk Prema Jakub W. Jurek and Zhka Xu Absrac We use cross-seconal nformaon on he prces of G10 currency opons o calbrae a non-gaussan model of prcng kernel dynamcs and consruc esmaes
More informationA Neural Network Approach to Time Series Forecasting
A Neural Nework Approach o Tme Seres Forecasng Iffa A. Gheyas, Lesle S. Smh Absrac We propose a smple approach for forecasng unvarae me seres. The proposed algorhm s an ensemble learnng echnque ha combnes
More informationThe Proposed Mathematical Models for Decision- Making and Forecasting on Euro-Yen in Foreign Exchange Market
Iranan Economc Revew, Vol.6, No.30, Fall 20 The Proposed Mahemacal Models for Decson- Makng and Forecasng on Euro-Yen n Foregn Exchange Marke Abdorrahman Haer Masoud Rabban Al Habbna Receved: 20/07/24
More informationAdjusted-Productivity Growth for Resource Rents: Kuwait Oil Industry
Appled Economcs and Fnance Vol. 3, No. 2; May 2016 ISSN 2332-7294 E-ISSN 2332-7308 Publshed by Redfame Publshng URL: hp://aef.redfame.com Adjused-Producvy Growh for Resource Rens: Kuwa Ol Indusry 1 Acng
More informationFinancial Innovation and Asset Price Volatility. Online Technical Appendix
Fnancal Innovaon and Asse Prce Volaly Onlne Techncal Aendx Felx Kubler and Karl Schmedders In hs echncal aendx we derve all numbered equaons dslayed n he aer Equaons For he wo models n he aer, he frs se
More informationModelling Inflation Rate Volatility in Kenya Using Arch -Type Model Family
Research Journal of Fnance and Accounng ISSN -697 (Paper) ISSN -847 (Onlne) Vol.7, No.3, 6 www.se.org Modellng Inflaon Rae Volaly n Kenya Usng Arch -Type Model Famly Johnson Okeyo Mwank Ivv Phlp Ngare.School
More informationSTOCHASTIC LOCAL VOLATILITY
STOCHASTIC OCA VOATIITY Carol Alexander Char of Rsk Managemen and Drecor of Research ISMA Cenre, Busness School, The Unversy of Readng PO Box 4, Readng, RG6 6BA Uned Kngdom c.alexander@smacenre.rdg.ac.uk
More informationVolatility Modeling for Forecasting Stock Index with Fixed Parameter Distributional Assumption
Journal of Appled Fnance & Banng, vol. 3, no. 1, 13, 19-1 ISSN: 179-5 (prn verson), 179-599 (onlne) Scenpress Ld, 13 Volaly Modelng for Forecasng Soc Index wh Fxed Parameer Dsrbuonal Assumpon Md. Mosafzur
More informationExplaining Product Release Planning Results Using Concept Analysis
Explanng Produc Release Plannng Resuls Usng Concep Analyss Gengshen Du, Thomas Zmmermann, Guenher Ruhe Deparmen of Compuer Scence, Unversy of Calgary 2500 Unversy Drve NW, Calgary, Albera T2N 1N4, Canada
More informationTHE TYRANNY OF THE IDENTITY: GROWTH ACCOUNTING REVISITED
THE TYRANNY OF THE IDENTITY: GROWTH ACCOUNTING REVISITED Jesus Felpe Economcs and Research Deparmen Asan Developmen Bank Manla (Phlppnes) e-mal: jfelpe@adb.org JSL McCombe Cenre for Economc and Publc Polcy
More informationMeasurement of EU agricultural total factor productivity growth
Measuremen of EU agrculural oal facor producvy groh Koen Mondelaers DG for Agrculure and Rural Developmen Un E.2 Economc modellng and Oulook General objecves for he CAP,2 : Relevancy Promong a vable food
More informationProperty of stocks and wealth effects on consumption
Propery of socks and wealh effecs on consumpon RICARDO M. SOUSA Unversy of Mnho Deparmen of Economcs Campus of Gualar, 470-057 - BRAGA PORTUGAL E-mal: rjsousa@eeg.umnho.p March 2003 Absrac Recen flucuaons
More informationA New N-factor Affine Term Structure Model of Futures Price for CO 2 Emissions Allowances: Empirical Evidence from the EU ETS
WSEAS RASACIOS on BUSIESS and ECOOMICS Ka Chang, Su-Sheng Wang, Je-Mn Huang A ew -facor Affne erm Srucure Model of Fuures Prce for CO Emssons Allowances: Emprcal Evdence from he EU ES KAI CHAG, SU-SHEG
More informationConvertible Bonds and Stock Liquidity. Author. Published. Journal Title DOI. Copyright Statement. Downloaded from. Griffith Research Online
Converble Bonds and Sock Lqudy Auhor Wes, Jason Publshed 2012 Journal Tle Asa-Pacfc Fnancal Markes DOI hps://do.org/10.1007/s10690-011-9139-3 Copyrgh Saemen 2011 Sprnger Japan. Ths s an elecronc verson
More informationThe Selection Ability of Italian Mutual Fund. By Valter Lazzari and Marco Navone
The Selecon Ably of Ialan Muual Fund By Valer Lazzar and Marco Navone Workng Paper N. 1/3 Ocober 23 THE SELECTION ABILITY OF ITALIAN MUTUAL FUND MANAGERS By Valer Lazzar Professor of Bankng and Fnance
More informationComparing Sharpe and Tint Surplus Optimization to the Capital Budgeting Approach with Multiple Investments in the Froot and Stein Framework.
Comparng Sharpe and Tn Surplus Opmzaon o he Capal Budgeng pproach wh Mulple Invesmens n he Froo and Sen Framework Harald Bogner Frs Draf: Sepember 9 h 015 Ths Draf: Ocober 1 h 015 bsrac Below s shown ha
More informationTHE APPLICATION OF REGRESSION ANALYSIS IN TESTING UNCOVERED INTEREST RATE PARITY
QUANTITATIVE METHOD IN ECONOMIC Vol. XIV, No., 03, pp. 3 4 THE APPLICATION OF REGREION ANALYI IN TETING UNCOVERED INTERET RATE PARITY Joanna Kselńsa, Kaarzyna Czech Faculy of Economcs cences Warsaw Unversy
More informationReal-Time Forecasting Revisited: Letting the Data Decide
Real-Tme Forecasng Revsed: Leng he Daa Decde Jackson Kchen John Kchen Revsed Sepember 2012 Orgnal verson June 2012 Paper prepared for he Naonal Assocaon for Busness Economcs 2012 Menns Award for presenaon
More informationCryptographic techniques used to provide integrity of digital content in long-term storage
RB/3/2011 Crypographc echnques used o provde negry of dgal conen n long-erm sorage REPORT ON THE PROBLEM Problem presened by Marn Šmka Paweł Wojcechowsk Polsh Secury Prnng Works (PWPW) 1 Repor auhors Małgorzaa
More informationDEPARTMENT OF ECONOMETRICS AND BUSINESS STATISTICS
ISSN 440-77X AUSTRALIA DEPARTMENT OF ECONOMETRICS AND BUSINESS STATISTICS Assocaon beween Markov regme-swchng marke volaly and bea rsk: Evdence from Dow Jones ndusral secures Don U.A. Galagedera and Roland
More informationThe Empirical Research of Price Fluctuation Rules and Influence Factors with Fresh Produce Sequential Auction Limei Cui
6h Inernaonal Conference on Sensor Nework and Compuer Engneerng (ICSNCE 016) The Emprcal Research of Prce Flucuaon Rules and Influence Facors wh Fresh Produce Sequenal Aucon Lme Cu Qujng Normal Unversy,
More informationForecasting Value-at-Risk under Different Distributional Assumptions
Arcle Forecasng Value-a-Rsk under Dfferen Dsrbuonal Assumpons Manuela Braone 1, and Ncolas K. Scholes 1,,, Receved: 1 July 1; Acceped: 1 December 1; Publshed: 11 January 16 Academc Edors: Fredj Jawad,
More informationImproving Earnings per Share: An Illusory Motive in Stock Repurchases
Inernaonal Journal of Busness and Economcs, 2009, Vol. 8, No. 3, 243-247 Improvng Earnngs per Share: An Illusory Move n Sock Repurchases Jong-Shn We Deparmen of Inernaonal Busness Admnsraon, Wenzao Ursulne
More informationConvexity Adjustments in Inflation linked Derivatives using a multi-factor version of the Jarrow and Yildirim (2003) Model
Imperal College of Scence, echnology and edcne Deparmen of ahemacs Convexy Adjusmens n Inflaon lnked Dervaves usng a mul-facor verson of he Jarrow and Yldrm (003 odel Hongyun L Sepember 007 Submed o Imperal
More information1 Purpose of the paper
Moneary Economics 2 F.C. Bagliano - Sepember 2017 Noes on: F.X. Diebold and C. Li, Forecasing he erm srucure of governmen bond yields, Journal of Economerics, 2006 1 Purpose of he paper The paper presens
More informationThe Net Benefit to Government of Higher Education: A Balance Sheet Approach
The Ne Benef o Governmen of Hgher Educaon: A Balance Shee Approach Davd Johnson and Roger Wlkns Melbourne Insue of Appled Economc and Socal Research The Unversy of Melbourne Melbourne Insue Workng Paper
More informationBusiness cycle, credit risk and economic capital determination by commercial banks
Busness cycle, cred rsk and economc capal deermnaon by commercal banks Alexs Dervz and Narcsa Kadlčáková 1 Czech Naonal Bank 1. Inroducon Regular assessmens of he defaul rsk of bank clens and esmaons of
More information