Advanced time-series analysis (University of Lund, Economic History Department)
|
|
- Donna Rice
- 5 years ago
- Views:
Transcription
1 Advanced time-series analysis (University of Lund, Economic History Department) 30 Jan-3 February and March 2012 Lecture 3 Monte Carlo simulations and Bootstrapping. 3.a. What is a Monte Carlo simulation? Imagine you have a truffle-seeking pig. You would like to know how good she is. You cannot just let her loose in the forest and count the number of truffles he found, since you do not know how rich in truffles the area was. You need a controlled experiment. So you choose an area, rail it off, and you randomly bury truffles. Now, you can let the pig loose, and you can count how many of the truffles were found. But one experience is not enough: after all, what if your pig was just too lucky or was disturbed by a squirrel? So you should repeat the experiment as many times as possible. As the number of experiments goes toward infinity, the hit ratio of your pig that you calculated from the experiments will get closer and closer to the real hit probability of the pig (law of large numbers ). This is a Monte Carlo simulation. Monte Carlo experiments are based on repeated random sampling. It was first used in Los Alamos National Laboratory by John von Neuman, Stanislaw Ulam and Nicholas Metropolis, the first two were mathematicians and the third one was nuclear physicist, to find out how far neutrons would travel in different materials (important for radiation shielding). The codename for the then secret method was Monte Carlo. We use Monte Carlo experiments in cases when we know the random process behind a phenomenon, but an analytical solution would be far too difficult to achieve. 3.b A simple example: the critical values for a DF unit-root test We know that if the null-hypothesis of the DF test is correct, the distribution of the coefficients will not follow the student s t distribution. As a consequence, you will need different critical values. One way to have an adjusted table is to design an MC experiment! We are going to assume that the DGP is the following: y 1 y where WN(0,1) and y0 1 t t 1 t So what we are going to do is the following: t We choose sample size, let us say 100. Then we will use a random number generator to create a series based on the above DGP. This is basically a sample. Then we run a basic DF test: yt 0 1yt 1 ut and we calculate the test statistics as α 1 divided by the standard error of it. WE store this statistics. Then we redo this experiment 9999 times (so we have experiments in total) and store all the test statistics. Our estimate for the critical value of the unit root test statistics is going to be the lower fifth percentile (in case of a one-tailed test). 1
2 the code in Eviews wfcreate dfsimu u 'we create the workfile called dfsimu with observations series testresults 'this is the variable where we store our testresults for final analysis scalar fivepercentcrit 'this is where we store our estimate of the 5% critical value for!i=1 to 'we will repeat this 1000 times series y=0 ' we create the series y with the starting value 0 smpl 'no we need to reset the sample to y=1+y(-1)+nrnd ' we generate the random series according our assumptions equation eq.ls d(y) c y(-1) 'we run the test regression on the generated 100 observations 'we reset the sample to the testresults(!i)=eq.@tstats(2) 'save the test stat 'return to the start of the for loop fivepercentcrit=@quantile(testresults,.05) By changing the parameters you can come to a table of your own: Simulated critical values of a DF test (10000 experiments) level of significance (one-tailed) 10% 5% 1% n=100 without constant with constant n=50 without constant with constant c. Storage and hunger in an ancient society There is a debate on the ability to have long-term storage of grain in pre-industrial societies. The historical evidence is sporadic and not decisive. McCloskey and Nash 1 utilized a simple probabilistic approach to find out whether the storage was. They assumed that grain production followed a normal probability distribution with mean 100 and standard deviation 35 and found out that this leads to very long waiting times between famines. Van Leeuwen et al (2011) 2 continued in their 1 McCloskey, D. N. and Nash, J. (1984) Corn at interest: the extent and cost of grain storage in medieval England, American Economic Review, 74, 1, van Leeuwen, B., Földvári, P. and Pirngruber, R. (2011) Markets in pre-industrial societies: storage in Hellenistic Babylonia in the medieval English mirror. Journal of Global History (2011) 6, pp
3 footsteps and after revising the parameters of the distribution of production and used a Monte Carlo simulation to find out which type of inter-annual grain storage was feasible. Obviously, without storage you can easily calculate the chance of a hunger. If the production is normally distributed with mean 100 and std. dev. 35, and the average yield is just 40% higher than the famine limit (so the famine limit is 60). All you need is to look it up in a standard normal probability distribution table how likely it is that you obtain a value equal or lower than Since z 1.142, we are interested in the p-value of the standard normal distribution 35 at z, which is about So with this the chance of a hunger would 12.7%. The average waiting time between two famines is so roughly 7 years. This is not far from the picture suggested by historical records. If we assume that there was storage we cannot do this easily, so we need a simulation. We assume the following rules: 1. If the production is above the famine limit, 0<c<1 share of the output is stored. This is our carryover parameter. The storage cannot be negative. Also you will not store if storage would reduce your consumption below the limit. 2. If the production is below or equal the famine limit, agents will consume as much from the storage that their consumption just exceeds the famine limit. 3. Storage is assumed not to decay over time (this is a simplification). Obviously, the process is serially correlated so any outcome depends on the history of the process. What we do here is that we simulate 1000 periods 1000 times. 'wfcreate storage u series q 'output series cons 'consumption series storage=0 'storage starting level scalar co=0.05 'carryover series famine 'here we store if there was a famine, 1 means famine 0 otherwise scalar limit=70 'famine limit series chfamine 'chance of famine scalar waiting 'number of periods between two famines for!j=1 to 1000 for!i=2 to 1001 q(!i)=35*@nrnd+100 'output in period i is simulated as a normally distributed proces ' if there is no famine, you will store some of the harvest and consume the rest if q(!i)>limit*(1+co) then storage(!i)=storage(!i-1)+co*q(!i) cons(!i)=q(!i)*(1-co) ' if you have just enough you will not store this obviously never happens since we have anormal random variable and the chance that is takes the value 60 is null. Still, for theoretical consistency I include this here. if q(!i)=limit*(1+co) then storage(!i)=storage(!i-1) cons(!i)=q(!i) ' if your output was not enough then you will use up the storage provided it is enough 3
4 if q(!i)<limit*(1+co) and storage(!i-1)>=(limit-q(!i)) then storage(!i)=storage(!i-1)-(limit-q(!i)) cons(!i)=limit ' if storage is not enough you use up as much as possible if q(!i)<limit*(1+co) and storage(!i-1)>0 and storage(!i-1)<(limit-q(!i)) then storage(!i)=0 cons(!i)=q(!i)+storage(!i-1) 'if you do not have storage at all, you obviously cannot do anything if q(!i)<limit*(1+co) and storage(!i-1)=0 then storage(!i)=0 cons(!i)=q(!i) smpl famine=cons<limit smpl chfamine(!j)=@mean(famine) scalar waiting=1/@mean(chfamine)-1 I am sure that you could find an even more elegant (=parsimonious) way to write the code but this will do for now. Again, you can create a table: Average waiting time between famines under different carryovers carryover 0% 1% 2% 5% famine limit= famine limit= Note: output is assume to be a normally distributed random variable with mean 100 and standard deviation 35 Obviously, the 5% carryover can be rejected, as it would have basically eradicated society wide hunger. Depending on what we believe about the productivity of pre-industrial agriculture, however, some minor storage could have been possible. In the paper we used a production with much lower standard deviation, based on new data. This enabled us to reject the possibility of any systematic storage. 3.d. Bootstrapping Unlike in case of an MC experiment, we use the bootstrapping when we do not know the Data Generating Process, i.e., we have no idea about the underlying distribution of the data. Since we do not know the population you will create samples using your sample only by drawing with replacement. For this reason we call this procedure resampling. It is important that the new samples should have the same size and the original sample! For an example we take the data of Austrian conscripts height by Komlos. 4
5 Looking at the sample, we estimate the mean at cm. it is quite obvious that the heights are not normally distributed. This is not surprising, the data is truncated: recruitment boards anly registered those with a minimum posture. You can also run a regression: Alternatively you can have an estimate by bootstrap: matrix height=@convert(heightcm) 'here I have to convert the heights to a vector series mean for!i=1 to 1000 matrix bs=@resample(height) 'resampling with replacement mtos(bs,h) 'conversion of the vector to a variable mean(!i)=@mean(h) mean.hist Series: MEAN Sample Observations 1000 Mean Median Maximum Minimum Std. Dev Skewness Kurtosis Jarque-Bera Probability We get about the same mean and standard deviation as in case of a regression with constant. This is not surprising: whatever distribution the data has the sample mean is going to be an unbiased estimator of the expected value of the population. The distribution of the parameter estimate is going to converge to normal. So we did not gain anything from this exercise except that you now can do a bootstrap (actually bootstrapping is rarely very useful). 5
6 Let us use bootstrapping in a regression. We regress height on age: So, on average if a recruit was one year older it was 16,5 mm taller in the sample. Our 95% confidence interval for this coefficient would be: Z I used standard normal distribution instead of a Student s t because at high degrees of freedom (it is now ) it converges to standard normal. How would you do the same with bootstrapping? 1. You run the regression 2. You save the residual and the fitted value 3. You resample the residual and add it to the fitted value. This will be you dependent variable in step You rerun the regression with the variable and store the coefficient 5. You repeat steps 3 and 4 many hundred or thousand times. 6. You take the 2.5 and 97.5 percentiles of the bootstrapped coefficients. These are the limits of your confidence interval. The Eviews code is: equation eq.ls heightcm c age 'I run the regression eq.makeresids res 'I store the residuals series yfit=heightcm-res 'I store fitted values vector resm=@convert(res) 'here I have to convert the residuals to a vector series coeff for!i=1 to matrix bsres=@resample(resm) 'resampling with replacement mtos(bsres,residual) 'conversion of the vector to a variable equation eq.ls yfit+residual(1) c age 'I rerun the regression on the bootstrapped dependent var coeff(!i)=eq.@coefs(2) ' I save the coefficient 6
7 scalar scalar 0.975) After 1000 repeated resampling, I got as confidence interval. After repetitions I obtain , so having much more resampling did not improve the results much. 7
Midterm Exam 2 Answers
Economics 31 Menzie D. Chinn Fall 4 Social Sciences 7418 University of Wisconsin-Madison Midterm Exam Answers This exam is 6 minutes long. There are 8 questions on the exam be sure to check that you answer
Lectures 5 & 6: Hypothesis Testing
Lectures 5 & 6: Hypothesis Testing in which you learn to apply the concept of statistical significance to OLS estimates, learn the concept of t values, how to use them in regression work and come across
Introduction to Econometrics. Review of Probability & Statistics
1 Introduction to Econometrics Review of Probability & Statistics Peerapat Wongchaiwat, Ph.D. wongchaiwat@hotmail.com Introduction 2 What is Econometrics? Econometrics consists of the application of mathematical
Volatility. Gerald P. Dwyer. February Clemson University
Volatility Gerald P. Dwyer Clemson University February 2016 Outline 1 Volatility Characteristics of Time Series Heteroskedasticity Simpler Estimation Strategies Exponentially Weighted Moving Average Use
The Empirical Rule, z-scores, and the Rare Event Approach
Overview The Empirical Rule, z-scores, and the Rare Event Approach Look at Chebyshev s Rule and the Empirical Rule Explore some applications of the Empirical Rule How to calculate and use z-scores Introducing
Chapter 23: Inferences About Means
Chapter 3: Inferences About Means Sample of Means: number of observations in one sample the population mean (theoretical mean) sample mean (observed mean) is the theoretical standard deviation of the population
Business Analytics and Data Mining Modeling Using R Prof. Gaurav Dixit Department of Management Studies Indian Institute of Technology, Roorkee
Business Analytics and Data Mining Modeling Using R Prof. Gaurav Dixit Department of Management Studies Indian Institute of Technology, Roorkee Lecture - 04 Basic Statistics Part-1 (Refer Slide Time: 00:33)
Chapter 8 Handout: Interval Estimates and Hypothesis Testing
Chapter 8 Handout: Interval Estimates and Hypothesis esting Preview Clint s Assignment: aking Stock General Properties of the Ordinary Least Squares (OLS) Estimation Procedure Estimate Reliability: Interval
Topics for Today. Sampling Distribution. The Central Limit Theorem. Stat203 Page 1 of 28 Fall 2011 Week 5 Lecture 2
Topics for Today Sampling Distribution The Central Limit Theorem Stat203 Page 1 of 28 Law of Large Numbers Draw at from any population with mean μ. As the number of individuals, the mean,, of the sample
Lecture 30. DATA 8 Summer Regression Inference
DATA 8 Summer 2018 Lecture 30 Regression Inference Slides created by John DeNero (denero@berkeley.edu) and Ani Adhikari (adhikari@berkeley.edu) Contributions by Fahad Kamran (fhdkmrn@berkeley.edu) and
The Central Limit Theorem
The Central Limit Theorem Patrick Breheny September 27 Patrick Breheny University of Iowa Biostatistical Methods I (BIOS 5710) 1 / 31 Kerrich s experiment Introduction 10,000 coin flips Expectation and
Lecture 14. Analysis of Variance * Correlation and Regression. The McGraw-Hill Companies, Inc., 2000
Lecture 14 Analysis of Variance * Correlation and Regression Outline Analysis of Variance (ANOVA) 11-1 Introduction 11-2 Scatter Plots 11-3 Correlation 11-4 Regression Outline 11-5 Coefficient of Determination
Lecture 14. Outline. Outline. Analysis of Variance * Correlation and Regression Analysis of Variance (ANOVA)
Outline Lecture 14 Analysis of Variance * Correlation and Regression Analysis of Variance (ANOVA) 11-1 Introduction 11- Scatter Plots 11-3 Correlation 11-4 Regression Outline 11-5 Coefficient of Determination
HYPOTHESIS TESTING. Hypothesis Testing
MBA 605 Business Analytics Don Conant, PhD. HYPOTHESIS TESTING Hypothesis testing involves making inferences about the nature of the population on the basis of observations of a sample drawn from the population.
23. MORE HYPOTHESIS TESTING
23. MORE HYPOTHESIS TESTING The Logic Behind Hypothesis Testing For simplicity, consider testing H 0 : µ = µ 0 against the two-sided alternative H A : µ µ 0. Even if H 0 is true (so that the expectation
5.2 Tests of Significance
5.2 Tests of Significance Example 5.7. Diet colas use artificial sweeteners to avoid sugar. Colas with artificial sweeteners gradually lose their sweetness over time. Manufacturers therefore test new colas
Least Absolute Value vs. Least Squares Estimation and Inference Procedures in Regression Models with Asymmetric Error Distributions
Journal of Modern Applied Statistical Methods Volume 8 Issue 1 Article 13 5-1-2009 Least Absolute Value vs. Least Squares Estimation and Inference Procedures in Regression Models with Asymmetric Error
10/4/2013. Hypothesis Testing & z-test. Hypothesis Testing. Hypothesis Testing
& z-test Lecture Set 11 We have a coin and are trying to determine if it is biased or unbiased What should we assume? Why? Flip coin n = 100 times E(Heads) = 50 Why? Assume we count 53 Heads... What could
9/2/2010. Wildlife Management is a very quantitative field of study. throughout this course and throughout your career.
Introduction to Data and Analysis Wildlife Management is a very quantitative field of study Results from studies will be used throughout this course and throughout your career. Sampling design influences
CHAPTER 6: SPECIFICATION VARIABLES
Recall, we had the following six assumptions required for the Gauss-Markov Theorem: 1. The regression model is linear, correctly specified, and has an additive error term. 2. The error term has a zero
Lecture 8. Using the CLR Model
Lecture 8. Using the CLR Model Example of regression analysis. Relation between patent applications and R&D spending Variables PATENTS = No. of patents (in 1000) filed RDEXP = Expenditure on research&development
MA Advanced Econometrics: Applying Least Squares to Time Series
MA Advanced Econometrics: Applying Least Squares to Time Series Karl Whelan School of Economics, UCD February 15, 2011 Karl Whelan (UCD) Time Series February 15, 2011 1 / 24 Part I Time Series: Standard
Chapter 12 - Lecture 2 Inferences about regression coefficient
Chapter 12 - Lecture 2 Inferences about regression coefficient April 19th, 2010 Facts about slope Test Statistic Confidence interval Hypothesis testing Test using ANOVA Table Facts about slope In previous
Business Statistics. Lecture 10: Course Review
Business Statistics Lecture 10: Course Review 1 Descriptive Statistics for Continuous Data Numerical Summaries Location: mean, median Spread or variability: variance, standard deviation, range, percentiles,
The Prediction of Monthly Inflation Rate in Romania 1
Economic Insights Trends and Challenges Vol.III (LXVI) No. 2/2014 75-84 The Prediction of Monthly Inflation Rate in Romania 1 Mihaela Simionescu Institute for Economic Forecasting of the Romanian Academy,
Permutation Tests. Noa Haas Statistics M.Sc. Seminar, Spring 2017 Bootstrap and Resampling Methods
Permutation Tests Noa Haas Statistics M.Sc. Seminar, Spring 2017 Bootstrap and Resampling Methods The Two-Sample Problem We observe two independent random samples: F z = z 1, z 2,, z n independently of
Probability Methods in Civil Engineering Prof. Dr. Rajib Maity Department of Civil Engineering Indian Institution of Technology, Kharagpur
Probability Methods in Civil Engineering Prof. Dr. Rajib Maity Department of Civil Engineering Indian Institution of Technology, Kharagpur Lecture No. # 36 Sampling Distribution and Parameter Estimation
MS&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
Finite-sample quantiles of the Jarque-Bera test
Finite-sample quantiles of the Jarque-Bera test Steve Lawford Department of Economics and Finance, Brunel University First draft: February 2004. Abstract The nite-sample null distribution of the Jarque-Bera
Chapter 7 Comparison of two independent samples
Chapter 7 Comparison of two independent samples 7.1 Introduction Population 1 µ σ 1 1 N 1 Sample 1 y s 1 1 n 1 Population µ σ N Sample y s n 1, : population means 1, : population standard deviations N
The bootstrap. Patrick Breheny. December 6. The empirical distribution function The bootstrap
Patrick Breheny December 6 Patrick Breheny BST 764: Applied Statistical Modeling 1/21 The empirical distribution function Suppose X F, where F (x) = Pr(X x) is a distribution function, and we wish to estimate
Economics 471: Econometrics Department of Economics, Finance and Legal Studies University of Alabama
Economics 471: Econometrics Department of Economics, Finance and Legal Studies University of Alabama Course Packet The purpose of this packet is to show you one particular dataset and how it is used in
Chapter 24. Comparing Means
Chapter 4 Comparing Means!1 /34 Homework p579, 5, 7, 8, 10, 11, 17, 31, 3! /34 !3 /34 Objective Students test null and alternate hypothesis about two!4 /34 Plot the Data The intuitive display for comparing
+ Specify 1 tail / 2 tail
Week 2: Null hypothesis Aeroplane seat designer wonders how wide to make the plane seats. He assumes population average hip size μ = 43.2cm Sample size n = 50 Question : Is the assumption μ = 43.2cm reasonable?
ECO220Y Simple Regression: Testing the Slope
ECO220Y Simple Regression: Testing the Slope Readings: Chapter 18 (Sections 18.3-18.5) Winter 2012 Lecture 19 (Winter 2012) Simple Regression Lecture 19 1 / 32 Simple Regression Model y i = β 0 + β 1 x
Functional Form. So far considered models written in linear form. Y = b 0 + b 1 X + u (1) Implies a straight line relationship between y and X
Functional Form So far considered models written in linear form Y = b 0 + b 1 X + u (1) Implies a straight line relationship between y and X Functional Form So far considered models written in linear form
T.I.H.E. IT 233 Statistics and Probability: Sem. 1: 2013 ESTIMATION AND HYPOTHESIS TESTING OF TWO POPULATIONS
ESTIMATION AND HYPOTHESIS TESTING OF TWO POPULATIONS In our work on hypothesis testing, we used the value of a sample statistic to challenge an accepted value of a population parameter. We focused only
Hypothesis Tests and Estimation for Population Variances. Copyright 2014 Pearson Education, Inc.
Hypothesis Tests and Estimation for Population Variances 11-1 Learning Outcomes Outcome 1. Formulate and carry out hypothesis tests for a single population variance. Outcome 2. Develop and interpret confidence
LECTURE 5. Introduction to Econometrics. Hypothesis testing
LECTURE 5 Introduction to Econometrics Hypothesis testing October 18, 2016 1 / 26 ON TODAY S LECTURE We are going to discuss how hypotheses about coefficients can be tested in regression models We will
TESTING FOR NORMALITY IN THE LINEAR REGRESSION MODEL: AN EMPIRICAL LIKELIHOOD RATIO TEST
Econometrics Working Paper EWP0402 ISSN 1485-6441 Department of Economics TESTING FOR NORMALITY IN THE LINEAR REGRESSION MODEL: AN EMPIRICAL LIKELIHOOD RATIO TEST Lauren Bin Dong & David E. A. Giles Department
CHAPTER 5 FUNCTIONAL FORMS OF REGRESSION MODELS
CHAPTER 5 FUNCTIONAL FORMS OF REGRESSION MODELS QUESTIONS 5.1. (a) In a log-log model the dependent and all explanatory variables are in the logarithmic form. (b) In the log-lin model the dependent variable
Review of Statistics
Review of Statistics Topics Descriptive Statistics Mean, Variance Probability Union event, joint event Random Variables Discrete and Continuous Distributions, Moments Two Random Variables Covariance and
Solutions: Monday, October 15
Amherst College Department of Economics Economics 360 Fall 2012 1. Consider Nebraska petroleum consumption. Solutions: Monday, October 15 Petroleum Consumption Data for Nebraska: Annual time series data
Hypothesis testing. Data to decisions
Hypothesis testing Data to decisions The idea Null hypothesis: H 0 : the DGP/population has property P Under the null, a sample statistic has a known distribution If, under that that distribution, the
Last week: Sample, population and sampling distributions finished with estimation & confidence intervals
Past weeks: Measures of central tendency (mean, mode, median) Measures of dispersion (standard deviation, variance, range, etc). Working with the normal curve Last week: Sample, population and sampling
1/24/2008. Review of Statistical Inference. C.1 A Sample of Data. C.2 An Econometric Model. C.4 Estimating the Population Variance and Other Moments
/4/008 Review of Statistical Inference Prepared by Vera Tabakova, East Carolina University C. A Sample of Data C. An Econometric Model C.3 Estimating the Mean of a Population C.4 Estimating the Population
Section 9.4. Notation. Requirements. Definition. Inferences About Two Means (Matched Pairs) Examples
Objective Section 9.4 Inferences About Two Means (Matched Pairs) Compare of two matched-paired means using two samples from each population. Hypothesis Tests and Confidence Intervals of two dependent means
STAT440/840: Statistical Computing
First Prev Next Last STAT440/840: Statistical Computing Paul Marriott pmarriott@math.uwaterloo.ca MC 6096 February 2, 2005 Page 1 of 41 First Prev Next Last Page 2 of 41 Chapter 3: Data resampling: the
ECON Introductory Econometrics. Lecture 2: Review of Statistics
ECON415 - Introductory Econometrics Lecture 2: Review of Statistics Monique de Haan (moniqued@econ.uio.no) Stock and Watson Chapter 2-3 Lecture outline 2 Simple random sampling Distribution of the sample
Hypothesis testing I. - In particular, we are talking about statistical hypotheses. [get everyone s finger length!] n =
Hypothesis testing I I. What is hypothesis testing? [Note we re temporarily bouncing around in the book a lot! Things will settle down again in a week or so] - Exactly what it says. We develop a hypothesis,
Review of Statistics 101
Review of Statistics 101 We review some important themes from the course 1. Introduction Statistics- Set of methods for collecting/analyzing data (the art and science of learning from data). Provides methods
9. Linear Regression and Correlation
9. Linear Regression and Correlation Data: y a quantitative response variable x a quantitative explanatory variable (Chap. 8: Recall that both variables were categorical) For example, y = annual income,
Quantitative Methods for Economics, Finance and Management (A86050 F86050)
Quantitative Methods for Economics, Finance and Management (A86050 F86050) Matteo Manera matteo.manera@unimib.it Marzio Galeotti marzio.galeotti@unimi.it 1 This material is taken and adapted from Guy Judge
Last two weeks: Sample, population and sampling distributions finished with estimation & confidence intervals
Past weeks: Measures of central tendency (mean, mode, median) Measures of dispersion (standard deviation, variance, range, etc). Working with the normal curve Last two weeks: Sample, population and sampling
Mathematics for Economics MA course
Mathematics for Economics MA course Simple Linear Regression Dr. Seetha Bandara Simple Regression Simple linear regression is a statistical method that allows us to summarize and study relationships between
4. Examples. Results: Example 4.1 Implementation of the Example 3.1 in SAS. In SAS we can use the Proc Model procedure.
4. Examples Example 4.1 Implementation of the Example 3.1 in SAS. In SAS we can use the Proc Model procedure. Simulate data from t-distribution with ν = 6. SAS: data tdist; do i = 1 to 500; y = tinv(ranuni(158),6);
Keppel, G. & Wickens, T.D. Design and Analysis Chapter 2: Sources of Variability and Sums of Squares
Keppel, G. & Wickens, T.D. Design and Analysis Chapter 2: Sources of Variability and Sums of Squares K&W introduce the notion of a simple experiment with two conditions. Note that the raw data (p. 16)
Longest runs in coin tossing. Teaching recursive formulae, asymptotic theorems and computer simulations
tmcs-karacsony 2011/11/20 19:51 page 261 #1 9/2 (2011), 261 274 Longest runs in coin tossing. Teaching recursive formulae, asymptotic theorems and computer simulations Zsolt Karácsony and Józsefné Libor
Bootstrapping, Permutations, and Monte Carlo Testing
Bootstrapping, Permutations, and Monte Carlo Testing Problem: Population of interest is extremely rare spatially and you are interested in using a 95% CI to estimate total abundance. The sampling design
Ch Inference for Linear Regression
Ch. 12-1 Inference for Linear Regression ACT = 6.71 + 5.17(GPA) For every increase of 1 in GPA, we predict the ACT score to increase by 5.17. population regression line β (true slope) μ y = α + βx mean
t-test for b Copyright 2000 Tom Malloy. All rights reserved. Regression
t-test for b Copyright 2000 Tom Malloy. All rights reserved. Regression Recall, back some time ago, we used a descriptive statistic which allowed us to draw the best fit line through a scatter plot. We
Multivariate Regression Model Results
Updated: August, 0 Page of Multivariate Regression Model Results 4 5 6 7 8 This exhibit provides the results of the load model forecast discussed in Schedule. Included is the forecast of short term system
Lecture 8. Using the CLR Model. Relation between patent applications and R&D spending. Variables
Lecture 8. Using the CLR Model Relation between patent applications and R&D spending Variables PATENTS = No. of patents (in 000) filed RDEP = Expenditure on research&development (in billions of 99 $) The
Eco517 Fall 2004 C. Sims MIDTERM EXAM
Eco517 Fall 2004 C. Sims MIDTERM EXAM Answer all four questions. Each is worth 23 points. Do not devote disproportionate time to any one question unless you have answered all the others. (1) We are considering
13.7 ANOTHER TEST FOR TREND: KENDALL S TAU
13.7 ANOTHER TEST FOR TREND: KENDALL S TAU In 1969 the U.S. government instituted a draft lottery for choosing young men to be drafted into the military. Numbers from 1 to 366 were randomly assigned to
Lecture 8: Functional Form
Lecture 8: Functional Form What we know now OLS - fitting a straight line y = b 0 + b 1 X through the data using the principle of choosing the straight line that minimises the sum of squared residuals
Mathematical Notation Math Introduction to Applied Statistics
Mathematical Notation Math 113 - Introduction to Applied Statistics Name : Use Word or WordPerfect to recreate the following documents. Each article is worth 10 points and should be emailed to the instructor
L21: Chapter 12: Linear regression
L21: Chapter 12: Linear regression Department of Statistics, University of South Carolina Stat 205: Elementary Statistics for the Biological and Life Sciences 1 / 37 So far... 12.1 Introduction One sample
The regression model with one stochastic regressor.
The regression model with one stochastic regressor. 3150/4150 Lecture 6 Ragnar Nymoen 30 January 2012 We are now on Lecture topic 4 The main goal in this lecture is to extend the results of the regression
MONT 105Q Mathematical Journeys Lecture Notes on Statistical Inference and Hypothesis Testing March-April, 2016
MONT 105Q Mathematical Journeys Lecture Notes on Statistical Inference and Hypothesis Testing March-April, 2016 Background For the purposes of this introductory discussion, imagine first that you are a
Robustness of Simultaneous Estimation Methods to Varying Degrees of Correlation Between Pairs of Random Deviates
Global Journal of Mathematical Sciences: Theory and Practical. Volume, Number 3 (00), pp. 5--3 International Research Publication House http://www.irphouse.com Robustness of Simultaneous Estimation Methods
Inference for the mean of a population. Testing hypotheses about a single mean (the one sample t-test). The sign test for matched pairs
Stat 528 (Autumn 2008) Inference for the mean of a population (One sample t procedures) Reading: Section 7.1. Inference for the mean of a population. The t distribution for a normal population. Small sample
Hypothesis testing Goodness of fit Multicollinearity Prediction. Applied Statistics. Lecturer: Serena Arima
Applied Statistics Lecturer: Serena Arima Hypothesis testing for the linear model Under the Gauss-Markov assumptions and the normality of the error terms, we saw that β N(β, σ 2 (X X ) 1 ) and hence s
Section 3: Simple Linear Regression
Section 3: Simple Linear Regression Carlos M. Carvalho The University of Texas at Austin McCombs School of Business http://faculty.mccombs.utexas.edu/carlos.carvalho/teaching/ 1 Regression: General Introduction
Econometrics Homework 1
Econometrics Homework Due Date: March, 24. by This problem set includes questions for Lecture -4 covered before midterm exam. Question Let z be a random column vector of size 3 : z = @ (a) Write out z
determine whether or not this relationship is.
Section 9-1 Correlation A correlation is a between two. The data can be represented by ordered pairs (x,y) where x is the (or ) variable and y is the (or ) variable. There are several types of correlations
-However, this definition can be expanded to include: biology (biometrics), environmental science (environmetrics), economics (econometrics).
Chemometrics Application of mathematical, statistical, graphical or symbolic methods to maximize chemical information. -However, this definition can be expanded to include: biology (biometrics), environmental
Physics 509: Bootstrap and Robust Parameter Estimation
Physics 509: Bootstrap and Robust Parameter Estimation Scott Oser Lecture #20 Physics 509 1 Nonparametric parameter estimation Question: what error estimate should you assign to the slope and intercept
Hint: The following equation converts Celsius to Fahrenheit: F = C where C = degrees Celsius F = degrees Fahrenheit
Amherst College Department of Economics Economics 360 Fall 2014 Exam 1: Solutions 1. (10 points) The following table in reports the summary statistics for high and low temperatures in Key West, FL from
Lecture 6a: Unit Root and ARIMA Models
Lecture 6a: Unit Root and ARIMA Models 1 2 Big Picture A time series is non-stationary if it contains a unit root unit root nonstationary The reverse is not true. For example, y t = cos(t) + u t has no
Outline. PubH 5450 Biostatistics I Prof. Carlin. Confidence Interval for the Mean. Part I. Reviews
Outline Outline PubH 5450 Biostatistics I Prof. Carlin Lecture 11 Confidence Interval for the Mean Known σ (population standard deviation): Part I Reviews σ x ± z 1 α/2 n Small n, normal population. Large
Chapter 23. Inference About Means
Chapter 23 Inference About Means 1 /57 Homework p554 2, 4, 9, 10, 13, 15, 17, 33, 34 2 /57 Objective Students test null and alternate hypotheses about a population mean. 3 /57 Here We Go Again Now that
Inferences for Regression
Inferences for Regression An Example: Body Fat and Waist Size Looking at the relationship between % body fat and waist size (in inches). Here is a scatterplot of our data set: Remembering Regression In
Institute of Actuaries of India
Institute of Actuaries of India Subject CT3 Probability & Mathematical Statistics May 2011 Examinations INDICATIVE SOLUTION Introduction The indicative solution has been written by the Examiners with the
Section 1 : Introduction to the Potential Outcomes Framework. Andrew Bertoli 4 September 2013
Section 1 : Introduction to the Potential Outcomes Framework Andrew Bertoli 4 September 2013 Roadmap 1. Preview 2. Helpful Tips 3. Potential Outcomes Framework 4. Experiments vs. Observational Studies
Inferential statistics
Inferential statistics Inference involves making a Generalization about a larger group of individuals on the basis of a subset or sample. Ahmed-Refat-ZU Null and alternative hypotheses In hypotheses testing,
The t-statistic. Student s t Test
The t-statistic 1 Student s t Test When the population standard deviation is not known, you cannot use a z score hypothesis test Use Student s t test instead Student s t, or t test is, conceptually, very
LECTURE 15: SIMPLE LINEAR REGRESSION I
David Youngberg BSAD 20 Montgomery College LECTURE 5: SIMPLE LINEAR REGRESSION I I. From Correlation to Regression a. Recall last class when we discussed two basic types of correlation (positive and negative).
Biostatistics Quantitative Data
Biostatistics Quantitative Data Descriptive Statistics Statistical Models One-sample and Two-Sample Tests Introduction to SAS-ANALYST T- and Rank-Tests using ANALYST Thomas Scheike Quantitative Data This
Econometrics. 4) Statistical inference
30C00200 Econometrics 4) Statistical inference Timo Kuosmanen Professor, Ph.D. http://nomepre.net/index.php/timokuosmanen Today s topics Confidence intervals of parameter estimates Student s t-distribution
Exercise sheet 3 The Multiple Regression Model
Exercise sheet 3 The Multiple Regression Model Note: In those problems that include estimations and have a reference to a data set the students should check the outputs obtained with Gretl. 1. Let the
If we want to analyze experimental or simulated data we might encounter the following tasks:
Chapter 1 Introduction If we want to analyze experimental or simulated data we might encounter the following tasks: Characterization of the source of the signal and diagnosis Studying dependencies Prediction
Semester 2, 2015/2016
ECN 3202 APPLIED ECONOMETRICS 2. Simple linear regression B Mr. Sydney Armstrong Lecturer 1 The University of Guyana 1 Semester 2, 2015/2016 PREDICTION The true value of y when x takes some particular
1 Quantitative Techniques in Practice
1 Quantitative Techniques in Practice 1.1 Lecture 2: Stationarity, spurious regression, etc. 1.1.1 Overview In the rst part we shall look at some issues in time series economics. In the second part we
Multiple Regression Theory 2006 Samuel L. Baker
MULTIPLE REGRESSION THEORY 1 Multiple Regression Theory 2006 Samuel L. Baker Multiple regression is regression with two or more independent variables on the right-hand side of the equation. Use multiple
y ˆ i = ˆ " T u i ( i th fitted value or i th fit)
1 2 INFERENCE FOR MULTIPLE LINEAR REGRESSION Recall Terminology: p predictors x 1, x 2,, x p Some might be indicator variables for categorical variables) k-1 non-constant terms u 1, u 2,, u k-1 Each u
Sampling Distributions: Central Limit Theorem
Review for Exam 2 Sampling Distributions: Central Limit Theorem Conceptually, we can break up the theorem into three parts: 1. The mean (µ M ) of a population of sample means (M) is equal to the mean (µ)
MORE ON SIMPLE REGRESSION: OVERVIEW
FI=NOT0106 NOTICE. Unless otherwise indicated, all materials on this page and linked pages at the blue.temple.edu address and at the astro.temple.edu address are the sole property of Ralph B. Taylor and
Chapter 5: HYPOTHESIS TESTING
MATH411: Applied Statistics Dr. YU, Chi Wai Chapter 5: HYPOTHESIS TESTING 1 WHAT IS HYPOTHESIS TESTING? As its name indicates, it is about a test of hypothesis. To be more precise, we would first translate
Chapter 9 Inferences from Two Samples
Chapter 9 Inferences from Two Samples 9-1 Review and Preview 9-2 Two Proportions 9-3 Two Means: Independent Samples 9-4 Two Dependent Samples (Matched Pairs) 9-5 Two Variances or Standard Deviations Review