Application of ARIMA Model for River Discharges Analysis
|
|
- Cleopatra Houston
- 6 years ago
- Views:
Transcription
1 Alcaon of ARIMA Model for Rver Dscharges Analyss Bhola NS Ghmre Journal of Neal Physcal Socey Volume 4, Issue 1, February 17 ISSN: X Edors: Dr. Go Chandra Kahle Dr. Devendra Adhkar Mr. Deeendra Parajul JNPS, 4 (1), 7-3 (17) Publshed by: Neal Physcal Socey P.O. Box : 934 Tr-Chandra Camus Kahmandu, Neal Emal: nsedor@gmal.com
2 JNPS 4 (1), 7-3 (17) ISSN: X Neal Physcal Socey Research Arcle Alcaon of ARIMA Model for Rver Dscharges Analyss Bhola NS Ghmre Dearmen of Cvl Engneerng, Pulchowk Camus, Insue of Engneerng, T. U., Lalur, Neal Corresondng Emal: bholag@oe.edu.n ABSTRACT Tme seres daa ofen arse when monorng hydrologcal rocesses. Mos of he hydrologcal daa are me relaed and drecly or ndrecly her analyss relaed wh me comonen. Tme seres analyss accouns for he fac ha daa ons aken over me may have an nernal srucure (such as auocorrelaon, rend or seasonal varaon) ha should be accouned for. Many mehods and aroaches for formulang me seres forecasng models are avalable n leraure. Ths sudy wll gve a bref overvew of auo-regressve negraed movng average (ARIMA) rocess and s alcaon o forecas he rver dscharges for a rver. The develoed ARIMA model s esed successfully for wo hydrologcal saons for a rver n US. Keywords: Tme Seres Analyss, ARIMA Model, Hydrologcal Process, Auocorrelaon, Seasonal Varaon. INTRODUCTION Auo-Regressve Inegraed Movng Average (ARIMA) mehod s wdely used n feld of me seres modelng and analyss. These models were descrbed by Box and Jenkns (1976) and furher dscussed by Waler (Chafeld, 1996). The Box- Jenkns aroach n hydrologcal modelng s used by several researchers. Chew e al. (1993) conduced a comarson of sx ranfall-runoff modelng aroaches o smulae daly, monhly and annual flows n egh unregulaed cachmens. Langu (1993) used me seres analyss o deec he changes n ranfall and runoff aerns. Kuo and Sun (1993) used he me seres model for en days sream flow forecas and generae synhess hydrograh caused by yhoons n Tanshu Rver n Tawan. Nall and Moman (9) used he me seres analyss for ranfall daa n Jordan. Ths aer s amed o show he usefulness of hs oular echnue ARIMA for a ycal case sudy. ARIMA MODEL Trend and redcon of me seres can be comued by usng ARIMA model. ARIMA (,d,) model s a comlex lnear model. In sascs, normally n me seres analyss, ARIMA model s generalzaon of auoregressve movng average (ARMA) models, somemes called Box-Jenkns models afer he erave Box-Jenkns mehodology. Gven a me seres of daa X, he ARMA model s a ool for undersandng and, erhas, redcng fuure values n hs seres. The model consss of wo ars, an auoregressve (AR) ar and a movng average (MA) ar. And ha of n ARIMA model he hrd ar negraed (I) ncluded. The model s usually hen referred o as he ARIMA (,d,) model where s he order of he auoregressve ar, d s he order of non seasonal dfferences and s he order of he movng average ar. The noaon AR() refers o he auoregressve model of order, whch can be wren as: (1) X c X 1 where φ 1,.., φ are he arameers of he model, c s a consan and ε s whe nose. The consan erm s omed by many auhors for smlcy. Smlarly, he noaon MA () refers o he movng average model of order. Ths can be wren as: X 1 () Where he θ 1,..,θ are he arameers of he model, μ s he execaon of X (ofen assumed eual o ), and he ε, ε -1,, ε - are he whe nose error erms. The movng average model s a fne mulse resonse fler wh some addonal nerreaon laced on. Whle combng hese wo models, he ARMA (,) s obaned. 7
3 Alcaon of ARIMA Model for Rver Dscharges Analyss X c X 1 1 (3) The error erm ε are generally assumed o be ndeenden dencally-dsrbued random varables samled from a normal dsrbuon wh zero mean: ε ~ N(,σ ) where, σ s he varance. Some researchers have used he euaon n a lag oeraor form. In lag oeraor form, AR () model s gven by- 1 L X X (4) 1 And MA() model s gven by- X 1 L (5) 1 Where, φ and θ are defned by he arameers conanng nsde he arenhess of each model. Combnng hese models we can manulae and wre n he followng form. 1 L X 1 L (6) 1 1 Assume ha he olynomal of frs erm of above euaon has a unary roo of mullcy d. Then hs euaon can be udaed ncludng he dfference erm, whch can be exressed as, 1 d d L ) 1 L (1 L) (7) 1 1 An ARIMA (, d, ) rocess exresses hs olynomal facorzaon roery, and s fnally wren as: d 1 L ) (1 L) 1 L (8) 1 1 More recsely, he ARIMA (, d, ) model can be wren as: d ( B)(1 B) y ( B) (9) The model whle used wh seasonal flucuaon, wh seasonal lengh s, he rocess s called SARIMA (, d, )(P, D, Q) s, where, d, reresens he order of rocess AR, order of dfference (I) and order of rocess (MA) for non seasonal ar and P, D, Q, reresens he order of seasonal rocess AR, order of seasonal dfference and order of seasonal MA and s s he lengh of seasonal erod. The general euaon of SARIMA model s: d s s D s ( B)(1 B) P( B )(1 B ) y ( B) Q( B ) (1) Where, φ (B) s auo regressve oeraor, θ(b) s he oeraor of movng average; ΦP(Bs) s seasonal auoregressve oeraor, ΘQ(Bs) s seasonal oeraor of movng averages, ε s whe nose. STATISTICAL TESTS FOR MODEL PERFORMANCE There are several sascal ess for model erformance. In hs sudy some of hese ess are used whch are easy o undersand and use. Coeffcen of Correlaon A very moran ar of sascs s descrbng he relaonsh beween wo (or more) varables. One of he mos fundamenal conces n research s he conce of correlaon. If wo varables are correlaed, hs means ha can use nformaon abou one varable o redc he values of he oher varable. The coeffcen of correlaon s gven by followng euaon. r ( ( x x) ( y y) (11) x x) ( y y) Where, r s correlaon coeffcen; x and y are ndeenden (observed) and deenden (redced) varables and x and y are her corresondng means. Roo Mean Suare Error The roo mean suare error (RMSE)) (also roo mean suare devaon (RMSD)) s a freuenlyused measure of he dfferences beween values redced by a model or an esmaor and he values acually observed from he hng beng modeled or esmaed. RMSE s a good measure of accuracy. These ndvdual dfferences are also called resduals, and he RMSE serves o aggregae hem no a sngle measure of redcve ower. The mahemacal form of RMSE s gven by: x y RMSE (1) N Where N s oal number of daa se and oher varables are same as earler euaon. These wo 8
4 Bhola NS Ghmre euaons are used for he comarson of observed and redced values. AIC and BIC crera The wo mos commonly used enalzed model selecon crera, he Akake s nformaon creron (AIC) and he Bayesan nformaon creron (BIC), are examned and comared for ARIMA model selecon. AIC - In general case, AIC k nln( SSE / n) (13) Where k s he number of arameers n he sascal model, n s he number of observaons and SSE s suare sum of error gven by - SSE n (14) 1 BIC - In general BIC s gven by- BIC k ln( n) nln( SSE / n) (15) The mnmum values of hese AIC and BIC crera gve he beer model erformance. CASE STUDY DESCRIPTION Dscharge Daa For he alcaon demonsraon of ARIMA model, he me seres daly dscharge daa wo saons n Schuylkll Rver a Berne (Saon no: 1475, La. 4 º 31'1'' and Long. 75 º 59'55'') and Phladelha (Saon no: 14745, La. 39 º 58'4'' and Long. 75 º 11'''), USA are aken. The cachmens area of Berne saon s abou km and ha of Phladelha saon s km. Ths nformaon was obaned from USGS webse. The daa from he erod Ocober 1, o Seember 3, 6 were aken for boh of he saons. Inal all sx years daa were aken for ARIMA model develomen and fnally usng model las one year daa (Ocober 1, 6 o Seember 3, 7) were redced for boh he saons. Some of he sascal arameers for hese ses are shown n Table 1 of he dscharge daa. The arameers μ, σ, σ/μ, C sk, C kr, X max, X mn are mean, sandard devaon, varance, skew-ness, kuross, maxmum and mnmum values resecvely. The dscharge lms of Berne saon are.13 o 97.1 m 3 /s and ha of Phladelha saon are.4 o m 3 /s. Table 1. The daly sascal arameers for Schuylkll Rver. Saon Basn Area (Km ) μ σ σ/μ C sk C kr X max X mn Berne 1475 Phladelha Develomen of ARIMA Models From he me seres lo for he gven daa (fgure 1), can be observed ha here s no seasonaly for daly daa. In fac, s very dffcul o fx he seasonaly for daly daa and due o he large san of me (365 days), s unrelable oo. So, he model formulaon has done whou seasonaly. The ARIMA models for he boh saons are develoed by usng SPSS. Inally, he several models were esed based on he AIC and BIC creron. The AIC and BIC values for few models for hese are gven n he followng Table. Table. The AIC and BIC values for some esng ARIMA models for Berne daa and Phladelha daa. ARIMA - Berne Daa (1,,) (1,1,) (,,) (1,,1) (1,1,1) (1,1,) (1,,1) (1,,) (,,) AIC BIC ARIMA Phladelha Daa (,1,) (,,) (1,,) (1,1,) (1,1,1) (,11) (,,1) AIC BIC
5 Alcaon of ARIMA Model for Rver Dscharges Analyss From Table, can be judged ha on he rncle of AIC and BIC es, ARIMA (1, 1, ) s suable for Berne saon and ARIMA (1, 1, ) s suable for Phladelha saon. However whle we observed he correlaon marx of he ARIMA arameers for Berne Saon from Table 3, he arameers MA(1) and MA() has very hgh correlaon aroachng o uny. So ha her effecs n ARIMA model are neglgble and we can reduce hs MA arameer. Then he ARIMA (1, 1, 1) s roosed for he furher analyss even hough ARIMA (1, 1, ) has farly less AIC and BIC values. The correlaon marx and arameers for fnal model for he Berne saon are gven n he followng Table 4. Table 3. The correlaon marx for arameers of Berne daa. ARIMA (1,1,) Non-Seasonal φ 1 θ 1 θ Consan Non-Seasonal φ θ θ Consan 1. Table 4. The correlaon marx and fnal arameers for Berne Saon. a. The correlaon marx b. The fnal arameers ARIMA (1,1,1) Non-Seasonal Non-Seasonal φ 1 θ 1 φ θ Consan 1. Consan ARIMA (1, 1, 1) Esmaes Sd Error Non-Seasonal φ θ Consan Smlarly, for Phladelha saon, AIC and BIC values for some models, correlaon marx for seleced model and fnal ARIMA model arameers are gven n Table 5. Table 5. The correlaon marx and fnal arameers for Phladelha daa. a. The correlaon marx b. The fnal arameers ARIMA (,1,1) Non- Seasonal Non-Seasonal φ 1 φ θ 1 φ φ θ Consan 1. Consan ARIMA (,1,1) Esmaes Sd Error Non-Seasonal φ φ θ Consan As dscussed n earler, for ARIMA (1, 1, 1), he general euaon can be reduced as: ( B)(1 B) y ( B) (16) Whle subsung he model arameers n he above euaon and smlfy, we ge he fnal model for Berne saon as: y y.698 y.991e.1 (17) 1 1 Smlarly, for ARIMA (, 1, 1) model, he general euaon s gven by: ( B)(1 B) y ( B) (18) 1 1 3
6 Predced Dscharge (Mm 3 ) Predced Dscharge (Mm 3 ) Dscharge (Mm 3 ) Dscharge (Mm 3 ) Bhola NS Ghmre Based on he model arameers, he fnal model for Phladelha saon s gven by: y y.98 y.141 y.968 e.9 (19) RESULT AND DISCUSSIONS The me seres lo of orgnal daa and redced daa from he fnal models for boh he saons are gven n fgure 1. Ths shows ha he general rend s followed by redced daa o ha of observed daa. In he lo, he sx years observed daa followed by las one year redced daa for he boh saons. 1 5 Observed Predced Berne Observed Predced Phladelha Tme n Days Tme n Days Fg. 1. Tme seres lo of observed and redced daly dscharge a Berne and Phladelha Saons. The scaer los of hese wo saons are also gven n fgure. Ths shows ha hey are ue good models. The R for Berne saon and Phladelha saons are:.991 and.9688 resecvely. 7 Berne 1 Phladelha 6 R =.99 1 R = Observed Dscharge (Mm 3 ) Observed Dscharge (Mm 3 ) Fg.. Scaer lo of observed and redced daly dscharge a Berne and Phladelha Saons. CONCLUSIONS Tme seres analyss for he rver dscharge shows ha s an moran ool for modelng and forecasng. ARIMA (1, 1, 1) model s fed for Berne saon and ARIMA (, 1, 1) s fed for Phladelha saon. Boh he saons are les n he same rvers bu hey have dfferen cachmens coverage. So should be noed ha even he rver s same, deendng uon he cachmens characerscs, alcable models are dfferen 31
7 Alcaon of ARIMA Model for Rver Dscharges Analyss ndvdual ses. The coeffcen of deermnaons (.99 for Berne and.969 for Phladelha) shows ha he model s useful for runoff forecasng. REFERENCES: Box, G. E. P., and Jenkns, G. M. (1976). Tme seres analyss: forecasng and conrol. Revsed Edn. Holden-Day, San Francsco.... Chafeld, C. (1996). Analyss of me seres: an nroducon. 5 h ed. Chaman and hall, Boca Raon. Chew, F. H. S.; Sewardson, M. J., and McMahon, T. A. (1993). Comarson of sx ranfallrunoff modellng aroaches. J. of Hydrology, 147(1-4): Langu, E. M. (1993). Deecon of changes n ranfall and runoff aerns. J. of Hydrology, 147 (1-4): Kuo, J. T., and Sun, Y. H. (1993). An nervenon model for average 1 day sream flow forecas and synhess. J. of Hydrology, 151(1): Nal, P. E., and Moman, M. (9). Tme seres analyss model for ranfall daa n Jrdan: case sudy for usng Tme seres analyss. Amercan J. of Envronmenal Scences, 5(5):
Advanced time-series analysis (University of Lund, Economic History Department)
Advanced me-seres analss (Unvers of Lund, Economc Hsor Dearmen) 3 Jan-3 Februar and 6-3 March Lecure 4 Economerc echnues for saonar seres : Unvarae sochasc models wh Box- Jenns mehodolog, smle forecasng
More informationNPTEL Project. Econometric Modelling. Module23: Granger Causality Test. Lecture35: Granger Causality Test. Vinod Gupta School of Management
P age NPTEL Proec Economerc Modellng Vnod Gua School of Managemen Module23: Granger Causaly Tes Lecure35: Granger Causaly Tes Rudra P. Pradhan Vnod Gua School of Managemen Indan Insue of Technology Kharagur,
More informationTHEORETICAL AUTOCORRELATIONS. ) if often denoted by γ. Note that
THEORETICAL AUTOCORRELATIONS Cov( y, y ) E( y E( y))( y E( y)) ρ = = Var( y) E( y E( y)) =,, L ρ = and Cov( y, y ) s ofen denoed by whle Var( y ) f ofen denoed by γ. Noe ha γ = γ and ρ = ρ and because
More informationRELATIONSHIP BETWEEN VOLATILITY AND TRADING VOLUME: THE CASE OF HSI STOCK RETURNS DATA
RELATIONSHIP BETWEEN VOLATILITY AND TRADING VOLUME: THE CASE OF HSI STOCK RETURNS DATA Mchaela Chocholaá Unversy of Economcs Braslava, Slovaka Inroducon (1) one of he characersc feaures of sock reurns
More informationDynamic Regressions with Variables Observed at Different Frequencies
Dynamc Regressons wh Varables Observed a Dfferen Frequences Tlak Abeysnghe and Anhony S. Tay Dearmen of Economcs Naonal Unversy of Sngaore Ken Rdge Crescen Sngaore 96 January Absrac: We consder he roblem
More informationLecture 18: The Laplace Transform (See Sections and 14.7 in Boas)
Lecure 8: The Lalace Transform (See Secons 88- and 47 n Boas) Recall ha our bg-cure goal s he analyss of he dfferenal equaon, ax bx cx F, where we emloy varous exansons for he drvng funcon F deendng on
More informationEcon107 Applied Econometrics Topic 5: Specification: Choosing Independent Variables (Studenmund, Chapter 6)
Econ7 Appled Economercs Topc 5: Specfcaon: Choosng Independen Varables (Sudenmund, Chaper 6 Specfcaon errors ha we wll deal wh: wrong ndependen varable; wrong funconal form. Ths lecure deals wh wrong ndependen
More informationJanuary Examinations 2012
Page of 5 EC79 January Examnaons No. of Pages: 5 No. of Quesons: 8 Subjec ECONOMICS (POSTGRADUATE) Tle of Paper EC79 QUANTITATIVE METHODS FOR BUSINESS AND FINANCE Tme Allowed Two Hours ( hours) Insrucons
More informationData Collection Definitions of Variables - Conceptualize vs Operationalize Sample Selection Criteria Source of Data Consistency of Data
Apply Sascs and Economercs n Fnancal Research Obj. of Sudy & Hypoheses Tesng From framework objecves of sudy are needed o clarfy, hen, n research mehodology he hypoheses esng are saed, ncludng esng mehods.
More information( ) () we define the interaction representation by the unitary transformation () = ()
Hgher Order Perurbaon Theory Mchael Fowler 3/7/6 The neracon Represenaon Recall ha n he frs par of hs course sequence, we dscussed he chrödnger and Hesenberg represenaons of quanum mechancs here n he chrödnger
More informationEP2200 Queuing theory and teletraffic systems. 3rd lecture Markov chains Birth-death process - Poisson process. Viktoria Fodor KTH EES
EP Queung heory and eleraffc sysems 3rd lecure Marov chans Brh-deah rocess - Posson rocess Vora Fodor KTH EES Oulne for oday Marov rocesses Connuous-me Marov-chans Grah and marx reresenaon Transen and
More informationF-Tests and Analysis of Variance (ANOVA) in the Simple Linear Regression Model. 1. Introduction
ECOOMICS 35* -- OTE 9 ECO 35* -- OTE 9 F-Tess and Analyss of Varance (AOVA n he Smple Lnear Regresson Model Inroducon The smple lnear regresson model s gven by he followng populaon regresson equaon, or
More informationFI 3103 Quantum Physics
/9/4 FI 33 Quanum Physcs Aleander A. Iskandar Physcs of Magnesm and Phooncs Research Grou Insu Teknolog Bandung Basc Conces n Quanum Physcs Probably and Eecaon Value Hesenberg Uncerany Prncle Wave Funcon
More informationDepartment of Economics University of Toronto
Deparmen of Economcs Unversy of Torono ECO408F M.A. Economercs Lecure Noes on Heeroskedascy Heeroskedascy o Ths lecure nvolves lookng a modfcaons we need o make o deal wh he regresson model when some of
More informationRobustness Experiments with Two Variance Components
Naonal Insue of Sandards and Technology (NIST) Informaon Technology Laboraory (ITL) Sascal Engneerng Dvson (SED) Robusness Expermens wh Two Varance Componens by Ana Ivelsse Avlés avles@ns.gov Conference
More information[ ] 2. [ ]3 + (Δx i + Δx i 1 ) / 2. Δx i-1 Δx i Δx i+1. TPG4160 Reservoir Simulation 2018 Lecture note 3. page 1 of 5
TPG460 Reservor Smulaon 08 page of 5 DISCRETIZATIO OF THE FOW EQUATIOS As we already have seen, fne dfference appromaons of he paral dervaves appearng n he flow equaons may be obaned from Taylor seres
More informationSolution in semi infinite diffusion couples (error function analysis)
Soluon n sem nfne dffuson couples (error funcon analyss) Le us consder now he sem nfne dffuson couple of wo blocks wh concenraon of and I means ha, n a A- bnary sysem, s bondng beween wo blocks made of
More informationAdditive Outliers (AO) and Innovative Outliers (IO) in GARCH (1, 1) Processes
Addve Oulers (AO) and Innovave Oulers (IO) n GARCH (, ) Processes MOHAMMAD SAID ZAINOL, SITI MERIAM ZAHARI, KAMARULZAMMAN IBRAHIM AZAMI ZAHARIM, K. SOPIAN Cener of Sudes for Decson Scences, FSKM, Unvers
More informationV.Abramov - FURTHER ANALYSIS OF CONFIDENCE INTERVALS FOR LARGE CLIENT/SERVER COMPUTER NETWORKS
R&RATA # Vol.) 8, March FURTHER AALYSIS OF COFIDECE ITERVALS FOR LARGE CLIET/SERVER COMPUTER ETWORKS Vyacheslav Abramov School of Mahemacal Scences, Monash Unversy, Buldng 8, Level 4, Clayon Campus, Wellngon
More informationEndogeneity. Is the term given to the situation when one or more of the regressors in the model are correlated with the error term such that
s row Endogeney Is he erm gven o he suaon when one or more of he regressors n he model are correlaed wh he error erm such ha E( u 0 The 3 man causes of endogeney are: Measuremen error n he rgh hand sde
More informationAnalysis And Evaluation of Econometric Time Series Models: Dynamic Transfer Function Approach
1 Appeared n Proceedng of he 62 h Annual Sesson of he SLAAS (2006) pp 96. Analyss And Evaluaon of Economerc Tme Seres Models: Dynamc Transfer Funcon Approach T.M.J.A.COORAY Deparmen of Mahemacs Unversy
More informationTSS = SST + SSE An orthogonal partition of the total SS
ANOVA: Topc 4. Orhogonal conrass [ST&D p. 183] H 0 : µ 1 = µ =... = µ H 1 : The mean of a leas one reamen group s dfferen To es hs hypohess, a basc ANOVA allocaes he varaon among reamen means (SST) equally
More informationVariants of Pegasos. December 11, 2009
Inroducon Varans of Pegasos SooWoong Ryu bshboy@sanford.edu December, 009 Youngsoo Cho yc344@sanford.edu Developng a new SVM algorhm s ongong research opc. Among many exng SVM algorhms, we wll focus on
More informationOn One Analytic Method of. Constructing Program Controls
Appled Mahemacal Scences, Vol. 9, 05, no. 8, 409-407 HIKARI Ld, www.m-hkar.com hp://dx.do.org/0.988/ams.05.54349 On One Analyc Mehod of Consrucng Program Conrols A. N. Kvko, S. V. Chsyakov and Yu. E. Balyna
More information( t) Outline of program: BGC1: Survival and event history analysis Oslo, March-May Recapitulation. The additive regression model
BGC1: Survval and even hsory analyss Oslo, March-May 212 Monday May 7h and Tuesday May 8h The addve regresson model Ørnulf Borgan Deparmen of Mahemacs Unversy of Oslo Oulne of program: Recapulaon Counng
More informationApproximate Analytic Solution of (2+1) - Dimensional Zakharov-Kuznetsov(Zk) Equations Using Homotopy
Arcle Inernaonal Journal of Modern Mahemacal Scences, 4, (): - Inernaonal Journal of Modern Mahemacal Scences Journal homepage: www.modernscenfcpress.com/journals/jmms.aspx ISSN: 66-86X Florda, USA Approxmae
More informationEconometric Modelling of. Selected Approaches. Michaela Chocholatá University of Economics Bratislava
Economerc Modellng of Fnancal Tme Seres: Seleced Aroaches Mchaela Chocholaá Unversy of Economcs Braslava The man am of he resenaon a bref nroducon no he economerc modellng of he fnancal me seres ars: volaly
More informationA New Generalized Gronwall-Bellman Type Inequality
22 Inernaonal Conference on Image, Vson and Comung (ICIVC 22) IPCSIT vol. 5 (22) (22) IACSIT Press, Sngaore DOI:.7763/IPCSIT.22.V5.46 A New Generalzed Gronwall-Bellman Tye Ineualy Qnghua Feng School of
More informationImprovement in Estimating Population Mean using Two Auxiliary Variables in Two-Phase Sampling
Rajesh ngh Deparmen of ascs, Banaras Hndu Unvers(U.P.), Inda Pankaj Chauhan, Nrmala awan chool of ascs, DAVV, Indore (M.P.), Inda Florenn marandache Deparmen of Mahemacs, Unvers of New Meco, Gallup, UA
More informationForecast of Stock Index Volatility Using Grey GARCH-Type Models
Send Orders for Rerns o rerns@benhamscence.ae he Oen Cybernecs & Sysemcs Journal, 015, 9, 93-98 93 Oen Access Forecas of Sock Index Volaly Usng Grey GARCH-ye Models L-Yan Geng 1, and Zhan-Fu Zhang 1 School
More informationA New Method for Computing EM Algorithm Parameters in Speaker Identification Using Gaussian Mixture Models
0 IACSI Hong Kong Conferences IPCSI vol. 9 (0) (0) IACSI Press, Sngaore A New ehod for Comung E Algorhm Parameers n Seaker Idenfcaon Usng Gaussan xure odels ohsen Bazyar +, Ahmad Keshavarz, and Khaoon
More informationNew M-Estimator Objective Function. in Simultaneous Equations Model. (A Comparative Study)
Inernaonal Mahemacal Forum, Vol. 8, 3, no., 7 - HIKARI Ld, www.m-hkar.com hp://dx.do.org/.988/mf.3.3488 New M-Esmaor Objecve Funcon n Smulaneous Equaons Model (A Comparave Sudy) Ahmed H. Youssef Professor
More information- The whole joint distribution is independent of the date at which it is measured and depends only on the lag.
Saionary Processes Sricly saionary - The whole join disribuion is indeenden of he dae a which i is measured and deends only on he lag. - E y ) is a finie consan. ( - V y ) is a finie consan. ( ( y, y s
More informationIn the complete model, these slopes are ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL. (! i+1 -! i ) + [(!") i+1,q - [(!
ANALYSIS OF VARIANCE FOR THE COMPLETE TWO-WAY MODEL The frs hng o es n wo-way ANOVA: Is here neracon? "No neracon" means: The man effecs model would f. Ths n urn means: In he neracon plo (wh A on he horzonal
More informationWe are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors
We are InechOen, he world s leadng ublsher of Oen Access books Bul by scenss, for scenss 4,100 116,000 120M Oen access books avalable Inernaonal auhors and edors Downloads Our auhors are among he 154 Counres
More informationRobustness of DEWMA versus EWMA Control Charts to Non-Normal Processes
Journal of Modern Appled Sascal Mehods Volume Issue Arcle 8 5--3 Robusness of D versus Conrol Chars o Non- Processes Saad Saeed Alkahan Performance Measuremen Cener of Governmen Agences, Insue of Publc
More informationJohn Geweke a and Gianni Amisano b a Departments of Economics and Statistics, University of Iowa, USA b European Central Bank, Frankfurt, Germany
Herarchcal Markov Normal Mxure models wh Applcaons o Fnancal Asse Reurns Appendx: Proofs of Theorems and Condonal Poseror Dsrbuons John Geweke a and Gann Amsano b a Deparmens of Economcs and Sascs, Unversy
More informationModelling Abrupt Shift in Time Series Using Indicator Variable: Evidence of Nigerian Insurance Stock
Inernaonal Journal of Fnance and Accounng 5, 4(): 9-3 DOI:.593/.fa.54. Modellng Abru Shf n Tme Seres Usng Indcaor Varable: Evdence of Ngeran Insurance Sock H. G. Dkko, O. E. Asrbo, A. Samson,* Dearmen
More informationImprovement in Estimating Population Mean using Two Auxiliary Variables in Two-Phase Sampling
Improvemen n Esmang Populaon Mean usng Two Auxlar Varables n Two-Phase amplng Rajesh ngh Deparmen of ascs, Banaras Hndu Unvers(U.P.), Inda (rsnghsa@ahoo.com) Pankaj Chauhan and Nrmala awan chool of ascs,
More information2. SPATIALLY LAGGED DEPENDENT VARIABLES
2. SPATIALLY LAGGED DEPENDENT VARIABLES In hs chaper, we descrbe a sascal model ha ncorporaes spaal dependence explcly by addng a spaally lagged dependen varable y on he rgh-hand sde of he regresson equaon.
More informationCHAPTER FOUR REPEATED MEASURES IN TOXICITY TESTING
CHAPTER FOUR REPEATED MEASURES IN TOXICITY TESTING 4. Inroducon The repeaed measures sudy s a very commonly used expermenal desgn n oxcy esng because no only allows one o nvesgae he effecs of he oxcans,
More informationNotes on the stability of dynamic systems and the use of Eigen Values.
Noes on he sabl of dnamc ssems and he use of Egen Values. Source: Macro II course noes, Dr. Davd Bessler s Tme Seres course noes, zarads (999) Ineremporal Macroeconomcs chaper 4 & Techncal ppend, and Hamlon
More informationChapter Lagrangian Interpolation
Chaper 5.4 agrangan Inerpolaon Afer readng hs chaper you should be able o:. dere agrangan mehod of nerpolaon. sole problems usng agrangan mehod of nerpolaon and. use agrangan nerpolans o fnd deraes and
More informationLet s treat the problem of the response of a system to an applied external force. Again,
Page 33 QUANTUM LNEAR RESPONSE FUNCTON Le s rea he problem of he response of a sysem o an appled exernal force. Agan, H() H f () A H + V () Exernal agen acng on nernal varable Hamlonan for equlbrum sysem
More informationCubic Bezier Homotopy Function for Solving Exponential Equations
Penerb Journal of Advanced Research n Compung and Applcaons ISSN (onlne: 46-97 Vol. 4, No.. Pages -8, 6 omoopy Funcon for Solvng Eponenal Equaons S. S. Raml *,,. Mohamad Nor,a, N. S. Saharzan,b and M.
More informationKeywords: Hedonic regressions; hedonic indexes; consumer price indexes; superlative indexes.
Hedonc Imuaon versus Tme Dummy Hedonc Indexes Erwn Dewer, Saeed Herav and Mck Slver December 5, 27 (wh a commenary by Jan de Haan) Dscusson Paer 7-7, Dearmen of Economcs, Unversy of Brsh Columba, 997-873
More informationKayode Ayinde Department of Pure and Applied Mathematics, Ladoke Akintola University of Technology P. M. B. 4000, Ogbomoso, Oyo State, Nigeria
Journal of Mahemacs and Sascs 3 (4): 96-, 7 ISSN 549-3644 7 Scence Publcaons A Comparave Sudy of he Performances of he OLS and some GLS Esmaors when Sochasc egressors are boh Collnear and Correlaed wh
More informationDynamic Team Decision Theory. EECS 558 Project Shrutivandana Sharma and David Shuman December 10, 2005
Dynamc Team Decson Theory EECS 558 Proec Shruvandana Sharma and Davd Shuman December 0, 005 Oulne Inroducon o Team Decson Theory Decomposon of he Dynamc Team Decson Problem Equvalence of Sac and Dynamc
More informationTime Scale Evaluation of Economic Forecasts
CENTRAL BANK OF CYPRUS EUROSYSTEM WORKING PAPER SERIES Tme Scale Evaluaon of Economc Forecass Anons Mchs February 2014 Worng Paper 2014-01 Cenral Ban of Cyprus Worng Papers presen wor n progress by cenral
More informationLecture 6: Learning for Control (Generalised Linear Regression)
Lecure 6: Learnng for Conrol (Generalsed Lnear Regresson) Conens: Lnear Mehods for Regresson Leas Squares, Gauss Markov heorem Recursve Leas Squares Lecure 6: RLSC - Prof. Sehu Vjayakumar Lnear Regresson
More informationPHYS 705: Classical Mechanics. Canonical Transformation
PHYS 705: Classcal Mechancs Canoncal Transformaon Canoncal Varables and Hamlonan Formalsm As we have seen, n he Hamlonan Formulaon of Mechancs,, are ndeenden varables n hase sace on eual foong The Hamlon
More informationOrdinary Differential Equations in Neuroscience with Matlab examples. Aim 1- Gain understanding of how to set up and solve ODE s
Ordnary Dfferenal Equaons n Neuroscence wh Malab eamples. Am - Gan undersandng of how o se up and solve ODE s Am Undersand how o se up an solve a smple eample of he Hebb rule n D Our goal a end of class
More information5th International Conference on Advanced Design and Manufacturing Engineering (ICADME 2015)
5h Inernaonal onference on Advanced Desgn and Manufacurng Engneerng (IADME 5 The Falure Rae Expermenal Sudy of Specal N Machne Tool hunshan He, a, *, La Pan,b and Bng Hu 3,c,,3 ollege of Mechancal and
More informationFORECASTS GENERATING FOR ARCH-GARCH PROCESSES USING THE MATLAB PROCEDURES
FORECASS GENERAING FOR ARCH-GARCH PROCESSES USING HE MALAB PROCEDURES Dušan Marček, Insiue of Comuer Science, Faculy of Philosohy and Science, he Silesian Universiy Oava he Faculy of Managemen Science
More informationInverse Joint Moments of Multivariate. Random Variables
In J Conem Mah Scences Vol 7 0 no 46 45-5 Inverse Jon Momens of Mulvarae Rom Varables M A Hussan Dearmen of Mahemacal Sascs Insue of Sascal Sudes Research ISSR Caro Unversy Egy Curren address: Kng Saud
More informationComb Filters. Comb Filters
The smple flers dscussed so far are characered eher by a sngle passband and/or a sngle sopband There are applcaons where flers wh mulple passbands and sopbands are requred Thecomb fler s an example of
More informationHEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD
Journal of Appled Mahemacs and Compuaonal Mechancs 3, (), 45-5 HEAT CONDUCTION PROBLEM IN A TWO-LAYERED HOLLOW CYLINDER BY USING THE GREEN S FUNCTION METHOD Sansław Kukla, Urszula Sedlecka Insue of Mahemacs,
More informationIs it necessary to seasonally adjust business and consumer surveys. Emmanuelle Guidetti
Is necessar o seasonall adjs bsness and consmer srves Emmanelle Gde Olne 1 BTS feares 2 Smlaon eercse 3 Seasonal ARIMA modellng 4 Conclsons Jan-85 Jan-87 Jan-89 Jan-91 Jan-93 Jan-95 Jan-97 Jan-99 Jan-01
More informationMethods for the estimation of missing values in time series
Edh Cowan Unversy Research Onlne Theses: Docoraes and Masers Theses 6 Mehods for he esmaon of mssng values n me seres Davd S. Fung Edh Cowan Unversy Recommended Caon Fung, D. S. (6). Mehods for he esmaon
More informationThe Development of a Three-Dimensional Material Point Method Computer Simulation Algorithm for Bullet Impact Studies
MERIN JOURN OF UNDERGRDUTE RESERH VO. 9, NOS. & 3 ( The Develomen of a Three-Dmensonal Maeral on Mehod omuer Smulaon lgorhm for ulle Imac Sudes M.J. onnolly, E. Maldonado and M.W. Roh Dearmen of hyscs
More informationTime-interval analysis of β decay. V. Horvat and J. C. Hardy
Tme-nerval analyss of β decay V. Horva and J. C. Hardy Work on he even analyss of β decay [1] connued and resuled n he developmen of a novel mehod of bea-decay me-nerval analyss ha produces hghly accurae
More informationCS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 4
CS434a/54a: Paern Recognon Prof. Olga Veksler Lecure 4 Oulne Normal Random Varable Properes Dscrmnan funcons Why Normal Random Varables? Analycally racable Works well when observaon comes form a corruped
More informationBOX-JENKINS MODEL NOTATION. The Box-Jenkins ARMA(p,q) model is denoted by the equation. pwhile the moving average (MA) part of the model is θ1at
BOX-JENKINS MODEL NOAION he Box-Jenkins ARMA(,q) model is denoed b he equaion + + L+ + a θ a L θ a 0 q q. () he auoregressive (AR) ar of he model is + L+ while he moving average (MA) ar of he model is
More informationCS286.2 Lecture 14: Quantum de Finetti Theorems II
CS286.2 Lecure 14: Quanum de Fne Theorems II Scrbe: Mara Okounkova 1 Saemen of he heorem Recall he las saemen of he quanum de Fne heorem from he prevous lecure. Theorem 1 Quanum de Fne). Le ρ Dens C 2
More informationFall 2009 Social Sciences 7418 University of Wisconsin-Madison. Problem Set 2 Answers (4) (6) di = D (10)
Publc Affars 974 Menze D. Chnn Fall 2009 Socal Scences 7418 Unversy of Wsconsn-Madson Problem Se 2 Answers Due n lecure on Thursday, November 12. " Box n" your answers o he algebrac quesons. 1. Consder
More informationChapter 3: Signed-rank charts
Chaer : gned-ran chars.. The hewhar-ye conrol char... Inroducon As menoned n Chaer, samles of fxed sze are aen a regular nervals and he long sasc s hen loed. The queson s: Whch qualy arameer should be
More informationPanel Data Regression Models
Panel Daa Regresson Models Wha s Panel Daa? () Mulple dmensoned Dmensons, e.g., cross-secon and me node-o-node (c) Pongsa Pornchawseskul, Faculy of Economcs, Chulalongkorn Unversy (c) Pongsa Pornchawseskul,
More informationAppendix H: Rarefaction and extrapolation of Hill numbers for incidence data
Anne Chao Ncholas J Goell C seh lzabeh L ander K Ma Rober K Colwell and Aaron M llson 03 Rarefacon and erapolaon wh ll numbers: a framewor for samplng and esmaon n speces dversy sudes cology Monographs
More informationTHE PREDICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS
THE PREICTION OF COMPETITIVE ENVIRONMENT IN BUSINESS INTROUCTION The wo dmensonal paral dfferenal equaons of second order can be used for he smulaon of compeve envronmen n busness The arcle presens he
More informationLecture VI Regression
Lecure VI Regresson (Lnear Mehods for Regresson) Conens: Lnear Mehods for Regresson Leas Squares, Gauss Markov heorem Recursve Leas Squares Lecure VI: MLSC - Dr. Sehu Vjayakumar Lnear Regresson Model M
More informationMachine Learning Linear Regression
Machne Learnng Lnear Regresson Lesson 3 Lnear Regresson Bascs of Regresson Leas Squares esmaon Polynomal Regresson Bass funcons Regresson model Regularzed Regresson Sascal Regresson Mamum Lkelhood (ML)
More informationMALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES. Institute for Mathematical Research, Universiti Putra Malaysia, UPM Serdang, Selangor, Malaysia
Malaysan Journal of Mahemacal Scences 9(2): 277-300 (2015) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Journal homeage: h://ensemumedumy/journal A Mehod for Deermnng -Adc Orders of Facorals 1* Rafka Zulkal,
More informationForecasting Post-War Tourist Arrivals to Sri Lanka Using Dynamic Transfer Function Modeling Method
IJMS 0 vol. 3 (): 111-1 Inernaonal Journal of Muldscplnary Sudes (IJMS) Volume 3, Issue, 0 Forecasng Pos-War Tours Arrvals o Sr Lanka Usng Dynamc Transfer Funcon Modelng Mehod Gnanapragasam SR 1 and Cooray
More informationArticle A Hybrid Method for Short-Term Wind Speed Forecasting
Arcle A Hybrd Mehod for Shor-Term Wnd Speed Forecasng Jnlang Zhang 1,2, *, YMng We 2,3, Zhong-fu Tan 1, Wang Ke 2,3 and We Tan 4 1 School of Economcs and Managemen, Norh Chna Elecrc Power Unversy, Bejng
More informationBernoulli process with 282 ky periodicity is detected in the R-N reversals of the earth s magnetic field
Submed o: Suden Essay Awards n Magnecs Bernoull process wh 8 ky perodcy s deeced n he R-N reversals of he earh s magnec feld Jozsef Gara Deparmen of Earh Scences Florda Inernaonal Unversy Unversy Park,
More information2/20/2013. EE 101 Midterm 2 Review
//3 EE Mderm eew //3 Volage-mplfer Model The npu ressance s he equalen ressance see when lookng no he npu ermnals of he amplfer. o s he oupu ressance. I causes he oupu olage o decrease as he load ressance
More informationSurvival Analysis and Reliability. A Note on the Mean Residual Life Function of a Parallel System
Communcaons n Sascs Theory and Mehods, 34: 475 484, 2005 Copyrgh Taylor & Francs, Inc. ISSN: 0361-0926 prn/1532-415x onlne DOI: 10.1081/STA-200047430 Survval Analyss and Relably A Noe on he Mean Resdual
More informationChapter 6: AC Circuits
Chaper 6: AC Crcus Chaper 6: Oulne Phasors and he AC Seady Sae AC Crcus A sable, lnear crcu operang n he seady sae wh snusodal excaon (.e., snusodal seady sae. Complee response forced response naural response.
More informationJ i-1 i. J i i+1. Numerical integration of the diffusion equation (I) Finite difference method. Spatial Discretization. Internal nodes.
umercal negraon of he dffuson equaon (I) Fne dfference mehod. Spaal screaon. Inernal nodes. R L V For hermal conducon le s dscree he spaal doman no small fne spans, =,,: Balance of parcles for an nernal
More informationグラフィカルモデルによる推論 確率伝搬法 (2) Kenji Fukumizu The Institute of Statistical Mathematics 計算推論科学概論 II (2010 年度, 後期 )
グラフィカルモデルによる推論 確率伝搬法 Kenj Fukuzu he Insue of Sascal Maheacs 計算推論科学概論 II 年度 後期 Inference on Hdden Markov Model Inference on Hdden Markov Model Revew: HMM odel : hdden sae fne Inference Coue... for any Naïve
More informationVECTOR AUTOREGRESSIVE MOVING AVERAGE MODEL
VECOR AUOREGRESSIVE MOVING AVERAGE MODEL Ran Kumar Paul Indan Agrculural Sac Reearch Inue Lbrar Avenue, New Delh rana@gmal.com Inroducon In man area of Agrculure relaed feld forecang roducon ver much eenal
More informationGENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS. Youngwoo Ahn and Kitae Kim
Korean J. Mah. 19 (2011), No. 3, pp. 263 272 GENERATING CERTAIN QUINTIC IRREDUCIBLE POLYNOMIALS OVER FINITE FIELDS Youngwoo Ahn and Kae Km Absrac. In he paper [1], an explc correspondence beween ceran
More informationfor regression Y ˆ ˆ corr Z, X 0 and
Lecure Noes 9 Econ B2000, MA Economercs Kevn R Foser, CCNY Fall 202 Insrumenal Varables Regresson vald nsrumen, some Z for regresson Y ˆ ˆ 0 X u corr Z, X 0 and o relevance: o exogeney: corr Z, u 0 o nsrumen
More informationLinear Response Theory: The connection between QFT and experiments
Phys540.nb 39 3 Lnear Response Theory: The connecon beween QFT and expermens 3.1. Basc conceps and deas Q: ow do we measure he conducvy of a meal? A: we frs nroduce a weak elecrc feld E, and hen measure
More information10. A.C CIRCUITS. Theoretically current grows to maximum value after infinite time. But practically it grows to maximum after 5τ. Decay of current :
. A. IUITS Synopss : GOWTH OF UNT IN IUIT : d. When swch S s closed a =; = d. A me, curren = e 3. The consan / has dmensons of me and s called he nducve me consan ( τ ) of he crcu. 4. = τ; =.63, n one
More informationA NEW TECHNIQUE FOR SOLVING THE 1-D BURGERS EQUATION
S19 A NEW TECHNIQUE FOR SOLVING THE 1-D BURGERS EQUATION by Xaojun YANG a,b, Yugu YANG a*, Carlo CATTANI c, and Mngzheng ZHU b a Sae Key Laboraory for Geomechancs and Deep Underground Engneerng, Chna Unversy
More informationPerformance Analysis for a Network having Standby Redundant Unit with Waiting in Repair
TECHNI Inernaonal Journal of Compung Scence Communcaon Technologes VOL.5 NO. July 22 (ISSN 974-3375 erformance nalyss for a Nework havng Sby edundan Un wh ang n epar Jendra Sngh 2 abns orwal 2 Deparmen
More informationMODELING TIME-VARYING TRADING-DAY EFFECTS IN MONTHLY TIME SERIES
MODELING TIME-VARYING TRADING-DAY EFFECTS IN MONTHLY TIME SERIES Wllam R. Bell, Census Bureau and Donald E. K. Marn, Howard Unversy and Census Bureau Donald E. K. Marn, Howard Unversy, Washngon DC 0059
More informationChapter 8 Dynamic Models
Chaper 8 Dnamc odels 8. Inroducon 8. Seral correlaon models 8.3 Cross-seconal correlaons and me-seres crosssecon models 8.4 me-varng coeffcens 8.5 Kalman fler approach 8. Inroducon When s mporan o consder
More informationChapter 6 DETECTION AND ESTIMATION: Model of digital communication system. Fundamental issues in digital communications are
Chaper 6 DCIO AD IMAIO: Fndaenal sses n dgal concaons are. Deecon and. saon Deecon heory: I deals wh he desgn and evalaon of decson ang processor ha observes he receved sgnal and gesses whch parclar sybol
More informationChickens vs. Eggs: Replicating Thurman and Fisher (1988) by Arianto A. Patunru Department of Economics, University of Indonesia 2004
Chicens vs. Eggs: Relicaing Thurman and Fisher (988) by Ariano A. Paunru Dearmen of Economics, Universiy of Indonesia 2004. Inroducion This exercise lays ou he rocedure for esing Granger Causaliy as discussed
More informationApplication of Vector Error Correction Model (VECM) and Impulse Response Function for Analysis Data Index of Farmers Terms of Trade
Indan Journal of Scence and echnology, Vol 0(9), DOI: 0.7485/js/07/v09/58, May 07 ISSN (Prn) : 0974-6846 ISSN (Onlne) : 0974-5645 Applcaon of Vecor Error Correcon Model (VECM) and Impulse Response Funcon
More informationDepartment of Economics University of Warsaw Warsaw, Poland Długa Str. 44/50.
MIGRATIOS OF HETEROGEEOUS POPULATIO OF DRIVERS ACROSS CLASSES OF A BOUS-MALUS SYSTEM BY WOJCIECH OTTO Dearmen of Economcs Unversy of Warsaw 00-24 Warsaw Poland Długa Sr. 44/50 woo@wne.uw.edu.l . ITRODUCTIO
More informationExistence and Uniqueness Results for Random Impulsive Integro-Differential Equation
Global Journal of Pure and Appled Mahemacs. ISSN 973-768 Volume 4, Number 6 (8), pp. 89-87 Research Inda Publcaons hp://www.rpublcaon.com Exsence and Unqueness Resuls for Random Impulsve Inegro-Dfferenal
More informationA New Generalisation of Sam-Solai s Multivariate symmetric Arcsine Distribution of Kind-1*
IOSR Journal o Mahemacs IOSRJM ISSN: 78-578 Volume, Issue May-June 0, PP 4-48 www.osrournals.org A New Generalsaon o Sam-Sola s Mulvarae symmerc Arcsne Dsrbuon o Knd-* Dr. G.S. Davd Sam Jayaumar. Dr.A.Solarau.
More informationAdvanced Machine Learning & Perception
Advanced Machne Learnng & Percepon Insrucor: Tony Jebara SVM Feaure & Kernel Selecon SVM Eensons Feaure Selecon (Flerng and Wrappng) SVM Feaure Selecon SVM Kernel Selecon SVM Eensons Classfcaon Feaure/Kernel
More information3. OVERVIEW OF NUMERICAL METHODS
3 OVERVIEW OF NUMERICAL METHODS 3 Inroducory remarks Ths chaper summarzes hose numercal echnques whose knowledge s ndspensable for he undersandng of he dfferen dscree elemen mehods: he Newon-Raphson-mehod,
More informationVolatility Interpolation
Volaly Inerpolaon Prelmnary Verson March 00 Jesper Andreasen and Bran Huge Danse Mares, Copenhagen wan.daddy@danseban.com brno@danseban.com Elecronc copy avalable a: hp://ssrn.com/absrac=69497 Inro Local
More information[Link to MIT-Lab 6P.1 goes here.] After completing the lab, fill in the following blanks: Numerical. Simulation s Calculations
Chaper 6: Ordnary Leas Squares Esmaon Procedure he Properes Chaper 6 Oulne Cln s Assgnmen: Assess he Effec of Sudyng on Quz Scores Revew o Regresson Model o Ordnary Leas Squares () Esmaon Procedure o he
More informationLaser Interferometer Space Antenna (LISA)
aser nerferomeer Sace Anenna SA Tme-elay nerferomery wh Movng Sacecraf Arrays Massmo Tno Je Proulson aboraory, Calforna nsue of Technology GSFC JP 8 h GWAW, ec 7-0, 00, Mlwaukee, Wsconsn WM Folkner e al,
More information