STAT 5302 Applied Regression Analysis. Hawkins

Size: px
Start display at page:

Download "STAT 5302 Applied Regression Analysis. Hawkins"

Transcription

1 Homework 3 sample solution 1. MinnLand data STAT 5302 Applied Regression Analysis. Hawkins newdata <- subset(minnland, year == 2010) fit1 <- with(newdata, lm(acreprice ~ region-1)) summary(fit1) regionnorthwest <2e-16 *** regionwest Central <2e-16 *** regioncentral <2e-16 *** regionsouth West <2e-16 *** regionsouth Central <2e-16 *** regionsouth East <2e-16 *** Residual standard error: 1330 on 1817 degrees of freedom Multiple R-squared: , Adjusted R-squared: F-statistic: 2496 on 6 and 1817 DF, p-value: < 2.2e-16 This initial no-intercept call gives the 2010 mean land price in each of the six regions. These range from $1509 in the Northwest to $4791 in the South Central. The regression output also gives a standard error for each of these means, something that we did not draw attention to in class. The fit with the intercept gives fit2 <- with(newdata, lm(acreprice ~ region)) summary(fit2) (Intercept) <2e-16 *** regionwest Central <2e-16 *** regioncentral <2e-16 *** regionsouth West <2e-16 *** regionsouth Central <2e-16 *** regionsouth East <2e-16 *** Residual standard error: 1330 on 1817 degrees of freedom Multiple R-squared: , Adjusted R-squared: F-statistic: on 5 and 1817 DF, p-value: < 2.2e-16 In this output, the intercept, $1509, is the mean price in the reference group Northwest. The coefficients of the subsequent terms give the difference between the mean prince in each region and this reference group So, for example, the Central region had a mean of $4074 in fit1. This is ( ) = $2565 higher than the reference group, a value that we see as the coefficient of Central in the second fit. R also reports a standard error, a t and a P value. These reflect the test of equality of the means of Central and Northwest; the enormous t value, 25.87, reflects that the two regions are highly significantly different. Changing the reference group using STAT 5302 Homework 3 Page 1

2 region2 <- with(newdata, relevel(region, "Central")) fit3 <- lm(newdata$acreprice ~ region2) summary(fit3) gives (Intercept) < 2e-16 *** region2northwest < 2e-16 *** region2west Central < 2e-16 *** region2south West e-05 *** region2south Central e-11 *** region2south East e-05 *** Residual standard error: 1330 on 1817 degrees of freedom Multiple R-squared: , Adjusted R-squared: F-statistic: on 5 and 1817 DF, p-value: < 2.2e-16 This fit and the preceding one have much in common the last three summary lines of output are identical, reflecting that this is the same fit expressed differently. The two sets of coefficients are also equivalent. Taking three groups, Northwest, with a mean of $1509, Central, with a mean of 4075 and say South West with a mean of $4554, in the earlier fit the coefficient was $ ) = $3045, but changing the reference group to Central changes the coefficient to $479. Note that in other regions all differ significantly from Central. fit4 <- with(newdata, lm(acreprice ~ region+productivity)) summary(fit4) (Intercept) regionwest Central e-11 *** regioncentral < 2e-16 *** regionsouth West < 2e-16 *** regionsouth Central < 2e-16 *** regionsouth East < 2e-16 *** productivity < 2e-16 *** --- Signif. codes: 0 *** ** 0.01 * Residual standard error: 1087 on 921 degrees of freedom (895 observations deleted due to missingness) Multiple R-squared: , Adjusted R-squared: F-statistic: on 6 and 921 DF, p-value: < 2.2e-16 This regression shows that, on aggregate across the regions, productivity has a highly significant effect on the land price, but the significance of the coefficients of the regions shows that it is not the root cause of the pricing differences. fit5 <- with(newdata, lm(acreprice ~region+productivity+region:productivity)) summary(fit5) Use fit5 to get the regression of acreprice on productivity separately in each region. Adding up the coefficients, we get: STAT 5302 Homework 3 Page 2

3 Region Intercept Slope Northwest West Central Central South West South Central South East (As a crude direct check, on this, running chekker <- subset(newdata, region=="central") fit6 <- with(chekker, lm(acreprice ~ productivity)) summary(fit6) gives (Intercept) ** productivity ***) The Type I analysis of variance of fit5 using anova(fit5) gives Analysis of Variance Table Response: acreprice Df Sum Sq Mean Sq F value Pr(>F) region < 2.2e-16 *** productivity < 2.2e-16 *** region:productivity e-11 *** Residuals In this problem, it is not particularly helpful; it tells us that the effect of productivity differs enormously between regions. STAT 5302 Homework 3 Page 3

4 2. Rolling the two steps into one, we get The plot before adding the smooth did not look wonderful, and this is reinforced when we add the smooth. The relationship seems to decelerate as we get to the higher ages. fit6 <- with(lakemary, lm(length ~ Age )) summary(fit6) (Intercept) <2e-16 *** Age <2e-16 *** Residual standard error: on 76 degrees of freedom Multiple R-squared: , Adjusted R-squared: F-statistic: on 1 and 76 DF, p-value: < 2.2e-16 Quite a strong regression fit7 <- update(fit6, ~.+factor(age)) summary(fit7) (1 not defined because of singularities) STAT 5302 Homework 3 Page 4

