The Hodrick-Prescott Filter
|
|
- Lily Pope
- 6 years ago
- Views:
Transcription
1 The Hodrick-Prescott Filter A Special Case of Penalized Spline Smoothing Alex Trindade Dept. of Mathematics & Statistics, Texas Tech University Joint work with Rob Paige, Missouri University of Science and Technology Funded in part by the National Security Agency Grants: H (Paige) & H (Trindade) June 2012 SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University JuneJoint 2012work with 1 / 38 Rob
2 Outline 1 Motivation Examples: the need for a semiparametric smoothing approach... 2 Semiparametric Models Partially linear models & penalized spline models Penalized splines as linear mixed models (LMMs) 3 Result #2: Estimation of Smoothing Parameter Estimators as roots of quadratic estimating equations (QEEs) Saddlepoint-based bootstrap (SPBB) inference for QEEs Exact ML & REML Inference in LMMs 4 Result #1: The Hodrick-Prescott Filter (HPF) HPF solution coincides with linear penalized spline 5 Examples and Simulations U.S. Gross National Product (GNP) Returns of GNP SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University JuneJoint 2012work with 2 / 38 Rob
3 The LIDAR Data (Ruppert, Wand, & Carroll, 2003) Model: y = µ(x) + error. Goal: estimate mean function µ(x), i.e. smooth data. log ratio range SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University JuneJoint 2012work with 3 / 38 Rob
4 Standard Parametric Approaches Fail Polynomial regression: µ(x) = β 0 + β 1 x + + β p x p p = 3, 4: don t capture sharp drop (x 600) p = 10: µ(x) too wiggly; µ (x) too variable (and changes sign!) Nonlinear regression, e.g. µ(x) = β 0 e β 1x, requires thinking about correct (nonlinear) relationship... Can often transform x and/or y such that get linear relationship, e.g. take logs above: log µ(x) = log(β 0 ) + β 1 x No transform here seems to help in linearizing (and stabilizing variance of) relationship. SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University JuneJoint 2012work with 4 / 38 Rob
5 The Fossil Data (Ruppert, Wand, & Carroll, 2003) strontium ratio dip here? obvious dip! age (millions of years) SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University JuneJoint 2012work with 5 / 38 Rob
6 U.S.A. Detrended Log(GNP) Quarterly Data (Hodrick & Prescott, 1997) A cycle of long duration? Time SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University JuneJoint 2012work with 6 / 38 Rob
7 A Simple Spline Model (LIDAR Data) Linear spline model with knots at κ 1 = 575 and κ 2 = 600 is promising µ(x) = β 0 + β 1 x + u 1 (x 575) + + u 2 (x 600) + Piecewise linear function, called spline (thin strip of flexible wood). A spline model is a flexible model. log ratio range SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University JuneJoint 2012work with 7 / 38 Rob
8 Polynomial Spline Models In general, have polynomial spline model of degree p with K knots µ(x) = β 0 + β 1 x + + β p x p + K k=1 u k (x κ k ) p + Given pairs of obs (x i, y i ), i = 1,..., n, can write mean in matrix form µ = X β + Zu Bθ Can also add vector of covariates w entering mean linearly µ = W γ }{{} + }{{} Bθ linear part nonlinear part Model can allow for autocorrelation, R, in residuals (e.g. time series): y = µ + ε, ε (0, σ 2 ε R) SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University JuneJoint 2012work with 8 / 38 Rob
9 Penalized Spline Models Polynomial spline model Estimate θ by minimizing µ = X β + Zu Bθ ˆθ PS = arg min θ { (y Bθ) R 1 (y Bθ) + αθ Dθ } D: a non-negative definite penalty matrix (usually θ Dθ = u u) α is a smoothing parameter controlling balance between: LIDAR: linear spline fits with max and min smoothing (24 knots) fidelity to data (α = 0) smoothness of fit (α = ) log ratio α = 0 α = range SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University JuneJoint 2012work with 9 / 38 Rob
10 Linear Mixed Model (LMM) Formulation With assumptions: [ ] u E = ε [ ] 0, Cov 0 [ ] [ ] u σ 2 = u I K 0 ε 0 σε 2 R Penalized spline can be recast as LMM with one variance component (Brumback, Ruppert, & Wand, 1999) y = X β }{{} + Zu }{{ + ε } fixed effects random effects BLUP of y in this context is ỹ = B θ, where θ = arg min θ { (y Bθ) R 1 (y Bθ) + σ2 ε σu 2 u u }. SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work10 with / 38 Rob
11 Smoothing Parameter is BLUP in LMM Formulation Define: α = σε 2 /σu 2 Let D be zero with I K as lower diagonal block: [ ] 0 0 D = J 0 I p+1,k = θ Dθ = u u K We obtain that: θ = ˆθ PS Fitted penalized spline for mean of y coincides with BLUP of y in LMM defined above. Implies BLUP-optimal value for α is: α = σ 2 ε /σ 2 u SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work11 with / 38 Rob
12 Estimation of Smoothing Parameter Since α is ratio of variance components in LMM many methods available in this context; typically require normal errors. Call these parametric. (SAS: proc mixed; R: lme.) Methods that do not require LMM representation and specification of error distribution: nonparametric. Examples (Parametric) Maximum Likelihood (ML) REstricted Maximum Likelihood (REML) Examples (Nonparametric) Akaike s Information Criterion (AIC) Generalized Cross-Validation (GCV) SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work12 with / 38 Rob
13 A Unified View of Smoothing Parameter Estimators (New) Above estimators can be viewed as roots of a quadratic estimating equation (QEE) in normal random variables Q(α) = y A α y The n n matrix A α has a (complicated, but) closed-form expression in each case... Example (A α for GCV estimator) Define: S α = B ( B R 1 B + αd ) 1 B R 1, and Ṡ α = S α / α. Then: A α = (I n S α )Ṡ α {1 Tr(S α )/n} (I n S α ) 2 Tr(S α )/n SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work13 with / 38 Rob
14 Example of QEE Based Inference: REML Example (A α for REML estimator) Assume: y u N(X β + Zu, σε 2 R), u N(0, σui 2 K ) Differentiation of restricted log-likelihood in β and σε 2 leads to profile REML(α, R) criterion, and A α =... too long to fit here... Krivobokova & Kauermann (2007): REML less sensitive to misspecification of residual correlation than AIC or GCV. Inspect ACF/PACF plots for ideas on what R should be... SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work14 with / 38 Rob
15 Saddlepoint-Based Bootstrap (SPBB) Inference for QEEs Technique pioneered by Paige, Trindade, & Fernando (2009): Assume y N(µ, Σ); want confidence interval (CI) for scalar α from Q(α). Estimator ˆα is (appropriate) root of Q(α). Q(α) usually obtained by differentiation of (profile) log likelihood (nuisance parameters replaced by conditional MLEs). Can express MGF of Q(α) in closed-form. Use MGF to saddlepoint approximate CDF of Q(α). Under monotonicity of Q(α) in α, relate saddlepoint approximations: CDF of Q(α) to CDF of ˆα. CI (α L, α U ) obtained by pivoting CDF of ˆα (inverting 2-sided test). SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work15 with / 38 Rob
16 SPBB Details: Slide 1 Approx (1 ζ)100% CI obtained by solving: Fˆα (ˆαobs α L, ˆθ(α L ), ˆσ 2 ε (α L ) ) = 1 ζ/2 Fˆα (ˆαobs α U, ˆθ(α U ), ˆσ 2 ε (α U ) ) = ζ/2 2nd order accurate: coverage prob is (1 ζ) + O(n 1 ) (DiCiccio & Romano, 1995). Fˆα ( ) is difficult to compute..., but if Q(α) is monotone (e.g. decreasing) in α: (Daniels, 1983). Fˆα (α) = P(ˆα α) = P(Q(α) 0) = F Q (0) F Q(α) ( ) also difficult to compute..., but have closed-form expression for its MGF (linear combo of indep noncentral χ 2 r.v. s). SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work16 with / 38 Rob
17 SPBB Details: Slide 2 Saddlepoint approximate CDF of Q(α): F Q(α) ( ) ˆF Q(α) ( ) (Lugannani & Rice, 1980). Saddlepoint approximations extremely accurate (3rd order accurate over sets of bounded central tendency), e.g. Butler (2007). (1 ζ)100% CI obtained by solving (numerically over a grid): ( ˆF Q(ˆαobs ) 0 αl, ˆθ(α L ), ˆσ ε 2 (α L ) ) = 1 ζ/2 ( ˆF Q(ˆαobs ) 0 αu, ˆθ(α U ), ˆσ ε 2 (α U ) ) = ζ/2 Still retain 2nd order accuracy in coverage prob (Paige, Trindade, & Fernando, 2009). SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work17 with / 38 Rob
18 SPBB Summary Method is similar to a parametric bootstrap: Use plug-in estimator ˆα; Usually sample from bootstrap distribution Fˆα ( ) via Monte Carlo simulation; BUT, replace simulation with saddlepoint approximation of Fˆα ( ). ( Bootstrap is in the name, but there s no re-sampling...) What is a saddlepoint approximation? Accurate method to approx CDF & PDF of X from its MGF, M X ( ). Alternative to Fourier inversion of characteristic function (FFT). Most computationally expensive step: for each x X, solve for ŝ log M X (s) s = x s=ŝ Original paper: Daniels (1954). Recent book: Butler (2007). SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work18 with / 38 Rob
19 Graphical Overview of SPBB (Result #2) Intractable! (Except via parametric bootstrap... slow.) Fˆα (ˆα obs ) (α L, α U ) ˆα solves pivot Q(α) = 0 ˆF Q(ˆαobs )(0) Q(α) monotone Fˆα (ˆα obs ) = F Q(ˆαobs )(0) saddlepoint approx via MGF of Q(α) SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work19 with / 38 Rob
20 Exact ML & REML Inference for α Exact finite sample inference for α = σ 2 ε /σ 2 u in LMMs with one variance component (Crainiceanu, Ruppert, Claeskens, & Wand, 2005): Note: asymptotic χ 2 dist is poor approx in finite samples due to substantial point mass at 0 (Crainiceanu & Ruppert, 2004). Assume y N(µ, Σ). Invert (restricted) likelihood ratio test corresponding to (REML-based) ML-based inference. Grid search needed to locate endpoints of CI (α L, α U ): At each grid value of α, use fast algorithm to draw large sample from H 0 distribution of RLRT/LRT statistic in order to (empirically) calculate p-value(α) p-value(α L ) = 0.95 and p-value(α U ) = 0.95 SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work20 with / 38 Rob
21 Comparison of SPBB and Exact Method Exact y normal. α: ratio of vars in LMM. Dataset-specific (needs manual tuning). LIDAR: 4 hours/ci. Performance: exact. Estimators: ML & REML only. SPBB y normal. α: root of QEE. QEE-specific (needs derivatives of QEE). LIDAR: 10 mins/ci. Performance: near-exact. Estimators: many (ML, REML, AIC, GCV, etc.) Fact SPBB is able to furnish a CI when no other computationally feasible method, exact or otherwise, exists, e.g. AIC & GCV. SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work21 with / 38 Rob
22 The Hodrick-Prescott Filter (Hodrick & Prescott, 1997) Smoothing technique for series y t : viewed as sum of trend µ t and residual component ε t y t = µ t + ε t, t = 1,..., n For choice of smoothing parameter α, trend fitted by minimizing weighted sum of squares ŷ HP = arg min µ R n n t=1 balance (y t µ t ) 2 {}}{ + α }{{} fidelity to data n t=1 (µ t 2µ t 1 + µ t 2 ) 2 }{{} smoothness of fit = (I n + α 2 2 ) 1 y, 2 = function({0, 1, 2}) SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work22 with / 38 Rob
23 HPF is Linear Penalized Spline (Result #1) Theorem Consider penalized spline model of degree p = 1, explanatory variable time x t = t, K = n 2 equi-spaced knots at 2,..., n 1, and uncorrelated residuals R = I n : n 1 y t = β 0 + β 1 t + u k (t k) + + ε t, t = 1,..., n. k=2 Then HPF solution coincides with penalized spline solution: ŷ PS = B PS (B PS B PS + αj 2,n 2 ) 1 B PS y = (I n + α 2 2 ) 1 y = ŷ HP SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work23 with / 38 Rob
24 HPF Coincides With Linear Penalized Spline: Proof Proof. HPF can be formulated as LMM (Schlicht, 2005): y = X HP β + Z HP u + ε B HP θ + ε where Z HP = 2 ( 2 2 ) 1, and X HP is any n 2 matrix satisfying: 2 X HP = 0, and X HP X HP = I 2 LMM formulation implies HPF solution is BLUP of y: ŷ HP = B HP (B HP B HP + αj 2,n 2 ) 1 B HP y SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work24 with / 38 Rob
25 Proof (continued...) Defining L = B 1 PS B HP, transform basis functions (but preserve degree and knot locations): ŷ HP = B HP [B HP B HP + αj 2,n 2 ] 1 B HP [ ( y = B PS L L B PS B PS + αl T J 2,n 2 L 1) 1 L] L B PS y = B PS [B PS B PS + αl T J 2,n 2 L 1] 1 B PS y =... = B PS (B PS B PS + αj 2,n 2 ) 1 B PS y = ŷ PS Equivalence follows if we can show in step above that L T J 2,n 2 L 1 = J 2,n 2 Decomposing L 1 into block structure, result follows...
26 Remarks Schlicht (2005): calls LMM representation of HPF a stochastic model, unaware that it is in fact a LMM; proposes maximum (Gaussian) likelihood estimator (MLE) for α; proposes also a method of moments estimate. Successive econometrics papers have built on this; propose new estimators e.g. Dermoune, Djehiche, and Rahmania (2008). But, there is already a vast existing body of work on the estimation of fixed and random components in LMMs... e.g. Christensen (1996). Fact LMM formulation of penalized spline endows it, and therefore HPF, with rich variety of existing estimation methods, e.g. REML (a parametric choice) and GCV (a nonparametric choice). SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work26 with / 38 Rob
27 Example 1: U.S. Gross National Product (GNP) Quarterly (seasonally adjusted) GNP for period to (billions of dollars). Have n = 189 pairs of (t, y t ) values: y t = log GNP, t = time (quarter since ). Strong linearly increasing trend, remove it! Hodrick & Prescott (1997): smooth via HPF with choice α = 1, 600 advocated for quarterly data (the HPF smooth). Theorem: HPF equivalent to 187-knot (K = 187) linear (p = 1) penalized spline model with uncorrelated residuals (R = I n )... Alternatives to HPF α = 1, 600 choice: REML and GCV. SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work27 with / 38 Rob
28 Example 2: Log Returns of GNP Analyze also first differences, {y t y t 1 } of log GNP (log returns), corresponds to a growth rate. Have n = 188 pairs of (t, y t ) values. Again compare smooths: HPF, REML, GCV. Theorem: HPF equivalent to 186-knot (K = 186) linear (p = 1) penalized spline modell with uncorrelated residuals (R = I n ). SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work28 with / 38 Rob
29 Log GNP and Log Returns of GNP Data Detrended Log GNP Log Returns of GNP Time Time ACF of Detrended Log GNP ACF of Log Returns of GNP Lag Lag SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work29 with / 38 Rob
30 Results: Fitted Models Log GNP: ACF/PACF of REML-IID & GCV-IID (R = I n ) residuals suggests AR(3) autocorrelation structure... so fit REML-AR(3). Log Returns GNP: REML-IID & GCV-IID (R = I n ) residuals white. Dataset Method α 95% Estimate (lo, pt, hi) Log GNP Log Returns GNP HPF (?, 1600,?) Exact-REML-IID (0.07, 0.16, 0.32) SPBB-GCV-IID (0.04, 0.15, 0.59) SPBB-REML-AR(3) (97.6, 337.3, 871.4) HPF (?, 1600,?) Exact-REML-IID (8, 595, 65, 302, 330, 436) SPBB-REML-IID (5, 572, 65, 302, 463, 185) SPBB-GCV-IID (6, 022, 79, 164, ) SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work30 with / 38 Rob
31 Smoothed Log GNP: REML & GCV with R = I n underfit... Detrended Log GNP HPF REML IID & GCV IID Quarter since 1947 SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work31 with / 38 Rob
32 Smoothed Log GNP: REML-AR(3) coincides with HPF! Detrended Log GNP HPF REML AR(3) Quarter since 1947 SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work32 with / 38 Rob
33 Smoothed Returns of Log GNP: HPF underfits? Log Returns of GNP REML IID & GCV IID HPF Quarter since 1947 SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work33 with / 38 Rob
34 Simulations: Performance of Exact vs. SPBB CI s for α 1,000 reps from REML-fitted model to Log Returns of GNP. Performance of SPBB is near-exact! SPBB: only method to get CI in case of GCV or REML (non-iid). Estimation Median Underage Coverage Overage Method CI Length Prob. Prob. Prob. Exact-REML 374, SPBB-REML 466, SPBB-GCV SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work34 with / 38 Rob
35 Summary Avoid blind usage of HPF: equivalent to linear penalized spline with equi-spaced knots and uncorrelated residuals. Advocate structured penalized spline approach: pay attention to degree of polynomial (p), number of knots (K), residual autocorrelaton (R). Use data-driven method for estimation: REML is less sensitive to misspecification of R than GCV or AIC. Use SPBB for α CI s: if desired; gives range of plausible smooths. SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work35 with / 38 Rob
36 Application: The Fossil Data strontium ratio Point and 95% Interval Smoothing Estimates REML & GCV point estimate (solid) Exact REML lower bound (dash) Exact REML upper bound (dot) SPBB GCV lower bound (dot dash) SPBB GCV upper bound (long dash) age (millions of years) SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work36 with / 38 Rob
37 Acknowledgements Reviewers: valuable suggestions. Ciprian Crainiceanu: exact ML & REML code. High Performance Computing Center, Texas Tech University: supercomputing. National Security Agency: not having to teach in summer. SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work37 with / 38 Rob
38 Key References Crainiceanu, C., Ruppert, D., Claeskens, G., and Wand, M. (2005), Exact likelihood ratio tests for penalized splines, Biometrika, 92, Hodrick, R.J., and Prescott, E.C. (1997), Postwar U.S. business cycles: an empirical investigation, J. Money, Credit, and Banking, 29, Krivobokova, T., and Kauermann, G. (2007), A Note on Penalized Spline Smoothing with Correlated Errors, J. Amer. Statist. Assoc., 102, Paige, R.L., Trindade, A.A. and Fernando, P.H. (2009), Saddlepoint-based bootstrap inference for quadratic estimating equations, Scand. J. Stat., 36, Paige, R.L. and Trindade, A.A. (2010), The Hodrick-Prescott Filter: A Special Case of Penalized Spline Smoothing, Electronic J. Statist., 4, Paige, R.L., and Trindade, A.A. (2012), Fast and Accurate Inference for the Smoothing Parameter in Semiparametric Models, Austral. N. Zeal. J. Statist., (to appear). Ruppert, D., Wand, M.P. and Carroll, R.J. (2003), Semiparametric Regression, London: Cambridge. Schlicht, E. (2005), Estimating the smoothing parameter in the so-called Hodrick-Prescott filter, J. Japan Statist. Soc., 35, SPBB (Dept. Inference of Mathematics for HPF & Statistics, Texas Tech University June 2012 Joint work38 with / 38 Rob
Saddlepoint-Based Bootstrap Inference in Dependent Data Settings
Saddlepoint-Based Bootstrap Inference in Dependent Data Settings Alex Trindade Dept. of Mathematics & Statistics, Texas Tech University Rob Paige, Missouri University of Science and Technology Indika Wickramasinghe,
More informationLikelihood Ratio Tests. that Certain Variance Components Are Zero. Ciprian M. Crainiceanu. Department of Statistical Science
1 Likelihood Ratio Tests that Certain Variance Components Are Zero Ciprian M. Crainiceanu Department of Statistical Science www.people.cornell.edu/pages/cmc59 Work done jointly with David Ruppert, School
More informationExact Likelihood Ratio Tests for Penalized Splines
Exact Likelihood Ratio Tests for Penalized Splines By CIPRIAN CRAINICEANU, DAVID RUPPERT, GERDA CLAESKENS, M.P. WAND Department of Biostatistics, Johns Hopkins University, 615 N. Wolfe Street, Baltimore,
More informationAndreas Blöchl: Trend Estimation with Penalized Splines as Mixed Models for Series with Structural Breaks
Andreas Blöchl: Trend Estimation with Penalized Splines as Mixed Models for Series with Structural Breaks Munich Discussion Paper No. 2014-2 Department of Economics University of Munich Volkswirtschaftliche
More informationNonparametric Small Area Estimation Using Penalized Spline Regression
Nonparametric Small Area Estimation Using Penalized Spline Regression 0verview Spline-based nonparametric regression Nonparametric small area estimation Prediction mean squared error Bootstrapping small
More informationMath 494: Mathematical Statistics
Math 494: Mathematical Statistics Instructor: Jimin Ding jmding@wustl.edu Department of Mathematics Washington University in St. Louis Class materials are available on course website (www.math.wustl.edu/
More informationRestricted Likelihood Ratio Tests in Nonparametric Longitudinal Models
Restricted Likelihood Ratio Tests in Nonparametric Longitudinal Models Short title: Restricted LR Tests in Longitudinal Models Ciprian M. Crainiceanu David Ruppert May 5, 2004 Abstract We assume that repeated
More informationSome properties of Likelihood Ratio Tests in Linear Mixed Models
Some properties of Likelihood Ratio Tests in Linear Mixed Models Ciprian M. Crainiceanu David Ruppert Timothy J. Vogelsang September 19, 2003 Abstract We calculate the finite sample probability mass-at-zero
More informationPenalized Splines, Mixed Models, and Recent Large-Sample Results
Penalized Splines, Mixed Models, and Recent Large-Sample Results David Ruppert Operations Research & Information Engineering, Cornell University Feb 4, 2011 Collaborators Matt Wand, University of Wollongong
More informationOn the Behavior of Marginal and Conditional Akaike Information Criteria in Linear Mixed Models
On the Behavior of Marginal and Conditional Akaike Information Criteria in Linear Mixed Models Thomas Kneib Institute of Statistics and Econometrics Georg-August-University Göttingen Department of Statistics
More informationFall 2017 STAT 532 Homework Peter Hoff. 1. Let P be a probability measure on a collection of sets A.
1. Let P be a probability measure on a collection of sets A. (a) For each n N, let H n be a set in A such that H n H n+1. Show that P (H n ) monotonically converges to P ( k=1 H k) as n. (b) For each n
More informationStatistics - Lecture One. Outline. Charlotte Wickham 1. Basic ideas about estimation
Statistics - Lecture One Charlotte Wickham wickham@stat.berkeley.edu http://www.stat.berkeley.edu/~wickham/ Outline 1. Basic ideas about estimation 2. Method of Moments 3. Maximum Likelihood 4. Confidence
More informationModel Selection and Geometry
Model Selection and Geometry Pascal Massart Université Paris-Sud, Orsay Leipzig, February Purpose of the talk! Concentration of measure plays a fundamental role in the theory of model selection! Model
More informationBootstrap. Director of Center for Astrostatistics. G. Jogesh Babu. Penn State University babu.
Bootstrap G. Jogesh Babu Penn State University http://www.stat.psu.edu/ babu Director of Center for Astrostatistics http://astrostatistics.psu.edu Outline 1 Motivation 2 Simple statistical problem 3 Resampling
More informationBetter Bootstrap Confidence Intervals
by Bradley Efron University of Washington, Department of Statistics April 12, 2012 An example Suppose we wish to make inference on some parameter θ T (F ) (e.g. θ = E F X ), based on data We might suppose
More informationHigh Dimensional Empirical Likelihood for Generalized Estimating Equations with Dependent Data
High Dimensional Empirical Likelihood for Generalized Estimating Equations with Dependent Data Song Xi CHEN Guanghua School of Management and Center for Statistical Science, Peking University Department
More informationNow consider the case where E(Y) = µ = Xβ and V (Y) = σ 2 G, where G is diagonal, but unknown.
Weighting We have seen that if E(Y) = Xβ and V (Y) = σ 2 G, where G is known, the model can be rewritten as a linear model. This is known as generalized least squares or, if G is diagonal, with trace(g)
More informationIntegrated Likelihood Estimation in Semiparametric Regression Models. Thomas A. Severini Department of Statistics Northwestern University
Integrated Likelihood Estimation in Semiparametric Regression Models Thomas A. Severini Department of Statistics Northwestern University Joint work with Heping He, University of York Introduction Let Y
More informationIndependent and conditionally independent counterfactual distributions
Independent and conditionally independent counterfactual distributions Marcin Wolski European Investment Bank M.Wolski@eib.org Society for Nonlinear Dynamics and Econometrics Tokyo March 19, 2018 Views
More informationTwo Applications of Nonparametric Regression in Survey Estimation
Two Applications of Nonparametric Regression in Survey Estimation 1/56 Jean Opsomer Iowa State University Joint work with Jay Breidt, Colorado State University Gerda Claeskens, Université Catholique de
More informationBayesian estimation of the discrepancy with misspecified parametric models
Bayesian estimation of the discrepancy with misspecified parametric models Pierpaolo De Blasi University of Torino & Collegio Carlo Alberto Bayesian Nonparametrics workshop ICERM, 17-21 September 2012
More informationSome Theories about Backfitting Algorithm for Varying Coefficient Partially Linear Model
Some Theories about Backfitting Algorithm for Varying Coefficient Partially Linear Model 1. Introduction Varying-coefficient partially linear model (Zhang, Lee, and Song, 2002; Xia, Zhang, and Tong, 2004;
More informationTREND ESTIMATION AND THE HODRICK-PRESCOTT FILTER
J. Japan Statist. Soc. Vol. 38 No. 1 2008 41 49 TREND ESTIMATION AND THE HODRICK-PRESCOTT FILTER Andrew Harvey* and Thomas Trimbur** The article analyses the relationship between unobserved component trend-cycle
More information9. Model Selection. statistical models. overview of model selection. information criteria. goodness-of-fit measures
FE661 - Statistical Methods for Financial Engineering 9. Model Selection Jitkomut Songsiri statistical models overview of model selection information criteria goodness-of-fit measures 9-1 Statistical models
More informationINFERENCE APPROACHES FOR INSTRUMENTAL VARIABLE QUANTILE REGRESSION. 1. Introduction
INFERENCE APPROACHES FOR INSTRUMENTAL VARIABLE QUANTILE REGRESSION VICTOR CHERNOZHUKOV CHRISTIAN HANSEN MICHAEL JANSSON Abstract. We consider asymptotic and finite-sample confidence bounds in instrumental
More informationTAKEHOME FINAL EXAM e iω e 2iω e iω e 2iω
ECO 513 Spring 2015 TAKEHOME FINAL EXAM (1) Suppose the univariate stochastic process y is ARMA(2,2) of the following form: y t = 1.6974y t 1.9604y t 2 + ε t 1.6628ε t 1 +.9216ε t 2, (1) where ε is i.i.d.
More informationSimultaneous Confidence Bands for the Coefficient Function in Functional Regression
University of Haifa From the SelectedWorks of Philip T. Reiss August 7, 2008 Simultaneous Confidence Bands for the Coefficient Function in Functional Regression Philip T. Reiss, New York University Available
More informationinterval forecasting
Interval Forecasting Based on Chapter 7 of the Time Series Forecasting by Chatfield Econometric Forecasting, January 2008 Outline 1 2 3 4 5 Terminology Interval Forecasts Density Forecast Fan Chart Most
More informationSpatially Adaptive Smoothing Splines
Spatially Adaptive Smoothing Splines Paul Speckman University of Missouri-Columbia speckman@statmissouriedu September 11, 23 Banff 9/7/3 Ordinary Simple Spline Smoothing Observe y i = f(t i ) + ε i, =
More informationHypothesis Testing in Smoothing Spline Models
Hypothesis Testing in Smoothing Spline Models Anna Liu and Yuedong Wang October 10, 2002 Abstract This article provides a unified and comparative review of some existing test methods for the hypothesis
More informationLikelihood ratio testing for zero variance components in linear mixed models
Likelihood ratio testing for zero variance components in linear mixed models Sonja Greven 1,3, Ciprian Crainiceanu 2, Annette Peters 3 and Helmut Küchenhoff 1 1 Department of Statistics, LMU Munich University,
More informationDiscussion of the paper Inference for Semiparametric Models: Some Questions and an Answer by Bickel and Kwon
Discussion of the paper Inference for Semiparametric Models: Some Questions and an Answer by Bickel and Kwon Jianqing Fan Department of Statistics Chinese University of Hong Kong AND Department of Statistics
More informationTesting Restrictions and Comparing Models
Econ. 513, Time Series Econometrics Fall 00 Chris Sims Testing Restrictions and Comparing Models 1. THE PROBLEM We consider here the problem of comparing two parametric models for the data X, defined by
More informationRIDGE REGRESSION REPRESENTATIONS OF THE GENERALIZED HODRICK-PRESCOTT FILTER
J. Japan Statist. Soc. Vol. 45 No. 2 2015 121 128 RIDGE REGRESSION REPRESENTATIONS OF THE GENERALIZED HODRICK-PRESCOTT FILTER Hiroshi Yamada* The Hodrick-Prescott (HP) filter is a popular econometric tool
More informationChapter 3: Regression Methods for Trends
Chapter 3: Regression Methods for Trends Time series exhibiting trends over time have a mean function that is some simple function (not necessarily constant) of time. The example random walk graph from
More informationIssues on quantile autoregression
Issues on quantile autoregression Jianqing Fan and Yingying Fan We congratulate Koenker and Xiao on their interesting and important contribution to the quantile autoregression (QAR). The paper provides
More informationAR, MA and ARMA models
AR, MA and AR by Hedibert Lopes P Based on Tsay s Analysis of Financial Time Series (3rd edition) P 1 Stationarity 2 3 4 5 6 7 P 8 9 10 11 Outline P Linear Time Series Analysis and Its Applications For
More information1 General problem. 2 Terminalogy. Estimation. Estimate θ. (Pick a plausible distribution from family. ) Or estimate τ = τ(θ).
Estimation February 3, 206 Debdeep Pati General problem Model: {P θ : θ Θ}. Observe X P θ, θ Θ unknown. Estimate θ. (Pick a plausible distribution from family. ) Or estimate τ = τ(θ). Examples: θ = (µ,
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 informationREGRESSION WITH SPATIALLY MISALIGNED DATA. Lisa Madsen Oregon State University David Ruppert Cornell University
REGRESSION ITH SPATIALL MISALIGNED DATA Lisa Madsen Oregon State University David Ruppert Cornell University SPATIALL MISALIGNED DATA 10 X X X X X X X X 5 X X X X X 0 X 0 5 10 OUTLINE 1. Introduction 2.
More informationModel comparison and selection
BS2 Statistical Inference, Lectures 9 and 10, Hilary Term 2008 March 2, 2008 Hypothesis testing Consider two alternative models M 1 = {f (x; θ), θ Θ 1 } and M 2 = {f (x; θ), θ Θ 2 } for a sample (X = x)
More informationST495: Survival Analysis: Hypothesis testing and confidence intervals
ST495: Survival Analysis: Hypothesis testing and confidence intervals Eric B. Laber Department of Statistics, North Carolina State University April 3, 2014 I remember that one fateful day when Coach took
More informationResearch Projects. Hanxiang Peng. March 4, Department of Mathematical Sciences Indiana University-Purdue University at Indianapolis
Hanxiang Department of Mathematical Sciences Indiana University-Purdue University at Indianapolis March 4, 2009 Outline Project I: Free Knot Spline Cox Model Project I: Free Knot Spline Cox Model Consider
More informationReview of Classical Least Squares. James L. Powell Department of Economics University of California, Berkeley
Review of Classical Least Squares James L. Powell Department of Economics University of California, Berkeley The Classical Linear Model The object of least squares regression methods is to model and estimate
More informationMISCELLANEOUS TOPICS RELATED TO LIKELIHOOD. Copyright c 2012 (Iowa State University) Statistics / 30
MISCELLANEOUS TOPICS RELATED TO LIKELIHOOD Copyright c 2012 (Iowa State University) Statistics 511 1 / 30 INFORMATION CRITERIA Akaike s Information criterion is given by AIC = 2l(ˆθ) + 2k, where l(ˆθ)
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 informationStatistics of stochastic processes
Introduction Statistics of stochastic processes Generally statistics is performed on observations y 1,..., y n assumed to be realizations of independent random variables Y 1,..., Y n. 14 settembre 2014
More informationTime Series and Forecasting Lecture 4 NonLinear Time Series
Time Series and Forecasting Lecture 4 NonLinear Time Series Bruce E. Hansen Summer School in Economics and Econometrics University of Crete July 23-27, 2012 Bruce Hansen (University of Wisconsin) Foundations
More informationCalibration Estimation of Semiparametric Copula Models with Data Missing at Random
Calibration Estimation of Semiparametric Copula Models with Data Missing at Random Shigeyuki Hamori 1 Kaiji Motegi 1 Zheng Zhang 2 1 Kobe University 2 Renmin University of China Institute of Statistics
More informationOn the Behavior of Marginal and Conditional Akaike Information Criteria in Linear Mixed Models
On the Behavior of Marginal and Conditional Akaike Information Criteria in Linear Mixed Models Thomas Kneib Department of Mathematics Carl von Ossietzky University Oldenburg Sonja Greven Department of
More informationEstimation of cumulative distribution function with spline functions
INTERNATIONAL JOURNAL OF ECONOMICS AND STATISTICS Volume 5, 017 Estimation of cumulative distribution function with functions Akhlitdin Nizamitdinov, Aladdin Shamilov Abstract The estimation of the cumulative
More informationTerminology Suppose we have N observations {x(n)} N 1. Estimators as Random Variables. {x(n)} N 1
Estimation Theory Overview Properties Bias, Variance, and Mean Square Error Cramér-Rao lower bound Maximum likelihood Consistency Confidence intervals Properties of the mean estimator Properties of the
More informationCan we do statistical inference in a non-asymptotic way? 1
Can we do statistical inference in a non-asymptotic way? 1 Guang Cheng 2 Statistics@Purdue www.science.purdue.edu/bigdata/ ONR Review Meeting@Duke Oct 11, 2017 1 Acknowledge NSF, ONR and Simons Foundation.
More information11. Bootstrap Methods
11. Bootstrap Methods c A. Colin Cameron & Pravin K. Trivedi 2006 These transparencies were prepared in 20043. They can be used as an adjunct to Chapter 11 of our subsequent book Microeconometrics: Methods
More informationStatistics 203: Introduction to Regression and Analysis of Variance Course review
Statistics 203: Introduction to Regression and Analysis of Variance Course review Jonathan Taylor - p. 1/?? Today Review / overview of what we learned. - p. 2/?? General themes in regression models Specifying
More informationWavelet Methods for Time Series Analysis. Part IV: Wavelet-Based Decorrelation of Time Series
Wavelet Methods for Time Series Analysis Part IV: Wavelet-Based Decorrelation of Time Series DWT well-suited for decorrelating certain time series, including ones generated from a fractionally differenced
More informationOne-stage dose-response meta-analysis
One-stage dose-response meta-analysis Nicola Orsini, Alessio Crippa Biostatistics Team Department of Public Health Sciences Karolinska Institutet http://ki.se/en/phs/biostatistics-team 2017 Nordic and
More informationLecture 3: Statistical sampling uncertainty
Lecture 3: Statistical sampling uncertainty c Christopher S. Bretherton Winter 2015 3.1 Central limit theorem (CLT) Let X 1,..., X N be a sequence of N independent identically-distributed (IID) random
More informationEconometric Forecasting
Robert M. Kunst robert.kunst@univie.ac.at University of Vienna and Institute for Advanced Studies Vienna October 1, 2014 Outline Introduction Model-free extrapolation Univariate time-series models Trend
More information401 Review. 6. Power analysis for one/two-sample hypothesis tests and for correlation analysis.
401 Review Major topics of the course 1. Univariate analysis 2. Bivariate analysis 3. Simple linear regression 4. Linear algebra 5. Multiple regression analysis Major analysis methods 1. Graphical analysis
More informationRLRsim: Testing for Random Effects or Nonparametric Regression Functions in Additive Mixed Models
RLRsim: Testing for Random Effects or Nonparametric Regression Functions in Additive Mixed Models Fabian Scheipl 1 joint work with Sonja Greven 1,2 and Helmut Küchenhoff 1 1 Department of Statistics, LMU
More informationLikelihood-based inference with missing data under missing-at-random
Likelihood-based inference with missing data under missing-at-random Jae-kwang Kim Joint work with Shu Yang Department of Statistics, Iowa State University May 4, 014 Outline 1. Introduction. Parametric
More informationNonparametric Small Area Estimation Using Penalized Spline Regression
Nonparametric Small Area Estimation Using Penalized Spline Regression J. D. Opsomer Iowa State University G. Claeskens Katholieke Universiteit Leuven M. G. Ranalli Universita degli Studi di Perugia G.
More informationRegularization in Cox Frailty Models
Regularization in Cox Frailty Models Andreas Groll 1, Trevor Hastie 2, Gerhard Tutz 3 1 Ludwig-Maximilians-Universität Munich, Department of Mathematics, Theresienstraße 39, 80333 Munich, Germany 2 University
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 informationSTAT 512 sp 2018 Summary Sheet
STAT 5 sp 08 Summary Sheet Karl B. Gregory Spring 08. Transformations of a random variable Let X be a rv with support X and let g be a function mapping X to Y with inverse mapping g (A = {x X : g(x A}
More informationAFT Models and Empirical Likelihood
AFT Models and Empirical Likelihood Mai Zhou Department of Statistics, University of Kentucky Collaborators: Gang Li (UCLA); A. Bathke; M. Kim (Kentucky) Accelerated Failure Time (AFT) models: Y = log(t
More informationEconometrics of Panel Data
Econometrics of Panel Data Jakub Mućk Meeting # 6 Jakub Mućk Econometrics of Panel Data Meeting # 6 1 / 36 Outline 1 The First-Difference (FD) estimator 2 Dynamic panel data models 3 The Anderson and Hsiao
More information1 Mixed effect models and longitudinal data analysis
1 Mixed effect models and longitudinal data analysis Mixed effects models provide a flexible approach to any situation where data have a grouping structure which introduces some kind of correlation between
More informationSupplement to Quantile-Based Nonparametric Inference for First-Price Auctions
Supplement to Quantile-Based Nonparametric Inference for First-Price Auctions Vadim Marmer University of British Columbia Artyom Shneyerov CIRANO, CIREQ, and Concordia University August 30, 2010 Abstract
More informationEconometrics of Panel Data
Econometrics of Panel Data Jakub Mućk Meeting # 2 Jakub Mućk Econometrics of Panel Data Meeting # 2 1 / 26 Outline 1 Fixed effects model The Least Squares Dummy Variable Estimator The Fixed Effect (Within
More informationSTATS 200: Introduction to Statistical Inference. Lecture 29: Course review
STATS 200: Introduction to Statistical Inference Lecture 29: Course review Course review We started in Lecture 1 with a fundamental assumption: Data is a realization of a random process. The goal throughout
More informationEconometrics I KS. Module 2: Multivariate Linear Regression. Alexander Ahammer. This version: April 16, 2018
Econometrics I KS Module 2: Multivariate Linear Regression Alexander Ahammer Department of Economics Johannes Kepler University of Linz This version: April 16, 2018 Alexander Ahammer (JKU) Module 2: Multivariate
More informationWhat s New in Econometrics. Lecture 13
What s New in Econometrics Lecture 13 Weak Instruments and Many Instruments Guido Imbens NBER Summer Institute, 2007 Outline 1. Introduction 2. Motivation 3. Weak Instruments 4. Many Weak) Instruments
More information10. Time series regression and forecasting
10. Time series regression and forecasting Key feature of this section: Analysis of data on a single entity observed at multiple points in time (time series data) Typical research questions: What is the
More informationA NOTE ON A NONPARAMETRIC REGRESSION TEST THROUGH PENALIZED SPLINES
Statistica Sinica 24 (2014), 1143-1160 doi:http://dx.doi.org/10.5705/ss.2012.230 A NOTE ON A NONPARAMETRIC REGRESSION TEST THROUGH PENALIZED SPLINES Huaihou Chen 1, Yuanjia Wang 2, Runze Li 3 and Katherine
More informationA Bayesian Nonparametric Approach to Monotone Missing Data in Longitudinal Studies with Informative Missingness
A Bayesian Nonparametric Approach to Monotone Missing Data in Longitudinal Studies with Informative Missingness A. Linero and M. Daniels UF, UT-Austin SRC 2014, Galveston, TX 1 Background 2 Working model
More informationCalibration Estimation of Semiparametric Copula Models with Data Missing at Random
Calibration Estimation of Semiparametric Copula Models with Data Missing at Random Shigeyuki Hamori 1 Kaiji Motegi 1 Zheng Zhang 2 1 Kobe University 2 Renmin University of China Econometrics Workshop UNC
More informationIllustration of the Varying Coefficient Model for Analyses the Tree Growth from the Age and Space Perspectives
TR-No. 14-06, Hiroshima Statistical Research Group, 1 11 Illustration of the Varying Coefficient Model for Analyses the Tree Growth from the Age and Space Perspectives Mariko Yamamura 1, Keisuke Fukui
More informationSTAT Financial Time Series
STAT 6104 - Financial Time Series Chapter 4 - Estimation in the time Domain Chun Yip Yau (CUHK) STAT 6104:Financial Time Series 1 / 46 Agenda 1 Introduction 2 Moment Estimates 3 Autoregressive Models (AR
More informationThe Restricted Likelihood Ratio Test at the Boundary in Autoregressive Series
The Restricted Likelihood Ratio Test at the Boundary in Autoregressive Series Willa W. Chen Rohit S. Deo July 6, 009 Abstract. The restricted likelihood ratio test, RLRT, for the autoregressive coefficient
More informationUncertainty Quantification for Inverse Problems. November 7, 2011
Uncertainty Quantification for Inverse Problems November 7, 2011 Outline UQ and inverse problems Review: least-squares Review: Gaussian Bayesian linear model Parametric reductions for IP Bias, variance
More informationWeb Appendix for Hierarchical Adaptive Regression Kernels for Regression with Functional Predictors by D. B. Woodard, C. Crainiceanu, and D.
Web Appendix for Hierarchical Adaptive Regression Kernels for Regression with Functional Predictors by D. B. Woodard, C. Crainiceanu, and D. Ruppert A. EMPIRICAL ESTIMATE OF THE KERNEL MIXTURE Here we
More informationBootstrap and Parametric Inference: Successes and Challenges
Bootstrap and Parametric Inference: Successes and Challenges G. Alastair Young Department of Mathematics Imperial College London Newton Institute, January 2008 Overview Overview Review key aspects of frequentist
More informationProf. Dr. Roland Füss Lecture Series in Applied Econometrics Summer Term Introduction to Time Series Analysis
Introduction to Time Series Analysis 1 Contents: I. Basics of Time Series Analysis... 4 I.1 Stationarity... 5 I.2 Autocorrelation Function... 9 I.3 Partial Autocorrelation Function (PACF)... 14 I.4 Transformation
More informationCOMS 4721: Machine Learning for Data Science Lecture 1, 1/17/2017
COMS 4721: Machine Learning for Data Science Lecture 1, 1/17/2017 Prof. John Paisley Department of Electrical Engineering & Data Science Institute Columbia University OVERVIEW This class will cover model-based
More informationCOMPARING PARAMETRIC AND SEMIPARAMETRIC ERROR CORRECTION MODELS FOR ESTIMATION OF LONG RUN EQUILIBRIUM BETWEEN EXPORTS AND IMPORTS
Applied Studies in Agribusiness and Commerce APSTRACT Center-Print Publishing House, Debrecen DOI: 10.19041/APSTRACT/2017/1-2/3 SCIENTIFIC PAPER COMPARING PARAMETRIC AND SEMIPARAMETRIC ERROR CORRECTION
More informationMasters Comprehensive Examination Department of Statistics, University of Florida
Masters Comprehensive Examination Department of Statistics, University of Florida May 6, 003, 8:00 am - :00 noon Instructions: You have four hours to answer questions in this examination You must show
More informationSTAT 443 Final Exam Review. 1 Basic Definitions. 2 Statistical Tests. L A TEXer: W. Kong
STAT 443 Final Exam Review L A TEXer: W Kong 1 Basic Definitions Definition 11 The time series {X t } with E[X 2 t ] < is said to be weakly stationary if: 1 µ X (t) = E[X t ] is independent of t 2 γ X
More informationGeneralized Linear Models. Kurt Hornik
Generalized Linear Models Kurt Hornik Motivation Assuming normality, the linear model y = Xβ + e has y = β + ε, ε N(0, σ 2 ) such that y N(μ, σ 2 ), E(y ) = μ = β. Various generalizations, including general
More informationNonparametric Tests for Multi-parameter M-estimators
Nonparametric Tests for Multi-parameter M-estimators John Robinson School of Mathematics and Statistics University of Sydney The talk is based on joint work with John Kolassa. It follows from work over
More informationA Diagnostic for Seasonality Based Upon Autoregressive Roots
A Diagnostic for Seasonality Based Upon Autoregressive Roots Tucker McElroy (U.S. Census Bureau) 2018 Seasonal Adjustment Practitioners Workshop April 26, 2018 1 / 33 Disclaimer This presentation is released
More informationDo Markov-Switching Models Capture Nonlinearities in the Data? Tests using Nonparametric Methods
Do Markov-Switching Models Capture Nonlinearities in the Data? Tests using Nonparametric Methods Robert V. Breunig Centre for Economic Policy Research, Research School of Social Sciences and School of
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 informationEcon 623 Econometrics II Topic 2: Stationary Time Series
1 Introduction Econ 623 Econometrics II Topic 2: Stationary Time Series In the regression model we can model the error term as an autoregression AR(1) process. That is, we can use the past value of the
More informationRobustness to Parametric Assumptions in Missing Data Models
Robustness to Parametric Assumptions in Missing Data Models Bryan Graham NYU Keisuke Hirano University of Arizona April 2011 Motivation Motivation We consider the classic missing data problem. In practice
More informationAlternative implementations of Monte Carlo EM algorithms for likelihood inferences
Genet. Sel. Evol. 33 001) 443 45 443 INRA, EDP Sciences, 001 Alternative implementations of Monte Carlo EM algorithms for likelihood inferences Louis Alberto GARCÍA-CORTÉS a, Daniel SORENSEN b, Note a
More informationNonparametric Bayesian Methods - Lecture I
Nonparametric Bayesian Methods - Lecture I Harry van Zanten Korteweg-de Vries Institute for Mathematics CRiSM Masterclass, April 4-6, 2016 Overview of the lectures I Intro to nonparametric Bayesian statistics
More informationMIT Spring 2015
Regression Analysis MIT 18.472 Dr. Kempthorne Spring 2015 1 Outline Regression Analysis 1 Regression Analysis 2 Multiple Linear Regression: Setup Data Set n cases i = 1, 2,..., n 1 Response (dependent)
More informationReview of Statistics
Review of Statistics Topics Descriptive Statistics Mean, Variance Probability Union event, joint event Random Variables Discrete and Continuous Distributions, Moments Two Random Variables Covariance and
More information