Space-time modelling of air pollution with array methods
|
|
- Alan Berry
- 6 years ago
- Views:
Transcription
1 Space-time modelling of air pollution with array methods Dae-Jin Lee Royal Statistical Society Conference Edinburgh 2009 D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 1
2 Motivation Increasing research on modelling spatio-temporal data Wide variety of approaches and research communities E.g. : Environmental data, epidemiologic studies, disease mapping applications,... European Environmental Agency (EEA): Monitoring networks EMEP project (European Monitoring and Evaluation Programme) Ozone (O 3 ) is currently one of the air pollutants of most concern in Europe. D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 2
3 Monitoring stations across Europe Sweden Latitude UK Austria Spain sample of 45 monitoring stations Longitude Monitoring station D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 3
4 O 3 time series plot for selected locations Seasonal pattern: O Spain Sweden Austria UK year D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 4
5 Motivation Spatio-temporal data Response variable, y ijt measured over geographical locations, s = (xi, x j ), with i, j = 1,.., n and over time periods, xt, for t = 1,..., T ISSUE: huge amount of data available Smoothing techniques: Study spatial and temporal trends. Space and time interactions. 3-dimensional smoothing: P-splines and GLAM. D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 5
6 Outline 1 P-splines in spatio-temporal smoothing context 2 ANOVA-Type Interaction Models 3 Application to air pollution data 4 Concluding remarks D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 6
7 Outline 1 P-splines in spatio-temporal smoothing context 2 ANOVA-Type Interaction Models 3 Application to air pollution data 4 Concluding remarks D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 7
8 Penalized splines GLAM in 3d 3d-case: f (x 1, x 2, x 3 ) = Bθ θ can be expressed as a 3d-array Θ = {θ} ijk of dim. c 1 c 2 c 3 θ (1,1,c3 ) θ (1,c2,c 3 ) layer 1,...,c 3 columns θ (1,c2,1) θ (1,1,1) 1,...,c 2 rows 1,...,c 1 θ (c1,1,c 3 ) θ (c1,c 2,c 3 ) θ (c1,1,1) θ (c1,c 2,1) D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 8
9 3d-Penalty matrix: Set penalties over the 3d-array Θ: P = λ 1 D 1D 1 I c2 I c3 +λ 2 I c1 D 2D 2 I c3 +λ t I c1 I c2 D td t {z } {z } {z } row-wise column-wise layer-wise For spatio-temporal data: f ( longitude, latitude, time) }{{} Space Spatial anisotropy (λ 1 λ 2 ), different amount of smoothing for latitude and longitude. Temporal smoothing (λ t) Space-time interaction. D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 9
10 Penalized splines Spatio-Temporal data smoothing For spatio-temporal data, we propose: Spatio-temporal B-splines Basis: B = B s B t, of dim. nt c 1 c 2 c 3 where B s is the spatial B-spline basis (B 1 B 2 ) and B t is the B-spline basis for time of dim. t c 3. Note that: GLAM framework E[Y] = B t ΘB s Mixed model D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 10
11 Penalized splines Spatio-Temporal data smoothing Most of the space-time approaches consider Simplest additive model (ignores space-time interaction): f s (x 1, x 2 ) + f t (x t ) We proposed space-time interaction models 3d: f st (x 1, x 2, x t ) Space-time ANOVA: f s (x 1, x 2 ) + f t (x t ) + f st (x 1, x 2, x t ) D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 11
12 Penalized splines Spatio-Temporal data smoothing Most of the space-time approaches consider Simplest additive model (ignores space-time interaction): f s (x 1, x 2 ) + f t (x t ) We proposed space-time interaction models 3d: f st (x 1, x 2, x t ) Space-time ANOVA: f s (x 1, x 2 ) + f t (x t ) + f st (x 1, x 2, x t ) D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 11
13 Penalized splines Spatio-Temporal data smoothing Most of the space-time approaches consider Simplest additive model (ignores space-time interaction): f s (x 1, x 2 ) + f t (x t ) We proposed space-time interaction models 3d: f st (x 1, x 2, x t ) Space-time ANOVA: f s (x 1, x 2 ) + f t (x t ) + f st (x 1, x 2, x t ) D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 11
14 Outline 1 P-splines in spatio-temporal smoothing context 2 ANOVA-Type Interaction Models 3 Application to air pollution data 4 Concluding remarks D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 12
15 ANOVA-Type Interaction Models Chen (1993), Gu (2002): Smoothing-Spline ANOVA (SS-ANOVA). Interpretation as main effects and interactions. Models of type: ŷ = f (x 1 ) + f (x 2 ) + f (x t ) Main/additive effects + f (x 1, x 2 ) + f (x 1, x t ) + f (x 2, x t ) 2-way interactions + f (x 1, x 2, x t ) 3-way interactions PROBLEMS: identifiability, and basis dimension ( curse of dimensionality ) D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 13
16 P-spline ANOVA model for spatio-temporal smoothing Lee and Durbán (2009b), consider: where y = γ + f s (x 1, x 2 ) + f t (time) + f st (x 1, x 2, time) + ɛ, f s (x 1, x 2 ) Spatial 2d smooth surface f t (time) Smooth time trend f st (x 1, x 2, time) Space-time interaction We need to construct an identifiable model. Our approach is based on: low-rank basis (P-splines) the mixed model representation and SVD properties. D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 14
17 P-spline ANOVA model for spatio-temporal smoothing Model: by = Intercept + f s(x 1, x 2 ) + f t(time) + f st(x 1, x 2, time) with Basis and Coefficients: ( ) Bθ = 1 nt B s 1 t 1 n B t B s B t and Penalty: 0 P = where τ 1 D 1D 1 I c2 + τ 2 I c1 D 2D 2 τ t D td t γ θ (s) θ (t) θ (st) P st = λ 1 D 1 D 1 I c2 I c3 + λ 2 I c1 D 2 D 2 I c3 + λ t I c1 I c2 D t D t P st D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 15
18 P-spline ANOVA model for spatio-temporal smoothing Model: by = Intercept + f s(x 1, x 2 ) + f s(time) + f st(x 1, x 2, time) with Basis and Coefficients: ( ) Bθ = 1 nt B s 1 t 1 n B t B s B t and Penalty: 0 P = where τ 1 D 1D 1 I c2 + τ 2 I c1 D 2D 2 τ t D td t γ θ (s) θ (t) θ (st) P st = λ 1 D 1 D 1 I c2 I c3 + λ 2 I c1 D 2 D 2 I c3 + λ t I c1 I c2 D t D t P st D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 15
19 P-spline ANOVA model for spatio-temporal smoothing Model: by = Intercept + f s(x 1, x 2 ) + f s(time) + f st(x 1, x 2, time) with Basis and Coefficients: ( ) Bθ = 1 nt B s 1 t 1 n B t B s B t and Penalty: 0 P = where τ 1 D 1D 1 I c2 + τ 2 I c1 D 2D 2 τ t D td t γ θ (s) θ (t) θ (st) P st = λ 1 D 1 D 1 I c2 I c3 + λ 2 I c1 D 2 D 2 I c3 + λ t I c1 I c2 D t D t P st D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 15
20 P-spline ANOVA model for spatio-temporal smoothing Model: by = Intercept + f s(x 1, x 2 ) + f s(time) + f st(x 1, x 2, time) with Basis and Coefficients: ( ) Bθ = 1 nt B s 1 t 1 n B t B s B t and Penalty: 0 P = where τ 1 D 1D 1 I c2 + τ 2 I c1 D 2D 2 τ t D td t γ θ (s) θ (t) θ (st) P st = λ 1 D 1 D 1 I c2 I c3 + λ 2 I c1 D 2 D 2 I c3 + λ t I c1 I c2 D t D t P st D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 15
21 P-spline ANOVA model for spatio-temporal smoothing Model: by = Intercept + f s(x 1, x 2 ) + f s(time) + f st(x 1, x 2, time) with Basis and Coefficients: ( ) Bθ = 1 nt B s 1 t 1 n B t B s B t and Penalty: 0 P = where τ 1 D 1D 1 I c2 + τ 2 I c1 D 2D 2 τ t D td t γ θ (s) θ (t) θ (st) P st = λ 1 D 1 D 1 I c2 I c3 + λ 2 I c1 D 2 D 2 I c3 + λ t I c1 I c2 D t D t P st D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 15
22 P-spline ANOVA model for spatio-temporal smoothing Model: by = Intercept + f s(x 1, x 2 ) + f s(time) + f st(x 1, x 2, time) with Basis and Coefficients: ( ) Bθ = 1 nt B s 1 t 1 n B t B s B t and Penalty: 0 P = where τ 1 D 1D 1 I c2 + τ 2 I c1 D 2D 2 τ t D td t γ θ (s) θ (t) θ (st) P st = λ 1 D 1 D 1 I c2 I c3 + λ 2 I c1 D 2 D 2 I c3 + λ t I c1 I c2 D t D t P st D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 15
23 In GLAM notation: ( B s 1 t ) θ (s) 1 t Θ s B s ( 1 n B t ) θ (t) B t Θ t 1 n ( B s B t ) θ (s,t) B t Θ st B s D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 16
24 P-spline ANOVA model for spatio-temporal smoothing We avoid identifiability problems using Mixed model reparameterization and SVD properties For each term we have: Basis [ X : Z ] f s (x 1, x 2 ) x 1 : x 2 (1) f t (x t ) x t (2) f st (x 1, x 2, x t ) x 1 : x 2 : x t (3) Some terms in (1) and (2) also appear in (3). Remove linearly dependent columns in the basis» Constraints D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 17
25 P-spline ANOVA model for spatio-temporal smoothing We avoid identifiability problems using Mixed model reparameterization and SVD properties For each term we have: Basis [ X : Z ] f s (x 1, x 2 ) x 1 : x 2 (1) f t (x t ) x t (2) f st (x 1, x 2, x t ) x 1 : x 2 : x t (3) Some terms in (1) and (2) also appear in (3). Remove linearly dependent columns in the basis» Constraints D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 17
26 P-spline ANOVA model for spatio-temporal smoothing Equivalent to apply constraints over the coefficients: ct t=1 θ(t) t = 0 (time) c 1 i θ (st) t,ij = c 2 j θ (st) t,ij = c 1 i c 2 j θ (st) t,ij = 0 (space-time) or in array form over Θ ijt» Array
27 c1 i θ (st) t,ij = c 2 j θ (st) t,ij = c 1 c2 i j θ (st) t,ij = 0 Centering and scaling matrix: ( I c 11 /c ) θ (1,1,c3 ) θ (1,c2,c 3 ) layer 1,...,c columns 3 θ (1,1,1) 1,...,c 2 θ (1,c2,1) rows 1,...,c 1 θ (c1,1,c 3 ) θ (c1,c 2,c 3 ) θ (c1,1,1) θ (c1,c 2,1) D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 19
28 c1 i θ (st) t,ij = c 2 j θ (st) t,ij = c 1 c2 i j θ (st) t,ij = 0 Centering and scaling matrix: ( I c 11 /c ) θ (1,1,c3 ) θ (1,c2,c 3 ) layer 1,...,c columns 3 θ (1,1,1) 1,...,c 2 θ (1,c2,1) rows 1,...,c 1 θ (c1,1,c 3 ) θ (c1,c 2,c 3 ) θ (c1,1,1) θ (c1,c 2,1) D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 19
29 c1 i θ (st) t,ij = c 2 j θ (st) t,ij = c 1 c2 i j θ (st) t,ij = 0 Centering and scaling matrix: ( I c 11 /c ) θ (1,1,c3 ) θ (1,c2,c 3 ) layer 1,...,c columns 3 θ (1,1,1) 1,...,c 2 θ (1,c2,1) rows 1,...,c 1 θ (c1,1,c 3 ) θ (c1,c 2,c 3 ) θ (c1,1,1) θ (c1,c 2,1) D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 19
30 In practice We only need to construct the matrices X, Z and penalty F f s (x 1, x 2 ) f t (x t ) f st (x 1, x 2, x t ) X by columns x 1 : x 2 x t (x 1, x 2, x t ) Z by blocks F blockdiagonal F s F t F st (λ 1, λ 2, λ t ) λ t (τ 1, τ 2, τ t ) And estimate by REML D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 20
31 Outline 1 P-splines in spatio-temporal smoothing context 2 ANOVA-Type Interaction Models 3 Application to air pollution data 4 Concluding remarks D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 21
32 Ozone pollution in Europe Lee and Durbán (2009b) Sample of 45 monitoring stations Monthly averages of O 3 levels (in µg/m 3 units) from january 1999 to december 2005 (t = 1,..., 84) Models: Additive: Spatio-temporal Interaction: 3d: f s (x 1, x 2 ) + f t (x t ) f st(x 1, x 2, x t) ANOVA: f s(x 1, x 2 ) + f t(x t) + f st(x 1, x 2, x t) D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 22
33 Spatial 2d + time f s (x 1, x 2 ) + f t (x t ) Longitude f(time) Latitude year Space-time interaction is not considered time smooth trend is additive D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 23
34 Spatio-temporal ANOVA model Play animation ^ y = f(space) + f(space,time) + f(time) 1999 : 1 D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 24
35 Comparison of fitted values Additive VS ANOVA Additive model fit ANOVA model fit f s(x 1, x 2 ) + f s(x t) f s(x 1, x 2 ) + f t(x t) + f st(x 1, x 2, x t) O Spain Sweden Austria UK O Spain Sweden Austria UK year year Additive model assumes a spatial smooth surface over all monitoring stations that remains constant over time. ANOVA model captures individual characteristics of the stations throughout time. D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 25
36 Comparison of Models ANOVA, 3d and Additive Model AIC df ANOVA 14280,73 366,03 3d interaction 14537,22 765,05 Additive 16506,28 65,98 Observations: Best overall performance of ANOVA in terms of AIC also with less d.f. than 3d ANOVA model is more realistic than Additive, and easier to decompose and interpret in terms of the fit. D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 26
37 Outline 1 P-splines in spatio-temporal smoothing context 2 ANOVA-Type Interaction Models 3 Application to air pollution data 4 Concluding remarks D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 27
38 Concluding remarks Use of Array methods in space-time data speeds up calculations data and model matrices are stored efficiently Spatio-temporal ANOVA model Flexible formulation and interpretability of smoothing Computational efficiency: nested B-spline bases consider a smaller basis B t for the interaction B s B t... such that rank( B t ) < rank(b t ) The size of the full basis B is reduced and model is nested f s (x 1, x 2 ) + f t (x t ) + f st(x 1, x 2, x t ) f s (x 1, x 2 ) + f t (x t ) + f st (x 1, x 2, x t ) D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 28
39 Concluding remarks Use of Array methods in space-time data speeds up calculations data and model matrices are stored efficiently Spatio-temporal ANOVA model Flexible formulation and interpretability of smoothing Computational efficiency: nested B-spline bases consider a smaller basis B t for the interaction B s B t... such that rank( B t ) < rank(b t ) The size of the full basis B is reduced and model is nested f s (x 1, x 2 ) + f t (x t ) + f st(x 1, x 2, x t ) f s (x 1, x 2 ) + f t (x t ) + f st (x 1, x 2, x t ) D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 28
40 Concluding remarks Use of Array methods in space-time data speeds up calculations data and model matrices are stored efficiently Spatio-temporal ANOVA model Flexible formulation and interpretability of smoothing Computational efficiency: nested B-spline bases consider a smaller basis B t for the interaction B s B t... such that rank( B t ) < rank(b t ) The size of the full basis B is reduced and model is nested f s (x 1, x 2 ) + f t (x t ) + f st(x 1, x 2, x t ) f s (x 1, x 2 ) + f t (x t ) + f st (x 1, x 2, x t ) D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 28
41 Concluding remarks Use of Array methods in space-time data speeds up calculations data and model matrices are stored efficiently Spatio-temporal ANOVA model Flexible formulation and interpretability of smoothing Computational efficiency: nested B-spline bases consider a smaller basis B t for the interaction B s B t... such that rank( B t ) < rank(b t ) The size of the full basis B is reduced and model is nested f s (x 1, x 2 ) + f t (x t ) + f st(x 1, x 2, x t ) f s (x 1, x 2 ) + f t (x t ) + f st (x 1, x 2, x t ) D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 28
42 THANKS FOR YOUR ATTENTION!!!
43 Bibliography Spatial and Spatio-temporal smoothing with P-splines: Lee, D.-J. and Durbán, M. (2009a) Smooth-CAR mixed models for spatial count data. CSDA 53(8): Lee, D.-J. and Durbán, M. (2009b) P-spline ANOVA-Type interaction models for spatio-temporal smoothing Submitted P-splines: Eilers, PHC. and Marx, BD. Flexible smoothing with B-splines and penalties Stat. Sci. (1996), 11: Currie, ID., Durbán M. and Eilers, PHC. (2006) Generalized linear array models with applications to multidimensional smoothing JRSSB, 68:1-22 D.-J. Lee (Uc3m) GLAM: Array methods in Statistics RSS 09 - Edinburgh # 30
P -spline ANOVA-type interaction models for spatio-temporal smoothing
P -spline ANOVA-type interaction models for spatio-temporal smoothing Dae-Jin Lee 1 and María Durbán 1 1 Department of Statistics, Universidad Carlos III de Madrid, SPAIN. e-mail: dae-jin.lee@uc3m.es and
More informationFlexible Spatio-temporal smoothing with array methods
Int. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin (Session IPS046) p.849 Flexible Spatio-temporal smoothing with array methods Dae-Jin Lee CSIRO, Mathematics, Informatics and
More informationCurrie, Iain Heriot-Watt University, Department of Actuarial Mathematics & Statistics Edinburgh EH14 4AS, UK
An Introduction to Generalized Linear Array Models Currie, Iain Heriot-Watt University, Department of Actuarial Mathematics & Statistics Edinburgh EH14 4AS, UK E-mail: I.D.Currie@hw.ac.uk 1 Motivating
More informationGLAM An Introduction to Array Methods in Statistics
GLAM An Introduction to Array Methods in Statistics Iain Currie Heriot Watt University GLAM A Generalized Linear Array Model is a low-storage, high-speed, method for multidimensional smoothing, when data
More informationMultidimensional Density Smoothing with P-splines
Multidimensional Density Smoothing with P-splines Paul H.C. Eilers, Brian D. Marx 1 Department of Medical Statistics, Leiden University Medical Center, 300 RC, Leiden, The Netherlands (p.eilers@lumc.nl)
More informationMULTIDIMENSIONAL COVARIATE EFFECTS IN SPATIAL AND JOINT EXTREMES
MULTIDIMENSIONAL COVARIATE EFFECTS IN SPATIAL AND JOINT EXTREMES Philip Jonathan, Kevin Ewans, David Randell, Yanyun Wu philip.jonathan@shell.com www.lancs.ac.uk/ jonathan Wave Hindcasting & Forecasting
More informationA Hierarchical Perspective on Lee-Carter Models
A Hierarchical Perspective on Lee-Carter Models Paul Eilers Leiden University Medical Centre L-C Workshop, Edinburgh 24 The vantage point Previous presentation: Iain Currie led you upward From Glen Gumbel
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 informationUsing P-splines to smooth two-dimensional Poisson data
1 Using P-splines to smooth two-dimensional Poisson data Maria Durbán 1, Iain Currie 2, Paul Eilers 3 17th IWSM, July 2002. 1 Dept. Statistics and Econometrics, Universidad Carlos III de Madrid, Spain.
More informationAn Introduction to GAMs based on penalized regression splines. Simon Wood Mathematical Sciences, University of Bath, U.K.
An Introduction to GAMs based on penalied regression splines Simon Wood Mathematical Sciences, University of Bath, U.K. Generalied Additive Models (GAM) A GAM has a form something like: g{e(y i )} = η
More informationBayesian covariate models in extreme value analysis
Bayesian covariate models in extreme value analysis David Randell, Philip Jonathan, Kathryn Turnbull, Mathew Jones EVA 2015 Ann Arbor Copyright 2015 Shell Global Solutions (UK) EVA 2015 Ann Arbor June
More informationSTATISTICS 174: APPLIED STATISTICS FINAL EXAM DECEMBER 10, 2002
Time allowed: 3 HOURS. STATISTICS 174: APPLIED STATISTICS FINAL EXAM DECEMBER 10, 2002 This is an open book exam: all course notes and the text are allowed, and you are expected to use your own calculator.
More informationSpatio-Temporal Expectile Regression Models
Spatio-Temporal Expectile Regression Models Elmar Spiegel University of Goettingen Abstract Spatio-temporal models are becoming increasingly popular in recent regression research. However, they usually
More informationOne-Way ANOVA Model. group 1 y 11, y 12 y 13 y 1n1. group 2 y 21, y 22 y 2n2. group g y g1, y g2, y gng. g is # of groups,
One-Way ANOVA Model group 1 y 11, y 12 y 13 y 1n1 group 2 y 21, y 22 y 2n2 group g y g1, y g2, y gng g is # of groups, n i denotes # of obs in the i-th group, and the total sample size n = g i=1 n i. 1
More informationLinear Regression Models P8111
Linear Regression Models P8111 Lecture 25 Jeff Goldsmith April 26, 2016 1 of 37 Today s Lecture Logistic regression / GLMs Model framework Interpretation Estimation 2 of 37 Linear regression Course started
More informationA Modern Look at Classical Multivariate Techniques
A Modern Look at Classical Multivariate Techniques Yoonkyung Lee Department of Statistics The Ohio State University March 16-20, 2015 The 13th School of Probability and Statistics CIMAT, Guanajuato, Mexico
More informationEstimating the long-term health impact of air pollution using spatial ecological studies. Duncan Lee
Estimating the long-term health impact of air pollution using spatial ecological studies Duncan Lee EPSRC and RSS workshop 12th September 2014 Acknowledgements This is joint work with Alastair Rushworth
More informationStat 217 Final Exam. Name: May 1, 2002
Stat 217 Final Exam Name: May 1, 2002 Problem 1. Three brands of batteries are under study. It is suspected that the lives (in weeks) of the three brands are different. Five batteries of each brand are
More informationPartial factor modeling: predictor-dependent shrinkage for linear regression
modeling: predictor-dependent shrinkage for linear Richard Hahn, Carlos Carvalho and Sayan Mukherjee JASA 2013 Review by Esther Salazar Duke University December, 2013 Factor framework The factor framework
More informationVine Copulas. Spatial Copula Workshop 2014 September 22, Institute for Geoinformatics University of Münster.
Spatial Workshop 2014 September 22, 2014 Institute for Geoinformatics University of Münster http://ifgi.uni-muenster.de/graeler 1 spatio-temporal data Typically, spatio-temporal data is given at a set
More informationYORK UNIVERSITY. Faculty of Science Department of Mathematics and Statistics MATH M Test #1. July 11, 2013 Solutions
YORK UNIVERSITY Faculty of Science Department of Mathematics and Statistics MATH 222 3. M Test # July, 23 Solutions. For each statement indicate whether it is always TRUE or sometimes FALSE. Note: For
More informationRecovering Indirect Information in Demographic Applications
Recovering Indirect Information in Demographic Applications Jutta Gampe Abstract In many demographic applications the information of interest can only be estimated indirectly. Modelling events and rates
More informationMultiple linear regression: estimation and model fitting
Multiple linear regression: estimation and model fitting January 25 Introduction The goal of today s class is to set up a multiple regression model in terms of matrices and then solve for the regression
More informationTowards a Regression using Tensors
February 27, 2014 Outline Background 1 Background Linear Regression Tensorial Data Analysis 2 Definition Tensor Operation Tensor Decomposition 3 Model Attention Deficit Hyperactivity Disorder Data Analysis
More informationSpatial bias modeling with application to assessing remotely-sensed aerosol as a proxy for particulate matter
Spatial bias modeling with application to assessing remotely-sensed aerosol as a proxy for particulate matter Chris Paciorek Department of Biostatistics Harvard School of Public Health application joint
More informationSpatiotemporal smoothing and sulphur dioxide trends over Europe
Appl. Statist. (2009) 58, Part 5, pp. 737 752 Spatiotemporal smoothing and sulphur dioxide trends over Europe Adrian W. Bowman, University of Glasgow, UK Marco Giannitrapani Shell Global Solutions International,
More informationFunctional SVD for Big Data
Functional SVD for Big Data Pan Chao April 23, 2014 Pan Chao Functional SVD for Big Data April 23, 2014 1 / 24 Outline 1 One-Way Functional SVD a) Interpretation b) Robustness c) CV/GCV 2 Two-Way Problem
More informationSTAT 430 (Fall 2017): Tutorial 8
STAT 430 (Fall 2017): Tutorial 8 Balanced Incomplete Block Design Luyao Lin November 7th/9th, 2017 Department Statistics and Actuarial Science, Simon Fraser University Block Design Complete Random Complete
More informationObjectives Simple linear regression. Statistical model for linear regression. Estimating the regression parameters
Objectives 10.1 Simple linear regression Statistical model for linear regression Estimating the regression parameters Confidence interval for regression parameters Significance test for the slope Confidence
More informationMultivariate Calibration with Robust Signal Regression
Multivariate Calibration with Robust Signal Regression Bin Li and Brian Marx from Louisiana State University Somsubhra Chakraborty from Indian Institute of Technology Kharagpur David C Weindorf from Texas
More informationAn application of the GAM-PCA-VAR model to respiratory disease and air pollution data
An application of the GAM-PCA-VAR model to respiratory disease and air pollution data Márton Ispány 1 Faculty of Informatics, University of Debrecen Hungary Joint work with Juliana Bottoni de Souza, Valdério
More informationA short introduction to INLA and R-INLA
A short introduction to INLA and R-INLA Integrated Nested Laplace Approximation Thomas Opitz, BioSP, INRA Avignon Workshop: Theory and practice of INLA and SPDE November 7, 2018 2/21 Plan for this talk
More informationLecture 10. Factorial experiments (2-way ANOVA etc)
Lecture 10. Factorial experiments (2-way ANOVA etc) Jesper Rydén Matematiska institutionen, Uppsala universitet jesper@math.uu.se Regression and Analysis of Variance autumn 2014 A factorial experiment
More informationAndrew B. Lawson 2019 BMTRY 763
BMTRY 763 FMD Foot and mouth disease (FMD) is a viral disease of clovenfooted animals and being extremely contagious it spreads rapidly. It reduces animal production, and can sometimes be fatal in young
More informationBiostatistics 533 Classical Theory of Linear Models Spring 2007 Final Exam. Please choose ONE of the following options.
1 Biostatistics 533 Classical Theory of Linear Models Spring 2007 Final Exam Name: Problems do not have equal value and some problems will take more time than others. Spend your time wisely. You do not
More informationSymmetry and Separability In Spatial-Temporal Processes
Symmetry and Separability In Spatial-Temporal Processes Man Sik Park, Montserrat Fuentes Symmetry and Separability In Spatial-Temporal Processes 1 Motivation In general, environmental data have very complex
More informationANOVA Situation The F Statistic Multiple Comparisons. 1-Way ANOVA MATH 143. Department of Mathematics and Statistics Calvin College
1-Way ANOVA MATH 143 Department of Mathematics and Statistics Calvin College An example ANOVA situation Example (Treating Blisters) Subjects: 25 patients with blisters Treatments: Treatment A, Treatment
More informationPaper Review: NONSTATIONARY COVARIANCE MODELS FOR GLOBAL DATA
Paper Review: NONSTATIONARY COVARIANCE MODELS FOR GLOBAL DATA BY MIKYOUNG JUN AND MICHAEL L. STEIN Presented by Sungkyu Jung April, 2009 Outline 1 Introduction 2 Covariance Models 3 Application: Level
More information3. For a given dataset and linear model, what do you think is true about least squares estimates? Is Ŷ always unique? Yes. Is ˆβ always unique? No.
7. LEAST SQUARES ESTIMATION 1 EXERCISE: Least-Squares Estimation and Uniqueness of Estimates 1. For n real numbers a 1,...,a n, what value of a minimizes the sum of squared distances from a to each of
More informationIntroduction to the Tensor Train Decomposition and Its Applications in Machine Learning
Introduction to the Tensor Train Decomposition and Its Applications in Machine Learning Anton Rodomanov Higher School of Economics, Russia Bayesian methods research group (http://bayesgroup.ru) 14 March
More informationLecture 1 From Continuous-Time to Discrete-Time
Lecture From Continuous-Time to Discrete-Time Outline. Continuous and Discrete-Time Signals and Systems................. What is a signal?................................2 What is a system?.............................
More informationModel checking overview. Checking & Selecting GAMs. Residual checking. Distribution checking
Model checking overview Checking & Selecting GAMs Simon Wood Mathematical Sciences, University of Bath, U.K. Since a GAM is just a penalized GLM, residual plots should be checked exactly as for a GLM.
More informationSpatio-temporal modelling of daily air temperature in Catalonia
Spatio-temporal modelling of daily air temperature in Catalonia M. Saez 1,, M.A. Barceló 1,, A. Tobias 3, D. Varga 1,4 and R. Ocaña-Riola 5 1 Research Group on Statistics, Applied Economics and Health
More informationLinear Algebra I Lecture 8
Linear Algebra I Lecture 8 Xi Chen 1 1 University of Alberta January 25, 2019 Outline 1 2 Gauss-Jordan Elimination Given a system of linear equations f 1 (x 1, x 2,..., x n ) = 0 f 2 (x 1, x 2,..., x n
More information1-Way ANOVA MATH 143. Spring Department of Mathematics and Statistics Calvin College
1-Way ANOVA MATH 143 Department of Mathematics and Statistics Calvin College Spring 2010 The basic ANOVA situation Two variables: 1 Categorical, 1 Quantitative Main Question: Do the (means of) the quantitative
More informationBayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes
Bayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes Alan Gelfand 1 and Andrew O. Finley 2 1 Department of Statistical Science, Duke University, Durham, North
More informationINTRODUCTION. In March 1998, the tender for project CT.98.EP.04 was awarded to the Department of Medicines Management, Keele University, UK.
INTRODUCTION In many areas of Europe patterns of drug use are changing. The mechanisms of diffusion are diverse: introduction of new practices by new users, tourism and migration, cross-border contact,
More informationAdvanced Statistics I : Gaussian Linear Model (and beyond)
Advanced Statistics I : Gaussian Linear Model (and beyond) Aurélien Garivier CNRS / Telecom ParisTech Centrale Outline One and Two-Sample Statistics Linear Gaussian Model Model Reduction and model Selection
More informationBayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes
Bayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes Sudipto Banerjee 1 and Andrew O. Finley 2 1 Biostatistics, School of Public Health, University of Minnesota,
More informationEstimating prediction error in mixed models
Estimating prediction error in mixed models benjamin saefken, thomas kneib georg-august university goettingen sonja greven ludwig-maximilians-university munich 1 / 12 GLMM - Generalized linear mixed models
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 informationVariable Selection and Model Choice in Survival Models with Time-Varying Effects
Variable Selection and Model Choice in Survival Models with Time-Varying Effects Boosting Survival Models Benjamin Hofner 1 Department of Medical Informatics, Biometry and Epidemiology (IMBE) Friedrich-Alexander-Universität
More informationHierarchical Nearest-Neighbor Gaussian Process Models for Large Geo-statistical Datasets
Hierarchical Nearest-Neighbor Gaussian Process Models for Large Geo-statistical Datasets Abhirup Datta 1 Sudipto Banerjee 1 Andrew O. Finley 2 Alan E. Gelfand 3 1 University of Minnesota, Minneapolis,
More informationTopic 17 - Single Factor Analysis of Variance. Outline. One-way ANOVA. The Data / Notation. One way ANOVA Cell means model Factor effects model
Topic 17 - Single Factor Analysis of Variance - Fall 2013 One way ANOVA Cell means model Factor effects model Outline Topic 17 2 One-way ANOVA Response variable Y is continuous Explanatory variable is
More informationCOMS 4721: Machine Learning for Data Science Lecture 18, 4/4/2017
COMS 4721: Machine Learning for Data Science Lecture 18, 4/4/2017 Prof. John Paisley Department of Electrical Engineering & Data Science Institute Columbia University TOPIC MODELING MODELS FOR TEXT DATA
More informationParameter Estimation in the Spatio-Temporal Mixed Effects Model Analysis of Massive Spatio-Temporal Data Sets
Parameter Estimation in the Spatio-Temporal Mixed Effects Model Analysis of Massive Spatio-Temporal Data Sets Matthias Katzfuß Advisor: Dr. Noel Cressie Department of Statistics The Ohio State University
More informationMEI Exam Review. June 7, 2002
MEI Exam Review June 7, 2002 1 Final Exam Revision Notes 1.1 Random Rules and Formulas Linear transformations of random variables. f y (Y ) = f x (X) dx. dg Inverse Proof. (AB)(AB) 1 = I. (B 1 A 1 )(AB)(AB)
More informationarxiv: v4 [stat.me] 14 Sep 2015
Does non-stationary spatial data always require non-stationary random fields? Geir-Arne Fuglstad 1, Daniel Simpson 1, Finn Lindgren 2, and Håvard Rue 1 1 Department of Mathematical Sciences, NTNU, Norway
More informationBiostatistics 533 Classical Theory of Linear Models Spring 2007 Final Exam. Please choose ONE of the following options.
1 Biostatistics 533 Classical Theory of Linear Models Spring 2007 Final Exam Name: KEY Problems do not have equal value and some problems will take more time than others. Spend your time wisely. You do
More informationBayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes
Bayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes Andrew O. Finley Department of Forestry & Department of Geography, Michigan State University, Lansing
More informationarxiv: v1 [stat.me] 27 Jul 2016
Spatial Models for Field Trials María Xosé Rodríguez - Álvarez1,2, Martin P. Boer 3, Fred A. van Eeuwijk 3, Paul H. C. Eilers 4 arxiv:167.8255v1 [stat.me] 27 Jul 216 1 Department of Statistics and Operations
More informationMultiple Predictor Variables: ANOVA
What if you manipulate two factors? Multiple Predictor Variables: ANOVA Block 1 Block 2 Block 3 Block 4 A B C D B C D A C D A B D A B C Randomized Controlled Blocked Design: Design where each treatment
More informationLinear Regression. In this problem sheet, we consider the problem of linear regression with p predictors and one intercept,
Linear Regression In this problem sheet, we consider the problem of linear regression with p predictors and one intercept, y = Xβ + ɛ, where y t = (y 1,..., y n ) is the column vector of target values,
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 informationThe influence of mid-latitude cyclones on European background surface ozone
1 The influence of mid-latitude cyclones on European background surface ozone K. Emma Knowland (USRA/GESTAR NASA/GMAO) In collaboration with: Ruth Doherty (U. Edinburgh) Kevin Hodges (U. Reading) Lesley
More informationGeneral Functional Concurrent Model
General Functional Concurrent Model Janet S. Kim Arnab Maity Ana-Maria Staicu June 17, 2016 Abstract We propose a flexible regression model to study the association between a functional response and multiple
More informationSTAT 360-Linear Models
STAT 360-Linear Models Instructor: Yogendra P. Chaubey Sample Test Questions Fall 004 Note: The following questions are from previous tests and exams. The final exam will be for three hours and will contain
More informationProperties of Matrices and Operations on Matrices
Properties of Matrices and Operations on Matrices A common data structure for statistical analysis is a rectangular array or matris. Rows represent individual observational units, or just observations,
More informationRestricted Maximum Likelihood in Linear Regression and Linear Mixed-Effects Model
Restricted Maximum Likelihood in Linear Regression and Linear Mixed-Effects Model Xiuming Zhang zhangxiuming@u.nus.edu A*STAR-NUS Clinical Imaging Research Center October, 015 Summary This report derives
More informationBayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes
Bayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes Andrew O. Finley 1 and Sudipto Banerjee 2 1 Department of Forestry & Department of Geography, Michigan
More informationAnalysing geoadditive regression data: a mixed model approach
Analysing geoadditive regression data: a mixed model approach Institut für Statistik, Ludwig-Maximilians-Universität München Joint work with Ludwig Fahrmeir & Stefan Lang 25.11.2005 Spatio-temporal regression
More informationLasso Maximum Likelihood Estimation of Parametric Models with Singular Information Matrices
Article Lasso Maximum Likelihood Estimation of Parametric Models with Singular Information Matrices Fei Jin 1,2 and Lung-fei Lee 3, * 1 School of Economics, Shanghai University of Finance and Economics,
More informationAn Introduction to Functional Data Analysis
An Introduction to Functional Data Analysis Chongzhi Di Fred Hutchinson Cancer Research Center cdi@fredhutch.org Biotat 578A: Special Topics in (Genetic) Epidemiology November 10, 2015 Textbook Ramsay
More informationGaussian processes for spatial modelling in environmental health: parameterizing for flexibility vs. computational efficiency
Gaussian processes for spatial modelling in environmental health: parameterizing for flexibility vs. computational efficiency Chris Paciorek March 11, 2005 Department of Biostatistics Harvard School of
More informationGeneralized Functional Concurrent Model
Biometrics, 1 25 2015 Generalized Functional Concurrent Model Janet S. Kim, Arnab Maity, and Ana-Maria Staicu Department of Statistics, North Carolina State University, Raleigh, North Carolina 27695, U.S.A.
More informationHigh-dimensional Statistics
High-dimensional Statistics Pradeep Ravikumar UT Austin Outline 1. High Dimensional Data : Large p, small n 2. Sparsity 3. Group Sparsity 4. Low Rank 1 Curse of Dimensionality Statistical Learning: Given
More informationOverview. Multilevel Models in R: Present & Future. What are multilevel models? What are multilevel models? (cont d)
Overview Multilevel Models in R: Present & Future What are multilevel models? Linear mixed models Symbolic and numeric phases Generalized and nonlinear mixed models Examples Douglas M. Bates U of Wisconsin
More informationMultiple Predictor Variables: ANOVA
Multiple Predictor Variables: ANOVA 1/32 Linear Models with Many Predictors Multiple regression has many predictors BUT - so did 1-way ANOVA if treatments had 2 levels What if there are multiple treatment
More informationSTAT 135 Lab 9 Multiple Testing, One-Way ANOVA and Kruskal-Wallis
STAT 135 Lab 9 Multiple Testing, One-Way ANOVA and Kruskal-Wallis Rebecca Barter April 6, 2015 Multiple Testing Multiple Testing Recall that when we were doing two sample t-tests, we were testing the equality
More informationGeneralized Additive Models
Generalized Additive Models The Model The GLM is: g( µ) = ß 0 + ß 1 x 1 + ß 2 x 2 +... + ß k x k The generalization to the GAM is: g(µ) = ß 0 + f 1 (x 1 ) + f 2 (x 2 ) +... + f k (x k ) where the functions
More informationSpacetime models in R-INLA. Elias T. Krainski
Spacetime models in R-INLA Elias T. Krainski 2 Outline Separable space-time models Infant mortality in Paraná PM-10 concentration in Piemonte, Italy 3 Multivariate dynamic regression model y t : n observations
More informationSmooth additive models for large datasets
Smooth additive models for large datasets Simon Wood, University of Bristol 1, U.K. in collaboration with Zheyuan Li, Gavin Shaddick, Nicole Augustin University of Bath Funded by EPSRC and the University
More informationPenalized Balanced Sampling. Jay Breidt
Penalized Balanced Sampling Jay Breidt Colorado State University Joint work with Guillaume Chauvet (ENSAI) February 4, 2010 1 / 44 Linear Mixed Models Let U = {1, 2,...,N}. Consider linear mixed models
More informationModeling Multi-Way Functional Data With Weak Separability
Modeling Multi-Way Functional Data With Weak Separability Kehui Chen Department of Statistics University of Pittsburgh, USA @CMStatistics, Seville, Spain December 09, 2016 Outline Introduction. Multi-way
More informationMa 3/103: Lecture 25 Linear Regression II: Hypothesis Testing and ANOVA
Ma 3/103: Lecture 25 Linear Regression II: Hypothesis Testing and ANOVA March 6, 2017 KC Border Linear Regression II March 6, 2017 1 / 44 1 OLS estimator 2 Restricted regression 3 Errors in variables 4
More informationTranscendence in Positive Characteristic
Transcendence in Positive Characteristic Galois Group Examples and Applications W. Dale Brownawell Matthew Papanikolas Penn State University Texas A&M University Arizona Winter School 2008 March 18, 2008
More informationFACTORIAL DESIGNS and NESTED DESIGNS
Experimental Design and Statistical Methods Workshop FACTORIAL DESIGNS and NESTED DESIGNS Jesús Piedrafita Arilla jesus.piedrafita@uab.cat Departament de Ciència Animal i dels Aliments Items Factorial
More informationECON 450 Development Economics
ECON 450 Development Economics Statistics Background University of Illinois at Urbana-Champaign Summer 2017 Outline 1 Introduction 2 3 4 5 Introduction Regression analysis is one of the most important
More informationBasis Penalty Smoothers. Simon Wood Mathematical Sciences, University of Bath, U.K.
Basis Penalty Smoothers Simon Wood Mathematical Sciences, University of Bath, U.K. Estimating functions It is sometimes useful to estimate smooth functions from data, without being too precise about the
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 informationForecasting in Hierarchical Models
Forecasting in Hierarchical Models Lucy Morgan Supervisor: Nikolaos Kourentzes 20 th February 2015 Introduction Forecasting is the process of making statements about events whose actual outcomes (typically)
More information20.1. Balanced One-Way Classification Cell means parametrization: ε 1. ε I. + ˆɛ 2 ij =
20. ONE-WAY ANALYSIS OF VARIANCE 1 20.1. Balanced One-Way Classification Cell means parametrization: Y ij = µ i + ε ij, i = 1,..., I; j = 1,..., J, ε ij N(0, σ 2 ), In matrix form, Y = Xβ + ε, or 1 Y J
More informationSTAT 525 Fall Final exam. Tuesday December 14, 2010
STAT 525 Fall 2010 Final exam Tuesday December 14, 2010 Time: 2 hours Name (please print): Show all your work and calculations. Partial credit will be given for work that is partially correct. Points will
More informationOutline. Statistical inference for linear mixed models. One-way ANOVA in matrix-vector form
Outline Statistical inference for linear mixed models Rasmus Waagepetersen Department of Mathematics Aalborg University Denmark general form of linear mixed models examples of analyses using linear mixed
More informationBNAD 276 Lecture 10 Simple Linear Regression Model
1 / 27 BNAD 276 Lecture 10 Simple Linear Regression Model Phuong Ho May 30, 2017 2 / 27 Outline 1 Introduction 2 3 / 27 Outline 1 Introduction 2 4 / 27 Simple Linear Regression Model Managerial decisions
More informationStat 315c: Transposable Data Rasch model and friends
Stat 315c: Transposable Data Rasch model and friends Art B. Owen Stanford Statistics Art B. Owen (Stanford Statistics) Rasch and friends 1 / 14 Categorical data analysis Anova has a problem with too much
More informationChapter 1. Linear Regression with One Predictor Variable
Chapter 1. Linear Regression with One Predictor Variable 1.1 Statistical Relation Between Two Variables To motivate statistical relationships, let us consider a mathematical relation between two mathematical
More informationSimulation Based Aircraft Trajectories Design for Air Traffic Management
Simulation Based Aircraft Trajectories Design for Air Traffic Management N. Kantas 1 A. Lecchini-Visintini 2 J. M. Maciejowski 1 1 Cambridge University Engineering Dept., UK 2 Dept. of Engineering, University
More informationData Integration Model for Air Quality: A Hierarchical Approach to the Global Estimation of Exposures to Ambient Air Pollution
Data Integration Model for Air Quality: A Hierarchical Approach to the Global Estimation of Exposures to Ambient Air Pollution Matthew Thomas 9 th January 07 / 0 OUTLINE Introduction Previous methods for
More informationDummy Variables. Susan Thomas IGIDR, Bombay. 24 November, 2008
IGIDR, Bombay 24 November, 2008 The problem of structural change Model: Y i = β 0 + β 1 X 1i + ɛ i Structural change, type 1: change in parameters in time. Y i = α 1 + β 1 X i + e 1i for period 1 Y i =
More information