5 (Intercept) e-05 *** Age e-11 *** factor(age) factor(age) *** factor(age) *** factor(age) factor(age)6 NA NA NA NA Residual standard error: on 72 degrees of freedom Multiple R-squared: , Adjusted R-squared: F-statistic: 58.9 on 5 and 72 DF, p-value: < 2.2e-16 This fit does not look particularly happy. If the regression were truly linear in Age, then adding factor(age) to the fit6 model should not produce anything significant. But actually, we see striking significances on both ages 3 and 4, in the middle of the age range. This is our indicator that the visual impression of a decelerating trend was right, and the linear model does not hold. Now adding the square and cube of age lakemary <- transform(lakemary, Agesq = Age^2, Agecu = Age^3) fit8 <- update(fit6, ~.+Agesq) (Intercept) Age e-12 *** Agesq e-06 *** Residual standard error: on 75 degrees of freedom Multiple R-squared: , Adjusted R-squared: F-statistic: on 2 and 75 DF, p-value: < 2.2e-16 Agesq is highly significant in this regression, confirming the lack of fit of the linear model fit9 <- update(fit8, ~.+Agecu) (Intercept) Age ** Agesq Agecu Signif. codes: 0 *** ** 0.01 * Residual standard error: on 74 degrees of freedom Multiple R-squared: , Adjusted R-squared: F-statistic: on 3 and 74 DF, p-value: < 2.2e-16 But Agecu is not significant in this model. So if we want to model with a polynomial nop higher than cubic, a quadratic does it. The sequential analysis of variance is our most compact and efficient way of deciding what degree polynomial to use. anova(fit9) STAT 5302 Homework 3 Page 5

6 Analysis of Variance Table Response: Length Df Sum Sq Mean Sq F value Pr(>F) Age < 2.2e-16 *** Agesq e-06 *** Agecu Residuals Trimming from the bottom, we don t need a cube, but we do need a square. So the model would be Length = Age 4.72 Age 2 3. First step was to recode the variables. This is not strictly necessary we could fit the QRSM using the original units, as the book does, but I believe using the coded units makes the interpretation easier. So: cakes <- transform(cakes, x1=(x1-35)/2, x2=(x2-350)/10) cakes <- transform(cakes, x1sq=x1^2, x2sq=x2^2, prod=x1*x2) The QRSM gives (Intercept) e-11 *** x ** x *** x1sq ** x2sq e-05 *** prod ** Residual standard error: on 8 degrees of freedom Multiple R-squared: , Adjusted R-squared: F-statistic: 29.6 on 5 and 8 DF, p-value: 5.864e-05 Turning to the questions: There is indeed significant curvature. Both squared terms are highly significant There is an interaction between the two. The flat spot is a maximum, as both squared term coefficients are negative. To get the estimate of the location of the flat spot, we need to solve x x x , giving x Plugging this in, we get the predicted y at the flat spot of For completeness, here is a contour plot of the fitted QRSM (done in R, but not discussed in the text.) STAT 5302 Homework 3 Page 6

7 STAT 5302 Homework 3 Page 7

PSYC 6140 November 16, 2005 ANOVA output in R

PSYC 6140 November 16, 2005 ANOVA output in R PSYC 6140 November 16, 2005 ANOVA output in R Type I, Type II and Type III Sums of Squares are displayed in ANOVA tables in a mumber of packages. The car library in R makes these available in R. This handout

More information

To: Professor Roger Bohn & Hyeonsu Kang Subject: Big Data, Assignment April 13th. From: xxxx (anonymized) Date: 4/11/2016

To: Professor Roger Bohn & Hyeonsu Kang Subject: Big Data, Assignment April 13th. From: xxxx (anonymized) Date: 4/11/2016 To: Professor Roger Bohn & Hyeonsu Kang Subject: Big Data, Assignment April 13th. From: xxxx (anonymized) Date: 4/11/2016 Data Preparation: 1. Separate trany variable into Manual which takes value of 1

More information

Wine-Tasting by Numbers: Using Binary Logistic Regression to Reveal the Preferences of Experts

Wine-Tasting by Numbers: Using Binary Logistic Regression to Reveal the Preferences of Experts Wine-Tasting by Numbers: Using Binary Logistic Regression to Reveal the Preferences of Experts When you need to understand situations that seem to defy data analysis, you may be able to use techniques

More information

The R survey package used in these examples is version 3.22 and was run under R v2.7 on a PC.

The R survey package used in these examples is version 3.22 and was run under R v2.7 on a PC. CHAPTER 7 ANALYSIS EXAMPLES REPLICATION-R SURVEY PACKAGE 3.22 GENERAL NOTES ABOUT ANALYSIS EXAMPLES REPLICATION These examples are intended to provide guidance on how to use the commands/procedures for

More information

Homework 1 - Solutions. Problem 2

Homework 1 - Solutions. Problem 2 Homework 1 - Solutions Problem 2 a) soma

More information

AJAE Appendix: Testing Household-Specific Explanations for the Inverse Productivity Relationship

AJAE Appendix: Testing Household-Specific Explanations for the Inverse Productivity Relationship AJAE Appendix: Testing Household-Specific Explanations for the Inverse Productivity Relationship Juliano Assunção Department of Economics PUC-Rio Luis H. B. Braido Graduate School of Economics Getulio

More information

Missing Data Treatments

Missing Data Treatments Missing Data Treatments Lindsey Perry EDU7312: Spring 2012 Presentation Outline Types of Missing Data Listwise Deletion Pairwise Deletion Single Imputation Methods Mean Imputation Hot Deck Imputation Multiple

More information

wine 1 wine 2 wine 3 person person person person person

wine 1 wine 2 wine 3 person person person person person 1. A trendy wine bar set up an experiment to evaluate the quality of 3 different wines. Five fine connoisseurs of wine were asked to taste each of the wine and give it a rating between 0 and 10. The order

More information

Analysis of Things (AoT)

Analysis of Things (AoT) Analysis of Things (AoT) Big Data & Machine Learning Applied to Brent Crude Executive Summary Data Selecting & Visualising Data We select historical, monthly, fundamental data We check for correlations

More information

BORDEAUX WINE VINTAGE QUALITY AND THE WEATHER ECONOMETRIC ANALYSIS

