Time Series Analysis With Aviation Data - George Mason University

26d ago
25 Views
0 Downloads
971.10 KB
49 Pages
Last View : Today
Last Download : n/a
Upload by : Julia Hutchens
Transcription

Time Series Analysis of AviationDataDr. Richard XieFebruary, 2012

What is a Time Series A time series is a sequence of observations inchorological order, such as– Daily closing price of stock MSFT in the past tenyears– Weekly unemployment claims in the past 2 years– Monthly airline revenue passenger miles in thepast ten years Time series analysis is useful when– No other data available– System too complicated to model in detail

Where to Get the Data?

What Information Are You Interested In? How the data changes from month to month, yearto year?––––Any trend?How fluctuated the curve is?Any seasonal effects?Any unusual years/months which have significantlysmall or large number? Can we forecast future value based on the timeseries?

Let’s Work on the Data But first, what tool will you use?––––––Pencil and quadrille pad (or back of an envelope)ExcelMatlab, Mathematica, MapleSAS, SPSS, STATA, RROOT, PAW, KNIME, Data Applied, etc.Others

Use R! R is free R is a language, not just a statistical tool R makes graphics and visualization of the bestquality A flexible statistical analysis toolkit Access to powerful, cutting-edge analytics A robust, vibrant community Unlimited possibilities

Where to Download R To download R– Go to http://www.r-project.org/– Choose a CRAN Mirror, such as http://cran.cnr.berkeley.edu/– Click the link to download R according to your operatingsystem (Linux, MacOS X, or Windows) Or Enhanced Version of R Distributed by 3rd Party– Revolution R cts/revolution-r.php)– Free academic version of Revolution R ucts/revolutionenterprise.php)CRAN: Comprehensive R Archive Network

R References An Introduction to R, W.N.Venables,D.M.Smith and R Development Core Team More Documents/Tutorials, go to– http://cran.cnr.berkeley.edu/other-docs.html

