DETECTION AND ATTRIBUTION
|
|
- Francis Park
- 5 years ago
- Views:
Transcription
1 DETECTION AND ATTRIBUTION Richard L. Smith Department of Statistics and Operations Research University of North Carolina Chapel Hill, NC , USA 41st Winter Conference in Statistics Statistics and Climate Storhonga, Sweden March
2 2
3 3
4 4
5 5
6 6
7 7
8 8
9 Basic idea: How can we measure the agreement between an artificially generated climate signal from a climate model, and real data as measured by surface observation stations or satellites? Both observations and model data aggregated into both the observational and model-generated data are aggregated into grid cells, but several thousand of them. Example: Use 5 o latitude and longitude grid cells, implying 36 latitude classifications (from 90 o N to 90 o S), and similarly 72 longitude increments, giving a total of 2,592 grid cells. 9
10 Use l to denote the number of grid cells. One common kind of study is to generate m possible response patterns x k, k = 1,..., m using the climate model, each one an l-dimensional time series corresponding to the response to one particular kind of forcing factor. For example, x 1 may be the response due to the increase in carbon dioxide, x 2 the response due to changes in sulfate aerosols, x 3 the response due to change in solar forcing, x 4 the response due to volcanic eruptions. Assuming that the observed climate signal is a linear combination of the components due to different forcing factors, this suggests a model of the form y = Xβ + u, where y is l 1, X = ( x 1... x m ) is a l m matrix consisting of all the modeled responses to forcing factors, and β is the vector of coefficients which we are trying to estimate. 10
11 Here β 1 could be the main component of interest (the greenhouse gas effect), while the full family of regression coefficients β 1, β 2,... is needed to determine the influence of all the components. We say that the greenhouse gas signal is detected if we can calculate an estimate and standard error for β 1, that lead to rejection of the null hypothesis H 0 : β 1 = 0. Given that we have detected a signal, we then estimate all the coefficients β 1, β 2,... to attribute the observed climate change to the different components. 11
12 If the residual vector u has covariance matrix C N, then the optimal generalized least squares (GLS) estimator of β is and ˆβ = (X T CN 1 X) 1 X T CN 1 y, (1) Var{ˆβ} = (X T C 1 N X) 1. (2) The crucial difficulty is that C N is unknown and therefore we need some reliable method of estimating it. One way of deriving (1) and (2) from the standard ordinary least squares (OLS) regression equations is to define a matrix P so that P C N P T = I l, and then to write the regression equation in the form P y = P Xβ + P u where the covariance matrix of P u is I l. The OLS estimator for β is then ˆβ = (X T P T P X) 1 X T P T P y with covariance (X T P T P X) 1. But if P and C N are invertible we will have P T P = CN 1, which leads to (1) and (2). 12
13 Usual approach in climate studies: take model observations from runs of the GCM in stationary conditions without forcing factors, typically 1,000 2,000 years long. With observations Y N from N replications of unforced model runs centered to 0 mean, estimate Ĉ N = 1 N Y NY T N. (3) Difficulty with (3): typically l > N and so Ĉ N is a singular matrix. Even though this in itself is not a fatal objection, e.g. valid versions of (1) and (2) can be given involving generalized inverses, the real difficulty is that with so many degrees of freedom in the covariance function, the low-amplitude components of covariance are not reliably estimated, and therefore, any direct attempt to apply (2) by substituting Ĉ N for C N will give poor estimates for ˆβ and innaccurate estimates of uncertainty. 13
14 Solution: define a κ l transformation matrix P (κ) corresponding to the κ largest EOFs (or principal components) of Y N, so that P (κ) Ĉ N P (κ)t = I κ, and then define the OLS regression equation which leads to an estimator P (κ) y = P (κ) Xβ + P (κ) u, β = (X T P (κ)t P (κ) X) 1 X T P (κ)t P (κ) y (4) and estimated covariance (X T P (κ)t P (κ) X) 1. (5) 14
15 Refinements (Allen and Tett 1999) 1. Estimate of covariance likely underestimate, alternative Var{ β} ˆ = (X T ĈN 1 X) 1 X T Ĉ 1 1 N C 1 N2 CN 1 X(X T Ĉ 1 1 N X) 1 1 Ĉ N1 estimate of C N from EOF construction, Ĉ N2 independent estimate from separate run of the climate model. 2. Correct for serial correlation 3. Effect of random variation in X? 4. How to choose κ? Proposed residual test r 2 = (y X β) T Ĉ 1 N (y X β) χ 2 κ m, Rejection of r 2 taken as indication that κ is too large. 15
16 16
17 Human influence on climate has been detected in surface air temperature, sea level pressure, tropopause height, ocean heat content, etc...but not, so far, precipitation. Models predict... Anthropogenic forcing implies increase in global mean precipitation Increasing at high latitudes Decreasing at sub-tropical latitudes We compare increases in land precipitation with changes in 14 climate models. 17
18 18
19 Trends over two periods Focus on 40 o S 70 o N Climate models: ANT (27 simulations from 8 models), anthropogenic signals only (GHGs+aerosols) ALL (50 simulations from 10 models), includes volcanic, solar NAT4 (15 simulations from 4 models), natural external forcing only Also look at ANT4, ALL4 (same as ANT, ALL but restricted to same 4 models as are in NAT4) 19
20 Methods Analyzed trends in annual zonal mean precipitation anomalies relative to Trends in observed and simulated precipitation compared using optimal fingerprint methodology. Linear trends from observations and averages of multi-model simulations exhibit important areas of consistency in spatial distribution of precipitation change. But Inconsistency on o N High level of uncertainty 20
21 21
22 22
23 Estimated combined effects of anthropogenic and natural forcings by regressing observed trends on anthropogenic trends Response to ALL or ANT was detected in observed trends, but not response to NAT alone. To separate contribution from ANT and NAT, need twosignal attribution analysis Response to anthropogenic forcing separable from that due to natural forcings and internal variability The two responses can be reliably separated, though the Allen-Tett for internal residual consistency failed (underestimate of internal variability?) 23
24 24
25 25
26 Detection of anthropogenic influence is robust... Little uncertainty about sign of trend Robust to different constructions of dataset (e.g. gridding, treatment of missing observations) Robust to use of signal patterns from subsets of models Consistent with our understanding of mechanisms but... Internal variability estimate is smaller than variability estimated from observations Discrepancy in magnitude of changes (generally, models underestimate observed trends) 26
27 27
28 28
29 Methods (Allen and Stott 2003) y = m i=1 C N = E{uu T }. x i β i + u = Xβ + u, Prewhitening matrix P (κ N), E{P uu T P T } = I κ. Estimate β minimizes r 2 ( β) = (P X β P y) T (P X β P y) = ũ T P T P ũ, β = (X T P T P X) 1 X T P T P y = F y. Rows of F are distinguishing fingerprints. Also r 2 min χ2 κ m. 29
30 If we ignore uncertainty in estimate of noise variance, β N[β, V ( β)], V ( β) = (X T P T P X) 1. Equivalently, ( β β) T (X T P P T X)( β β) χ 2 m. Use to construct joint confidence intervals for β. Also compute (one-dimensional) confidence interval for any linear function c T β. 30
31 Accounting for uncertainty in estimated noise patterns Estimation of P based on a set of noise realizations Ŷ 1. Take a second independent sample Ŷ 2. v columns of Ŷ 2 not independent DF based on Ŷ 2 is v 2 < v. Typically use overlapping segments, but set v 2 to be 1.5 times the number of nonoverlapping segments r 2 min = κ ˆV ( β) = F T Ŷ 2 Ŷ T 2 F v 2, ( β β) T ˆV ( β) 1 ( β β) mf m,v2, i=1 (P ũũ T P T ) i,i (P Ŷ 2 Ŷ T 2 P T ) i,i /v 2 (κ m)f κ m,v2. Final distributional result approximate. 31
32 Total Least Squares Early regression-based approaches to detection and attribution (e.g. Hegerl et al. (1996), Allen and Tett (1999)) relied on standard regression equations of the form Y = m j=1 β j X j + η (6) where Y is the observational record (e.g. a vector of trend in temperature means), X 1,..., X m are the signals from m climate models, and η is an error term. Instead of ordinary least squares, Allen and Stott (2003) proposed to fit (6) by total least squares, which allows for errors in the X j s as well as Y. Technique extended by Huntingford, Stott, Allen and Lambert (2006). I want to propose a (non-bayesian) alternative approach, stemming from Gleser (1981). 32
33 Gleser s formulation of errors in variables (EIV) x i = ( xi1 x i2 ) = ( ui1 u i2 ) + ( ei1 e i2 ), (7) x i1 and x i2 observations of dimensions p and r respectively, u i1 and u i2 are true unobserved signals, e i1 and e i2 noise with covariances σ 2 I p and σ 2 I r. Also MLE: choose B and u 11,..., u n1 to minimize Q = 1 2σ 2 (x i1 u i1 ) T (x i1 u i1 ) i u i2 = Bu i1. (8) i (x i2 Bu i1 ) T (x i2 Bu i1 ). (9) 33
34 Computing the MLE Choose B to minimize Q = 1 2σ 2 (x i2 Bx i1 ) T (I r + BB T ) 1 (x i2 Bx i1 ). (10) i (also called generalized least squares by Gleser) 34
35 Distribution Theory Assume estimator ˆB n based on n observations, U 1 the p n matrix whose columns are u 11,..., u n1. [Assumption A:] e i are i.i.d. random vectors with mean 0 and common covariance matrix σ 2 I p+r [Assumption C:] = lim n n 1 U 1 U1 T exists, positive definite. [Assumption E:] The cross-moments of the common distribution of the e i are identical, up to and including moments of order 4, to the corresponding moments of the multivariate normal distribution with the same mean and covariance matrix. Then the elements of the n 1/2 ( ˆB n B) have an asymptotic rpdimensional normal distribution with mean 0 and the covariance between the (i, j) and (i, j ) elements is given by σ 2 [σ 2 1 (I p + B T B) ] jj [I r + BB T ] ii. 35
36 Application to Climate I Identify x i2 with Y (single observation, dimension r) Identify x i1 with (X 1,..., X m ) T (dimension rm). Also B = ( β 1 I r β 2 I r... β m I r ). (11) GLSE chooses β 1,..., β m to minimize S = (Y j β j X j ) T (Y j β j X j ) 1 + βj 2. (12) Equivalent to Allen and Stott (2003). In principle, we could use Gleser s theory to approximate the asymptotic (co-)variances of the estimators. However, here n =
37 Application to Climate II Apply Gleser s formulation to individual years. Assume an EOF rotation applied so that x i2 is r-dimensional observation for year i, x i1 is rm-dimensional vector of model runs for year i. Assume: x i2 = u i2 + e i2 (r 1), x i1 = u i1 + e i1 (p 1), e i2 independent with mean 0 and covariance matrix σ 2 I r, e i1 independent with mean 0 and covariance matrix σ 2 I p, u i2 = Bu i1 where B = ( ) β 1 I r β 2 I r... β m I r. GLSE chooses β 1,..., β p to minimize Q = 1 2σ 2 i (Y i j β j X ij ) T (Y i j β j X ij ) 1 + βj 2. (13) or let B be completely arbitrary??? (Needs rm < n, so we really would have to use a low-order EOF decomposition for this.) 37
38 Comments / Caveats Gleser s theorem not directly applicable for various reasons, but could adapt. Covariance matrix C now needed for individual years, not linear trends. Simpler than earlier discussion. Could also examine whether forced model runs are consistent with estimated covariances. As presently formulated, doesn t allow for (temporal) autocorrelation. 38
39 Detection and Attribution Analysis for Precipitation We re-analyze the data from the paper of Zhang et al., Nature 2007, discussed earlier. 39
40 Datasets Latitude-band averages 75 years 11 ten-degree latitude bands (from 40 o S to 70 o N) One observational dataset Several forced model runs under different forcings (ANT, ALL, NAT) Control runs typically 6 control runs of 1000 years broken up into 75-year segments 40
41 Traditional D&A Analysis 1. For each of 11 latitude bands, linear trend is calculated by OLS, including control runs 2. Covariance matrix from control run used to define the weights for EOF decomposition 3. Truncate the EOFs, e.g. by Allen-Tett test (anticipate 5 EOFs retained 4. Same formation of linear trends and EOFs applied to forced model runs to develop a set of comparison vectors. 5. Optimal weights for D&A defined by regressing observational data EOFs on forcing model EOFs. Calculate standard errors, tests of significance, etc. 41
42 Estimation of the Covariance Matrix The data we have available consists of 17 runs of 75 years data from control runs of HADCM3. The data are 11-dimensional because of the 11 latitude bands. The covariance matrix C is defined to be the covariance matrix of estimators of linear trends over the 75 years. Possible methods of estimating C include: 42
43 Method 1 Calculate linear trends for each of the 11 latitude bands and each of the 17 runs; then estimate C (say, as Ĉ 1 ) as the sample covariance matrix Ĉ 1 = (B B) T (B B) (14) 16 where B is the matrix of linear trend estimates and B is the mean of B over all rows. 43
44 Method 2 Suppose Y ijk is observation in year i from latitude band j and control model run k. Suppose that, for each i and k, {Y ijk, j = 1,..., 11} is an independent vector from N[0, V ], where V is the covariance matrix associated with a single year of control data. Linear trend estimator for region j and model run k given by B jk = 75 i=1 for weights w i. Then for any two regions j, j, Cov(B jk, B j k ) = 75 w i Y ijk (15) i=1 If ˆV is sample covariance matrix, estimate wi 2 V jj. (16) Ĉ 2 = ( w 2 i )ˆV. (17) 44
45 Method 3 Correct for autocorrelation in method 2. Assume h-order autocovariances negligible for h > H some given H. Then we can modify (16) to Cov(B jk, B j k ) = H h= H 75 i=1 w i w i+h Cov(Y ijk, Y (i+h)j k ) (18) where w i+h is defined to be 0 if i + h < 1 or i + h > 75. If ˆv hjj denotes sample covariance of Y ijk and Y (i+h)j k (computed from all control years of data) then this implies a covariance estimator Ĉ 3 with entries Ĉ 3 (j, j ) = H h= H 75 i=1 w i w i+hˆv hjj. (19) 45
46 Looking for autocorrelations in the control runs Sample autocorrelations and partial autocorrelations were calculated from the control data. They were considered statistically significant if outside the range ± 2 n These results suggest that significant autocorrelations do not persist beyond lag 2. However, as a test, the estimator Ĉ 3 was computed for each of H = 1, 2, 3, 4, 5. 46
47 Testing one estimated covariance matrix against another Advantage of Ĉ 1 makes no assumption at all about the autocorrelation structure, but it s based on a small sample (17). Ĉ 2 or Ĉ 3 far larger sample size, but maybe incorrect assumptions. Suggests a hypothesis test that of either H 0 : C = Ĉ 2 or H 0 : C = Ĉ 3, treating Ĉ 1 as a sample covariance matrix. Korin (1968) proposed likelihood ratio test statistic { ( ) } Σ0 L = N log p + tr(sσ 1 S 0 ) (20) where S is a p p sample covariance matrix with N degrees of freedom, and Σ 0 is the null hypothesis value of the true covariance matrix. Sampling distributions by Bartlett correction to the χ 2 test, an alternative approximation proposed by Korin, or simulation. 47
48 Application Ĉ 1 against Ĉ 2 (as null hypothesis): reject under Bartlett s approximation, accept under Korin s or simulation (p-values.047,.077,.11) but maybe not a clear-cut result. However for Ĉ 1 against Ĉ 3 (as null hypothesis), with any of H = 1, 2, 3, 4, 5, accept null hypothesis. As a specific example, for H = 2, Bartlett s, Korin s and the simulation test yield p- values 0.33, 0.42 and Conclusion: matter. Ĉ 3 seems fine, and the actual value of H doesn t 48
49 On the other hand, where we see real differences among the estimators is in computing their inverses. As an example, the diagonal entries of Ĉ1 1 are: The corresponding values for Ĉ 1 3 with H = 0 5 are: H=0: H=1: H=2: H=3: H=4: H=5: Evidently, there is much more stability in any of the versions of Ĉ3 1 (raising the question that maybe no dimension reduction is actually needed??) 49
50 Results Applied to D&A 50
51 Analysis using full data matrix (11 PCs) Analyses using Ĉ 1, Ĉ 2, Ĉ 3 (H = 2) Models ˆβ SE t statistic P-value ALL ANT NAT Models ˆβ SE t statistic P-value ALL ANT NAT Models ˆβ SE t statistic P-value ALL ANT NAT
52 Analysis using 5 PCs Analyses using Ĉ 1, Ĉ 2, Ĉ 3 (H = 2) Models ˆβ SE t statistic P-value ALL ANT NAT Models ˆβ SE t statistic P-value ALL ANT NAT Models ˆβ SE t statistic P-value ALL ANT NAT
53 Conclusions (preliminary) Using the full data matrix, there is a lot of discrepancy between the results based on Ĉ 1 and either of Ĉ 2 or Ĉ 3, but the results comparing Ĉ 2 and Ĉ 3 are quite close Using 5 PCs, there is much better agreement between the three covariance models, but results are quite different from those using all 11 PCs Based on all 11 PCs with either Ĉ 2 or Ĉ 3, we may claim to have attribution for the all forcings model, but not for either ANT or NAT. 53
54 More on Detection and Attribution Bayesian approaches: Levine and Berliner (J. Clim 12, , 1999) Levine, Berliner and Shea (J. Clim 13, , 2000) Berliner and Kim (J. Clim 21, , 2008) Detection and attribution applied to multi-model ensembles: Lopez, Tebaldi, New, Stainforth, Allen and Kettleborough (J. Clim 19, , 2006) 54
55 References Allen MR and Stott PA (2003), Estimating signal amplitudes in optimal fingerprinting. Part I: Theory. Clim. Dyn. 21, Allen MR and Tett SFB (1999), Checking for model consistency in optimal fingerprinting. Climate Dynamics 15, Berliner LM, Levine RA and Shea DJ (2000), Bayesian Climate Change Assessment. Journal of Climate 13, Berliner LM and Kim Y (2008), Bayesian design and analysis for superensemble-based climate forecasting. Journal of Climate 21, Gleser LJ (1981), Estimation of a multivariate Errors in Variables regression model: Large sample results. Annals of Statistics 9, Hasselmann K (1997), Multi-pattern fingerprint method for detection and attribution of climate change. Climate Dynamics 13,
56 Hegerl GC and North GR (1997), Comparison of statistical optimal approaches to detecting anthropogenic climate change. Journal of Climate 10, Hegerl G, von Storch H, Hasselmann K, Santer BD, Cubasch U and Jones PD (1996), Detecting greenhouse-gas-induced climate change with an optimal fingerprint method. Journal of Climate 9, Levine RA and Berliner LM (1999), Statistical Principles for Climate Change Studies. Journal of Climate 12, Lopez A, Tebaldi C, New M, Stainforth D, Allen M and Kettleborough J (2006). Two approaches to quantifying uncertainty in global temperature changes. Journal of Climate 19, Santer BD, Taylor KE, Wigley TML, Johns TC, Jones PD, Karoly DJ, Mitchell JFB, Oort AH, Penner JE, Ramaswamy V, Schwarzkopf MD, Stouffer RJ and Tett S. (1996), A search for human influences on the thermal structure of the atmosphere. Nature 382,
A HIERARCHICAL MODEL FOR REGRESSION-BASED CLIMATE CHANGE DETECTION AND ATTRIBUTION
A HIERARCHICAL MODEL FOR REGRESSION-BASED CLIMATE CHANGE DETECTION AND ATTRIBUTION Richard L Smith University of North Carolina and SAMSI Joint Statistical Meetings, Montreal, August 7, 2013 www.unc.edu/~rls
More informationROBUST MEASUREMENT OF THE DURATION OF
ROBUST MEASUREMENT OF THE DURATION OF THE GLOBAL WARMING HIATUS Ross McKitrick Department of Economics University of Guelph Revised version, July 3, 2014 Abstract: The IPCC has drawn attention to an apparent
More informationAnthropogenic warming of central England temperature
ATMOSPHERIC SCIENCE LETTERS Atmos. Sci. Let. 7: 81 85 (2006) Published online 18 September 2006 in Wiley InterScience (www.interscience.wiley.com).136 Anthropogenic warming of central England temperature
More informationEvidence of Decadal Climate Prediction Skill Resulting from Changes in Anthropogenic Forcing
15 OCTOBER 2006 L E E E T A L. 5305 Evidence of Decadal Climate Prediction Skill Resulting from Changes in Anthropogenic Forcing TERRY C. K. LEE Department of Mathematics and Statistics, University of
More informationSupplementary Information for:
Supplementary Information for: Linkage between global sea surface temperature and hydroclimatology of a major river basin of India before and after 1980 P. Sonali, Ravi S. Nanjundiah, & D. Nagesh Kumar
More informationA STATISTICAL APPROACH TO OPERATIONAL ATTRIBUTION
A STATISTICAL APPROACH TO OPERATIONAL ATTRIBUTION Richard L. Smith Department of Statistics and Operations Research University of North Carolina Chapel Hill, NC 27599-3260, USA rls@email.unc.edu IDAG Meeting
More informationAttribution of anthropogenic influence on seasonal sea level pressure
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 36, L23709, doi:10.1029/2009gl041269, 2009 Attribution of anthropogenic influence on seasonal sea level pressure N. P. Gillett 1 and P. A.
More informationESTIMATING THE MEAN LEVEL OF FINE PARTICULATE MATTER: AN APPLICATION OF SPATIAL STATISTICS
ESTIMATING THE MEAN LEVEL OF FINE PARTICULATE MATTER: AN APPLICATION OF SPATIAL STATISTICS Richard L. Smith Department of Statistics and Operations Research University of North Carolina Chapel Hill, N.C.,
More informationPeter Hoff Linear and multilinear models April 3, GLS for multivariate regression 5. 3 Covariance estimation for the GLM 8
Contents 1 Linear model 1 2 GLS for multivariate regression 5 3 Covariance estimation for the GLM 8 4 Testing the GLH 11 A reference for some of this material can be found somewhere. 1 Linear model Recall
More informationOn Sampling Errors in Empirical Orthogonal Functions
3704 J O U R N A L O F C L I M A T E VOLUME 18 On Sampling Errors in Empirical Orthogonal Functions ROBERTA QUADRELLI, CHRISTOPHER S. BRETHERTON, AND JOHN M. WALLACE University of Washington, Seattle,
More informationHeteroskedasticity. Part VII. Heteroskedasticity
Part VII Heteroskedasticity As of Oct 15, 2015 1 Heteroskedasticity Consequences Heteroskedasticity-robust inference Testing for Heteroskedasticity Weighted Least Squares (WLS) Feasible generalized Least
More informationMultiple Linear Regression
Multiple Linear Regression Simple linear regression tries to fit a simple line between two variables Y and X. If X is linearly related to Y this explains some of the variability in Y. In most cases, there
More informationA Comparison of Surface Air Temperature Variability in Three 1000-Yr Coupled Ocean Atmosphere Model Integrations
VOLUME 13 JOURNAL OF CLIMATE 1FEBRUARY 2000 A Comparison of Surface Air Temperature Variability in Three 1000-Yr Coupled Ocean Atmosphere Model Integrations RONALD J. STOUFFER NOAA/Geophysical Fluid Dynamics
More informationWeighted Least Squares
Weighted Least Squares The standard linear model assumes that Var(ε i ) = σ 2 for i = 1,..., n. As we have seen, however, there are instances where Var(Y X = x i ) = Var(ε i ) = σ2 w i. Here w 1,..., w
More informationMultivariate Regression Analysis
Matrices and vectors The model from the sample is: Y = Xβ +u with n individuals, l response variable, k regressors Y is a n 1 vector or a n l matrix with the notation Y T = (y 1,y 2,...,y n ) 1 x 11 x
More informationProbability and Statistics Notes
Probability and Statistics Notes Chapter Seven Jesse Crawford Department of Mathematics Tarleton State University Spring 2011 (Tarleton State University) Chapter Seven Notes Spring 2011 1 / 42 Outline
More informationThe use of marine data for attribution of climate change and constraining climate predictions
The use of marine data for attribution of climate change and constraining climate predictions Peter Stott, Climar III, Thursday 8 May 2008 The use of marine data to: Quantify the contribution of human
More informationUncertainties in the attribution of greenhouse gas warming and implications for climate prediction
Uncertainties in the attribution of greenhouse gas warming and implications for climate prediction Gareth S. Jones, Peter A. Stott and John F. B. Mitchell - Met Office Hadley Centre, Exeter, UK Paper accepted
More informationDetection and Attribution of Observed Changes in Northern Hemisphere Spring Snow Cover
6904 J O U R N A L O F C L I M A T E VOLUME 26 Detection and Attribution of Observed Changes in Northern Hemisphere Spring Snow Cover DAVID E. RUPP AND PHILIP W. MOTE Oregon Climate Change Research Institute,
More informationA test for improved forecasting performance at higher lead times
A test for improved forecasting performance at higher lead times John Haywood and Granville Tunnicliffe Wilson September 3 Abstract Tiao and Xu (1993) proposed a test of whether a time series model, estimated
More informationPoint-to-point response to reviewers comments
Point-to-point response to reviewers comments Reviewer #1 1) The authors analyze only one millennial reconstruction (Jones, 1998) with the argument that it is the only one available. This is incorrect.
More informationEconomics 582 Random Effects Estimation
Economics 582 Random Effects Estimation Eric Zivot May 29, 2013 Random Effects Model Hence, the model can be re-written as = x 0 β + + [x ] = 0 (no endogeneity) [ x ] = = + x 0 β + + [x ] = 0 [ x ] = 0
More informationMIT Joint Program on the Science and Policy of Global Change Constraining Climate Model Properties Using Optimal Fingerprint Detection Methods
MIT Joint Program on the Science and Policy of Global Change Constraining Climate Model Properties Using Optimal Fingerprint Detection Methods Chris E. Forest, Myles R. Allen, Andrei P. Sokolov and Peter
More informationAdvanced Econometrics
Based on the textbook by Verbeek: A Guide to Modern Econometrics Robert M. Kunst robert.kunst@univie.ac.at University of Vienna and Institute for Advanced Studies Vienna May 16, 2013 Outline Univariate
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 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 informationEuropean temperatures prior to the 20 th century fundamentally a consequence of unforced variability? (Bengtsson et al., 2006)
01-12-2000 1 European temperatures prior to the 20 th century fundamentally a consequence of unforced variability? (Bengtsson et al., 2006) Not really (Hegerl, Luterbacher et al. 2011 NGeo) Observations
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 informationBAYESIAN KRIGING AND BAYESIAN NETWORK DESIGN
BAYESIAN KRIGING AND BAYESIAN NETWORK DESIGN Richard L. Smith Department of Statistics and Operations Research University of North Carolina Chapel Hill, N.C., U.S.A. J. Stuart Hunter Lecture TIES 2004
More informationCircle the single best answer for each multiple choice question. Your choice should be made clearly.
TEST #1 STA 4853 March 6, 2017 Name: Please read the following directions. DO NOT TURN THE PAGE UNTIL INSTRUCTED TO DO SO Directions This exam is closed book and closed notes. There are 32 multiple choice
More informationLecture 15. Hypothesis testing in the linear model
14. Lecture 15. Hypothesis testing in the linear model Lecture 15. Hypothesis testing in the linear model 1 (1 1) Preliminary lemma 15. Hypothesis testing in the linear model 15.1. Preliminary lemma Lemma
More informationstatistical methods for tailoring seasonal climate forecasts Andrew W. Robertson, IRI
statistical methods for tailoring seasonal climate forecasts Andrew W. Robertson, IRI tailored seasonal forecasts why do we make probabilistic forecasts? to reduce our uncertainty about the (unknown) future
More informationAttributing northern high-latitude precipitation change over the period to human influence
Clim Dyn (2015) 45:1713 1726 DOI 10.1007/s00382-014-2423-y Attributing northern high-latitude precipitation change over the period 1966 2005 to human influence Hui Wan Xuebin Zhang Francis Zwiers Seung
More informationIntroduction to climate modelling: Evaluating climate models
Introduction to climate modelling: Evaluating climate models Why? How? Professor David Karoly School of Earth Sciences, University of Melbourne Experiment design Detection and attribution of climate change
More informationA Bivariate Time Series Approach to Anthropogenic Trend Detection in Hemispheric Mean Temperatures
18 JOURNAL OF CLIMATE A Bivariate Time Series Approach to Anthropogenic Trend Detection in Hemispheric Mean Temperatures RICHARD L. SMITH Department of Statistics, University of rth Carolina at Chapel
More informationLinear models and their mathematical foundations: Simple linear regression
Linear models and their mathematical foundations: Simple linear regression Steffen Unkel Department of Medical Statistics University Medical Center Göttingen, Germany Winter term 2018/19 1/21 Introduction
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 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 informationExplaining Changes in Extremes and Decadal Climate Fluctuations
Explaining Changes in Extremes and Decadal Climate Fluctuations Gerald A. Meehl Julie Arblaster, Claudia Tebaldi, Aixue Hu, Ben Santer Explaining changes implies attributing those changes to some cause
More information1 Introduction to Generalized Least Squares
ECONOMICS 7344, Spring 2017 Bent E. Sørensen April 12, 2017 1 Introduction to Generalized Least Squares Consider the model Y = Xβ + ɛ, where the N K matrix of regressors X is fixed, independent of the
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION doi:10.1038/nature11784 Methods The ECHO-G model and simulations The ECHO-G model 29 consists of the 19-level ECHAM4 atmospheric model and 20-level HOPE-G ocean circulation model.
More informationEmpirical Economic Research, Part II
Based on the text book by Ramanathan: Introductory Econometrics Robert M. Kunst robert.kunst@univie.ac.at University of Vienna and Institute for Advanced Studies Vienna December 7, 2011 Outline Introduction
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 informationHow unusual is the recent series of warm years?
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L24706, doi:10.1029/2008gl036228, 2008 How unusual is the recent series of warm years? E. Zorita, 1 T. F. Stocker, 2 and H. von Storch
More informationIntroductory Econometrics
Based on the textbook by Wooldridge: : A Modern Approach Robert M. Kunst robert.kunst@univie.ac.at University of Vienna and Institute for Advanced Studies Vienna December 11, 2012 Outline Heteroskedasticity
More informationLECTURE ON HAC COVARIANCE MATRIX ESTIMATION AND THE KVB APPROACH
LECURE ON HAC COVARIANCE MARIX ESIMAION AND HE KVB APPROACH CHUNG-MING KUAN Institute of Economics Academia Sinica October 20, 2006 ckuan@econ.sinica.edu.tw www.sinica.edu.tw/ ckuan Outline C.-M. Kuan,
More informationQuick Review on Linear Multiple Regression
Quick Review on Linear Multiple Regression Mei-Yuan Chen Department of Finance National Chung Hsing University March 6, 2007 Introduction for Conditional Mean Modeling Suppose random variables Y, X 1,
More informationStatistics for analyzing and modeling precipitation isotope ratios in IsoMAP
Statistics for analyzing and modeling precipitation isotope ratios in IsoMAP The IsoMAP uses the multiple linear regression and geostatistical methods to analyze isotope data Suppose the response variable
More informationSTATISTICS OF CLIMATE EXTREMES// TRENDS IN CLIMATE DATASETS
STATISTICS OF CLIMATE EXTREMES// TRENDS IN CLIMATE DATASETS Richard L Smith Departments of STOR and Biostatistics, University of North Carolina at Chapel Hill and Statistical and Applied Mathematical Sciences
More informationQuantitative Analysis of Financial Markets. Summary of Part II. Key Concepts & Formulas. Christopher Ting. November 11, 2017
Summary of Part II Key Concepts & Formulas Christopher Ting November 11, 2017 christopherting@smu.edu.sg http://www.mysmu.edu/faculty/christophert/ Christopher Ting 1 of 16 Why Regression Analysis? Understand
More informationCh 2: Simple Linear Regression
Ch 2: Simple Linear Regression 1. Simple Linear Regression Model A simple regression model with a single regressor x is y = β 0 + β 1 x + ɛ, where we assume that the error ɛ is independent random component
More informationHadley Centre for Climate Prediction and Research, Met Office, FitzRoy Road, Exeter, EX1 3PB, UK.
Temperature Extremes, the Past and the Future. S Brown, P Stott, and R Clark Hadley Centre for Climate Prediction and Research, Met Office, FitzRoy Road, Exeter, EX1 3PB, UK. Tel: +44 (0)1392 886471 Fax
More informationDetection and attribution, forced changes, natural variability, signal and noise, ensembles
ETH Zurich Reto Knutti Detection and attribution, forced changes, natural variability, signal and noise, ensembles Reto Knutti, IAC ETH What s wrong with this presentation? For the next two decades, a
More information1.5 Statistical significance of linear trend
CHAPTER 1. DATA ANALYSIS 42 1.5 Statistical significance of linear trend Temporal trends can be calculated for any time varying data. When we perform the linear regression, we may find non-zero linear
More informationChapter 6. Panel Data. Joan Llull. Quantitative Statistical Methods II Barcelona GSE
Chapter 6. Panel Data Joan Llull Quantitative Statistical Methods II Barcelona GSE Introduction Chapter 6. Panel Data 2 Panel data The term panel data refers to data sets with repeated observations over
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 informationAsymptotic Multivariate Kriging Using Estimated Parameters with Bayesian Prediction Methods for Non-linear Predictands
Asymptotic Multivariate Kriging Using Estimated Parameters with Bayesian Prediction Methods for Non-linear Predictands Elizabeth C. Mannshardt-Shamseldin Advisor: Richard L. Smith Duke University Department
More informationHeteroskedasticity and Autocorrelation
Lesson 7 Heteroskedasticity and Autocorrelation Pilar González and Susan Orbe Dpt. Applied Economics III (Econometrics and Statistics) Pilar González and Susan Orbe OCW 2014 Lesson 7. Heteroskedasticity
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 informationCLIMATE EXTREMES AND GLOBAL WARMING: A STATISTICIAN S PERSPECTIVE
CLIMATE EXTREMES AND GLOBAL WARMING: A STATISTICIAN S PERSPECTIVE Richard L. Smith Department of Statistics and Operations Research University of North Carolina, Chapel Hill rls@email.unc.edu Statistics
More informationRichard L. Smith Department of Statistics and Operations Research University of North Carolina Chapel Hill, NC
EXTREME VALUE THEORY Richard L. Smith Department of Statistics and Operations Research University of North Carolina Chapel Hill, NC 27599-3260 rls@email.unc.edu AMS Committee on Probability and Statistics
More informationTotal Least Squares Approach in Regression Methods
WDS'08 Proceedings of Contributed Papers, Part I, 88 93, 2008. ISBN 978-80-7378-065-4 MATFYZPRESS Total Least Squares Approach in Regression Methods M. Pešta Charles University, Faculty of Mathematics
More informationMS&E 226: Small Data
MS&E 226: Small Data Lecture 15: Examples of hypothesis tests (v5) Ramesh Johari ramesh.johari@stanford.edu 1 / 32 The recipe 2 / 32 The hypothesis testing recipe In this lecture we repeatedly apply the
More informationEconometrics. Week 8. Fall Institute of Economic Studies Faculty of Social Sciences Charles University in Prague
Econometrics Week 8 Institute of Economic Studies Faculty of Social Sciences Charles University in Prague Fall 2012 1 / 25 Recommended Reading For the today Instrumental Variables Estimation and Two Stage
More informationAsymptotic distribution of GMM Estimator
Asymptotic distribution of GMM Estimator Eduardo Rossi University of Pavia Econometria finanziaria 2010 Rossi (2010) GMM 2010 1 / 45 Outline 1 Asymptotic Normality of the GMM Estimator 2 Long Run Covariance
More information[y i α βx i ] 2 (2) Q = i=1
Least squares fits This section has no probability in it. There are no random variables. We are given n points (x i, y i ) and want to find the equation of the line that best fits them. We take the equation
More informationSpecification testing in panel data models estimated by fixed effects with instrumental variables
Specification testing in panel data models estimated by fixed effects wh instrumental variables Carrie Falls Department of Economics Michigan State Universy Abstract I show that a handful of the regressions
More informationA Practitioner s Guide to Cluster-Robust Inference
A Practitioner s Guide to Cluster-Robust Inference A. C. Cameron and D. L. Miller presented by Federico Curci March 4, 2015 Cameron Miller Cluster Clinic II March 4, 2015 1 / 20 In the previous episode
More informationDetection and attribution of Atlantic salinity changes
GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L21702, doi:10.1029/2008gl035874, 2008 Detection and attribution of Atlantic salinity changes Peter A. Stott, 1 Rowan T. Sutton, 2 and Doug M. Smith 1 Received 1
More informationThe importance of including variability in climate
D. M. H. Sexton and G. R. Harris SUPPLEMENTARY INFORMATION The importance of including variability in climate change projections used for adaptation Modification to time scaling to allow for short-term
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 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 informationStatistical downscaling daily rainfall statistics from seasonal forecasts using canonical correlation analysis or a hidden Markov model?
Statistical downscaling daily rainfall statistics from seasonal forecasts using canonical correlation analysis or a hidden Markov model? Andrew W. Robertson International Research Institute for Climate
More informationUncertainty in Models of the Ocean and Atmosphere
Uncertainty in Models of the Ocean and Atmosphere Robert N. Miller College of Oceanic and Atmospheric Sciences Oregon State University Corvallis, OR 97330 Uncertainty in Models of the Ocean and Atmosphere
More informationCointegrating Regressions with Messy Regressors: J. Isaac Miller
NASMES 2008 June 21, 2008 Carnegie Mellon U. Cointegrating Regressions with Messy Regressors: Missingness, Mixed Frequency, and Measurement Error J. Isaac Miller University of Missouri 1 Messy Data Example
More informationProblems. Suppose both models are fitted to the same data. Show that SS Res, A SS Res, B
Simple Linear Regression 35 Problems 1 Consider a set of data (x i, y i ), i =1, 2,,n, and the following two regression models: y i = β 0 + β 1 x i + ε, (i =1, 2,,n), Model A y i = γ 0 + γ 1 x i + γ 2
More informationMultiple Linear Regression
Multiple Linear Regression University of California, San Diego Instructor: Ery Arias-Castro http://math.ucsd.edu/~eariasca/teaching.html 1 / 42 Passenger car mileage Consider the carmpg dataset taken from
More informationIntermediate Econometrics
Intermediate Econometrics Heteroskedasticity Text: Wooldridge, 8 July 17, 2011 Heteroskedasticity Assumption of homoskedasticity, Var(u i x i1,..., x ik ) = E(u 2 i x i1,..., x ik ) = σ 2. That is, the
More informationSensitivity of climate change signals deduced from multi-model Monte Carlo experiments
CLIMATE RESEARCH Vol. 22: 189 24, 22 Published November 4 Clim Res Sensitivity of climate change signals deduced from multi-model Monte Carlo experiments Heiko Paeth*, Andreas Hense Meteorologisches Institut,
More informationSTAT 100C: Linear models
STAT 100C: Linear models Arash A. Amini June 9, 2018 1 / 56 Table of Contents Multiple linear regression Linear model setup Estimation of β Geometric interpretation Estimation of σ 2 Hat matrix Gram matrix
More informationFractional Hot Deck Imputation for Robust Inference Under Item Nonresponse in Survey Sampling
Fractional Hot Deck Imputation for Robust Inference Under Item Nonresponse in Survey Sampling Jae-Kwang Kim 1 Iowa State University June 26, 2013 1 Joint work with Shu Yang Introduction 1 Introduction
More informationSTA 114: Statistics. Notes 21. Linear Regression
STA 114: Statistics Notes 1. Linear Regression Introduction A most widely used statistical analysis is regression, where one tries to explain a response variable Y by an explanatory variable X based on
More informationEconometrics 2, Class 1
Econometrics 2, Class Problem Set #2 September 9, 25 Remember! Send an email to let me know that you are following these classes: paul.sharp@econ.ku.dk That way I can contact you e.g. if I need to cancel
More informationCasuality and Programme Evaluation
Casuality and Programme Evaluation Lecture V: Difference-in-Differences II Dr Martin Karlsson University of Duisburg-Essen Summer Semester 2017 M Karlsson (University of Duisburg-Essen) Casuality and Programme
More informationTime Series Analysis. James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY
Time Series Analysis James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY & Contents PREFACE xiii 1 1.1. 1.2. Difference Equations First-Order Difference Equations 1 /?th-order Difference
More informationRepeated observations on the same cross-section of individual units. Important advantages relative to pure cross-section data
Panel data Repeated observations on the same cross-section of individual units. Important advantages relative to pure cross-section data - possible to control for some unobserved heterogeneity - possible
More information(x t. x t +1. TIME SERIES (Chapter 8 of Wilks)
45 TIME SERIES (Chapter 8 of Wilks) In meteorology, the order of a time series matters! We will assume stationarity of the statistics of the time series. If there is non-stationarity (e.g., there is a
More informationA Non-Parametric Approach of Heteroskedasticity Robust Estimation of Vector-Autoregressive (VAR) Models
Journal of Finance and Investment Analysis, vol.1, no.1, 2012, 55-67 ISSN: 2241-0988 (print version), 2241-0996 (online) International Scientific Press, 2012 A Non-Parametric Approach of Heteroskedasticity
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 informationAutoregressive Moving Average (ARMA) Models and their Practical Applications
Autoregressive Moving Average (ARMA) Models and their Practical Applications Massimo Guidolin February 2018 1 Essential Concepts in Time Series Analysis 1.1 Time Series and Their Properties Time series:
More informationEvaluation of proxy-based millennial reconstruction methods
Clim Dyn DOI 10.1007/s0038-007-0351-9 Evaluation of proxy-based millennial reconstruction methods Terry C. K. Lee Æ Francis W. Zwiers Æ Min Tsao Received: 3 July 007 / Accepted: 9 November 007 Ó Springer-Verlag
More informationConstraints on climate change from a multi-thousand member ensemble of simulations
GEOPHYSICAL RESEARCH LETTERS, VOL. 32, L23825, doi:10.1029/2005gl024452, 2005 Constraints on climate change from a multi-thousand member ensemble of simulations C. Piani, D. J. Frame, D. A. Stainforth,
More informationBootstrapping Heteroskedasticity Consistent Covariance Matrix Estimator
Bootstrapping Heteroskedasticity Consistent Covariance Matrix Estimator by Emmanuel Flachaire Eurequa, University Paris I Panthéon-Sorbonne December 2001 Abstract Recent results of Cribari-Neto and Zarkos
More informationAnalysis of Climate Model Outputs Using Linear and Nonlinear Statistical Techniques
WDS'09 Proceedings of Contributed Papers, Part III, 235 240, 2009. ISBN 978-80-7378-103-3 MATFYZPRESS Analysis of Climate Model Outputs Using Linear and Nonlinear Statistical Techniques Z. Chladova Faculty
More informationChapter 14. Linear least squares
Serik Sagitov, Chalmers and GU, March 5, 2018 Chapter 14 Linear least squares 1 Simple linear regression model A linear model for the random response Y = Y (x) to an independent variable X = x For a given
More informationCorrection notice. Nature Climate Change 2, (2012)
Correction notice Nature Climate Change 2, 524 549 (2012) Human-induced global ocean warming on multidecadal timescales P. J. Gleckler, B. D. Santer, C. M. Domingues, D.W. Pierce, T. P. Barnett, J. A.
More informationUnit roots in vector time series. Scalar autoregression True model: y t 1 y t1 2 y t2 p y tp t Estimated model: y t c y t1 1 y t1 2 y t2
Unit roots in vector time series A. Vector autoregressions with unit roots Scalar autoregression True model: y t y t y t p y tp t Estimated model: y t c y t y t y t p y tp t Results: T j j is asymptotically
More informationLecture 2: Linear Models. Bruce Walsh lecture notes Seattle SISG -Mixed Model Course version 23 June 2011
Lecture 2: Linear Models Bruce Walsh lecture notes Seattle SISG -Mixed Model Course version 23 June 2011 1 Quick Review of the Major Points The general linear model can be written as y = X! + e y = vector
More informationThe Standard Linear Model: Hypothesis Testing
Department of Mathematics Ma 3/103 KC Border Introduction to Probability and Statistics Winter 2017 Lecture 25: The Standard Linear Model: Hypothesis Testing Relevant textbook passages: Larsen Marx [4]:
More informationTopic 4 Unit Roots. Gerald P. Dwyer. February Clemson University
Topic 4 Unit Roots Gerald P. Dwyer Clemson University February 2016 Outline 1 Unit Roots Introduction Trend and Difference Stationary Autocorrelations of Series That Have Deterministic or Stochastic Trends
More information