ECAS Summer Course. Quantile Regression for Longitudinal Data. Roger Koenker University of Illinois at Urbana-Champaign
|
|
- Daniel Copeland
- 5 years ago
- Views:
Transcription
1 ECAS Summer Course 1 Quantile Regression for Longitudinal Data Roger Koenker University of Illinois at Urbana-Champaign La Roche-en-Ardennes: September 2005 Part I: Penalty Methods for Random Effects Part II: Estimation of Reference Growth Charts Part II is based on joint work with: Ying Wei, Columbia, Biostatistics Anneli Pere, Children s Hospital, Helsinki and Xuming He, UIUC.
2 Penalty Methods and Random Effects 2 Classical Models with Independent Random Effects Additive random effects have nice conditional mean functions and nasty conditional quantile functions. Transformations to annihilate fixed effects fail for quantile regression models
3 Penalty Methods and Random Effects 2 Classical Models with Independent Random Effects Additive random effects have nice conditional mean functions and nasty conditional quantile functions. Transformations to annihilate fixed effects fail for quantile regression models Classical Models with Fixed Effects Fixed effect models are computationally tractable Penalty methods for fixed effect models offer an alternative to classical random effects estimation
4 A Classical Result 3 Model: or, y ij = x ijβ + α i + u ij j = 1,...m i, i = 1,..., n, y = Xβ + Zα + u. If u and α are independent Gaussian vectors with u N (0, R) and α N (0, Q), the optimal weighted least squares estimator is, ˆβ = (X (R + ZQZ ) 1 X) 1 X (R + ZQZ ) 1 y.
5 A Classical Result 3 Model: or, y ij = x ijβ + α i + u ij j = 1,...m i, i = 1,..., n, y = Xβ + Zα + u. If u and α are independent Gaussian vectors with u N (0, R) and α N (0, Q), the optimal weighted least squares estimator is, ˆβ = (X (R + ZQZ ) 1 X) 1 X (R + ZQZ ) 1 y. Theorem ˆβ = argmin (α,β) y Xβ Zα 2 R 1 + α 2 Q 1, where x 2 A = x Ax.
6 Commentary 4 History Henderson (1950) on dairy cows, Goldberger (1962) on BLUP, Lindley and Smith (1972) on Bayesian hierarchical models, and Robinson (1991) for a survey.
7 Commentary 4 History Henderson (1950) on dairy cows, Goldberger (1962) on BLUP, Lindley and Smith (1972) on Bayesian hierarchical models, and Robinson (1991) for a survey. Interpretation Maximum likelihood estimations for the Gaussian random effects models is identical to penalized least squares estimation of the fixed effects model. Random effects is just shrinkage estimation.
8 Commentary 4 History Henderson (1950) on dairy cows, Goldberger (1962) on BLUP, Lindley and Smith (1972) on Bayesian hierarchical models, and Robinson (1991) for a survey. Interpretation Maximum likelihood estimations for the Gaussian random effects models is identical to penalized least squares estimation of the fixed effects model. Random effects is just shrinkage estimation. Ideology Our focus is on design of effective algorithms for data analysis. Models are useful as an interpretative device, but they shouldn t be allowed to become sacred objects.
9 Quantile Regression with Fixed Effects 5 Model: Q yij (τ x ij ) = α i + x ijβ(τ) j = 1,...m i, i = 1,..., n.
10 Quantile Regression with Fixed Effects 5 Model: Q yij (τ x ij ) = α i + x ijβ(τ) j = 1,...m i, i = 1,..., n. Estimator: With (the inevitable) ρ τ (u) = u(τ I(u < 0)), min (α,β) q k=1 n j=1 m i i=1 w k ρ τk (y ij α i x ijβ(τ k ))
11 Quantile Regression with Fixed Effects 5 Model: Q yij (τ x ij ) = α i + x ijβ(τ) j = 1,...m i, i = 1,..., n. Estimator: With (the inevitable) ρ τ (u) = u(τ I(u < 0)), min (α,β) q k=1 n j=1 m i i=1 w k ρ τk (y ij α i x ijβ(τ k )) Implementation: Linear programming formulation Interior point algorithms exploiting sparse linear algebra
12 Quantile Regression with Penalized Fixed Effects 6 Model: (as before) Q yij (τ x ij ) = α i + x ijβ(τ) j = 1,...m i, i = 1,..., n. Estimator: min (α,β) q k=1 Limiting Forms: n j=1 m i i=1 w k ρ τk (y ij α i x ijβ(τ k )) + λ q α i i=1 as λ 0 we obtain fixed effects, as λ we force all α i to 0,
13 L 1 versus L 2 Shrinkage 7 α α λ λ Figure 1: Shrinkage of the fixed effect parameter estimates, ˆα i. The left panel illustrates an example of the l 1 shrinkage effect. The right panel illustrates an example of the l 2 shrinkage effect. There are m = 5 observations on each of n = 50 cross-sectional units.
14 Sparsity 8 row column Figure 2: Sparsity of the design matrix for the penalized fixed effects quantile regression estimator. Grey areas indicate portions of the matrix that are nonzero, white space indicates zeros. Here, n = 50 and m i = 5 for all i. Three quantiles are being simultaneously estimated with the fixed effects restricted to be identical across quantiles.
15 Regularity Conditions 9 A1. The y ij are independent with conditional distribution functions, F ij, given x ij, and differentiable conditional densities, 0 < f ij <, with bounded derivatives f ij, at ξ ij(τ), for j = 1,..., m, i = 1,..., n, A2. Let ω = τ(1 τ) and denote Φ = diag(f ij (ξ ij (τ))). The limiting forms of the following matrices are positive definite: D 0 = lim m n ω m ( ( D 1 = lim m 1 m n A3. max 1 i n x ij < M. 1 j m Z Z X Z/ n Z ΦZ X ΦZ/ n Z X/ n X X/n ), Z ΦX/ n X ΦX/n ).
16 A Theorem 10 Theorem 1. Under conditions A1-3, with n a /m 0 for some a > 0, the ˆδ 1 component of the minimizer, ˆδ, converges in distribution to a Gaussian random vector with mean zero and covariance matrix given by the lower p by p block of the matrix D1 1 D 0D1 1.
17 More Regularity Conditions 11 B1. The y ij are independent with conditional distribution functions, F ij, given x ij, and differentiable conditional densities, 0 < f ij <, with bounded derivatives f ij, at ξ ij(τ), for j = 1,..., m, i = 1,..., n, B2. Let Ω denote the q by q matrix with typical element τ k τ l τ k τ l and Φ k = diag(f ij (ξ ij (τ k ))). The limiting forms of the following matrices are positive definite: ( D 0 = lim m 1 m n D 1 = lim m 1 m n w ΩwZ Z W Ωw X Z/ n w ΩW Z X/ n W ΩW X X/n ), wk Z Φ k Z w 1 Z Φ 1 X/ n w q Z Φ q X/ n w 1 X Φ 1 Z/ n w 1 X Φ 1 X/n w q X Φ q Z/ n 0 w q X Φ q X/n.
18 B3. max 1 i n x ij < M. 1 j m 12
19 Another Theorem 13 Theorem 2. Under conditions B1-3, provided that λ m / m λ 0, and n a /m 0 for some a > 0, the first component ˆδ 1 minimizing V mn has the same limiting distribution as the as the first component of the minimizer of, V 0 (δ) = δ B δ D 1 δ + λ 0 δ s where B is a zero mean Gaussian vector with covariance matrix D 0, and s = (s 0 0 pq) and s 0 = (sgn(α i )).
20 Monte-Carlo Models 14 We consider two very simple models: and where y ij = α i + x ij β + u ij (1) y ij = α i + x ij β + (1 + x ij γ)u ij. (2) x ij = δ i + v ij (3) with γ i and v ij independent and identically distributed over i and i, j respectively, then the interclass correlation coefficient, ρ x = σ 2 δ/(σ 2 δ + σ 2 u),
21 Monte Carlo: Location Shift Model 15 LS PLS LSFE QR PQR QRFE N Bias RMSE t 3 Bias RMSE χ 2 3 Bias RMSE
22 Monte Carlo: Location-Scale Shift Model 16 LS PLS LSFE QR PQR QRFE N Bias RMSE t 3 Bias RMSE χ 2 3 Bias RMSE
23 Conclusions 17 Classical random effects in Gaussian models is just penalized least squares estimation in which fixed effects are shrunken toward their prior mean of zero. Other penalties may be useful in other settings.
24 Conclusions 17 Classical random effects in Gaussian models is just penalized least squares estimation in which fixed effects are shrunken toward their prior mean of zero. Other penalties may be useful in other settings. L 1 penalty methods have attractive properties in a variety of model selection, image analysis and function estimation problems. Interior point algorithms exploiting inherently sparse linear algebra facilate efficient computations.
25 Conclusions 17 Classical random effects in Gaussian models is just penalized least squares estimation in which fixed effects are shrunken toward their prior mean of zero. Other penalties may be useful in other settings. L 1 penalty methods have attractive properties in a variety of model selection, image analysis and function estimation problems. Interior point algorithms exploiting inherently sparse linear algebra facilate efficient computations. There is more to longitudinal data analysis than the estimation of location shift effects. Many interesting questions can only be addressed with broader models that incorporate distributional effects for treatments and other covariates.
26 Quantile Regression Methods for Reference Growth Charts 18 Roger Koenker University of Illinois at Urbana-Champaign Based on joint work with: Ying Wei, Columbia, Biostatistics Anneli Pere, Children s Hospital, Helsinki Xuming He, UIUC.
27 Quetelet s (1871) Growth Chart 19
28 Penalized Maximum Likelihood Estimation 20 References: Cole (1988), Cole and Green (1992), and Carey(2002) Data: {Y i (t i,j ) : j = 1,..., J i, i = 1,..., n.} Model: Z(t) = (Y (t)/µ(t))λ(t) 1 λ(t)σ(t) N (0, 1) Estimation: max l(λ, µ, σ) ν λ (λ (t)) 2 dt ν µ (µ (t)) 2 dt ν σ (σ (t)) 2 dt, l(λ, µ, σ) = n [λ(t i ) log(y (t i )/µ(t i )) log σ(t i ) 1 2 Z2 (t i )], i=1
29 Quantiles as Argmins 21 The τ th quantile of a random variable Y having distribution function F is: α(τ) = argmin ρ τ (y α)df (y) where ρ τ (u) = u (τ I(u < 0)).
30 Quantiles as Argmins 21 The τ th quantile of a random variable Y having distribution function F is: α(τ) = argmin ρ τ (y α)df (y) where ρ τ (u) = u (τ I(u < 0)). The τth sample quantile is thus: ˆα(τ) = argmin ρ τ (y α)df n (y) n = argmin n 1 ρ τ (y i α) i=1
31 Quantile Regression 22 The τth conditional quantile function of Y X = x is g(τ x) = argmin g G ρ τ (y g(x))df A natural estimator of g(τ x) is ĝ(τ x) = argmin g G n i=1 ρ τ (y i g(x i )) with G chosen as a finite dimensional linear space, g(x) = p ϕ j (x)β j. j=1
32 Choice of Basis 23 There are many possible choices for the basis expansion {ϕ j }. We opt for the (very conventional) cubic B-spline functions: Age In R these quantile regression models can be estimated with the command. fit <- rq(y bs(x,knots=knots),tau = 1:9/10) Similar functionality in SAS is coming real soon now.
33 Data 24 Longitudinal measurements on height for 2514 Finnish children,
34 Data 24 Longitudinal measurements on height for 2514 Finnish children, 1143 boys, 1162 girls all healthy, full-term, singleton births,
35 Data 24 Longitudinal measurements on height for 2514 Finnish children, 1143 boys, 1162 girls all healthy, full-term, singleton births, About 20 measurements per child,
36 Data 24 Longitudinal measurements on height for 2514 Finnish children, 1143 boys, 1162 girls all healthy, full-term, singleton births, About 20 measurements per child, Two cohorts: 1096 born between , 1209 born between
37 Data 24 Longitudinal measurements on height for 2514 Finnish children, 1143 boys, 1162 girls all healthy, full-term, singleton births, About 20 measurements per child, Two cohorts: 1096 born between , 1209 born between Sample constitutes 0.5 percent of Finns born in these periods.
38 Unconditional Reference Quantiles Boys Years Box Cox Parameter Functions λ(t) Height (cm) µ(t) LMS edf = (7,10,7) LMS edf = (22,25,22) QR edf = 16 σ(t) Age (years) Age (years)
39 26
40 Unconditional Reference Quantiles Boys 2 18 Years Box Cox Parameter Functions λ(t) µ(t) Height (cm) LMS edf = (7,10,7) σ(t) LMS edf = (22,25,22) QR edf = Age (years) Age (years)
41 28
42 Unconditional Reference Quantiles Girls Years Box Cox Parameter Functions λ(t) µ(t) 90 Height (cm) LMS edf = (7,10,7) σ(t) LMS edf = (22,25,22) QR edf = Age (years) Age (years)
43 30
44 31 Unconditional Reference Quantiles Girls 2 18 Years Box Cox Parameter Functions λ(t) Height (cm) µ(t) LMS edf = (7,10,7) σ(t) LMS edf = (22,25,22) QR edf = Age (years) Age (years)
45 32
46 33 Growth Velocity Curves Boys years Girls years τ Boys 2 18 years Girls 2 18 years LMS (7,10,7) LMS (22,25,22) QR Velocity (cm/year) Age (years)
47 34
48 35 Estimated Age Specific Density Functions Age = 0.5 Age = 1 Age = 1.5 N= 235 N= 204 N= Boys Data QR LMS N= 215 N= 191 N= Girls Height (cm)
49 36
50 37 Height Density at Age Boys Pooled Born in 1960 Born in Girls Height(cm)
51 Conditioning on Prior Growth 38 It is often important to condition not only on age, but also on prior growth and possibly on other covariates. Autoregressive models are natural, but complicated due to the irregular spacing of typical longitudinal measurements. Data: {Y i (t i,j ) : j = 1,..., J i, i = 1,..., n.} Model: Q Yi (t i,j )(τ t i,j, Y i (t i,j 1 ), x i ) = g τ (t i,j ) + [α(τ) + β(τ)(t i,j t i,j 1 )]Y i (t i,j 1 ) + x i γ(τ).
52 AR Components of the Infant Conditional Growth Model 39 τ Boys Girls ˆα(τ) ˆβ(τ) ˆγ(τ) ˆα(τ) ˆβ(τ) ˆγ(τ) (0.020) (0.011) (0.011) (0.024) (0.011) (0.020) (0.019) (0.025) (0.029) (0.024) (0.024) (0.007) (0.006) (0.009) (0.009) (0.016) (0.016) (0.007) (0.009) (0.013) (0.017) (0.017) (0.018) (0.022) (0.021) (0.027) (0.027) (0.030) (0.038) (0.007) (0.006) (0.008) (0.012) (0.015) (0.024) (0.010) (0.009) (0.008) (0.009) (0.015) (0.018) (0.027)
53 AR Components of the Childrens Conditional Growth Model 40 τ Boys Girls ˆα(τ) ˆβ(τ) ˆγ(τ) ˆα(τ) ˆβ(τ) ˆγ(τ) (0.010) (0.002) (0.013) (0.012) (0.002) (0.005) (0.006) (0.004) (0.004) (0.009) (0.014) (0.001) (0.001) (0.001) (0.001) (0.001) (0.002) (0.007) (0.006) (0.004) (0.006) (0.009) (0.015) (0.006) (0.005) (0.007) (0.007) (0.008) (0.013) (0.001) (0.001) (0.001) (0.001) (0.001) (0.001) (0.015) (0.007) (0.006) (0.006) (0.006) (0.012) (0.018)
54 Transformation to Normality 41 A presumed advantage of univariate (age-specific) transformation to normality is that once observations are transformed to univariate Z-scores they are automatically prepared to longitudinal autoregression: Z t = α 0 + α 1 Z t 1 + U t Premise: Marginal Normality Joint Normality
55 Transformation to Normality 41 A presumed advantage of univariate (age-specific) transformation to normality is that once observations are transformed to univariate Z-scores they are automatically prepared to longitudinal autoregression: Z t = α 0 + α 1 Z t 1 + U t Premise: Marginal Normality Joint Normality Of course we know it isn t true, but we also think we know that counterexamples are pathological, and don t occur in nature.
56 Boys age 6.5 Theoretical Quantiles Sample Quantiles Boys age 6 Theoretical Quantiles Sample Quantiles Boys AR(1) residuals Theoretical Quantiles Sample Quantiles Girls age 6.5 Theoretical Quantiles Sample Quantiles Girls age 6 Theoretical Quantiles Sample Quantiles Girls AR(1) residuals Theoretical Quantiles Sample Quantiles
57 43
58 44 Height (cm) Subject 89 Subject Age (years)
59 45 1 Subject 89 Subject Quantile Unconditional Conditional Observed Height (cm) Height (cm) 0
60 Conclusions 46 Nonparametric quantile regression using B-splines offers a reasonable alternative to parametric methods for constructing reference growth charts.
61 Conclusions 46 Nonparametric quantile regression using B-splines offers a reasonable alternative to parametric methods for constructing reference growth charts. The flexibility of quantile regression methods exposes features of the data that are not easily observable with conventional parametric methods. Even for height data.
62 Conclusions 46 Nonparametric quantile regression using B-splines offers a reasonable alternative to parametric methods for constructing reference growth charts. The flexibility of quantile regression methods exposes features of the data that are not easily observable with conventional parametric methods. Even for height data. Longitudinal data can be easily accomodated into the quantile regression framework by adding covariates, including the use of autoregressive effects for unequally spaced measurements.
Quantile Regression Methods for Reference Growth Charts
Quantile Regression Methods for Reference Growth Charts 1 Roger Koenker University of Illinois at Urbana-Champaign ASA Workshop on Nonparametric Statistics Texas A&M, January 15, 2005 Based on joint work
More informationQuantile Regression for Longitudinal Data
Quantile Regression for Longitudinal Data Roger Koenker CEMMAP and University of Illinois, Urbana-Champaign LSE: 17 May 2011 Roger Koenker (CEMMAP & UIUC) Quantile Regression for Longitudinal Data LSE:
More informationNonparametric Quantile Regression
Nonparametric Quantile Regression Roger Koenker CEMMAP and University of Illinois, Urbana-Champaign LSE: 17 May 2011 Roger Koenker (CEMMAP & UIUC) Nonparametric QR LSE: 17.5.2010 1 / 25 In the Beginning,...
More informationQuantile Regression for Panel/Longitudinal Data
Quantile Regression for Panel/Longitudinal Data Roger Koenker University of Illinois, Urbana-Champaign University of Minho 12-14 June 2017 y it 0 5 10 15 20 25 i = 1 i = 2 i = 3 0 2 4 6 8 Roger Koenker
More informationQuantile Regression for Dynamic Panel Data
Quantile Regression for Dynamic Panel Data Antonio Galvao 1 1 Department of Economics University of Illinois NASM Econometric Society 2008 June 22nd 2008 Panel Data Panel data allows the possibility of
More informationECAS Summer Course. Fundamentals of Quantile Regression Roger Koenker University of Illinois at Urbana-Champaign
ECAS Summer Course 1 Fundamentals of Quantile Regression Roger Koenker University of Illinois at Urbana-Champaign La Roche-en-Ardennes: September, 2005 An Outline 2 Beyond Average Treatment Effects Computation
More informationA DYNAMIC QUANTILE REGRESSION TRANSFORMATION MODEL FOR LONGITUDINAL DATA
Statistica Sinica 19 (2009), 1137-1153 A DYNAMIC QUANTILE REGRESSION TRANSFORMATION MODEL FOR LONGITUDINAL DATA Yunming Mu and Ying Wei Portland State University and Columbia University Abstract: This
More informationQuantile Regression: Inference
Quantile Regression: Inference Roger Koenker University of Illinois, Urbana-Champaign Aarhus: 21 June 2010 Roger Koenker (UIUC) Introduction Aarhus: 21.6.2010 1 / 28 Inference for Quantile Regression Asymptotics
More informationQuantile Regression: Inference
Quantile Regression: Inference Roger Koenker University of Illinois, Urbana-Champaign 5th RMetrics Workshop, Meielisalp: 28 June 2011 Roger Koenker (UIUC) Introduction Meielisalp: 28.6.2011 1 / 29 Inference
More informationRegression based principal component analysis for sparse functional data with applications to screening pubertal growth paths.
Regression based principal component analysis for sparse functional data with applications to screening pubertal growth paths Wenfei Zhang Submitted in partial fulfillment of the requirements for the degree
More informationConfidence Intervals for Low-dimensional Parameters with High-dimensional Data
Confidence Intervals for Low-dimensional Parameters with High-dimensional Data Cun-Hui Zhang and Stephanie S. Zhang Rutgers University and Columbia University September 14, 2012 Outline Introduction Methodology
More informationPenalty Methods for Bivariate Smoothing and Chicago Land Values
Penalty Methods for Bivariate Smoothing and Chicago Land Values Roger Koenker University of Illinois, Urbana-Champaign Ivan Mizera University of Alberta, Edmonton Northwestern University: October 2001
More informationMixed Effects Multivariate Adaptive Splines Model for the Analysis of Longitudinal and Growth Curve Data. Heping Zhang.
Mixed Effects Multivariate Adaptive Splines Model for the Analysis of Longitudinal and Growth Curve Data Heping Zhang Yale University NIA Workshop on Statistical Methods for Longitudinal Data on Aging
More informationBeyond the Average Man: The Art of Unlikelihood
Beyond the Average Man: The Art of Unlikelihood Roger Koenker University of Illinois, Urbana-Champaign FEMES/Taipei: July, 2007 Roger Koenker (U. of Illinois) Beyond the Average Man FEMES/Taipei: July,
More informationSupporting Information for Estimating restricted mean. treatment effects with stacked survival models
Supporting Information for Estimating restricted mean treatment effects with stacked survival models Andrew Wey, David Vock, John Connett, and Kyle Rudser Section 1 presents several extensions to the simulation
More informationWEIGHTED QUANTILE REGRESSION THEORY AND ITS APPLICATION. Abstract
Journal of Data Science,17(1). P. 145-160,2019 DOI:10.6339/JDS.201901_17(1).0007 WEIGHTED QUANTILE REGRESSION THEORY AND ITS APPLICATION Wei Xiong *, Maozai Tian 2 1 School of Statistics, University of
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 informationQuantile Regression Computation: From Outside, Inside and the Proximal
Quantile Regression Computation: From Outside, Inside and the Proximal Roger Koenker University of Illinois, Urbana-Champaign University of Minho 12-14 June 2017 1.0 0.5 0.0 0.5 1.0 1.0 0.5 0.0 0.5 1.0
More informationMotivational Example
Motivational Example Data: Observational longitudinal study of obesity from birth to adulthood. Overall Goal: Build age-, gender-, height-specific growth charts (under 3 year) to diagnose growth abnomalities.
More informationModeling Real Estate Data using Quantile Regression
Modeling Real Estate Data using Semiparametric Quantile Regression Department of Statistics University of Innsbruck September 9th, 2011 Overview 1 Application: 2 3 4 Hedonic regression data for house prices
More informationStat 579: Generalized Linear Models and Extensions
Stat 579: Generalized Linear Models and Extensions Linear Mixed Models for Longitudinal Data Yan Lu April, 2018, week 15 1 / 38 Data structure t1 t2 tn i 1st subject y 11 y 12 y 1n1 Experimental 2nd subject
More informationA Short Introduction to the Lasso Methodology
A Short Introduction to the Lasso Methodology Michael Gutmann sites.google.com/site/michaelgutmann University of Helsinki Aalto University Helsinki Institute for Information Technology March 9, 2016 Michael
More informationAsymptotic standard errors of MLE
Asymptotic standard errors of MLE Suppose, in the previous example of Carbon and Nitrogen in soil data, that we get the parameter estimates For maximum likelihood estimation, we can use Hessian matrix
More informationReference growth charts for Saudi Arabian children and adolescents
Chapter 7 Reference growth charts for Saudi Arabian children and adolescents P. J. Foster and T. Kecojević Abstract The purpose of this study is to provide Saudi Arabian population reference growth standards
More informationQuantile Processes for Semi and Nonparametric Regression
Quantile Processes for Semi and Nonparametric Regression Shih-Kang Chao Department of Statistics Purdue University IMS-APRM 2016 A joint work with Stanislav Volgushev and Guang Cheng Quantile Response
More informationLecture 22 Survival Analysis: An Introduction
University of Illinois Department of Economics Spring 2017 Econ 574 Roger Koenker Lecture 22 Survival Analysis: An Introduction There is considerable interest among economists in models of durations, which
More informationINFERENCE ON QUANTILE REGRESSION FOR HETEROSCEDASTIC MIXED MODELS
Statistica Sinica 19 (2009), 1247-1261 INFERENCE ON QUANTILE REGRESSION FOR HETEROSCEDASTIC MIXED MODELS Huixia Judy Wang North Carolina State University Abstract: This paper develops two weighted quantile
More informationTechnical Vignette 5: Understanding intrinsic Gaussian Markov random field spatial models, including intrinsic conditional autoregressive models
Technical Vignette 5: Understanding intrinsic Gaussian Markov random field spatial models, including intrinsic conditional autoregressive models Christopher Paciorek, Department of Statistics, University
More informationSpring 2017 Econ 574 Roger Koenker. Lecture 14 GEE-GMM
University of Illinois Department of Economics Spring 2017 Econ 574 Roger Koenker Lecture 14 GEE-GMM Throughout the course we have emphasized methods of estimation and inference based on the principle
More informationCourse topics (tentative) The role of random effects
Course topics (tentative) random effects linear mixed models analysis of variance frequentist likelihood-based inference (MLE and REML) prediction Bayesian inference The role of random effects Rasmus Waagepetersen
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 informationReference Growth Charts for Saudi Arabian Children and Adolescents
Reference Growth Charts for Saudi Arabian Children and Adolescents P. J. Foster & T. Kecojević First version: 3 September 2009 Research Report No. 12, 2009, Probability and Statistics Group School of Mathematics,
More informationVariable Selection for Highly Correlated Predictors
Variable Selection for Highly Correlated Predictors Fei Xue and Annie Qu arxiv:1709.04840v1 [stat.me] 14 Sep 2017 Abstract Penalty-based variable selection methods are powerful in selecting relevant covariates
More informationA Blockwise Descent Algorithm for Group-penalized Multiresponse and Multinomial Regression
A Blockwise Descent Algorithm for Group-penalized Multiresponse and Multinomial Regression Noah Simon Jerome Friedman Trevor Hastie November 5, 013 Abstract In this paper we purpose a blockwise descent
More informationGeneralized Elastic Net Regression
Abstract Generalized Elastic Net Regression Geoffroy MOURET Jean-Jules BRAULT Vahid PARTOVINIA This work presents a variation of the elastic net penalization method. We propose applying a combined l 1
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 informationQuantile Regression Computation: From the Inside and the Outside
Quantile Regression Computation: From the Inside and the Outside Roger Koenker University of Illinois, Urbana-Champaign UI Statistics: 24 May 2008 Slides will be available on my webpages: click Software.
More informationModelling the Covariance
Modelling the Covariance Jamie Monogan Washington University in St Louis February 9, 2010 Jamie Monogan (WUStL) Modelling the Covariance February 9, 2010 1 / 13 Objectives By the end of this meeting, participants
More informationNonparametric Regression. Badr Missaoui
Badr Missaoui Outline Kernel and local polynomial regression. Penalized regression. We are given n pairs of observations (X 1, Y 1 ),...,(X n, Y n ) where Y i = r(x i ) + ε i, i = 1,..., n and r(x) = E(Y
More informationBuilding a Prognostic Biomarker
Building a Prognostic Biomarker Noah Simon and Richard Simon July 2016 1 / 44 Prognostic Biomarker for a Continuous Measure On each of n patients measure y i - single continuous outcome (eg. blood pressure,
More informationSample Size and Power Considerations for Longitudinal Studies
Sample Size and Power Considerations for Longitudinal Studies Outline Quantities required to determine the sample size in longitudinal studies Review of type I error, type II error, and power For continuous
More informationSingle Index Quantile Regression for Heteroscedastic Data
Single Index Quantile Regression for Heteroscedastic Data E. Christou M. G. Akritas Department of Statistics The Pennsylvania State University SMAC, November 6, 2015 E. Christou, M. G. Akritas (PSU) SIQR
More informationOutline. Mixed models in R using the lme4 package Part 3: Longitudinal data. Sleep deprivation data. Simple longitudinal data
Outline Mixed models in R using the lme4 package Part 3: Longitudinal data Douglas Bates Longitudinal data: sleepstudy A model with random effects for intercept and slope University of Wisconsin - Madison
More informationarxiv: v2 [stat.me] 4 Jun 2016
Variable Selection for Additive Partial Linear Quantile Regression with Missing Covariates 1 Variable Selection for Additive Partial Linear Quantile Regression with Missing Covariates Ben Sherwood arxiv:1510.00094v2
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 informationAn Encompassing Test for Non-Nested Quantile Regression Models
An Encompassing Test for Non-Nested Quantile Regression Models Chung-Ming Kuan Department of Finance National Taiwan University Hsin-Yi Lin Department of Economics National Chengchi University Abstract
More informationHierarchical Modeling for Univariate Spatial Data
Hierarchical Modeling for Univariate Spatial Data Geography 890, Hierarchical Bayesian Models for Environmental Spatial Data Analysis February 15, 2011 1 Spatial Domain 2 Geography 890 Spatial Domain This
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 informationQuantile regression for longitudinal data using the asymmetric Laplace distribution
Biostatistics (2007), 8, 1, pp. 140 154 doi:10.1093/biostatistics/kxj039 Advance Access publication on April 24, 2006 Quantile regression for longitudinal data using the asymmetric Laplace distribution
More informationGaussian Process Regression
Gaussian Process Regression 4F1 Pattern Recognition, 21 Carl Edward Rasmussen Department of Engineering, University of Cambridge November 11th - 16th, 21 Rasmussen (Engineering, Cambridge) Gaussian Process
More informationQuantile Regression: Past and Prospects
Quantile Regression: Past and Prospects Roger Koenker University of Illinois, Urbana-Champaign Northern Illinois University: 31 October 2014 Roger Koenker (UIUC) Quantile Regression NIU: 31.10.2014 1 /
More informationSupplementary Material to General Functional Concurrent Model
Supplementary Material to General Functional Concurrent Model Janet S. Kim Arnab Maity Ana-Maria Staicu June 17, 2016 This Supplementary Material contains six sections. Appendix A discusses modifications
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 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 informationInversion Base Height. Daggot Pressure Gradient Visibility (miles)
Stanford University June 2, 1998 Bayesian Backtting: 1 Bayesian Backtting Trevor Hastie Stanford University Rob Tibshirani University of Toronto Email: trevor@stat.stanford.edu Ftp: stat.stanford.edu:
More informationRecap. HW due Thursday by 5 pm Next HW coming on Thursday Logistic regression: Pr(G = k X) linear on the logit scale Linear discriminant analysis:
1 / 23 Recap HW due Thursday by 5 pm Next HW coming on Thursday Logistic regression: Pr(G = k X) linear on the logit scale Linear discriminant analysis: Pr(G = k X) Pr(X G = k)pr(g = k) Theory: LDA more
More informationBayesian Linear Regression
Bayesian Linear Regression Sudipto Banerjee 1 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A. September 15, 2010 1 Linear regression models: a Bayesian perspective
More informationThe Wishart distribution Scaled Wishart. Wishart Priors. Patrick Breheny. March 28. Patrick Breheny BST 701: Bayesian Modeling in Biostatistics 1/11
Wishart Priors Patrick Breheny March 28 Patrick Breheny BST 701: Bayesian Modeling in Biostatistics 1/11 Introduction When more than two coefficients vary, it becomes difficult to directly model each element
More informationRecovering Copulae from Conditional Quantiles
Wolfgang K. Härdle Chen Huang Alexander Ristig Ladislaus von Bortkiewicz Chair of Statistics C.A.S.E. Center for Applied Statistics and Economics HumboldtUniversität zu Berlin http://lvb.wiwi.hu-berlin.de
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 informationFor more information about how to cite these materials visit
Author(s): Kerby Shedden, Ph.D., 2010 License: Unless otherwise noted, this material is made available under the terms of the Creative Commons Attribution Share Alike 3.0 License: http://creativecommons.org/licenses/by-sa/3.0/
More informationRegression, Ridge Regression, Lasso
Regression, Ridge Regression, Lasso Fabio G. Cozman - fgcozman@usp.br October 2, 2018 A general definition Regression studies the relationship between a response variable Y and covariates X 1,..., X n.
More informationChapter 3. Linear Models for Regression
Chapter 3. Linear Models for Regression Wei Pan Division of Biostatistics, School of Public Health, University of Minnesota, Minneapolis, MN 55455 Email: weip@biostat.umn.edu PubH 7475/8475 c Wei Pan Linear
More informationPenalized Empirical Likelihood and Growing Dimensional General Estimating Equations
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 Biometrika (0000), 00, 0, pp. 1 15 C 0000 Biometrika Trust Printed
More informationFixed Effects, Invariance, and Spatial Variation in Intergenerational Mobility
American Economic Review: Papers & Proceedings 2016, 106(5): 400 404 http://dx.doi.org/10.1257/aer.p20161082 Fixed Effects, Invariance, and Spatial Variation in Intergenerational Mobility By Gary Chamberlain*
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 informationBayesian linear regression
Bayesian linear regression Linear regression is the basis of most statistical modeling. The model is Y i = X T i β + ε i, where Y i is the continuous response X i = (X i1,..., X ip ) T is the corresponding
More informationHierarchical Modelling for Univariate Spatial Data
Hierarchical Modelling for Univariate Spatial Data Sudipto Banerjee 1 and Andrew O. Finley 2 1 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A. 2 Department
More informationQuantile Regression: 40 Years On
Quantile Regression: 40 Years On Roger Koenker Department of Economics, University of Illinois, Urbana, IL 61801, USA; email: rkoenker@uiuc.edu xxxxxx 0000. 00:1 25 Copyright c 0000 by Annual Reviews.
More informationAsymmetric least squares estimation and testing
Asymmetric least squares estimation and testing Whitney Newey and James Powell Princeton University and University of Wisconsin-Madison January 27, 2012 Outline ALS estimators Large sample properties Asymptotic
More informationRegularization Parameter Selection for a Bayesian Multi-Level Group Lasso Regression Model with Application to Imaging Genomics
Regularization Parameter Selection for a Bayesian Multi-Level Group Lasso Regression Model with Application to Imaging Genomics arxiv:1603.08163v1 [stat.ml] 7 Mar 016 Farouk S. Nathoo, Keelin Greenlaw,
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 informationBayesian Inference. Chapter 4: Regression and Hierarchical Models
Bayesian Inference Chapter 4: Regression and Hierarchical Models Conchi Ausín and Mike Wiper Department of Statistics Universidad Carlos III de Madrid Master in Business Administration and Quantitative
More informationLecture 5 Models and methods for recurrent event data
Lecture 5 Models and methods for recurrent event data Recurrent and multiple events are commonly encountered in longitudinal studies. In this chapter we consider ordered recurrent and multiple events.
More informationBayesian Estimation of Input Output Tables for Russia
Bayesian Estimation of Input Output Tables for Russia Oleg Lugovoy (EDF, RANE) Andrey Polbin (RANE) Vladimir Potashnikov (RANE) WIOD Conference April 24, 2012 Groningen Outline Motivation Objectives Bayesian
More informationEconomics 536 Lecture 7. Introduction to Specification Testing in Dynamic Econometric Models
University of Illinois Fall 2016 Department of Economics Roger Koenker Economics 536 Lecture 7 Introduction to Specification Testing in Dynamic Econometric Models In this lecture I want to briefly describe
More informationmboost - Componentwise Boosting for Generalised Regression Models
mboost - Componentwise Boosting for Generalised Regression Models Thomas Kneib & Torsten Hothorn Department of Statistics Ludwig-Maximilians-University Munich 13.8.2008 Boosting in a Nutshell Boosting
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 informationPractical Bayesian Quantile Regression. Keming Yu University of Plymouth, UK
Practical Bayesian Quantile Regression Keming Yu University of Plymouth, UK (kyu@plymouth.ac.uk) A brief summary of some recent work of us (Keming Yu, Rana Moyeed and Julian Stander). Summary We develops
More informationHigh-dimensional Ordinary Least-squares Projection for Screening Variables
1 / 38 High-dimensional Ordinary Least-squares Projection for Screening Variables Chenlei Leng Joint with Xiangyu Wang (Duke) Conference on Nonparametric Statistics for Big Data and Celebration to Honor
More informationWU Weiterbildung. Linear Mixed Models
Linear Mixed Effects Models WU Weiterbildung SLIDE 1 Outline 1 Estimation: ML vs. REML 2 Special Models On Two Levels Mixed ANOVA Or Random ANOVA Random Intercept Model Random Coefficients Model Intercept-and-Slopes-as-Outcomes
More informationQuantile Regression: An Introductory Overview
Quantile Regression: An Introductory Overview Roger Koenker U. of Illinois Xuming He U. of Michigan Lan Wang U. of Minnesota JSM Baltimore: 2 August 2017 Koenker, He and Wang Quantile Regression JSM Baltimore:
More informationQuantile regression 40 years on
Quantile regression 40 years on Roger Koenker The Institute for Fiscal Studies Department of Economics, UCL cemmap working paper CWP36/17 QUANTILE REGRESSION 40 YEARS ON ROGER KOENKER Abstract. Since Quetelet
More informationFuncICA for time series pattern discovery
FuncICA for time series pattern discovery Nishant Mehta and Alexander Gray Georgia Institute of Technology The problem Given a set of inherently continuous time series (e.g. EEG) Find a set of patterns
More informationCALCULATION METHOD FOR NONLINEAR DYNAMIC LEAST-ABSOLUTE DEVIATIONS ESTIMATOR
J. Japan Statist. Soc. Vol. 3 No. 200 39 5 CALCULAION MEHOD FOR NONLINEAR DYNAMIC LEAS-ABSOLUE DEVIAIONS ESIMAOR Kohtaro Hitomi * and Masato Kagihara ** In a nonlinear dynamic model, the consistency and
More informationTesting Error Correction in Panel data
University of Vienna, Dept. of Economics Master in Economics Vienna 2010 The Model (1) Westerlund (2007) consider the following DGP: y it = φ 1i + φ 2i t + z it (1) x it = x it 1 + υ it (2) where the stochastic
More informationIntroduction to Regression
Introduction to Regression David E Jones (slides mostly by Chad M Schafer) June 1, 2016 1 / 102 Outline General Concepts of Regression, Bias-Variance Tradeoff Linear Regression Nonparametric Procedures
More informationSpatial Process Estimates as Smoothers: A Review
Spatial Process Estimates as Smoothers: A Review Soutir Bandyopadhyay 1 Basic Model The observational model considered here has the form Y i = f(x i ) + ɛ i, for 1 i n. (1.1) where Y i is the observed
More informationLecture 6: Methods for high-dimensional problems
Lecture 6: Methods for high-dimensional problems Hector Corrada Bravo and Rafael A. Irizarry March, 2010 In this Section we will discuss methods where data lies on high-dimensional spaces. In particular,
More informationCPSC 540: Machine Learning
CPSC 540: Machine Learning Empirical Bayes, Hierarchical Bayes Mark Schmidt University of British Columbia Winter 2017 Admin Assignment 5: Due April 10. Project description on Piazza. Final details coming
More informationMixed effects modeling with missing data using quantile regression and joint modeling
Graduate Theses and Dissertations Graduate College 2014 Mixed effects modeling with missing data using quantile regression and joint modeling Luke Karsten Fostvedt Iowa State University Follow this and
More informationStatistics 203: Introduction to Regression and Analysis of Variance Penalized models
Statistics 203: Introduction to Regression and Analysis of Variance Penalized models Jonathan Taylor - p. 1/15 Today s class Bias-Variance tradeoff. Penalized regression. Cross-validation. - p. 2/15 Bias-variance
More informationHigh-dimensional regression
High-dimensional regression Advanced Methods for Data Analysis 36-402/36-608) Spring 2014 1 Back to linear regression 1.1 Shortcomings Suppose that we are given outcome measurements y 1,... y n R, and
More informationPart 6: Multivariate Normal and Linear Models
Part 6: Multivariate Normal and Linear Models 1 Multiple measurements Up until now all of our statistical models have been univariate models models for a single measurement on each member of a sample of
More informationLog Gaussian Cox Processes. Chi Group Meeting February 23, 2016
Log Gaussian Cox Processes Chi Group Meeting February 23, 2016 Outline Typical motivating application Introduction to LGCP model Brief overview of inference Applications in my work just getting started
More informationA Penalized spline approach to functional mixed effects model analysis
Biometrics 63, 1?? December 2007 DOI: 10.1111/j.1541-0420.2005.00454.x A Penalized spline approach to functional mixed effects model analysis Huaihou Chen Department of Biostatistics, Mailman School of
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 informationHypothesis Testing Based on the Maximum of Two Statistics from Weighted and Unweighted Estimating Equations
Hypothesis Testing Based on the Maximum of Two Statistics from Weighted and Unweighted Estimating Equations Takeshi Emura and Hisayuki Tsukuma Abstract For testing the regression parameter in multivariate
More informationEcon 2148, fall 2017 Gaussian process priors, reproducing kernel Hilbert spaces, and Splines
Econ 2148, fall 2017 Gaussian process priors, reproducing kernel Hilbert spaces, and Splines Maximilian Kasy Department of Economics, Harvard University 1 / 37 Agenda 6 equivalent representations of the
More informationHierarchical Modelling for Univariate Spatial Data
Spatial omain Hierarchical Modelling for Univariate Spatial ata Sudipto Banerjee 1 and Andrew O. Finley 2 1 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A.
More information