Cross-Sectional Vs. Time Series Benchmarking in Small Area Estimation; Which Approach Should We Use? Danny Pfeffermann
|
|
- Mitchell Jenkins
- 6 years ago
- Views:
Transcription
1 Cross-Sectional Vs. Time Series Benchmarking in Small Area Estimation; Which Approach Should We Use? Danny Pfeffermann Joint work with Anna Sikov and Richard Tiller Graybill Conference on Modern Survey Statistics, Colorado State University, Fort Collins, 2013
2 What is benchmarking? Y - target characteristic in area d at time t, y- direct survey estimate, Y ˆ model - estimate obtained under a model. d 1,2,..., D Areas, t 1,2,... Time Benchmarking: modify model based estimates to satisfy: ˆ D model by d 1 B t ; 12 t =,,... ( B t known, e.g., B D b y ). t d 1 b fixed coefficients (relative size, scale factors, ). Requirement: B t sufficiently close to true value D by d 1. When B t function of { y} internal benchmarking. 2
3 Why benchmark? 1. The area model-dependent estimators may not agree with the direct survey estimator for an aggregate of the small areas, for which the direct estimator is deemed accurate. Guarantees consistency of publication between modeldependent small area estimators and design-based estimator in a larger area, which is required by statistical bureaus 2. Corrects for the bias resulting from model misspecification? 3
4 3. Problems considered in present presentation Consider e.g., a repeated monthly survey with sample estimates y in area d at month t. Denote by Yˆ time-series ( Yˆ cross-section ) the optimal unbenchmarked predictor under an appropriate times series (cross-sectional) model. 1- Should we benchmark, Yˆ time-series or Yˆ cross-section? 2- Does it matter which cross-sectional benchmarking method we use? Distinguish between correct model specification and when the model is misspecified. 4
5 Problems considered in present presentation (cont.) 3- Develop a time series two-stage benchmarking procedure for small area hierarchical estimates. First stage: benchmark concurrent model-based estimators at higher level of hierarchy to reliable aggregate of corresponding survey estimates. Second stage: benchmark concurrent model-based estimates at lower level of hierarchy to first stage benchmarked estimate of higher level to which they belong. 5
6 Example: Labour Force estimates in the U.S.A 6
7 Content of presentation 1. Review cross-sectional and times series benchmarking procedures for small area estimation; 2. Compare cross-sectional and times series benchmarking: 2.a Under correct model specification 2.b When models are misspecified 3. Review two-stage cross-sectional and times series benchmarking procedures; (if time allows!!) 4. Apply time series two-stage benchmarking to monthly unemployment series in the U.S.A. 7
8 Single- stage cross-sectional benchmarking procedures (Pfeffermann & Barnard 1991): No BLUP benchmarked predictor exists!! Pro-rata (ratio) benchmarking (in common use), Limitations: Y ˆ Y ˆ b y b Y ˆ ) ; bmk model D D model d, R d / d1 d d k1 k k 1- Adjusts all the small area model-based estimates the same way, irrespective of their precision, 2- Benchmarked estimates not consistent: if sample size in area d increases but sample sizes in other areas unchanged, Yˆ bmk d, R does not converge to true population value Y d. 8
9 Additive single-stage cross-sectional benchmarking Y ˆ Y ˆ ( b y by ˆ ); bmk model D D model d, d 1 k k A d k k1 k k D b d 1 d d 1. Coefficients { d } measure precision (next slide); distribute difference between benchmark and aggregate of model-based estimates between the areas. Unbiased under correct model. If d 0 nd model ˆbmk Y Yˆ A Y consistent. d, d d If ˆ model k Plim( Y y ) 0 Area k with accurate estimate n k k does not contribute to benchmarking in other areas. ˆ bmk Easy to estimate variance of Yd, A. 9
10 Examples of additive cross-sectional benchmarking Wang et al. (2008) minimize model s.t. by ˆ d d by d d, D D bmk d1 d1 ( ˆ ) D bmk φ EY d 1 d Y d d, A A. Sol.: d 2 under F-H 1 D 1 2 d bd / k 1 k bk λ. { φ d } represent precision of direct or model-based estimators. ˆ model 1 φ d [ Var( Y d )] Battese et al model D ˆ ˆmodel 1 φ bd[ Cov( Yd, d Y )] 1 k Pfeffermann & Barnard φd [ Var( y d )] Isaki et al k In practice, model parameters replaced by estimates. Bell et al. (2012) extend results to multiple constraints. 10
11 Self benchmarking and other benchmarking procedures You and Rao (2002) propose self-benchmarked estimators for unit-level model by modifying the estimator of β. Approach extended by Wang et al. (2008) to F-H model. Wang et al. (2008) consider alternative self-benchmarking, 2 x b σ. obtained by extending the covariates x i to (x, ) 11 i i i Di Ugarte et al. (2009) benchmark the BLUP under unit-level model to synthetic estimator for all areas under regression model with heterogeneous variances. Lahiri (1990), You et al. (2004), Datta et al. (2011), Nandram and Sayit (2011) and Steorts and Ghosh (2013) consider Bayesian & empirical Bayesian benchmarking.
12 Small area single-stage time series benchmarking Basic idea: Fit a time series model to direct estimators in all the areas, incorporate benchmark constraints into model equations. Variances obtained as part of model fitting. Pfeffermann & Tiller (2006) consider the following model for total unemployment in census divisions. yt ( y t,, ydt) = = sampling errors. Yt ( Y1 t,..., YDt) = true division totals, 1 direct estimates, et ( e1 t,, edt) 2 2 y Y e; E( e ) 0, E e e Σ Diag[,, ] t t t t t τt 1,, t D,, t. Division sampling errors independent between divisions but highly auto-correlated within a division, and heteroscedastic. Sampling rotation: (4 in, 8 out, 4 in). 12
13 Time series model for division d Totals Y assumed to evolve independently between divisions according to basic structural model (BSM, Harvey 1989). Model accounts for covariates, stochastic trend, stochastically varying seasonal effects and random irregular terms. Model written as: Y z; Tdd, t1 (state-space) Errors mutually independent white noise, E( ) Q. d ARIMA, regression with random coefficients and unit & area level models can all be expressed in state-space form. 13
14 Combining the separate division models y Y e Z e (measurement equation) ; t t t t t t (,..., ) 1, t t D t t t1 t t D D d T (state equation) ; Z z, T T ; - block diagonal E 0, E Q Q, E 0, t. t t t D d t Benchmark constraints D D MODEL D b y b z b Y d1 d1 d1 be =0, t = 1,2,... D d=1 (imposed in model equations). 14
15 Time series benchmarked predictor Denote by α ˆt,u the state predictor at time t without imposing D the constraint be = 0 at time t, but imposing the d=1 constraints in previous time points. Define, [ ( ˆ ft Var bz, u )] ; D d 1 Cov[( ˆ ), b z ( ˆ )]. dft, u d 1, u D ˆ bmk z 1 ˆ z D D Y ˆ, u dft b y b z ft d 1 d 1, u D = 1 b d 1 z dft ft 1. ( ). 15
16 Properties of time series benchmarked predictor Y ˆ z ˆ z ( b y b z ˆ ) bmk 1 D D, u dftft d 1 d 1, u. Y ˆ bmk member of cross-sectional benchmarked predictors, Y ˆ Y ˆ λ ( b y by ˆ ) (Wang et al. 2008). bmk model D D model d, d k 1 k k A d k1 k k ˆ model Y = z αˆ un-benchmarked predictor at time t;,u b ; λ b 1 D ˆmodel ˆmodel φ 1 2 = =b Cov Y D -Y, b k=1 kt Ykt Y kt zδdft k 1 kt b kt Pfeffermann & Barnard (1991). { [( ) ( )]} -1 16
17 Properties of time series benchmarked predictor (cont.) (a)- unbiasedness: if ˆ bmk EY ( d,t-1 Yd, t 1) 0 ( ˆ bmk EY Y ) 0. To warrant unbiasedness under the model, it suffices to initialize at time t= 1 with unbiased predictor. (b)- Consistency: Plim( y Y ) 0 & Plim( ˆ bmk Y y ) 0 n d Plim( ˆ bmk Y Y ) 0 (even if model misspecified). n d (c) ( ˆ bmk Var Y Y ) accounts for variances and covariances of sampling errors, variances and covariances of benchmark errors D d1 d1 D b y b z and their covariances with division sampling errors, and variances of model components. n d 17
18 Comparison: cross-sectional & time series benchmarking Generate series { y, t 1,2,...,75} for areas i 1,...,25 independently between the areas from model: ti y ti Y e ; e N(0, ) ti Y x L S ; ti ti ti ti 1 ti ti, ee tti E ti i t i, t, 1,2,...,75. L L R ; R R, ti t1, i t1, i Lti ti t1, i Rti N N 2 2 Lti (0, L), Rti (0, R) ; S cos S sin S ; ~ N(0, ) * 2 jti j jt, 1, i j jt, 1, i jti jti S S sin S cos S ; ~ N(0, ) * * * * 2 jti j j, t1, i j j, t1, i jti jti S 6 S S ; 2 j/12, j ti jti j j1 Used by BLS for production of Labor Force estimates. 18
19 Cross-sectional model operating at time t By repeated substitutions, for fixed starting values { L, R, S, S, j 1,...,6} (same in all areas), the cross-sectional * 0 0 j0 j0 model operating in area i at time t is the F-H model y x u e, ti 0 1 ij ti ti j 6 5 * 0 L0tR0( Sj0cos tsj0sin t), j1 6 j1 6 t t 6 t 5 t j * j ti Lli ( ) Rli jli cos ( ) jli sin ( ) l1 l1 j1 l1 6 j1 l1 6 u t l t l t l Eu ( ti ) 0, β 0 and 19 j t ( ti) ut L R ( ) S6 l1 Var u t t l t. σ constant across the areas for given t. 2 ut.
20 Simulation set up Divided the 25 areas into 5 equal groups. 2 Var( ) 0, Var( ) , 2 4 R Rti Var( ) S Sti Variances of similar magnitude to variances estimated by BLS for division total unemployment series (after scaling). L Lti Sampling errors e ti in first 4 groups generated from AR(1) models eti et 1, i ti, with different s and 2 e Var( eti) in 2 2 different groups: ( ρ 1,σ 1e ) = (0.15,3), ( ρ 2,σ 2e ) = (0.75,3), ( ) (0.15,5.5) 2 2 ρ 3,σ 3e =, 4 4e ( ρ,σ ) = (0.75,5.5). 20
21 Simulation set up (cont.) For group 5 used variances and autocorrelations estimated for scaled total unemployment series in division of West North Central. 4.4 tti First three autocorrelations are 0.38, 0.21, Autocorrelations at lags around β 1 = 1.5 and started with R 0, L 1, S S 0, j 1,...,6, * j 0 j such that E( Y ) = (1+1.5 x ) for all ( t,i ). ti i We compare the bias and RPMSE at t=25,50,75. 21
22 Simulation set up (cont.) For each method re-estimated model parameters at each of three time points with 2 u estimated by F-H method and BSM variances estimated by mle. True σ at three time points are σ , σ , σ u Variances & autocorrelations of sampling errors known. Same model parameters for five areas in each group. Used true U.S.A State x-values. Results in figures below are group averages of empirical bias and RPMSE over 400 simulated series. 22
23 Root Prediction MSE under correct model specification 23
24 Root Prediction MSE under correct model specification 24
25 Root Prediction MSE under correct model specification 25
26 Prediction bias when excluding covariate from the model 26
27 Prediction bias when excluding covariate from the model 27
28 Prediction bias when excluding covariate from the model 28
29 Mean R- Square when regressing y on x Group Time σ ut t= t= t=
30 Root PMSE when excluding covariate from the model 30
31 Root PMSE when excluding covariate from the model 31
32 Root PMSE when excluding covariate from the model 32
33 Conclusions 1- The cross-sectional procedures perform similarly and with 25 areas, there is no loss in efficiency from benchmarking. 2- Cross-sectional benchmarking does not correct for bias. 3- For sufficiently long time series, the extra efforts of fitting time series models and benchmarking the predictors are well rewarded in terms of the RPMSE. 4- The length of time series required for dominating the crosssectional benchmarked methods depends on variances of sampling errors and correlations between them. Time series predictors can be enhanced by smoothing. 33
34 Two-stage time series benchmarking First stage: benchmark concurrent model-based estimates at higher level of hierarchy to reliable aggregate of corresponding survey estimates. Second stage: benchmark concurrent model-based estimates at lower level of hierarchy to first stage benchmarked estimate of higher level to which they belong. Ghosh and Steorts (2013) consider Bayesian two-stage cross-sectional benchmarking, applied in a single run. 34
35 Why not benchmark second level areas in one step under time series benchmarking? 1- May not be feasible in real time production: For USA-CPS series, proposed procedure requires joint modeling of all the areas state-space model of order Delay in processing data for one second level area could hold up all the area estimates. 3- When 1 st level hierarchy composed of homogeneous 2 nd level areas, benchmarking more effectively tailored to 1 st level characteristics. 35
36 Second-stage time series benchmarking Suppose S states in Division d; Direct Target y Y e ; E( e ) 0, Cov( e, e ) 2 ds, t ds, t ds, t ds, t ds, ds*, t s, s* ds, t Y z ds, t ds, t ds, t ; T, E( ) 0, E( ) Q ds, t ds ds, t1 ds, t ds, t ds, t ds,* t t,* t ds, t S Benchmark: b y,, b z ˆ,,, s1 ds t ds t s1 ds t ds t ds t S Y ˆ bmk Guarantees that model-based estimates at lower level add up to published estimate at higher level to which they belong. bmk Benchmark error: ˆ bmk bmk S r ( Y Y ) ( z ˆ b z ds, t ds, tds, t ). s1 36
37 Benchmarking of second-level estimates (cont.) (,...,, ˆ bmk y y y Y ) d t d1, t ds, t d (, ˆ bmk y Y ) t bmk e ( e,..., e, r ) ( e, r ). d bmk d t d1, t ds, t t (,..., ), d t d1, t ds, t d S T, d ds Zt S z ds, t T, d d Z t Z t. bd1, tz d1, t,..., bds, tz ds, t Combined model: y Z e ; T d d d d d d d d t t t t t t1 t E( ) Q Q ; E( e e ) d d d d d t t S ds, t d t t. 37
38 Empirical results Total unemployment, CPS-USA, Jan Dec2009. First level- Census divisions, Second level- States. Illustrate: 1- Consistency of benchmarked predictors 2- Robustness to model failure 38
39 Consistency of benchmarked predictors S.E.( y ) /y (left) Y ˆ bmk / y (right) ds,t ds,t Massachusetts ds,t ds,t Jan-90 Jan-94 Jan-98 Jan-02 Jan Jan-90 Jan-94 Jan-98 Jan-02 Jan-06 39
40 Consistency of benchmarked predictors [S.E.( y ds,t ) /y ds,t] (left) [ Y ˆ bmk ds,t / y ds,t ] (right) New Hampshire Jan-90 Jan-94 Jan-98 Jan-02 Jan Jan-90 Jan-94 Jan-98 Jan-02 Jan-06 40
41 Direct, Benchmarked and Unbenchmarked estimates of Total Unemployment, Georgia (numbers in 000 s) CPS BMK UNBMK Jan-00 Jan-02 Jan-04 Jan-06 Jan-08 41
42 Direct, Benchmarked and Unbenchmarked estimates of Total Unemployment, Alabama (numbers in 000 s). 170 CPS BMK UNBMK Jan-00 Jan-02 Jan-04 Jan-06 Jan-08 42
43 Conclusions 1. Two-stage time series benchmarking consistent and protects against abrupt changes in time series model. 2. Due to bias correction, two-stage time series benchmarking occasionally produces predictors with lower RPMSE than unbenchmarked predictors (not illustrated). 3. Methodology can be extended to other time series models or to three- or higher-stage benchmarking. 4. Proposed method produces variances of benchmarked predictors as part of the model fitting. 43
Two-stage Benchmarking of Time-Series Models for. Small Area Estimation. Danny Pfeffermann, Richard Tiller
Two-sage Benchmarking of Time-Series Moels for Small Area Esimaion anny Pfeffermann, Souhampon Universiy, UK & Hebrew universiy, Israel Richar Tiller Bureau of Labor Saisics, U.S.A. Small Area Conference,
More informationContextual Effects in Modeling for Small Domains
University of Wollongong Research Online Applied Statistics Education and Research Collaboration (ASEARC) - Conference Papers Faculty of Engineering and Information Sciences 2011 Contextual Effects in
More informationMultiple Linear Regression
Multiple Linear Regression Asymptotics Asymptotics Multiple Linear Regression: Assumptions Assumption MLR. (Linearity in parameters) Assumption MLR. (Random Sampling from the population) We have a random
More informationA measurement error model approach to small area estimation
A measurement error model approach to small area estimation Jae-kwang Kim 1 Spring, 2015 1 Joint work with Seunghwan Park and Seoyoung Kim Ouline Introduction Basic Theory Application to Korean LFS Discussion
More informationStatistics 910, #5 1. Regression Methods
Statistics 910, #5 1 Overview Regression Methods 1. Idea: effects of dependence 2. Examples of estimation (in R) 3. Review of regression 4. Comparisons and relative efficiencies Idea Decomposition Well-known
More informationTime series models in the Frequency domain. The power spectrum, Spectral analysis
ime series models in the Frequency domain he power spectrum, Spectral analysis Relationship between the periodogram and the autocorrelations = + = ( ) ( ˆ α ˆ ) β I Yt cos t + Yt sin t t= t= ( ( ) ) cosλ
More informationEcon 424 Time Series Concepts
Econ 424 Time Series Concepts Eric Zivot January 20 2015 Time Series Processes Stochastic (Random) Process { 1 2 +1 } = { } = sequence of random variables indexed by time Observed time series of length
More informationF9 F10: Autocorrelation
F9 F10: Autocorrelation Feng Li Department of Statistics, Stockholm University Introduction In the classic regression model we assume cov(u i, u j x i, x k ) = E(u i, u j ) = 0 What if we break the assumption?
More informationØkonomisk Kandidateksamen 2004 (I) Econometrics 2. Rettevejledning
Økonomisk Kandidateksamen 2004 (I) Econometrics 2 Rettevejledning This is a closed-book exam (uden hjælpemidler). Answer all questions! The group of questions 1 to 4 have equal weight. Within each group,
More informationIntroduction to Regression Analysis. Dr. Devlina Chatterjee 11 th August, 2017
Introduction to Regression Analysis Dr. Devlina Chatterjee 11 th August, 2017 What is regression analysis? Regression analysis is a statistical technique for studying linear relationships. One dependent
More informationINDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -12 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.
INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY Lecture -12 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. Summary of the previous lecture Data Extension & Forecasting Moving
More informationSmall Area Estimation with Uncertain Random Effects
Small Area Estimation with Uncertain Random Effects Gauri Sankar Datta 1 2 3 and Abhyuday Mandal 4 Abstract Random effects models play an important role in model-based small area estimation. Random effects
More informationEstimating and Identifying Vector Autoregressions Under Diagonality and Block Exogeneity Restrictions
Estimating and Identifying Vector Autoregressions Under Diagonality and Block Exogeneity Restrictions William D. Lastrapes Department of Economics Terry College of Business University of Georgia Athens,
More informationEcon 423 Lecture Notes
Econ 423 Lecture Notes (hese notes are modified versions of lecture notes provided by Stock and Watson, 2007. hey are for instructional purposes only and are not to be distributed outside of the classroom.)
More informationIntroduction to Econometrics
Introduction to Econometrics STAT-S-301 Introduction to Time Series Regression and Forecasting (2016/2017) Lecturer: Yves Dominicy Teaching Assistant: Elise Petit 1 Introduction to Time Series Regression
More informationTIME SERIES MODELS FOR STATE LABOR FORCE ESTIMATES
TIME SERIES MODELS FOR STATE LABOR FORCE ESTIMATES Thomas D. Evans, Richard B. Tiller, and Tamara Sue Zimmerman, Bureau of Labor Statistics Tamara Zimmerman, Bureau of Labor Statistics, Room 4985, 2 Mass.
More information11.1 Gujarati(2003): Chapter 12
11.1 Gujarati(2003): Chapter 12 Time Series Data 11.2 Time series process of economic variables e.g., GDP, M1, interest rate, echange rate, imports, eports, inflation rate, etc. Realization An observed
More informationEmpirical Market Microstructure Analysis (EMMA)
Empirical Market Microstructure Analysis (EMMA) Lecture 3: Statistical Building Blocks and Econometric Basics Prof. Dr. Michael Stein michael.stein@vwl.uni-freiburg.de Albert-Ludwigs-University of Freiburg
More informationRegression of Time Series
Mahlerʼs Guide to Regression of Time Series CAS Exam S prepared by Howard C. Mahler, FCAS Copyright 2016 by Howard C. Mahler. Study Aid 2016F-S-9Supplement Howard Mahler hmahler@mac.com www.howardmahler.com/teaching
More informationTime-series small area estimation for unemployment based on a rotating panel survey
Discussion Paper Time-series small area estimation for unemployment based on a rotating panel survey The views expressed in this paper are those of the author and do not necessarily relect the policies
More informationOutline. Nature of the Problem. Nature of the Problem. Basic Econometrics in Transportation. Autocorrelation
1/30 Outline Basic Econometrics in Transportation Autocorrelation Amir Samimi What is the nature of autocorrelation? What are the theoretical and practical consequences of autocorrelation? Since the assumption
More informationHeteroscedasticity and Autocorrelation
Heteroscedasticity and Autocorrelation Carlo Favero Favero () Heteroscedasticity and Autocorrelation 1 / 17 Heteroscedasticity, Autocorrelation, and the GLS estimator Let us reconsider the single equation
More informationChapter 3 - Temporal processes
STK4150 - Intro 1 Chapter 3 - Temporal processes Odd Kolbjørnsen and Geir Storvik January 23 2017 STK4150 - Intro 2 Temporal processes Data collected over time Past, present, future, change Temporal aspect
More informationA time series is called strictly stationary if the joint distribution of every collection (Y t
5 Time series A time series is a set of observations recorded over time. You can think for example at the GDP of a country over the years (or quarters) or the hourly measurements of temperature over a
More informationClassic Time Series Analysis
Classic Time Series Analysis Concepts and Definitions Let Y be a random number with PDF f Y t ~f,t Define t =E[Y t ] m(t) is known as the trend Define the autocovariance t, s =COV [Y t,y s ] =E[ Y t t
More informationInverse of a Square Matrix. For an N N square matrix A, the inverse of A, 1
Inverse of a Square Matrix For an N N square matrix A, the inverse of A, 1 A, exists if and only if A is of full rank, i.e., if and only if no column of A is a linear combination 1 of the others. A is
More informationSTAT 520: Forecasting and Time Series. David B. Hitchcock University of South Carolina Department of Statistics
David B. University of South Carolina Department of Statistics What are Time Series Data? Time series data are collected sequentially over time. Some common examples include: 1. Meteorological data (temperatures,
More informationModel Mis-specification
Model Mis-specification Carlo Favero Favero () Model Mis-specification 1 / 28 Model Mis-specification Each specification can be interpreted of the result of a reduction process, what happens if the reduction
More informationIntroductory Econometrics
Based on the textbook by Wooldridge: : A Modern Approach Robert M. Kunst robert.kunst@univie.ac.at University of Vienna and Institute for Advanced Studies Vienna December 17, 2012 Outline Heteroskedasticity
More informationAn estimate of the long-run covariance matrix, Ω, is necessary to calculate asymptotic
Chapter 6 ESTIMATION OF THE LONG-RUN COVARIANCE MATRIX An estimate of the long-run covariance matrix, Ω, is necessary to calculate asymptotic standard errors for the OLS and linear IV estimators presented
More informationRegression #3: Properties of OLS Estimator
Regression #3: Properties of OLS Estimator Econ 671 Purdue University Justin L. Tobias (Purdue) Regression #3 1 / 20 Introduction In this lecture, we establish some desirable properties associated with
More informationEconometrics Summary Algebraic and Statistical Preliminaries
Econometrics Summary Algebraic and Statistical Preliminaries Elasticity: The point elasticity of Y with respect to L is given by α = ( Y/ L)/(Y/L). The arc elasticity is given by ( Y/ L)/(Y/L), when L
More informationEconometrics of Panel Data
Econometrics of Panel Data Jakub Mućk Meeting # 9 Jakub Mućk Econometrics of Panel Data Meeting # 9 1 / 22 Outline 1 Time series analysis Stationarity Unit Root Tests for Nonstationarity 2 Panel Unit Root
More informationTIME SERIES ANALYSIS AND FORECASTING USING THE STATISTICAL MODEL ARIMA
CHAPTER 6 TIME SERIES ANALYSIS AND FORECASTING USING THE STATISTICAL MODEL ARIMA 6.1. Introduction A time series is a sequence of observations ordered in time. A basic assumption in the time series analysis
More informationarxiv: v2 [math.st] 20 Jun 2014
A solution in small area estimation problems Andrius Čiginas and Tomas Rudys Vilnius University Institute of Mathematics and Informatics, LT-08663 Vilnius, Lithuania arxiv:1306.2814v2 [math.st] 20 Jun
More informationGLS. Miguel Sarzosa. Econ626: Empirical Microeconomics, Department of Economics University of Maryland
GLS Miguel Sarzosa Department of Economics University of Maryland Econ626: Empirical Microeconomics, 2012 1 When any of the i s fail 2 Feasibility 3 Now we go to Stata! GLS Fixes i s Failure Remember that
More informationAutocorrelation. Think of autocorrelation as signifying a systematic relationship between the residuals measured at different points in time
Autocorrelation Given the model Y t = b 0 + b 1 X t + u t Think of autocorrelation as signifying a systematic relationship between the residuals measured at different points in time This could be caused
More informationThe Ins and Outs of UK Unemployment
University of Warwick 9 July 2009 I Do in ows or out ows drive unemployment dynamics? I Spate of recent US work inspired by Hall (2005) and Shimer (2005, 2007) s claim that job nding dominates and separations
More informationFractional Hot Deck Imputation for Robust Inference Under Item Nonresponse in Survey Sampling
Fractional Hot Deck Imputation for Robust Inference Under Item Nonresponse in Survey Sampling Jae-Kwang Kim 1 Iowa State University June 26, 2013 1 Joint work with Shu Yang Introduction 1 Introduction
More informationReliability and Risk Analysis. Time Series, Types of Trend Functions and Estimates of Trends
Reliability and Risk Analysis Stochastic process The sequence of random variables {Y t, t = 0, ±1, ±2 } is called the stochastic process The mean function of a stochastic process {Y t} is the function
More informationUnivariate ARIMA Models
Univariate ARIMA Models ARIMA Model Building Steps: Identification: Using graphs, statistics, ACFs and PACFs, transformations, etc. to achieve stationary and tentatively identify patterns and model components.
More informationOn Modifications to Linking Variance Estimators in the Fay-Herriot Model that Induce Robustness
Statistics and Applications {ISSN 2452-7395 (online)} Volume 16 No. 1, 2018 (New Series), pp 289-303 On Modifications to Linking Variance Estimators in the Fay-Herriot Model that Induce Robustness Snigdhansu
More informationA Data-Driven Model for Software Reliability Prediction
A Data-Driven Model for Software Reliability Prediction Author: Jung-Hua Lo IEEE International Conference on Granular Computing (2012) Young Taek Kim KAIST SE Lab. 9/4/2013 Contents Introduction Background
More informationfor Sparsely Observed Functional Data University of Glasgow
Spatial Correlation Estimation for Sparsely Observed Functional Data Workshop on High dimensional and dependent functional data, Bristol, Sep 2012 Surajit Ray University of Glasgow Contributors t Chong
More informationEcon 300/QAC 201: Quantitative Methods in Economics/Applied Data Analysis. 17th Class 7/1/10
Econ 300/QAC 201: Quantitative Methods in Economics/Applied Data Analysis 17th Class 7/1/10 The only function of economic forecasting is to make astrology look respectable. --John Kenneth Galbraith show
More information1. You have data on years of work experience, EXPER, its square, EXPER2, years of education, EDUC, and the log of hourly wages, LWAGE
1. You have data on years of work experience, EXPER, its square, EXPER, years of education, EDUC, and the log of hourly wages, LWAGE You estimate the following regressions: (1) LWAGE =.00 + 0.05*EDUC +
More informationAustrian Inflation Rate
Austrian Inflation Rate Course of Econometric Forecasting Nadir Shahzad Virkun Tomas Sedliacik Goal and Data Selection Our goal is to find a relatively accurate procedure in order to forecast the Austrian
More informationINDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -18 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.
INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY Lecture -18 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. Summary of the previous lecture Model selection Mean square error
More informationNonparametric Small Area Estimation via M-quantile Regression using Penalized Splines
Nonparametric Small Estimation via M-quantile Regression using Penalized Splines Monica Pratesi 10 August 2008 Abstract The demand of reliable statistics for small areas, when only reduced sizes of the
More informationHB Methods for Combining Estimates from Multiple Surveys
Hierarchical Bayesian Methods for Combining Estimates from Multiple Surveys Adrijo Chakraborty NORC at the University of Chicago January 30, 2015 Joint work with Gauri Sankar Datta and Yang Cheng Outline
More informationMixed-Models. version 30 October 2011
Mixed-Models version 30 October 2011 Mixed models Mixed models estimate a vector! of fixed effects and one (or more) vectors u of random effects Both fixed and random effects models always include a vector
More informationMGR-815. Notes for the MGR-815 course. 12 June School of Superior Technology. Professor Zbigniew Dziong
Modeling, Estimation and Control, for Telecommunication Networks Notes for the MGR-815 course 12 June 2010 School of Superior Technology Professor Zbigniew Dziong 1 Table of Contents Preface 5 1. Example
More informationINDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -20 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.
INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY Lecture -20 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. Summary of the previous lecture Case study -3: Monthly streamflows
More informationParametric Bootstrap Methods for Bias Correction in Linear Mixed Models
CIRJE-F-801 Parametric Bootstrap Methods for Bias Correction in Linear Mixed Models Tatsuya Kubokawa University of Tokyo Bui Nagashima Graduate School of Economics, University of Tokyo April 2011 CIRJE
More informationSmall Area Confidence Bounds on Small Cell Proportions in Survey Populations
Small Area Confidence Bounds on Small Cell Proportions in Survey Populations Aaron Gilary, Jerry Maples, U.S. Census Bureau U.S. Census Bureau Eric V. Slud, U.S. Census Bureau Univ. Maryland College Park
More informationEmpirical Macroeconomics
Empirical Macroeconomics Francesco Franco Nova SBE April 21, 2015 Francesco Franco Empirical Macroeconomics 1/33 Growth and Fluctuations Supply and Demand Figure : US dynamics Francesco Franco Empirical
More informationJerry A. Hausinan. Working Paper No. 1133
NBER WORKING PAPER SERIES SEASONAL ADJUSTMENT WITH MEASUREMENT ERROR PRESENT Jerry A. Hausinan Mark Watson Working Paper No. 1133 NATIONAL BUREAU OF ECONOMIC RESEARCH 1050 Massachusetts Avenue Cambridge
More informationLinear Regression with Time Series Data
Econometrics 2 Linear Regression with Time Series Data Heino Bohn Nielsen 1of21 Outline (1) The linear regression model, identification and estimation. (2) Assumptions and results: (a) Consistency. (b)
More informationMixed-Model Estimation of genetic variances. Bruce Walsh lecture notes Uppsala EQG 2012 course version 28 Jan 2012
Mixed-Model Estimation of genetic variances Bruce Walsh lecture notes Uppsala EQG 01 course version 8 Jan 01 Estimation of Var(A) and Breeding Values in General Pedigrees The above designs (ANOVA, P-O
More informationSmall area solutions for the analysis of pollution data TIES Daniela Cocchi, Enrico Fabrizi, Carlo Trivisano. Università di Bologna
Small area solutions for the analysis of pollution data Daniela Cocchi, Enrico Fabrizi, Carlo Trivisano Università di Bologna TIES 00 Genova, June nd 00 Summary The small area problem in the environmental
More informationLeast Squares Estimation-Finite-Sample Properties
Least Squares Estimation-Finite-Sample Properties Ping Yu School of Economics and Finance The University of Hong Kong Ping Yu (HKU) Finite-Sample 1 / 29 Terminology and Assumptions 1 Terminology and Assumptions
More informationEstimation and Prediction Scenarios
Recursive BLUE BLUP and the Kalman filter: Estimation and Prediction Scenarios Amir Khodabandeh GNSS Research Centre, Curtin University of Technology, Perth, Australia IUGG 2011, Recursive 28 June BLUE-BLUP
More informationForecasting with R A practical workshop
Forecasting with R A practical workshop International Symposium on Forecasting 2016 19 th June 2016 Nikolaos Kourentzes nikolaos@kourentzes.com http://nikolaos.kourentzes.com Fotios Petropoulos fotpetr@gmail.com
More informationSome Time-Series Models
Some Time-Series Models Outline 1. Stochastic processes and their properties 2. Stationary processes 3. Some properties of the autocorrelation function 4. Some useful models Purely random processes, random
More informationEstimating σ 2. We can do simple prediction of Y and estimation of the mean of Y at any value of X.
Estimating σ 2 We can do simple prediction of Y and estimation of the mean of Y at any value of X. To perform inferences about our regression line, we must estimate σ 2, the variance of the error term.
More informationEmpirical Macroeconomics
Empirical Macroeconomics Francesco Franco Nova SBE April 5, 2016 Francesco Franco Empirical Macroeconomics 1/39 Growth and Fluctuations Supply and Demand Figure : US dynamics Francesco Franco Empirical
More informationStatistical signal processing
Statistical signal processing Short overview of the fundamentals Outline Random variables Random processes Stationarity Ergodicity Spectral analysis Random variable and processes Intuition: A random variable
More informationCovers Chapter 10-12, some of 16, some of 18 in Wooldridge. Regression Analysis with Time Series Data
Covers Chapter 10-12, some of 16, some of 18 in Wooldridge Regression Analysis with Time Series Data Obviously time series data different from cross section in terms of source of variation in x and y temporal
More informationAUTOCORRELATION. Phung Thanh Binh
AUTOCORRELATION Phung Thanh Binh OUTLINE Time series Gauss-Markov conditions The nature of autocorrelation Causes of autocorrelation Consequences of autocorrelation Detecting autocorrelation Remedial measures
More informationComments on Design-Based Prediction Using Auxilliary Information under Random Permutation Models (by Wenjun Li (5/21/03) Ed Stanek
Comments on Design-Based Prediction Using Auxilliary Information under Random Permutation Models (by Wenjun Li (5/2/03) Ed Stanek Here are comments on the Draft Manuscript. They are all suggestions that
More informationEconometrics of Panel Data
Econometrics of Panel Data Jakub Mućk Meeting # 1 Jakub Mućk Econometrics of Panel Data Meeting # 1 1 / 31 Outline 1 Course outline 2 Panel data Advantages of Panel Data Limitations of Panel Data 3 Pooled
More informationEconomics 536 Lecture 7. Introduction to Specification Testing in Dynamic Econometric Models
University of Illinois Fall 2016 Department of Economics Roger Koenker Economics 536 Lecture 7 Introduction to Specification Testing in Dynamic Econometric Models In this lecture I want to briefly describe
More informationChapter 9: Forecasting
Chapter 9: Forecasting One of the critical goals of time series analysis is to forecast (predict) the values of the time series at times in the future. When forecasting, we ideally should evaluate the
More informationThe Statistical Property of Ordinary Least Squares
The Statistical Property of Ordinary Least Squares The linear equation, on which we apply the OLS is y t = X t β + u t Then, as we have derived, the OLS estimator is ˆβ = [ X T X] 1 X T y Then, substituting
More informationQuestions and Answers on Heteroskedasticity, Autocorrelation and Generalized Least Squares
Questions and Answers on Heteroskedasticity, Autocorrelation and Generalized Least Squares L Magee Fall, 2008 1 Consider a regression model y = Xβ +ɛ, where it is assumed that E(ɛ X) = 0 and E(ɛɛ X) =
More informationEconometrics of Panel Data
Econometrics of Panel Data Jakub Mućk Meeting # 4 Jakub Mućk Econometrics of Panel Data Meeting # 4 1 / 30 Outline 1 Two-way Error Component Model Fixed effects model Random effects model 2 Non-spherical
More informationBayesian variable selection via. Penalized credible regions. Brian Reich, NCSU. Joint work with. Howard Bondell and Ander Wilson
Bayesian variable selection via penalized credible regions Brian Reich, NC State Joint work with Howard Bondell and Ander Wilson Brian Reich, NCSU Penalized credible regions 1 Motivation big p, small n
More informationWLS and BLUE (prelude to BLUP) Prediction
WLS and BLUE (prelude to BLUP) Prediction Rasmus Waagepetersen Department of Mathematics Aalborg University Denmark April 21, 2018 Suppose that Y has mean X β and known covariance matrix V (but Y need
More informationTime Series and Forecasting
Time Series and Forecasting Introduction to Forecasting n What is forecasting? n Primary Function is to Predict the Future using (time series related or other) data we have in hand n Why are we interested?
More informationAnalysis. Components of a Time Series
Module 8: Time Series Analysis 8.2 Components of a Time Series, Detection of Change Points and Trends, Time Series Models Components of a Time Series There can be several things happening simultaneously
More informationStat 579: Generalized Linear Models and Extensions
Stat 579: Generalized Linear Models and Extensions Linear Mixed Models for Longitudinal Data Yan Lu April, 2018, week 15 1 / 38 Data structure t1 t2 tn i 1st subject y 11 y 12 y 1n1 Experimental 2nd subject
More information{ } Stochastic processes. Models for time series. Specification of a process. Specification of a process. , X t3. ,...X tn }
Stochastic processes Time series are an example of a stochastic or random process Models for time series A stochastic process is 'a statistical phenomenon that evolves in time according to probabilistic
More informationA Review of the Problem of. Seasonal Adjustment Variances
A Review of the Problem of Seasonal Adjustment Variances William R. Bell U.S. Census Bureau William.R.Bell@ 1 Disclaimer: The views expressed here are those of the author and not necessarily those oftheu.s.censusbureau.
More informationLinear Model Under General Variance Structure: Autocorrelation
Linear Model Under General Variance Structure: Autocorrelation A Definition of Autocorrelation In this section, we consider another special case of the model Y = X β + e, or y t = x t β + e t, t = 1,..,.
More informationForecasting Unemployment Rates in the UK and EU
Forecasting Unemployment Rates in the UK and EU Team B9 Rajesh Manivannan (61710602) Kartick B Muthiah (61710764) Debayan Das (61710492) Devesh Kumar (61710353) Chirag Bhardwaj (61710812) Sreeharsha Konga
More informationmultilevel modeling: concepts, applications and interpretations
multilevel modeling: concepts, applications and interpretations lynne c. messer 27 october 2010 warning social and reproductive / perinatal epidemiologist concepts why context matters multilevel models
More informationBusiness Economics BUSINESS ECONOMICS. PAPER No. : 8, FUNDAMENTALS OF ECONOMETRICS MODULE No. : 3, GAUSS MARKOV THEOREM
Subject Business Economics Paper No and Title Module No and Title Module Tag 8, Fundamentals of Econometrics 3, The gauss Markov theorem BSE_P8_M3 1 TABLE OF CONTENTS 1. INTRODUCTION 2. ASSUMPTIONS OF
More informationAccounting for Hyperparameter Uncertainty in SAE Based on a State-Space Model: the Dutch Labour Force Survey Oksana Bollineni-Balabay, Jan van den
Accounting for Hyperparameter Uncertainty in SAE Based on a State-Space Model: the Dutch Labour Force Survey Oksana Bollineni-Balabay, Jan van den Brakel, Franz Palm The Dutch LFS monthly estimates for
More informationProblem Set #6: OLS. Economics 835: Econometrics. Fall 2012
Problem Set #6: OLS Economics 835: Econometrics Fall 202 A preliminary result Suppose we have a random sample of size n on the scalar random variables (x, y) with finite means, variances, and covariance.
More informationVector spaces. DS-GA 1013 / MATH-GA 2824 Optimization-based Data Analysis.
Vector spaces DS-GA 1013 / MATH-GA 2824 Optimization-based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_fall17/index.html Carlos Fernandez-Granda Vector space Consists of: A set V A scalar
More informationEstimation of Complex Small Area Parameters with Application to Poverty Indicators
1 Estimation of Complex Small Area Parameters with Application to Poverty Indicators J.N.K. Rao School of Mathematics and Statistics, Carleton University (Joint work with Isabel Molina from Universidad
More informationLab: Box-Jenkins Methodology - US Wholesale Price Indicator
Lab: Box-Jenkins Methodology - US Wholesale Price Indicator In this lab we explore the Box-Jenkins methodology by applying it to a time-series data set comprising quarterly observations of the US Wholesale
More information11. Further Issues in Using OLS with TS Data
11. Further Issues in Using OLS with TS Data With TS, including lags of the dependent variable often allow us to fit much better the variation in y Exact distribution theory is rarely available in TS applications,
More informationSection 2 NABE ASTEF 65
Section 2 NABE ASTEF 65 Econometric (Structural) Models 66 67 The Multiple Regression Model 68 69 Assumptions 70 Components of Model Endogenous variables -- Dependent variables, values of which are determined
More informationENSC327 Communications Systems 19: Random Processes. Jie Liang School of Engineering Science Simon Fraser University
ENSC327 Communications Systems 19: Random Processes Jie Liang School of Engineering Science Simon Fraser University 1 Outline Random processes Stationary random processes Autocorrelation of random processes
More informationNon-Stationary Time Series and Unit Root Testing
Econometrics II Non-Stationary Time Series and Unit Root Testing Morten Nyboe Tabor Course Outline: Non-Stationary Time Series and Unit Root Testing 1 Stationarity and Deviation from Stationarity Trend-Stationarity
More informationLinear Regression with Time Series Data
u n i v e r s i t y o f c o p e n h a g e n d e p a r t m e n t o f e c o n o m i c s Econometrics II Linear Regression with Time Series Data Morten Nyboe Tabor u n i v e r s i t y o f c o p e n h a g
More informationKalman Filter. Predict: Update: x k k 1 = F k x k 1 k 1 + B k u k P k k 1 = F k P k 1 k 1 F T k + Q
Kalman Filter Kalman Filter Predict: x k k 1 = F k x k 1 k 1 + B k u k P k k 1 = F k P k 1 k 1 F T k + Q Update: K = P k k 1 Hk T (H k P k k 1 Hk T + R) 1 x k k = x k k 1 + K(z k H k x k k 1 ) P k k =(I
More informationMultivariate Time Series
Multivariate Time Series Notation: I do not use boldface (or anything else) to distinguish vectors from scalars. Tsay (and many other writers) do. I denote a multivariate stochastic process in the form
More informationFORECASTING FLUCTUATIONS OF ASPHALT CEMENT PRICE INDEX IN GEORGIA
FORECASTING FLUCTUATIONS OF ASPHALT CEMENT PRICE INDEX IN GEORGIA Mohammad Ilbeigi, Baabak Ashuri, Ph.D., and Yang Hui Economics of the Sustainable Built Environment (ESBE) Lab, School of Building Construction
More information