Start an R Project Recommend using RStudio as the console(http://rstudio.org/download/) Create a project folder for storing R scripts, data, etc.– e.g. C:/Users/xie/Documents/SYST460/R projects/airline timeseries Open RStudio, navigate to the project folder1. Use getwd() tofind out the currentworking directory2. Click andselect theproject folder3. Click here4. Current workingdirectory is changed

Work on Data, Finally! Use Ctr Shft N to create a new script rm(list ls(all TRUE)) to clear the existingvariables in workspace, if any rpm read.csv("System Passenger Revenue Passenger Miles (Jan 1996 - Oct2011).csv") Use Ctr Enter to run the current line orselection

What the Data Looks Like? ls(rpm) plot(rpm), equivalent to plot(rpm Total rpm YYYYMM)

Plot It As A Time Series rpm.ts ts(as.numeric(rpm Total), start c(1996,01),freq 12) plot(rpm.ts,ylab 'Revenue PassengerMiles')Commands to checkproperties of rmp.ts

Trends and Seasonal Variation layout(1:2) plot(aggregate(rpm.ts)) boxplot(rpm.ts cycle(rpm.ts))Overall trendof increasingover yearsMax valueUpper QuartileMedianLower QuartileMin value

Window Function Extract a part of the time series betweenspecified start and end points rpm.Feb - window(rpm.ts, start c(1996,02), freq TRUE) rpm.Aug - window(rpm.ts, start c(1996,08), freq TRUE) mean(rpm.Feb)/mean(rpm.Aug)

Modeling Time Series - Notations

Modeling Time Series – Decomposition Models

Time Series Decomposition in R rpm.decom decompose(rpm.ts) plot(rpm.decom)

Auto-Correlation

AutoCorrelation and Correlogram rpm.acf acf(as.numeric(rpm.ts),lag.max 40)

Correlogram after Decomposition acf(as.numeric(rpm.decom random),na.action na.omit,lag.max 40)

Regression Trends: stochastic trends, deterministic trends Deterministic trends and seasonal variationcan be modeled using regression Deterministic trends are often used forprediction Time series regression differs from standardregression as time series tends to be seriallycorrelated

Linear Models

Fit A Linear Regression Model fit lm(rpm.ts time(rpm.ts)) plot(rpm.ts, type "o", ylab "RPM") abline(fit)

Diagnostic Plots hist(resid(fit)) acf(resid(fit)) AIC(fit)

Linear Model with Seasonal Variables

Linear Model with Seasonal Variables - R Seas cycle(rpm.ts)– Gives the positions in the cycle of each obsv. Time time(rpm.ts)– Creates the vector of times at which rpm.ts wassampled rpm.lm lm(rpm.ts 0 Time factor(Seas))– Fit rpm.ts to the linear model with seasonalvariables

Take a Look at rpm.lm summary(rpm.lm)ModelResiduals from fittingCoefficients

Correlogram of rpm.lm Residual acf(resid(rpm.lm),lag.max 40)It indicates strongpositive-autocorrelationResiduals are notpure randomnumbers, so itshould be furthermodeled

How Random is Random? - White Noise A time series {wt: t 1, 2, . . . , n} is discretewhite noise (DWN) if the variables w1,w2, . . .,wn are independent and identicallydistributed with a mean of zero. This implies that the variables all have thesame variance σ2 and Cor(wi,wj) 0 for all i j. If, in addition, the variables also follow anormal distribution (i.e., wt N(0, σ2)) theseries is called Gaussian white noise.

Simulate White Noise in R set.seed(1) w rnorm(100) plot(w, type "l") acf(w)

Random WalkLet {xt} be a time series. Then {xt} is a random walk if xt xt 1 wtwhere {wt} is a white noise series. Substituting xt 1 xt 2 wt 1 and then substituting for xt 2, followed by xt 3and so on gives: xt wt wt 1 wt 2 . . .In practice, the series will start at some timet 1. Hence, xt w1 w2 . . . wt

Simulate A Random Walk in R x - w - rnorm(1000) for (t in 2:1000) x[t] - x[t - 1] w[t] plot(x, type "l") acf(x)Your plot will bedifferent for sure!

Stationarity

Auto-Regressive (AR) Models The series {xt} is an autoregressive process oforder p, abbreviated to AR(p), ifxt α1xt 1 α2xt 2 . . . αpxt p wtwhere {wt} is white noise and the αi are themodel parameters with αp 0 for an order pprocess. The model is a regression of xt on past termsfrom the same series; hence the use of theterm ‘autoregressive’.

Fit an AR Model res.ar ar(resid(rpm.lm),method "ols") res.arxt 0.7153xt 1 0.1010xt 2 0.0579xt 3 wt acf(res.ar res[-(1:3)])OLS: Ordinary Least Square

Moving Average Model

Simulate an MA Process set.seed(1) b - c(0.8, 0.6, 0.4) x - w - rnorm(1000) for (t in 4:1000) {for (j in 1:3) x[t] - x[t] b[j] * w[t - j]} plot(x, type "l") acf(x)

Plot ResultsFirst 3 lags withsignificant correlationSimulated MA ProcessCorrelogram of the Simulated MA Process

ARMA Model

Fit An ARMA Model rpm.arma arima(rpm.ts,order c(1,0,1)) acf(as.numeric(rpm.arma resid))Strong seasonalinformation left inresiduals

ARIMA and SARIMA Model A time series {xt} follows an ARIMA(p, d, q)process if the dth differences of the {xt} seriesare an ARMA(p, q) process SARIMA is Seasonal ARIMA model whichextends ARIMA model with seasonal terms

Fit A SARIMA Model rpm.arima arima(rpm.ts,order c(1,1,1),seas list(order c(1,0,0),12))– First c(1,1,1): AR(1), first-order difference, MA(1)– Second c(1,0,0): seasonal terms on AR process, frequency12 acf(as.numeric(rpm.arima resid))

Forecast Predicting future values of a time series, xn m,using the set of present and past values of thetime series, x {xn, xn-1, , x1} The minimum mean square error predictor ofxn m is xnn m E ( xn m x) predict(model, newdata) method– model: a model object used for prediction– newdata: value of explanatory variables

Forecast in R new.t - seq(2011.750, len 2 * 12, by 1/12) new.dat - data.frame(Time new.t, Seas rep(1:12, 2)) rpm.pred ts(predict(rpm.lm,new.dat)[1:24],start c(2011,11),freq 12) ts.plot(rpm.ts,rpm.pred,lty 1:2)

Plot of Forecast

Forecast with SARIMA Model ts.plot( cbind( window(rpm.ts,start c(1996,1)), predict(rpm.arima,48) pred ), lty 1:2)

Homework Download and Install R with Rstudio Read An Introduction to R Download System Passenger - RevenueAircraft Miles Flown (000) (Jan 1996 - Oct2011) data from BTS Read the data into R using Rstudio Create a time series plot of the data, and plotits auto-correlation correlogram Decompose the time series and save the plot

Homework (Ctd.) Construct a linear regression model withoutseasonal factors, and plot the correlogram ofthe model’s residual data Construct a linear regression model withseasonal factors, and identifies thecharacteristics of the model residual. Fit an AR model to the model residual of theabove model Forecast the time series data into next 24months using the seasonal model

Exam Questions What are the data elements in a time series?What does auto-correlation mean?What are white noise and random walk?What are stationary models and non-stationarymodels?

Homework (Ctd.) Construct a linear regression model without seasonal factors, and plot the correlogram of the model's residual data Construct a linear regression model with seasonal factors, and identifies the characteristics of the model residual. Fit an AR model to the model residual of the above model

Related Documents:

SAF (Sustainable Aviation Fuel) a.k.a. aviation biofuel, biojet, alternative aviation fuel. Aviation Fuel: Maintains the certification basis of today’s aircraft and jet (gas turbine) engines by delivering the properties of ASTM D1655 – Aviation Turbine Fuel – enables drop-in approach – no changes to infrastructure or equipment,

Below are some Aviation Fun Facts to celebrate National Aviation History Month: National Aviation Day, August 19, is a United States national observation that celebrates the history and development of aviation. It was established in 1939 by Franklin Delano Roosevelt, who issued a presidential proclamation which designated the anniversary of .

AVIATION SAFETY Challenges and ways forward for a safe future 3 TABLE OF CONTENTS EXECUTIVE SUMMARY 5 CURRENT AVIATION-SAFETY CHALLENGES 7 The ‘big five’ aviation-safety challenges 9 . Safe. Research & Innovation Projects for Policy AVIATION SAFETY AT. AVIATION SAFETY 20.

1. Aviation Industry Context in the Digital Age 7 1.1 Adoption of emerging technologies in the aviation industry 7 1.2 What is at stake? 7 1.3 Need for global collaboration to enhance cyber resilience in aviation 9 2. Aviation Industry Analysis 10 2.1 Global risk insights 10 2.2 Key insights from industry stakeholders 10 2.3 Threats 11

The government is developing a long term Aviation Strategy to 2050 and beyond, the aim of which is to achieve a safe, secure and sustainable aviation sector that meets the needs of consumers and of a global, outward-looking Britain. The objectives of the strategy are to: help the aviation industry work for its customers

This is a workshop in practical aviation for teachers of the South Carolina public schools. Participants will be given a general orienta tion in aviation and its related fields. Emphasis will be placed on ways of introducing aviation into the classroom, the use of aviation

SMP Series page 73 FAKRA Connectors page 77 BNC Series page 79 TNC Series page 108 N Series page 133 7/16 Series page 149 UHF/MINI-UHF Series page 159 F Series page 167 Twin Series page 175 D-sub Series page 179 FME Series page 181 1.0/2.3 Series page 183 1.6/5.6 Series page 189 Filtered Series page 197

5 Series E39 5 Series E60/E61 5 Series F10/F11/F18 5 Series GT(F07) 6 Series E63/E64 6 Series F06/F12/F13 7 Series E38 7 Series E65/E66/E68 7 Series F01/F02/F03/F04 I Series I01/I12 X Series X3_E83 X Series X5_E

Forecasting economic time series using unobserved components time series models . time series analysis implies a speci c approach to the modelling of time series. It is somewhat di erent compared to the Box-Jenkins analysis. . and quarterly frequencies of time series), outliers, structural breaks and

