Air Traffic Forecast Empirical Research Based on the MCMC Method
|
|
- Irene Kelly
- 6 years ago
- Views:
Transcription
1 Compuer and Informaion Science; Vol. 5, No. 5; 0 ISSN E-ISSN Published by Canadian Cener of Science and Educaion Air Traffic Forecas Empirical Research Based on he MCMC Mehod Jian-bo Wang, Chong-jun Fan & Lei Bai Business School, Universiy of Shanghai for Science and Technology, Shanghai 00093, China Correspondence: Chong-jun Fan, Business School, Universiy of Shanghai for Science and Technology, Shanghai 00093, China. Tel: wangjianbai@homail.com Received: June, 0 Acceped: July, 0 Online Published: July 7, 0 doi:0.5539/cis.v5n5p50 URL: hp://dx.doi.org/0.5539/cis.v5n5p50 Absrac Airpor air raffic is one of he mos imporan and hardes one among all he airpor daa forecass. In his paper, he Markov Chain Mone Carlo (MCMC) mehod of applied heory of saisics has been inroduced ino he aviaion secor, and he discussion on airpor air raffic forecas has been conduced aking Shanghai airpor as he applicaion background. MCMC mehod is designed o solve he problem parameers expecaions where some baysian inference can no be direcly calculaed. This paper conducs he ieraive compuaion using WinBUGS, EViews is used o esimae he iniial ieraion parameer of he regression model, and i is on his basis ha he regression equaion under he MCMC mehod is obained. Keywords: MCMC, airpor, raffic daa, muliple linear regression. Inroducion In Bayesian saisical learning, he esimaion of overall parameer ofen needs high-dimensional probabiliy disribuion funcion o do he inegral. Under normal circumsances, here is no clear expression of his inegral, so i is very difficul o do i. In his case, he Markov Chain Mone Carlo (MCMC) mehod is a simple and effecive calculaion mehod. The basic idea of MCMC mehod is o consruc a Markov Chain, making is saionary disribuion he poserior disribuion of he parameers o be esimaed, hus generaing he samples of poserior disribuion hrough his Markov Chain, and doing he Mone Carlo inegraion of he samples of Markov Chain when hey reach saionary disribuion. For he reason ha he join disribuion of he parameers and variables in many ime series models are complex and high dimensional, i is impossible o exrac samples direcly he disribuion. MCMC can generae samples he join densiy of parameers, wihou having o know he specific forms of expression of densiy. The mos commonly used MCMC algorihm is Gibbs sampling, which is more accessible o gain he condiional disribuion densiy hrough simulaion, and i adops ieraive procedure o consruc he Markov Chain. Under he simulaion of large sample, i is expeced o converge o he arge disribuion (Lin & Han, 007; Sephen, 998; Lin & Han, 005). WinBUGS (Lunn, Spiegelhaler, Thomas, & Bes, 009) is a kind of specialized saisical sofware based on MCMC mehod, in which all unknown parameers are regarded as random variables in order o find ou he soluion of his kind of probabiliy model hen. In his paper, basing on he regression model, prereamen is done o he iniial value of he MCMC model, and empirical sudy is done regarding he relaionship beween air raffic, GDP growh rae and populaion growh (David, 00). The GDP means he money value of all he final goods and services produced in he domesic erriory of a counry in a year s ime. In his paper, we conduc hree pars, including MCMC mehod descripion, Empirical research on air raffic forecas and conclusions.. MCMC Mehod and Gibbs Sampling The MCMC mehod is a general simulaion mehod for sampling poserior disribuions and compuing poserior quaniies of ineres. MCMC mehods sample succeed a arge disribuion, and each sample depends on he previous one, hence he noion of he Markov Chain. A Markov Chain is a sequence of random variables,,,, in which he random variable depends on all s previous only hrough is immediae predecessor. The Markov Chain applied o sampling as a mechanism ha raverses randomly hrough a arge disribuion wihou having any memory of where i has been. Where i moves nex is enirely dependen on where i is now. 50
2 Compuer and Informaion Science Vol. 5, No. 5; 0 Mone Carlo, as in Mone Carlo inegraion, is mainly used o approximae an expecaion by using he Markov Chain samples. In he simples version n g( ) p( ) d g( ) () s n where g (.) is a funcion of ineres and are samples p( ) on is suppor S. This approximaes he expeced value of g( ). The Gibbs sampler is a special case of he Meropolis sampler in which he proposal disribuions exacly mach he poserior condiional disribuions and proposals are acceped 00% of he ime. The main feaures of he Gibbs sampling, a Markov Chain sample, is exraced a series of condiional disribuions (George & Edward, 99; Gelfand & Smih, 990). Suppose ( )' is he parameer vecor, py ( ) is he likelihood, and ( ) is he prior disribuion. The full poserior condiional disribuion of ( i j, i j is proporional o he join poserior densiy; ha is, ( i j, i j p( y ) ( ) () For insance, he one-dimensional condiional disribuion of given * j, j j k is compued as he following:, j k p( y (,,, ) (,,, )') (3) * * * * * j j k k The Gibbs sampler works as follows: (0) 0 (0) Se 0 and choose an arbirary iniial value of { }. Generae each componen of as follows: ( ) draw () () ( ) draw ( ) ( ) ( ), 3... ( ) draw ( ) ( ) k k,, k Se. If T he number of desired samples, reurn o sep. When Gibbs sampling ieraes enough number of seps, we can believe he sample is subjec o he arge disribuion, or converge o he arge disribuion. The choice of he iniial value will affec he convergence rae of MCMC. 3. Empirical Research on Air Traffic Forecas I is raher imporan o choose proper facors ahead of modeling and esimaion. The following pars include facors analysis and parameers esimaion. 3. Influencing Facors Analysis Air raffic forecas has always been a difficul and key poin of he airpor aviaion forecass. There are many aspecs of facors ha influence air raffic, in his paper, macro-economic facor and demographic impac facor are used o analyze he air raffic. For economic facor, GDP is an imporan sandard o measure he comprehensive level of he economic developmen; populaion growh rae index is used as demographic facor. A he same ime, because he daa of GDP and populaion growh rae are raher consecuive, he reference value of he daa is quie high. We conduc a linear regression of air raffic growh based on GDP growh rae and populaion growh rae. There a binary regression model is buil as follows: Y 0 X X (4) Where 0,, are parameers o be esimaed, X and X are know vecors sand for GDP growh and populaion growh. Y sands for Airpor raffic volume growh. This aricle uses Eviews simulaion, according o daa released by Naional Bureau of Saisic 990 o 009, and hen Equaion (4) would have he form: Y X 3.77X (5) 5
3 Compuer and Informaion Science Vol. 5, No. 5; 0 3. Parameers Esimaion Based on he MCMC Mehod WinBUGS is a special package based on MCMC mehod, which can easily conduc Gibbs sampling on many models and disribuions. Before implemening Bayesian analysis, programmers need o se prior disribuion and iniial values of variables. The analysis processes include modeling, daa loading, iniial values and ieraive analysis. In his case, he priors are defined as Regression equaion: muy[] i X [] i X [] i (6) Prior disribuion of he parameers: 0 0 ~ N (0,.0 E 6) ; ~ N (0,.0 E 6) ; ~ N (0,.0 E 6) ; ~ Gamma(0.00,0.00) ; Yi ~ N( muy[ i], ). The iniial parameer value: ; ; 3.77;. In order o ensure he convergence of he parameers, we firs conduc 5000 Gibbs pre-ieraions, and hen discard he iniial ieraion. We consruc he model using a graphical inerface Doodle (Figure ). Meanwhile, we se up hree Markov Chains in WinBUGS, and a he same ime ener he iniial values. This is because, for he convergence parameers, (and) hree ieraions resuls of he Markov Chain end o be overlapping. I is conducive o he observaion of he resuls. Figure. Bayesian modeling in WinBUGS Ieraions sar 500 o The saisical resuls show as in Table, which produces summary saisics for he variable, pooling over he chains seleced. In his case, we choose.5% and 97.5% as he confidence / inerval. The quaniy repored in he MC error column gives an esimae of / N, he Mone Carlo sandard error of he mean. During he 95% Confidence inervals, he average value of α0, α, α and τ are he inerval (-4.3, 0.4), (-0.383,.4), (-7.6, 4.69), (0.0044, ) respecively. Table. Poserior esimae saisics for moniored parameers wih ieraions node mean sd MC error.50% median 97.50% sar sample α α α τ E Before drawing our conclusion, we need o have a judgmen on densiy esimaion and convergence saisics of he parameers. The following plos a smoohed Poserior Kernel Densiy Esimae for he variable if i is coninuous or a hisogram if i is discree (Figure ). 5
4 Compuer and Informaion Science Vol. 5, No. 5; 0 Figure. Poserior Kernel Densiy Esimaion for moniored parameers wih ieraions In his case, we are running hree chains simulaneously. The race and hisory plos show each chain in a differen color. As he plos show (Figure 3), all he chains appear o be overlapping one anoher, herefor, we can be reasonably confiden ha convergence has been achieved. Figure 3. Trace plos for moniored parameers hree Chains ieraions 53
5 Compuer and Informaion Science Vol. 5, No. 5; 0 The following plos (Figure 4) show he Gelman-Rubin convergence saisics. The widh of he cenral 80% inerval of he pooled runs is green, he average widh of he 80% inervals wihin he individual runs is blue, and heir raio R (= pooled / wihin) is red. R would generally be expeced o be greaer han if he saring values are suiably over-dispersed. The proper one should be concerned boh wih convergence of R o, and wih convergence of boh he pooled and wihin inerval widhs o sabiliy. Figure 4. Parameers convergence diagnosics According o he daa of Table, as well as Poserior Kernel Densiy Esimaion, Trace Plos, Parameers Convergence Diagnosics, we believe ha parameers esimaion based on MCMC mehod is valid. The linear regression model equaion is shown as in Equaion (7). Y X 3.756X (7) 4. Conclusions Air raffic forecas has always been a problem roubling airpor senior decision-makers o make scienific decisions, bu he use of convenional forecasing mehods ofen canno ge saisfacory resuls. In his paper, MCMC simulaion mehods and WinBUGS sofware applicaions based on Gibbs sampling are used o conduc he volailiy analysis of he airpor air raffic daa. The analyical mehod proposed in his paper can provide new ideas and mehods for he scienific decision-making of he airpor senior decision-makers. We hope ha, based on he fundamenal predicion research of he MCMC mehod, airpors can have a deeper level of sudy. References David, P. M. S. (00). Acuarial modeling wih MCMC and BUGS. Norh American Acuarial Journal, 5(), Gelfand, A. E., & Smih, A. F. M. (990). Sampling-Based Approaches o Calculaing Marginal Densiies. J. American Saisical Associaion, 85, hp://dx.doi.org/0.080/ George, C., & Edward, I. G. (99). Explaining he Gibbs sampler. The American Saisician, 46, Lin, J., & Han, Y. (005). Exponenial Regression Model Based on MCMC Mehod and Is Applicaion. Operaions Research and Managemen Science, 4(4), Lin, J., & Han, Y. (007). Gamma Process Priors Model Based on MCMC Simulaion and Is Applicaion in Reliabiliy. Journal of Sysem Simulaion, 9(), Lunn, D., Spiegelhaler, D., Thomas, A., & Bes, N. (009). The BUGS projec: Evoluion, criique and fuure direcions. Saisics in Medicine, 8, hp://dx.doi.org/0.00/sim.3680 Sephen, P. B. (998). Series D (The Saisician): Chain Mone Carlo Mehod and Is Applicaion. Journal of he Royal Saisical Sociey, 47(),
Vehicle Arrival Models : Headway
Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where
More informationReliability of Technical Systems
eliabiliy of Technical Sysems Main Topics Inroducion, Key erms, framing he problem eliabiliy parameers: Failure ae, Failure Probabiliy, Availabiliy, ec. Some imporan reliabiliy disribuions Componen reliabiliy
More informationSimulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010
Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid
More informationGeorey E. Hinton. University oftoronto. Technical Report CRG-TR February 22, Abstract
Parameer Esimaion for Linear Dynamical Sysems Zoubin Ghahramani Georey E. Hinon Deparmen of Compuer Science Universiy oftorono 6 King's College Road Torono, Canada M5S A4 Email: zoubin@cs.orono.edu Technical
More informationBias in Conditional and Unconditional Fixed Effects Logit Estimation: a Correction * Tom Coupé
Bias in Condiional and Uncondiional Fixed Effecs Logi Esimaion: a Correcion * Tom Coupé Economics Educaion and Research Consorium, Naional Universiy of Kyiv Mohyla Academy Address: Vul Voloska 10, 04070
More informationTime series Decomposition method
Time series Decomposiion mehod A ime series is described using a mulifacor model such as = f (rend, cyclical, seasonal, error) = f (T, C, S, e) Long- Iner-mediaed Seasonal Irregular erm erm effec, effec,
More informationStochastic Model for Cancer Cell Growth through Single Forward Mutation
Journal of Modern Applied Saisical Mehods Volume 16 Issue 1 Aricle 31 5-1-2017 Sochasic Model for Cancer Cell Growh hrough Single Forward Muaion Jayabharahiraj Jayabalan Pondicherry Universiy, jayabharahi8@gmail.com
More informationExponential Weighted Moving Average (EWMA) Chart Under The Assumption of Moderateness And Its 3 Control Limits
DOI: 0.545/mjis.07.5009 Exponenial Weighed Moving Average (EWMA) Char Under The Assumpion of Moderaeness And Is 3 Conrol Limis KALPESH S TAILOR Assisan Professor, Deparmen of Saisics, M. K. Bhavnagar Universiy,
More informationDiebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles
Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance
More information) were both constant and we brought them from under the integral.
YIELD-PER-RECRUIT (coninued The yield-per-recrui model applies o a cohor, bu we saw in he Age Disribuions lecure ha he properies of a cohor do no apply in general o a collecion of cohors, which is wha
More informationExponential Smoothing
Exponenial moohing Inroducion A simple mehod for forecasing. Does no require long series. Enables o decompose he series ino a rend and seasonal effecs. Paricularly useful mehod when here is a need o forecas
More informationShiva Akhtarian MSc Student, Department of Computer Engineering and Information Technology, Payame Noor University, Iran
Curren Trends in Technology and Science ISSN : 79-055 8hSASTech 04 Symposium on Advances in Science & Technology-Commission-IV Mashhad, Iran A New for Sofware Reliabiliy Evaluaion Based on NHPP wih Imperfec
More information2017 3rd International Conference on E-commerce and Contemporary Economic Development (ECED 2017) ISBN:
7 3rd Inernaional Conference on E-commerce and Conemporary Economic Developmen (ECED 7) ISBN: 978--6595-446- Fuures Arbirage of Differen Varieies and based on he Coinegraion Which is under he Framework
More informationWhat Ties Return Volatilities to Price Valuations and Fundamentals? On-Line Appendix
Wha Ties Reurn Volailiies o Price Valuaions and Fundamenals? On-Line Appendix Alexander David Haskayne School of Business, Universiy of Calgary Piero Veronesi Universiy of Chicago Booh School of Business,
More information04. Kinetics of a second order reaction
4. Kineics of a second order reacion Imporan conceps Reacion rae, reacion exen, reacion rae equaion, order of a reacion, firs-order reacions, second-order reacions, differenial and inegraed rae laws, Arrhenius
More informationChapter 2. Models, Censoring, and Likelihood for Failure-Time Data
Chaper 2 Models, Censoring, and Likelihood for Failure-Time Daa William Q. Meeker and Luis A. Escobar Iowa Sae Universiy and Louisiana Sae Universiy Copyrigh 1998-2008 W. Q. Meeker and L. A. Escobar. Based
More informationOBJECTIVES OF TIME SERIES ANALYSIS
OBJECTIVES OF TIME SERIES ANALYSIS Undersanding he dynamic or imedependen srucure of he observaions of a single series (univariae analysis) Forecasing of fuure observaions Asceraining he leading, lagging
More informationDistribution of Estimates
Disribuion of Esimaes From Economerics (40) Linear Regression Model Assume (y,x ) is iid and E(x e )0 Esimaion Consisency y α + βx + he esimaes approach he rue values as he sample size increases Esimaion
More informationEvaluation of Mean Time to System Failure of a Repairable 3-out-of-4 System with Online Preventive Maintenance
American Journal of Applied Mahemaics and Saisics, 0, Vol., No., 9- Available online a hp://pubs.sciepub.com/ajams/// Science and Educaion Publishing DOI:0.69/ajams--- Evaluaion of Mean Time o Sysem Failure
More informationModal identification of structures from roving input data by means of maximum likelihood estimation of the state space model
Modal idenificaion of srucures from roving inpu daa by means of maximum likelihood esimaion of he sae space model J. Cara, J. Juan, E. Alarcón Absrac The usual way o perform a forced vibraion es is o fix
More informationNavneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi
Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec
More informationRobust estimation based on the first- and third-moment restrictions of the power transformation model
h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,
More informationNotes on Kalman Filtering
Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren
More informationNon-parametric techniques. Instance Based Learning. NN Decision Boundaries. Nearest Neighbor Algorithm. Distance metric important
on-parameric echniques Insance Based Learning AKA: neares neighbor mehods, non-parameric, lazy, memorybased, or case-based learning Copyrigh 2005 by David Helmbold 1 Do no fi a model (as do LDA, logisic
More informationZápadočeská Univerzita v Plzni, Czech Republic and Groupe ESIEE Paris, France
ADAPTIVE SIGNAL PROCESSING USING MAXIMUM ENTROPY ON THE MEAN METHOD AND MONTE CARLO ANALYSIS Pavla Holejšovsá, Ing. *), Z. Peroua, Ing. **), J.-F. Bercher, Prof. Assis. ***) Západočesá Univerzia v Plzni,
More informationTesting for a Single Factor Model in the Multivariate State Space Framework
esing for a Single Facor Model in he Mulivariae Sae Space Framework Chen C.-Y. M. Chiba and M. Kobayashi Inernaional Graduae School of Social Sciences Yokohama Naional Universiy Japan Faculy of Economics
More information20. Applications of the Genetic-Drift Model
0. Applicaions of he Geneic-Drif Model 1) Deermining he probabiliy of forming any paricular combinaion of genoypes in he nex generaion: Example: If he parenal allele frequencies are p 0 = 0.35 and q 0
More informationRecursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems
8 Froniers in Signal Processing, Vol. 1, No. 1, July 217 hps://dx.doi.org/1.2266/fsp.217.112 Recursive Leas-Squares Fixed-Inerval Smooher Using Covariance Informaion based on Innovaion Approach in Linear
More informationLecture 3: Exponential Smoothing
NATCOR: Forecasing & Predicive Analyics Lecure 3: Exponenial Smoohing John Boylan Lancaser Cenre for Forecasing Deparmen of Managemen Science Mehods and Models Forecasing Mehod A (numerical) procedure
More informationAn introduction to the theory of SDDP algorithm
An inroducion o he heory of SDDP algorihm V. Leclère (ENPC) Augus 1, 2014 V. Leclère Inroducion o SDDP Augus 1, 2014 1 / 21 Inroducion Large scale sochasic problem are hard o solve. Two ways of aacking
More informationIntroduction to Probability and Statistics Slides 4 Chapter 4
Inroducion o Probabiliy and Saisics Slides 4 Chaper 4 Ammar M. Sarhan, asarhan@mahsa.dal.ca Deparmen of Mahemaics and Saisics, Dalhousie Universiy Fall Semeser 8 Dr. Ammar Sarhan Chaper 4 Coninuous Random
More informationA Bayesian Approach to Spectral Analysis
Chirped Signals A Bayesian Approach o Specral Analysis Chirped signals are oscillaing signals wih ime variable frequencies, usually wih a linear variaion of frequency wih ime. E.g. f() = A cos(ω + α 2
More informationCHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK
175 CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK 10.1 INTRODUCTION Amongs he research work performed, he bes resuls of experimenal work are validaed wih Arificial Neural Nework. From he
More informationPhysics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle
Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,
More informationMulti-scale 2D acoustic full waveform inversion with high frequency impulsive source
Muli-scale D acousic full waveform inversion wih high frequency impulsive source Vladimir N Zubov*, Universiy of Calgary, Calgary AB vzubov@ucalgaryca and Michael P Lamoureux, Universiy of Calgary, Calgary
More informationEnsamble methods: Boosting
Lecure 21 Ensamble mehods: Boosing Milos Hauskrech milos@cs.pi.edu 5329 Senno Square Schedule Final exam: April 18: 1:00-2:15pm, in-class Term projecs April 23 & April 25: a 1:00-2:30pm in CS seminar room
More informationEstimation Uncertainty
Esimaion Uncerainy The sample mean is an esimae of β = E(y +h ) The esimaion error is = + = T h y T b ( ) = = + = + = = = T T h T h e T y T y T b β β β Esimaion Variance Under classical condiions, where
More informationMost Probable Phase Portraits of Stochastic Differential Equations and Its Numerical Simulation
Mos Probable Phase Porrais of Sochasic Differenial Equaions and Is Numerical Simulaion Bing Yang, Zhu Zeng and Ling Wang 3 School of Mahemaics and Saisics, Huazhong Universiy of Science and Technology,
More informationInventory Control of Perishable Items in a Two-Echelon Supply Chain
Journal of Indusrial Engineering, Universiy of ehran, Special Issue,, PP. 69-77 69 Invenory Conrol of Perishable Iems in a wo-echelon Supply Chain Fariborz Jolai *, Elmira Gheisariha and Farnaz Nojavan
More informationSTATE-SPACE MODELLING. A mass balance across the tank gives:
B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing
More informationEnsamble methods: Bagging and Boosting
Lecure 21 Ensamble mehods: Bagging and Boosing Milos Hauskrech milos@cs.pi.edu 5329 Senno Square Ensemble mehods Mixure of expers Muliple base models (classifiers, regressors), each covers a differen par
More information1 Differential Equation Investigations using Customizable
Differenial Equaion Invesigaions using Cusomizable Mahles Rober Decker The Universiy of Harford Absrac. The auhor has developed some plaform independen, freely available, ineracive programs (mahles) for
More informationThe Optimal Stopping Time for Selling an Asset When It Is Uncertain Whether the Price Process Is Increasing or Decreasing When the Horizon Is Infinite
American Journal of Operaions Research, 08, 8, 8-9 hp://wwwscirporg/journal/ajor ISSN Online: 60-8849 ISSN Prin: 60-8830 The Opimal Sopping Time for Selling an Asse When I Is Uncerain Wheher he Price Process
More informationIMPROVED HYPERBOLIC FUNCTION METHOD AND EXACT SOLUTIONS FOR VARIABLE COEFFICIENT BENJAMIN-BONA-MAHONY-BURGERS EQUATION
THERMAL SCIENCE, Year 015, Vol. 19, No. 4, pp. 1183-1187 1183 IMPROVED HYPERBOLIC FUNCTION METHOD AND EXACT SOLUTIONS FOR VARIABLE COEFFICIENT BENJAMIN-BONA-MAHONY-BURGERS EQUATION by Hong-Cai MA a,b*,
More informationNon-parametric techniques. Instance Based Learning. NN Decision Boundaries. Nearest Neighbor Algorithm. Distance metric important
on-parameric echniques Insance Based Learning AKA: neares neighbor mehods, non-parameric, lazy, memorybased, or case-based learning Copyrigh 2005 by David Helmbold 1 Do no fi a model (as do LTU, decision
More informationState-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter
Sae-Space Models Iniializaion, Esimaion and Smoohing of he Kalman Filer Iniializaion of he Kalman Filer The Kalman filer shows how o updae pas predicors and he corresponding predicion error variances when
More informationACE 562 Fall Lecture 8: The Simple Linear Regression Model: R 2, Reporting the Results and Prediction. by Professor Scott H.
ACE 56 Fall 5 Lecure 8: The Simple Linear Regression Model: R, Reporing he Resuls and Predicion by Professor Sco H. Irwin Required Readings: Griffihs, Hill and Judge. "Explaining Variaion in he Dependen
More informationHandout: Power Laws and Preferential Attachment
CS224W: Analysis of Neworks Fall 2017 Handou: Power Laws and Preferenial Aachmen 1 Preferenial Aachmen Empirical sudies of real world neworks revealed ha degree disribuion ofen follows a heavyailed disribuion,
More informationRandom Walk with Anti-Correlated Steps
Random Walk wih Ani-Correlaed Seps John Noga Dirk Wagner 2 Absrac We conjecure he expeced value of random walks wih ani-correlaed seps o be exacly. We suppor his conjecure wih 2 plausibiliy argumens and
More informationArticle from. Predictive Analytics and Futurism. July 2016 Issue 13
Aricle from Predicive Analyics and Fuurism July 6 Issue An Inroducion o Incremenal Learning By Qiang Wu and Dave Snell Machine learning provides useful ools for predicive analyics The ypical machine learning
More informationRecent Developments In Evolutionary Data Assimilation And Model Uncertainty Estimation For Hydrologic Forecasting Hamid Moradkhani
Feb 6-8, 208 Recen Developmens In Evoluionary Daa Assimilaion And Model Uncerainy Esimaion For Hydrologic Forecasing Hamid Moradkhani Cener for Complex Hydrosysems Research Deparmen of Civil, Consrucion
More informationT L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB
Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal
More informationA Robust Exponentially Weighted Moving Average Control Chart for the Process Mean
Journal of Modern Applied Saisical Mehods Volume 5 Issue Aricle --005 A Robus Exponenially Weighed Moving Average Conrol Char for he Process Mean Michael B. C. Khoo Universii Sains, Malaysia, mkbc@usm.my
More informationStability. Coefficients may change over time. Evolution of the economy Policy changes
Sabiliy Coefficiens may change over ime Evoluion of he economy Policy changes Time Varying Parameers y = α + x β + Coefficiens depend on he ime period If he coefficiens vary randomly and are unpredicable,
More informationTracking. Announcements
Tracking Tuesday, Nov 24 Krisen Grauman UT Ausin Announcemens Pse 5 ou onigh, due 12/4 Shorer assignmen Auo exension il 12/8 I will no hold office hours omorrow 5 6 pm due o Thanksgiving 1 Las ime: Moion
More informationFigure 1. Jaw RMS target-tracker difference for a.9hz sinusoidal target.
Approaches o Bayesian Smooh Unimodal Regression George Woodworh Dec 3, 999 (Draf - please send commens o george-woodworh@uiowa.edu). Background Speech Ariculaion Daa The daa in Figure were obained by asking
More informationSmoothing. Backward smoother: At any give T, replace the observation yt by a combination of observations at & before T
Smoohing Consan process Separae signal & noise Smooh he daa: Backward smooher: A an give, replace he observaion b a combinaion of observaions a & before Simple smooher : replace he curren observaion wih
More informationLecture 4 Kinetics of a particle Part 3: Impulse and Momentum
MEE Engineering Mechanics II Lecure 4 Lecure 4 Kineics of a paricle Par 3: Impulse and Momenum Linear impulse and momenum Saring from he equaion of moion for a paricle of mass m which is subjeced o an
More informationChapter 15. Time Series: Descriptive Analyses, Models, and Forecasting
Chaper 15 Time Series: Descripive Analyses, Models, and Forecasing Descripive Analysis: Index Numbers Index Number a number ha measures he change in a variable over ime relaive o he value of he variable
More informationA Generalized Poisson-Akash Distribution: Properties and Applications
Inernaional Journal of Saisics and Applicaions 08, 8(5): 49-58 DOI: 059/jsaisics0808050 A Generalized Poisson-Akash Disribuion: Properies and Applicaions Rama Shanker,*, Kamlesh Kumar Shukla, Tekie Asehun
More information3.1 More on model selection
3. More on Model selecion 3. Comparing models AIC, BIC, Adjused R squared. 3. Over Fiing problem. 3.3 Sample spliing. 3. More on model selecion crieria Ofen afer model fiing you are lef wih a handful of
More informationA Dynamic Model of Economic Fluctuations
CHAPTER 15 A Dynamic Model of Economic Flucuaions Modified for ECON 2204 by Bob Murphy 2016 Worh Publishers, all righs reserved IN THIS CHAPTER, OU WILL LEARN: how o incorporae dynamics ino he AD-AS model
More informationClass Meeting # 10: Introduction to the Wave Equation
MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion
More informationUSP. Surplus-Production Models
USP Surplus-Producion Models 2 Overview Purpose of slides: Inroducion o he producion model Overview of differen mehods of fiing Go over some criique of he mehod Source: Haddon 2001, Chaper 10 Hilborn and
More informationA new flexible Weibull distribution
Communicaions for Saisical Applicaions and Mehods 2016, Vol. 23, No. 5, 399 409 hp://dx.doi.org/10.5351/csam.2016.23.5.399 Prin ISSN 2287-7843 / Online ISSN 2383-4757 A new flexible Weibull disribuion
More informationTom Heskes and Onno Zoeter. Presented by Mark Buller
Tom Heskes and Onno Zoeer Presened by Mark Buller Dynamic Bayesian Neworks Direced graphical models of sochasic processes Represen hidden and observed variables wih differen dependencies Generalize Hidden
More informationWednesday, November 7 Handout: Heteroskedasticity
Amhers College Deparmen of Economics Economics 360 Fall 202 Wednesday, November 7 Handou: Heeroskedasiciy Preview Review o Regression Model o Sandard Ordinary Leas Squares (OLS) Premises o Esimaion Procedures
More informationIntroduction D P. r = constant discount rate, g = Gordon Model (1962): constant dividend growth rate.
Inroducion Gordon Model (1962): D P = r g r = consan discoun rae, g = consan dividend growh rae. If raional expecaions of fuure discoun raes and dividend growh vary over ime, so should he D/P raio. Since
More informationAppendix to Creating Work Breaks From Available Idleness
Appendix o Creaing Work Breaks From Available Idleness Xu Sun and Ward Whi Deparmen of Indusrial Engineering and Operaions Research, Columbia Universiy, New York, NY, 127; {xs2235,ww24}@columbia.edu Sepember
More informationA Shooting Method for A Node Generation Algorithm
A Shooing Mehod for A Node Generaion Algorihm Hiroaki Nishikawa W.M.Keck Foundaion Laboraory for Compuaional Fluid Dynamics Deparmen of Aerospace Engineering, Universiy of Michigan, Ann Arbor, Michigan
More informationPENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD
PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD HAN XIAO 1. Penalized Leas Squares Lasso solves he following opimizaion problem, ˆβ lasso = arg max β R p+1 1 N y i β 0 N x ij β j β j (1.1) for some 0.
More informationTwo Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017
Two Popular Bayesian Esimaors: Paricle and Kalman Filers McGill COMP 765 Sep 14 h, 2017 1 1 1, dx x Bel x u x P x z P Recall: Bayes Filers,,,,,,, 1 1 1 1 u z u x P u z u x z P Bayes z = observaion u =
More informationHPCFinance research project 8
HPCFinance research projec 8 Financial models, volailiy risk, and Bayesian algorihms Hanxue Yang Tampere Universiy of Technology March 14, 2016 Research projec 8 12/2012 11/2015, Tampere Universiy of Technology,
More information5.2. The Natural Logarithm. Solution
5.2 The Naural Logarihm The number e is an irraional number, similar in naure o π. Is non-erminaing, non-repeaing value is e 2.718 281 828 59. Like π, e also occurs frequenly in naural phenomena. In fac,
More informationSolutions to Odd Number Exercises in Chapter 6
1 Soluions o Odd Number Exercises in 6.1 R y eˆ 1.7151 y 6.3 From eˆ ( T K) ˆ R 1 1 SST SST SST (1 R ) 55.36(1.7911) we have, ˆ 6.414 T K ( ) 6.5 y ye ye y e 1 1 Consider he erms e and xe b b x e y e b
More informationOrdinary dierential equations
Chaper 5 Ordinary dierenial equaions Conens 5.1 Iniial value problem........................... 31 5. Forward Euler's mehod......................... 3 5.3 Runge-Kua mehods.......................... 36
More informationNature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time.
Supplemenary Figure 1 Spike-coun auocorrelaions in ime. Normalized auocorrelaion marices are shown for each area in a daase. The marix shows he mean correlaion of he spike coun in each ime bin wih he spike
More informationFinal Spring 2007
.615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o
More informationA Method for Setting the Artificial Boundary Conditions of Groundwater Model
Open Journal of Geology, 2013, 3, 50-54 doi:10.4236/ojg.2013.32b012 Published Online April 2013 (hp://www.scirp.org/journal/ojg) A Mehod for Seing he Arificial Boundary Condiions of Groundwaer Model Yipeng
More informationThe Potential Effectiveness of the Detection of Pulsed Signals in the Non-Uniform Sampling
The Poenial Effeciveness of he Deecion of Pulsed Signals in he Non-Uniform Sampling Arhur Smirnov, Sanislav Vorobiev and Ajih Abraham 3, 4 Deparmen of Compuer Science, Universiy of Illinois a Chicago,
More informationKriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Kriging Models Predicing Arazine Concenraions in Surface Waer Draining Agriculural Waersheds Paul L. Mosquin, Jeremy Aldworh, Wenlin Chen Supplemenal Maerial Number
More informationProbabilistic Models for Reliability Analysis of a System with Three Consecutive Stages of Deterioration
Yusuf I., Gaawa R.I. Volume, December 206 Probabilisic Models for Reliabiliy Analysis of a Sysem wih Three Consecuive Sages of Deerioraion Ibrahim Yusuf Deparmen of Mahemaical Sciences, Bayero Universiy,
More informationRobert Kollmann. 6 September 2017
Appendix: Supplemenary maerial for Tracable Likelihood-Based Esimaion of Non- Linear DSGE Models Economics Leers (available online 6 Sepember 207) hp://dx.doi.org/0.06/j.econle.207.08.027 Rober Kollmann
More informationA Specification Test for Linear Dynamic Stochastic General Equilibrium Models
Journal of Saisical and Economeric Mehods, vol.1, no.2, 2012, 65-70 ISSN: 2241-0384 (prin), 2241-0376 (online) Scienpress Ld, 2012 A Specificaion Tes for Linear Dynamic Sochasic General Equilibrium Models
More informationEnsemble Confidence Estimates Posterior Probability
Ensemble Esimaes Poserior Probabiliy Michael Muhlbaier, Aposolos Topalis, and Robi Polikar Rowan Universiy, Elecrical and Compuer Engineering, Mullica Hill Rd., Glassboro, NJ 88, USA {muhlba6, opali5}@sudens.rowan.edu
More information18 Biological models with discrete time
8 Biological models wih discree ime The mos imporan applicaions, however, may be pedagogical. The elegan body of mahemaical heory peraining o linear sysems (Fourier analysis, orhogonal funcions, and so
More informationComputer Simulates the Effect of Internal Restriction on Residuals in Linear Regression Model with First-order Autoregressive Procedures
MPRA Munich Personal RePEc Archive Compuer Simulaes he Effec of Inernal Resricion on Residuals in Linear Regression Model wih Firs-order Auoregressive Procedures Mei-Yu Lee Deparmen of Applied Finance,
More informationStochastic models and their distributions
Sochasic models and heir disribuions Couning cusomers Suppose ha n cusomers arrive a a grocery a imes, say T 1,, T n, each of which akes any real number in he inerval (, ) equally likely The values T 1,,
More informationACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H.
ACE 56 Fall 005 Lecure 5: he Simple Linear Regression Model: Sampling Properies of he Leas Squares Esimaors by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Inference in he Simple
More informationGMM - Generalized Method of Moments
GMM - Generalized Mehod of Momens Conens GMM esimaion, shor inroducion 2 GMM inuiion: Maching momens 2 3 General overview of GMM esimaion. 3 3. Weighing marix...........................................
More informationApplication of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing
Applicaion of a Sochasic-Fuzzy Approach o Modeling Opimal Discree Time Dynamical Sysems by Using Large Scale Daa Processing AA WALASZE-BABISZEWSA Deparmen of Compuer Engineering Opole Universiy of Technology
More informationImproved Approximate Solutions for Nonlinear Evolutions Equations in Mathematical Physics Using the Reduced Differential Transform Method
Journal of Applied Mahemaics & Bioinformaics, vol., no., 01, 1-14 ISSN: 179-660 (prin), 179-699 (online) Scienpress Ld, 01 Improved Approimae Soluions for Nonlinear Evoluions Equaions in Mahemaical Physics
More information0.1 MAXIMUM LIKELIHOOD ESTIMATION EXPLAINED
0.1 MAXIMUM LIKELIHOOD ESTIMATIO EXPLAIED Maximum likelihood esimaion is a bes-fi saisical mehod for he esimaion of he values of he parameers of a sysem, based on a se of observaions of a random variable
More informationBifurcation Analysis of a Stage-Structured Prey-Predator System with Discrete and Continuous Delays
Applied Mahemaics 4 59-64 hp://dx.doi.org/.46/am..4744 Published Online July (hp://www.scirp.org/ournal/am) Bifurcaion Analysis of a Sage-Srucured Prey-Predaor Sysem wih Discree and Coninuous Delays Shunyi
More informationThe complexity of climate model drifts
The complexiy of climae model drifs Davide Zanchein Angelo Rubino Maeregu Arisido Carlo Gaean Universiy of Venice, Dep. of Environmeal Sc., Informaics and Saisics A conribuion o PREFACE-WP10: (Saisical
More informationMath 10B: Mock Mid II. April 13, 2016
Name: Soluions Mah 10B: Mock Mid II April 13, 016 1. ( poins) Sae, wih jusificaion, wheher he following saemens are rue or false. (a) If a 3 3 marix A saisfies A 3 A = 0, hen i canno be inverible. True.
More informationWATER LEVEL TRACKING WITH CONDENSATION ALGORITHM
WATER LEVEL TRACKING WITH CONDENSATION ALGORITHM Shinsuke KOBAYASHI, Shogo MURAMATSU, Hisakazu KIKUCHI, Masahiro IWAHASHI Dep. of Elecrical and Elecronic Eng., Niigaa Universiy, 8050 2-no-cho Igarashi,
More informationLinear Gaussian State Space Models
Linear Gaussian Sae Space Models Srucural Time Series Models Level and Trend Models Basic Srucural Model (BSM Dynamic Linear Models Sae Space Model Represenaion Level, Trend, and Seasonal Models Time Varying
More informationZürich. ETH Master Course: L Autonomous Mobile Robots Localization II
Roland Siegwar Margaria Chli Paul Furgale Marco Huer Marin Rufli Davide Scaramuzza ETH Maser Course: 151-0854-00L Auonomous Mobile Robos Localizaion II ACT and SEE For all do, (predicion updae / ACT),
More informationSTRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN
Inernaional Journal of Applied Economerics and Quaniaive Sudies. Vol.1-3(004) STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN 001-004 OBARA, Takashi * Absrac The
More information