BORDEAUX WINE VINTAGE QUALITY AND THE WEATHER ECONOMETRIC ANALYSIS BORDEAUX WINE VINTAGE QUALITY AND THE WEATHER ECONOMETRIC ANALYSIS WINE PRICES OVER VINTAGES DATA The data sheet contains market prices for a collection of 13 high quality Bordeaux wines (not including

More information

Comparing R print-outs from LM, GLM, LMM and GLMM

Comparing R print-outs from LM, GLM, LMM and GLMM 3. Inference: interpretation of results, plotting results, confidence intervals, hypothesis tests (Wald,LRT). 4. Asymptotic distribution of maximum likelihood estimators and tests. 5. Checking the adequacy

More information

From VOC to IPA: This Beer s For You!

From VOC to IPA: This Beer s For You! From VOC to IPA: This Beer s For You! Joel Smith Statistician Minitab Inc. jsmith@minitab.com 2013 Minitab, Inc. Image courtesy of amazon.com The Data Online beer reviews Evaluated overall and: Appearance

More information

THE STATISTICAL SOMMELIER

THE STATISTICAL SOMMELIER THE STATISTICAL SOMMELIER An Introduction to Linear Regression 15.071 The Analytics Edge Bordeaux Wine Large differences in price and quality between years, although wine is produced in a similar way Meant

More information

Online Appendix to. Are Two heads Better Than One: Team versus Individual Play in Signaling Games. David C. Cooper and John H.

Online Appendix to. Are Two heads Better Than One: Team versus Individual Play in Signaling Games. David C. Cooper and John H. Online Appendix to Are Two heads Better Than One: Team versus Individual Play in Signaling Games David C. Cooper and John H. Kagel This appendix contains a discussion of the robustness of the regression

More information

Weather Sensitive Adjustment Using the WSA Factor Method

Weather Sensitive Adjustment Using the WSA Factor Method Weather Sensitive Adjustment Using the WSA Factor Method Alternate CBLs can be requested within elrs to compensate for the temperature differences between the CBL basis days and the temperature of the

More information

Missing value imputation in SAS: an intro to Proc MI and MIANALYZE

Missing value imputation in SAS: an intro to Proc MI and MIANALYZE Victoria SAS Users Group November 26, 2013 Missing value imputation in SAS: an intro to Proc MI and MIANALYZE Sylvain Tremblay SAS Canada Education Copyright 2010 SAS Institute Inc. All rights reserved.

More information

Curtis Miller MATH 3080 Final Project pg. 1. The first question asks for an analysis on car data. The data was collected from the Kelly

Curtis Miller MATH 3080 Final Project pg. 1. The first question asks for an analysis on car data. The data was collected from the Kelly Curtis Miller MATH 3080 Final Project pg. 1 Curtis Miller 4/10/14 MATH 3080 Final Project Problem 1: Car Data The first question asks for an analysis on car data. The data was collected from the Kelly

More information

INSTITUTE AND FACULTY OF ACTUARIES CURRICULUM 2019 SPECIMEN SOLUTIONS. Subject CS1B Actuarial Statistics

INSTITUTE AND FACULTY OF ACTUARIES CURRICULUM 2019 SPECIMEN SOLUTIONS. Subject CS1B Actuarial Statistics INSTITUTE AND FACULTY OF ACTUARIES CURRICULUM 2019 SPECIMEN SOLUTIONS Subject CS1B Actuarial Statistics Question 1 (i) # Data entry before

More information

Thought: The Great Coffee Experiment

Thought: The Great Coffee Experiment Thought: The Great Coffee Experiment 7/7/16 By Kevin DeLuca ThoughtBurner Opportunity Cost of Reading this ThoughtBurner post: $1.97 about 8.95 minutes I drink a lot of coffee. In fact, I m drinking a

More information

Gail E. Potter, Timo Smieszek, and Kerstin Sailer. April 24, 2015

Gail E. Potter, Timo Smieszek, and Kerstin Sailer. April 24, 2015 Supplementary Material to Modelling workplace contact networks: the effects of organizational structure, architecture, and reporting errors on epidemic predictions, published in Network Science Gail E.

More information

The Development of a Weather-based Crop Disaster Program

The Development of a Weather-based Crop Disaster Program The Development of a Weather-based Crop Disaster Program Eric Belasco Montana State University 2016 SCC-76 Conference Pensacola, FL March 19, 2016. Belasco March 2016 1 / 18 Motivation Recent efforts to

More information

Final Exam Financial Data Analysis (6 Credit points/imp Students) March 2, 2006

Final Exam Financial Data Analysis (6 Credit points/imp Students) March 2, 2006 Dr. Roland Füss Winter Term 2005/2006 Final Exam Financial Data Analysis (6 Credit points/imp Students) March 2, 2006 Note the following important information: 1. The total disposal time is 60 minutes.

More information

Credit Supply and Monetary Policy: Identifying the Bank Balance-Sheet Channel with Loan Applications. Web Appendix

Credit Supply and Monetary Policy: Identifying the Bank Balance-Sheet Channel with Loan Applications. Web Appendix Credit Supply and Monetary Policy: Identifying the Bank Balance-Sheet Channel with Loan Applications By GABRIEL JIMÉNEZ, STEVEN ONGENA, JOSÉ-LUIS PEYDRÓ, AND JESÚS SAURINA Web Appendix APPENDIX A -- NUMBER

More information

IT 403 Project Beer Advocate Analysis

IT 403 Project Beer Advocate Analysis 1. Exploratory Data Analysis (EDA) IT 403 Project Beer Advocate Analysis Beer Advocate is a membership-based reviews website where members rank different beers based on a wide number of categories. The

More information

Business Statistics /82 Spring 2011 Booth School of Business The University of Chicago Final Exam

Business Statistics /82 Spring 2011 Booth School of Business The University of Chicago Final Exam Business Statistics 41000-81/82 Spring 2011 Booth School of Business The University of Chicago Final Exam Name You may use a calculator and two cheat sheets. You have 3 hours. I pledge my honor that I

More information

Handling Missing Data. Ashley Parker EDU 7312

Handling Missing Data. Ashley Parker EDU 7312 Handling Missing Data Ashley Parker EDU 7312 Presentation Outline Types of Missing Data Treatments for Handling Missing Data Deletion Techniques Listwise Deletion Pairwise Deletion Single Imputation Techniques

More information

Relationships Among Wine Prices, Ratings, Advertising, and Production: Examining a Giffen Good

Relationships Among Wine Prices, Ratings, Advertising, and Production: Examining a Giffen Good Relationships Among Wine Prices, Ratings, Advertising, and Production: Examining a Giffen Good Carol Miu Massachusetts Institute of Technology Abstract It has become increasingly popular for statistics

More information

Valuation in the Life Settlements Market

Valuation in the Life Settlements Market Valuation in the Life Settlements Market New Empirical Evidence Jiahua (Java) Xu 1 1 Institute of Insurance Economics University of St.Gallen Western Risk and Insurance Association 2018 Annual Meeting

More information

The multivariate piecewise linear growth model for ZHeight and zbmi can be expressed as:

The multivariate piecewise linear growth model for ZHeight and zbmi can be expressed as: Bi-directional relationships between body mass index and height from three to seven years of age: an analysis of children in the United Kingdom Millennium Cohort Study Supplementary material The multivariate

More information

Summary of Main Points

Summary of Main Points 1 Model Selection in Logistic Regression Summary of Main Points Recall that the two main objectives of regression modeling are: Estimate the effect of one or more covariates while adjusting for the possible

More information

Regression Models for Saffron Yields in Iran

Regression Models for Saffron Yields in Iran Regression Models for Saffron ields in Iran Sanaeinejad, S.H., Hosseini, S.N 1 Faculty of Agriculture, Ferdowsi University of Mashhad, Iran sanaei_h@yahoo.co.uk, nasir_nbm@yahoo.com, Abstract: Saffron

More information

Zeitschrift für Soziologie, Jg., Heft 5, 2015, Online- Anhang

Zeitschrift für Soziologie, Jg., Heft 5, 2015, Online- Anhang I Are Joiners Trusters? A Panel Analysis of Participation and Generalized Trust Online Appendix Katrin Botzen University of Bern, Institute of Sociology, Fabrikstrasse 8, 3012 Bern, Switzerland; katrin.botzen@soz.unibe.ch

More information

Missing Data Methods (Part I): Multiple Imputation. Advanced Multivariate Statistical Methods Workshop

Missing Data Methods (Part I): Multiple Imputation. Advanced Multivariate Statistical Methods Workshop Missing Data Methods (Part I): Multiple Imputation Advanced Multivariate Statistical Methods Workshop University of Georgia: Institute for Interdisciplinary Research in Education and Human Development

More information

Gasoline Empirical Analysis: Competition Bureau March 2005

Gasoline Empirical Analysis: Competition Bureau March 2005 Gasoline Empirical Analysis: Update of Four Elements of the January 2001 Conference Board study: "The Final Fifteen Feet of Hose: The Canadian Gasoline Industry in the Year 2000" Competition Bureau March

More information

Activity 10. Coffee Break. Introduction. Equipment Required. Collecting the Data

Activity 10. Coffee Break. Introduction. Equipment Required. Collecting the Data . Activity 10 Coffee Break Economists often use math to analyze growth trends for a company. Based on past performance, a mathematical equation or formula can sometimes be developed to help make predictions

More information

Predicting Wine Quality

Predicting Wine Quality March 8, 2016 Ilker Karakasoglu Predicting Wine Quality Problem description: You have been retained as a statistical consultant for a wine co-operative, and have been asked to analyze these data. Each

More information

Climate change may alter human physical activity patterns

Climate change may alter human physical activity patterns In the format provided by the authors and unedited. SUPPLEMENTARY INFORMATION VOLUME: 1 ARTICLE NUMBER: 0097 Climate change may alter human physical activity patterns Nick Obradovich and James H. Fowler

More information

Investment Wines. - Risk Analysis. Prepared by: Michael Shortell & Adiam Woldetensae Date: 06/09/2015

Investment Wines. - Risk Analysis. Prepared by: Michael Shortell & Adiam Woldetensae Date: 06/09/2015 Investment Wines - Risk Analysis Prepared by: Michael Shortell & Adiam Woldetensae Date: 06/09/2015 Purpose Look at investment wines & examine factors that affect wine prices over time We will identify

More information

Panel A: Treated firm matched to one control firm. t + 1 t + 2 t + 3 Total CFO Compensation 5.03% 0.84% 10.27% [0.384] [0.892] [0.

Panel A: Treated firm matched to one control firm. t + 1 t + 2 t + 3 Total CFO Compensation 5.03% 0.84% 10.27% [0.384] [0.892] [0. Online Appendix 1 Table O1: Determinants of CMO Compensation: Selection based on both number of other firms in industry that have CMOs and number of other firms in industry with MBA educated executives

More information

A Hedonic Analysis of Retail Italian Vinegars. Summary. The Model. Vinegar. Methodology. Survey. Results. Concluding remarks.

A Hedonic Analysis of Retail Italian Vinegars. Summary. The Model. Vinegar. Methodology. Survey. Results. Concluding remarks. Vineyard Data Quantification Society "Economists at the service of Wine & Vine" Enometrics XX A Hedonic Analysis of Retail Italian Vinegars Luigi Galletto, Luca Rossetto Research Center for Viticulture

More information

Statistics 5303 Final Exam December 20, 2010 Gary W. Oehlert NAME ID#

Statistics 5303 Final Exam December 20, 2010 Gary W. Oehlert NAME ID# Statistics 5303 Final Exam December 20, 2010 Gary W. Oehlert NAME ID# This exam is open book, open notes; you may use a calculator. Do your own work! Use the back if more space is needed. There are nine

More information

Lollapalooza Did Not Attend (n = 800) Attended (n = 438)

Lollapalooza Did Not Attend (n = 800) Attended (n = 438) D SDS H F 1, 16 ( ) Warm-ups (A) Which bands come to ACL Fest? Is it true that if a band plays at Lollapalooza, then it is more likely to play at Austin City Limits (ACL) that year? To be able to provide

More information

Online Appendix to Voluntary Disclosure and Information Asymmetry: Evidence from the 2005 Securities Offering Reform

Online Appendix to Voluntary Disclosure and Information Asymmetry: Evidence from the 2005 Securities Offering Reform Online Appendix to Voluntary Disclosure and Information Asymmetry: Evidence from the 2005 Securities Offering Reform This document contains several additional results that are untabulated but referenced

More information

PEEL RIVER HEALTH ASSESSMENT

PEEL RIVER HEALTH ASSESSMENT PEEL RIVER HEALTH ASSESSMENT CONTENTS SUMMARY... 2 Overall River Health Scoring... 2 Overall Data Sufficiency Scoring... 2 HYDROLOGY... 3 Overall Hydrology River Health Scoring... 3 Hydrology Data Sufficiency...

More information

STACKING CUPS STEM CATEGORY TOPIC OVERVIEW STEM LESSON FOCUS OBJECTIVES MATERIALS. Math. Linear Equations

STACKING CUPS STEM CATEGORY TOPIC OVERVIEW STEM LESSON FOCUS OBJECTIVES MATERIALS. Math. Linear Equations STACKING CUPS STEM CATEGORY Math TOPIC Linear Equations OVERVIEW Students will work in small groups to stack Solo cups vs. Styrofoam cups to see how many of each it takes for the two stacks to be equal.

More information

Evaluating Population Forecast Accuracy: A Regression Approach Using County Data

Evaluating Population Forecast Accuracy: A Regression Approach Using County Data Evaluating Population Forecast Accuracy: A Regression Approach Using County Data Jeff Tayman, UC San Diego Stanley K. Smith, University of Florida Stefan Rayer, University of Florida Final formatted version

More information

February 26, The results below are generated from an R script.

February 26, The results below are generated from an R script. February 26, 2015 The results below are generated from an R script. weights = read.table("http://www.bio.ic.ac.uk/research/mjcraw/therbook/data/growth.txt", header = T) R functions: aov(y~a+b+a:b, data=mydata)

More information

Using Growing Degree Hours Accumulated Thirty Days after Bloom to Help Growers Predict Difficult Fruit Sizing Years

Using Growing Degree Hours Accumulated Thirty Days after Bloom to Help Growers Predict Difficult Fruit Sizing Years Using Growing Degree Hours Accumulated Thirty Days after Bloom to Help Growers Predict Difficult Fruit Sizing Years G. Lopez 1 and T. DeJong 2 1 Àrea de Tecnologia del Reg, IRTA, Lleida, Spain 2 Department

More information

Statistics & Agric.Economics Deptt., Tocklai Experimental Station, Tea Research Association, Jorhat , Assam. ABSTRACT

Statistics & Agric.Economics Deptt., Tocklai Experimental Station, Tea Research Association, Jorhat , Assam. ABSTRACT Two and a Bud 59(2):152-156, 2012 RESEARCH PAPER Global tea production and export trend with special reference to India Prasanna Kumar Bordoloi Statistics & Agric.Economics Deptt., Tocklai Experimental

More information

Appendix Table A1 Number of years since deregulation

Appendix Table A1 Number of years since deregulation Appendix Table A1 Number of years since deregulation This table presents the results of -in-s models incorporating the number of years since deregulation and using data for s with trade flows are above

More information

Napa Highway 29 Open Wineries

Napa Highway 29 Open Wineries 4 5 6 7 8 9 35 4 45 5 55 Sonoma State University Business 58-Business Intelligence Problem Set #6 Key Dr. Cuellar Trend Analysis-Analyzing Tasting Room Strategies 1. Graphical Analysis a. Show graphically

More information

Appendix A. Table A.1: Logit Estimates for Elasticities

Appendix A. Table A.1: Logit Estimates for Elasticities Estimates from historical sales data Appendix A Table A.1. reports the estimates from the discrete choice model for the historical sales data. Table A.1: Logit Estimates for Elasticities Dependent Variable:

More information

Decision making with incomplete information Some new developments. Rudolf Vetschera University of Vienna. Tamkang University May 15, 2017

Decision making with incomplete information Some new developments. Rudolf Vetschera University of Vienna. Tamkang University May 15, 2017 Decision making with incomplete information Some new developments Rudolf Vetschera University of Vienna Tamkang University May 15, 2017 Agenda Problem description Overview of methods Single parameter approaches

More information

Problem Set #3 Key. Forecasting

Problem Set #3 Key. Forecasting Problem Set #3 Key Sonoma State University Business 581E Dr. Cuellar The data set bus581e_ps3.dta is a Stata data set containing annual sales (cases) and revenue from December 18, 2004 to April 2 2011.

More information

Multiple Imputation for Missing Data in KLoSA

Multiple Imputation for Missing Data in KLoSA Multiple Imputation for Missing Data in KLoSA Juwon Song Korea University and UCLA Contents 1. Missing Data and Missing Data Mechanisms 2. Imputation 3. Missing Data and Multiple Imputation in Baseline

More information

Flexible Working Arrangements, Collaboration, ICT and Innovation

Flexible Working Arrangements, Collaboration, ICT and Innovation Flexible Working Arrangements, Collaboration, ICT and Innovation A Panel Data Analysis Cristian Rotaru and Franklin Soriano Analytical Services Unit Economic Measurement Group (EMG) Workshop, Sydney 28-29

More information

Economics 101 Spring 2016 Answers to Homework #1 Due Tuesday, February 9, 2016

Economics 101 Spring 2016 Answers to Homework #1 Due Tuesday, February 9, 2016 Economics 101 Spring 2016 Answers to Homework #1 Due Tuesday, February 9, 2016 Directions: The homework will be collected in a box before the large lecture. Please place your name, TA name and section

More information

Flexible Imputation of Missing Data

Flexible Imputation of Missing Data Chapman & Hall/CRC Interdisciplinary Statistics Series Flexible Imputation of Missing Data Stef van Buuren TNO Leiden, The Netherlands University of Utrecht The Netherlands crc pness Taylor &l Francis

More information

Cointegration Analysis of Commodity Prices: Much Ado about the Wrong Thing? Mindy L. Mallory and Sergio H. Lence September 17, 2010

Cointegration Analysis of Commodity Prices: Much Ado about the Wrong Thing? Mindy L. Mallory and Sergio H. Lence September 17, 2010 Cointegration Analysis of Commodity Prices: Much Ado about the Wrong Thing? Mindy L. Mallory and Sergio H. Lence September 17, 2010 Cointegration Analysis, Commodity Prices What is cointegration analysis?

More information

Online Appendix. for. Female Leadership and Gender Equity: Evidence from Plant Closure

Online Appendix. for. Female Leadership and Gender Equity: Evidence from Plant Closure Online Appendix for Female Leadership and Gender Equity: Evidence from Plant Closure Geoffrey Tate and Liu Yang In this appendix, we provide additional robustness checks to supplement the evidence in the

More information

After your yearly checkup, the doctor has bad news and good news.

After your yearly checkup, the doctor has bad news and good news. Modeling Belief How much do you believe it will rain? How strong is your belief in democracy? How much do you believe Candidate X? How much do you believe Car x is faster than Car y? How long do you think

More information

Growth in early yyears: statistical and clinical insights

Growth in early yyears: statistical and clinical insights Growth in early yyears: statistical and clinical insights Tim Cole Population, Policy and Practice Programme UCL Great Ormond Street Institute of Child Health London WC1N 1EH UK Child growth Growth is

More information

*p <.05. **p <.01. ***p <.001.

*p <.05. **p <.01. ***p <.001. Table 1 Weighted Descriptive Statistics and Zero-Order Correlations with Fatherhood Timing (N = 1114) Variables Mean SD Min Max Correlation Interaction time 280.70 225.47 0 1095 0.05 Interaction time with

More information

Biologist at Work! Experiment: Width across knuckles of: left hand. cm... right hand. cm. Analysis: Decision: /13 cm. Name

Biologist at Work! Experiment: Width across knuckles of: left hand. cm... right hand. cm. Analysis: Decision: /13 cm. Name wrong 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 right 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 score 100 98.6 97.2 95.8 94.4 93.1 91.7 90.3 88.9 87.5 86.1 84.7 83.3 81.9

More information

Notes on the Philadelphia Fed s Real-Time Data Set for Macroeconomists (RTDSM) Capacity Utilization. Last Updated: December 21, 2016

Notes on the Philadelphia Fed s Real-Time Data Set for Macroeconomists (RTDSM) Capacity Utilization. Last Updated: December 21, 2016 1 Notes on the Philadelphia Fed s Real-Time Data Set for Macroeconomists (RTDSM) Capacity Utilization Last Updated: December 21, 2016 I. General Comments This file provides documentation for the Philadelphia

More information

OF THE VARIOUS DECIDUOUS and

OF THE VARIOUS DECIDUOUS and (9) PLAXICO, JAMES S. 1955. PROBLEMS OF FACTOR-PRODUCT AGGRE- GATION IN COBB-DOUGLAS VALUE PRODUCTIVITY ANALYSIS. JOUR. FARM ECON. 37: 644-675, ILLUS. (10) SCHICKELE, RAINER. 1941. EFFECT OF TENURE SYSTEMS

More information

Measuring the extent of instability in foodgrains production in different districts of Karanataka INTRODUCTION. Research Paper

Measuring the extent of instability in foodgrains production in different districts of Karanataka INTRODUCTION. Research Paper Internationl Research Journal of Agricultural Economics and Statistics Volume 3 Issue 1 March, 2012 53-58 Research Paper Measuring the extent of instability in foodgrains production in different districts

More information

This appendix tabulates results summarized in Section IV of our paper, and also reports the results of additional tests.

This appendix tabulates results summarized in Section IV of our paper, and also reports the results of additional tests. Internet Appendix for Mutual Fund Trading Pressure: Firm-level Stock Price Impact and Timing of SEOs, by Mozaffar Khan, Leonid Kogan and George Serafeim. * This appendix tabulates results summarized in

More information

EFFECT OF TOMATO GENETIC VARIATION ON LYE PEELING EFFICACY TOMATO SOLUTIONS JIM AND ADAM DICK SUMMARY

EFFECT OF TOMATO GENETIC VARIATION ON LYE PEELING EFFICACY TOMATO SOLUTIONS JIM AND ADAM DICK SUMMARY EFFECT OF TOMATO GENETIC VARIATION ON LYE PEELING EFFICACY TOMATO SOLUTIONS JIM AND ADAM DICK 2013 SUMMARY Several breeding lines and hybrids were peeled in an 18% lye solution using an exposure time of

More information

Imputation of multivariate continuous data with non-ignorable missingness

Imputation of multivariate continuous data with non-ignorable missingness Imputation of multivariate continuous data with non-ignorable missingness Thais Paiva Jerry Reiter Department of Statistical Science Duke University NCRN Meeting Spring 2014 May 23, 2014 Thais Paiva, Jerry

More information

Wine Rating Prediction

Wine Rating Prediction CS 229 FALL 2017 1 Wine Rating Prediction Ke Xu (kexu@), Xixi Wang(xixiwang@) Abstract In this project, we want to predict rating points of wines based on the historical reviews from experts. The wine

More information

> Y=degre=="deces" > table(y) Y FALSE TRUE

> Y=degre==deces > table(y) Y FALSE TRUE - PARTIE 0 - > preambule=read.table( + "http://freakonometrics.free.fr/preambule.csv",header=true,sep=";") > table(preambule$y) 0 1 2 3 4 5 6 45 133 160 101 51 8 2 > reg0=glm(y/n~1,family="binomial",weights=n,data=preambule)

More information

Acetic acid dissociates immediately in solution. Reaction A does not react further following the sample taken at the end of

Acetic acid dissociates immediately in solution. Reaction A does not react further following the sample taken at the end of SUPPLEMENTAL Table S1. Assumptions made during fitting of the reaction model in Dynochem. ID Assumption A1 Acetic acid dissociates immediately in solution. A2 Reaction A does not react further following

More information

Modeling Wine Quality Using Classification and Regression. Mario Wijaya MGT 8803 November 28, 2017

Modeling Wine Quality Using Classification and Regression. Mario Wijaya MGT 8803 November 28, 2017 Modeling Wine Quality Using Classification and Mario Wijaya MGT 8803 November 28, 2017 Motivation 1 Quality How to assess it? What makes a good quality wine? Good or Bad Wine? Subjective? Wine taster Who

More information

COMPARISON OF CORE AND PEEL SAMPLING METHODS FOR DRY MATTER MEASUREMENT IN HASS AVOCADO FRUIT

COMPARISON OF CORE AND PEEL SAMPLING METHODS FOR DRY MATTER MEASUREMENT IN HASS AVOCADO FRUIT New Zealand Avocado Growers' Association Annual Research Report 2004. 4:36 46. COMPARISON OF CORE AND PEEL SAMPLING METHODS FOR DRY MATTER MEASUREMENT IN HASS AVOCADO FRUIT J. MANDEMAKER H. A. PAK T. A.

More information

Final Report to Delaware Soybean Board January 11, Delaware Soybean Board

Final Report to Delaware Soybean Board January 11, Delaware Soybean Board Final Report to Delaware Soybean Board January 11, 2017 Delaware Soybean Board (susanne@hammondmedia.com) Effect of Fertigation on Irrigated Full Season and Double Cropped Soybeans Cory Whaley, James Adkins,

More information

Can You Tell the Difference? A Study on the Preference of Bottled Water. [Anonymous Name 1], [Anonymous Name 2]

Can You Tell the Difference? A Study on the Preference of Bottled Water. [Anonymous Name 1], [Anonymous Name 2] Can You Tell the Difference? A Study on the Preference of Bottled Water [Anonymous Name 1], [Anonymous Name 2] Abstract Our study aims to discover if people will rate the taste of bottled water differently

More information

Northern Region Central Region Southern Region No. % of total No. % of total No. % of total Schools Da bomb

Northern Region Central Region Southern Region No. % of total No. % of total No. % of total Schools Da bomb Some Purr Words Laurie and Winifred Bauer A number of questions demanded answers which fell into the general category of purr words: words with favourable senses. Many of the terms supplied were given

More information

J. Best 1 A. Tepley 2

J. Best 1 A. Tepley 2 J. Best 1 2 1 Ecoinformatics Summer Institute, Clarkson University 2 Department of Geography Oregon State University Ecoinformatics Summer Institute Final Presentation Outline 1 Why are fires important?

More information

Internet Appendix to. The Price of Street Friends: Social Networks, Informed Trading, and Shareholder Costs. Jie Cai Ralph A.

Internet Appendix to. The Price of Street Friends: Social Networks, Informed Trading, and Shareholder Costs. Jie Cai Ralph A. Internet Appendix to The Price of Street Friends: Social Networks, Informed Trading, and Shareholder Costs Jie Cai Ralph A. Walkling Ke Yang October 2014 1 A11. Controlling for s Logically Associated with

More information

Bags not: avoiding the undesirable Laurie and Winifred Bauer

Bags not: avoiding the undesirable Laurie and Winifred Bauer Bags not: avoiding the undesirable Laurie and Winifred Bauer Question 10 asked how children claim the right not to do something: 10 Your class is waiting for the bus to arrive to take you on a trip. You

More information

Which of your fingernails comes closest to 1 cm in width? What is the length between your thumb tip and extended index finger tip? If no, why not?

Which of your fingernails comes closest to 1 cm in width? What is the length between your thumb tip and extended index finger tip? If no, why not? wrong 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 right 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 score 100 98.5 97.0 95.5 93.9 92.4 90.9 89.4 87.9 86.4 84.8 83.3 81.8 80.3 78.8 77.3 75.8 74.2

More information

Tips for Writing the RESULTS AND DISCUSSION:

Tips for Writing the RESULTS AND DISCUSSION: Tips for Writing the RESULTS AND DISCUSSION: 1. The contents of the R&D section depends on the sequence of procedures described in the Materials and Methods section of the paper. 2. Data should be presented

More information

MATERIALS AND METHODS

MATERIALS AND METHODS to yields of various sieved fractions and mean particle sizes (MPSs) from a micro hammer-cutter mill equipped with 2-mm and 6-mm screens (grinding time of this mill reported by other investigators was

More information

The R&D-patent relationship: An industry perspective

The R&D-patent relationship: An industry perspective Université Libre de Bruxelles (ULB) Solvay Brussels School of Economics and Management (SBS-EM) European Center for Advanced Research in Economics and Statistics (ECARES) The R&D-patent relationship: An

More information

Fleurieu zone (other)

Fleurieu zone (other) Fleurieu zone (other) Incorporating Southern Fleurieu and Kangaroo Island wine regions, as well as the remainder of the Fleurieu zone outside all GI regions Regional summary report 2006 South Australian

More information

Buying Filberts On a Sample Basis

Buying Filberts On a Sample Basis E 55 m ^7q Buying Filberts On a Sample Basis Special Report 279 September 1969 Cooperative Extension Service c, 789/0 ite IP") 0, i mi 1910 S R e, `g,,ttsoliktill:torvti EARs srin ITQ, E,6

More information

Archdiocese of New York Practice Items

Archdiocese of New York Practice Items Archdiocese of New York Practice Items Mathematics Grade 8 Teacher Sample Packet Unit 1 NY MATH_TE_G8_U1.indd 1 NY MATH_TE_G8_U1.indd 2 1. Which choice is equivalent to 52 5 4? A 1 5 4 B 25 1 C 2 1 D 25

More information

Hybrid ARIMA-ANN Modelling for Forecasting the Price of Robusta Coffee in India

Hybrid ARIMA-ANN Modelling for Forecasting the Price of Robusta Coffee in India International Journal of Current Microbiology and Applied Sciences ISSN: 2319-7706 Volume 6 Number 7 (2017) pp. 1721-1726 Journal homepage: http://www.ijcmas.com Original Research Article https://doi.org/10.20546/ijcmas.2017.607.207

More information

An application of cumulative prospect theory to travel time variability

An application of cumulative prospect theory to travel time variability Katrine Hjorth (DTU) Stefan Flügel, Farideh Ramjerdi (TØI) An application of cumulative prospect theory to travel time variability Sixth workshop on discrete choice models at EPFL August 19-21, 2010 Page

More information

USING STRUCTURAL TIME SERIES MODELS For Development of DEMAND FORECASTING FOR ELECTRICITY With Application to Resource Adequacy Analysis

USING STRUCTURAL TIME SERIES MODELS For Development of DEMAND FORECASTING FOR ELECTRICITY With Application to Resource Adequacy Analysis USING STRUCTURAL TIME SERIES MODELS For Development of DEMAND FORECASTING FOR ELECTRICITY With Application to Resource Adequacy Analysis December 31, 2014 INTRODUCTION In this paper we present the methodology,

More information

COMMUNICATION II Moisture Determination of Cocoa Beans by Microwave Oven

COMMUNICATION II Moisture Determination of Cocoa Beans by Microwave Oven Pertanika 13(1), 67-72 (10) COMMUNICATION II Moisture Determination of Cocoa Beans by Microwave Oven ABSTRAK Analisa statistik telah dijalankan antara kaedah menggunakan ketuhar mikrogelumbang dan kaedah

More information

Internet Appendix for Does Stock Liquidity Enhance or Impede Firm Innovation? *

Internet Appendix for Does Stock Liquidity Enhance or Impede Firm Innovation? * Internet Appendix for Does Stock Liquidity Enhance or Impede Firm Innovation? * This Internet Appendix provides supplemental analyses and robustness tests to the main results presented in Does Stock Liquidity

More information

Not to be published - available as an online Appendix only! 1.1 Discussion of Effects of Control Variables

Not to be published - available as an online Appendix only! 1.1 Discussion of Effects of Control Variables 1 Appendix Not to be published - available as an online Appendix only! 1.1 Discussion of Effects of Control Variables Table 1 in the main text includes a number of additional control variables. We find

More information

RELATIVE EFFICIENCY OF ESTIMATES BASED ON PERCENTAGES OF MISSINGNESS USING THREE IMPUTATION NUMBERS IN MULTIPLE IMPUTATION ANALYSIS ABSTRACT

RELATIVE EFFICIENCY OF ESTIMATES BASED ON PERCENTAGES OF MISSINGNESS USING THREE IMPUTATION NUMBERS IN MULTIPLE IMPUTATION ANALYSIS ABSTRACT RELATIVE EFFICIENCY OF ESTIMATES BASED ON PERCENTAGES OF MISSINGNESS USING THREE IMPUTATION NUMBERS IN MULTIPLE IMPUTATION ANALYSIS Nwakuya, M. T. (Ph.D) Department of Mathematics/Statistics University

More information

Senior poverty in Canada, : A decomposition analysis of income and poverty rates

Senior poverty in Canada, : A decomposition analysis of income and poverty rates Senior poverty in Canada, 1973-2006: A decomposition analysis of income and poverty rates Tammy Schirle Department of Economics Wilfrid Laurier University October 2010 Preliminary and Incomplete - Please

More information

PREDICTION MODEL FOR ESTIMATING PEACH FRUIT WEIGHT AND VOLUME ON THE BASIS OF FRUIT LINEAR MEASUREMENTS DURING GROWTH

PREDICTION MODEL FOR ESTIMATING PEACH FRUIT WEIGHT AND VOLUME ON THE BASIS OF FRUIT LINEAR MEASUREMENTS DURING GROWTH Journal of Fruit and Ornamental Plant Research Vol. 15, 2007: 65-69 PREDICTION MODEL FOR ESTIMATING PEACH FRUIT WEIGHT AND VOLUME ON THE BASIS OF FRUIT LINEAR MEASUREMENTS DURING GROWTH H ü s n ü D e m

More information

Relation between Grape Wine Quality and Related Physicochemical Indexes

Relation between Grape Wine Quality and Related Physicochemical Indexes Research Journal of Applied Sciences, Engineering and Technology 5(4): 557-5577, 013 ISSN: 040-7459; e-issn: 040-7467 Maxwell Scientific Organization, 013 Submitted: October 1, 01 Accepted: December 03,

More information

Oenometrics VII Conference Reims, May 11-13, Predicting Italian wines quality from weather data and experts ratings (DRAFT)

Oenometrics VII Conference Reims, May 11-13, Predicting Italian wines quality from weather data and experts ratings (DRAFT) Oenometrics VII Conference Reims, May 11-13, 2000 Predicting Italian wines quality from weather data and experts by Alessandro Corsi (U. Torino) - Orley Ashenfelter (U. Princeton) (DRAFT) Introduction

More information

Growth dynamics and forecasting of finger millet (Ragi) production in Karnataka

Growth dynamics and forecasting of finger millet (Ragi) production in Karnataka Growth dynamics and forecasting of finger millet (Ragi) production in Karnataka Veerabhadrappa Bellundagi*, K.B. Umesh and S.C. Ravi Department of Agricultural Economics, UAS, GKVK, Bengaluru, Karnataka,

More information