ATR-72/42/B-1900 Local Time Days of Operation Civil Aviation Authority of Nepal Civil Aviation Safety Regulation Directorate Special Domestic Winter Schedule 2021 Revised Operator : Buddha Air (P.) Ltd. Effective Date: 8 February to 31 March 2021 SN Flight No. A/C Type Sector Fare in NRS ATR-72/42/B-1900

The Initial Approach Fix Since July 1955, this magazine has been the voice of naval-aviation safety—your voice. Back then, it was called “The Approach, U.S. Naval Aviation Safety Review,” and was introduced to the fleet during a time when the aviation-Class

Police Aviation News 276 April 2019 2 EDITORIAL This short period experiment called Police Aviation News – PAN for short - has now been around for an unexpected 23 years, a time far longer than env

was accurate at the time of printing. However, this document is a guide and readers are reminded that the Defence Aviation Safety Assurance Manual (DASAMAN) and the Defence Aviation Safety Regulations (DASR) remain the authoritative documents for the management and regulation of aviation safety in the Defence. Reserved sections within

Caiado, J., Maharaj, E. A., and D’Urso, P. (2015) Time series clustering. In: Handbook of cluster analysis. Chapman and Hall/CRC. Andrés M. Alonso Time series clustering. Introduction Time series clustering by features Model based time series clustering Time series clustering by dependence Introduction to clustering

Aviation Transport of Ukraine Aviation Register of Interstate Aviation Committee (ARMAK) Certificates No. SPR-11, No. SPR-15, No. R-56, No. R-69 and others Bureau Veritas (France) Certificate EN 9100:2009 No. FR015515-1, EN 9110:2009 No. FR015516-1, ISO 9001:2008 No. UA227484 European Aviation (Germany) Certificate No. 216/2008, No. 1702/2003 .

Sustainable Aviation Fuels Aviation Kerosene Jet-A1 Aviation Gasoline AvGas Natural Gas Liquified Hydrogen . Catalytic hydrothermolysis jet fuel (CHJ), a type of synthetic kerosene 50% . (IATA/ICAO) Roadmap Jet market 2019 (Actual) Jet market 2050 (Projected) .

3 “Aviation is a small but important contributor to climate change. ICAO/CAEP’s initial estimate is that the total volume of aviation CO2 emissions in 2006 (both domestic and international) is in the range of 600 million tonnes. At present, aviation accounts for about 2% of total global CO2 emissions and about 12% of

The 39th Assembly is long anticipated to adopt a global measure to arrest gaseous emissions from international civil aviation. Further, the upcoming Assembly session is slated to discuss and adopt the Global Aviation Safety Plan and Global Air Navigation Plan: respectively ICAO, 2017-2019 Global Aviation Safety Plan, 2nd ed

Director of the European Aviation Safety Agency of 22.11.2007. 3 Regulation (EC) No 216/2008 of the European Parliament and of the Council of 20 February 2008 on common rules in the field of civil aviation and establishing a European Aviation Safety Agency, and

deserves a special place in our aviation history. The “Springfield Cadet” with original A-40 engine and . efforts to preserve and celebrate Oregon Aviation History. Fill in the blanks and return to: OR Aviation Historical Society . As far as the government and national media were concerned, the official starting point of the