PSYC 6140 November 16, 2005 ANOVA output in R
|
|
- Florence Collins
- 6 years ago
- Views:
Transcription
1 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 discusses how to use and interpret anova output for regression models with interaction. The following script shows output for a model with interaction using all three types of sum of squares. In summary: Note: 1. Type I and Type II SS compare models that obey the principle of marginality and, thus, these tables remain the same if the categorical variables is recoded or if the continuous variable is subject to an affine transformation. 2. Type III SS depend on the coding the library(car) source(" coursefun.r Last update: Oct. 27, 2005 An R script containing functions and some datasets for the courses PSYC6140 and MATH6630 in Help on the following functions is available by typing the name of the function. This text is available by typing coursefun A current version of this file can be sourced or downloaded from A copy is kept at: Please make corrections, changes and additions to the version on the wiki. They will be periodically transferred to the downloadable version. Functions: Tables: atotal: border an array with sums
2 ANOVA output in R 2 abind : glue two arrays together on a selected dimension Graphics: td : easy front end to trellis.device and related functions xqplot: extended quantile plots 3D graphics by John Fox: scatter3d identify3d ellipsoid Inference cell - a modified version of car::confidence.ellipse.lm that can add to a plot Graphics for linear algebra vplot - plots column vectors adding to current plot vell - ellipse as a 2 x n matrix vbox - box around unit circle orthog - 2 x 2 rotation orthog.prog 2 x 2 matrix of orthog projection To add functions, modify data( Prestige ) scatterplot.matrix(prestige) ## Let's just use those with non.missing type dd <- na.omit(prestige) dim(dd) [1] 98 6 dim(prestige) [1] attach(dd)
3 ANOVA output in R 3 The following object(s) are masked from package:datasets : women table(type) type bc prof wc income.log <- log(income ) scatterplot( income.log, prestige, groups = type)
4 ANOVA output in R 4 prestige bc prof wc income.log
5 ANOVA output in R 5 Model: E( Y ) = R Syntax P. of Marg. Comments 1 α + β X + γ 1D1+ γ 2D2 + δ1dx 1 + δ2dx Y ~ X*G 2 Yes Full model 2N α + γ 1D1+ γ 2D2 + δ1dx 1 + δ2dx Y ~ G + X:G 2 No Y ~ X*G X 3N α + β X + δ1dx 1 + δ2dx Y ~ X + X:G 2 No Y ~ X*G G 4N α + δ1dx 1 + δ2dx Y ~ X:G 2 No 5 α β X + γ 1D1 γ 2D2 + + Y ~ X + G Yes Additive model 6 α + γ 1D1+ γ 2D2 Y ~ G Yes 7 α + β X Y ~ X Yes 8 α Y ~ 1 Yes Intercept only model fit.add <- lm( prestige ~ income.log + type ) Additive model fit.int <- lm( prestige ~ income.log * type ) Full interaction model summary(fit.int) Call: lm(formula = prestige ~ income.log * type) Residuals: Min 1Q Median 3Q Max Coefficients: Estimate Std. Error t value Pr( t ) (Intercept) e-09 *** income.log e-10 *** typeprof e-05 *** typewc * income.log:typeprof *** income.log:typewc * Exercise: Draw a sketch of the fitted model and label it with each of the estimated values in this table
6 ANOVA output in R 6 Residual standard error: on 92 degrees of freedom Multiple R-Squared: , Adjusted R-squared: F-statistic: on 5 and 92 DF, p-value: < 2.2e-16 anova(fit.int) ANOVA for full model using anova Analysis of Variance Table Type I SS (sequential) Df Sum Sq Mean Sq F value Pr(F) income.log < 2.2e-16 *** (8) (7) / (1) M:ok type < 2.2e-16 *** (7) (5) / (1) M:ok income.log:type *** Residuals (5) (1) / (1) M:ok EXERCISE: What would you get if you specified the model as prestige ~ income.log * type anova( fit.add, fit.int) Analysis of Variance Table Model 1: prestige ~ income.log + type Model 2: prestige ~ income.log * type Res.Df RSS Df Sum of Sq F Pr(F) *** (5) (1) / (1) ## ## Interpreting coefficients ## contrasts(type) # bc is reference level prof wc bc 0 0 bc is reference level prof 1 0 wc 0 1
7 ANOVA output in R 7 fit.int$contrasts $type [1] "contr.treatment" anova( fit.int) Analysis of Variance Table Df Sum Sq Mean Sq F value Pr(F) income.log < 2.2e-16 *** type < 2.2e-16 *** income.log:type *** Residuals Anova( fit.int, type = "II") ANOVA using Anova Type II Anova Table (Type II tests) in library(car) NOTE: All satisfy PoM Sum Sq Df F value Pr(F) income.log e-11 *** (6) (5) / (1) [Different] type < 2.2e-16 *** (7) (5) / (1) [Same as Type I] income.log:type *** (5) (1) / (1) [Same as Type I] Residuals Anova( fit.int, type = "III") ANOVA using Anova Type III Anova Table (Type III tests) in library(car) Sum Sq Df F value Pr(F) (Intercept) e-09 *** income.log e-10 *** (2N) (1) / (1) [Depends on coding of G] type *** (3N) (1) / (1) [Depends on 0 of X] income.log:type *** (5) (1) / (1) [Same as I and II] Residuals
8 ANOVA output in R 8 fit.int.s <- lm( prestige ~ income.log * type, contrasts = list( type = contr.sum ) ) Sum to 0 coding fit.int.s$contrasts $type [,1] [,2] bc 1 0 Note: No type is reference level. prof 0 1 Reference level is at the mean wc -1-1 of the 3 types. summary(fit.int.s) Call: lm(formula = prestige ~ income.log * type, contrasts = list(type = contr.sum)) Residuals: Min 1Q Median 3Q Max Coefficients: Estimate Std. Error t value Pr( t ) (Intercept) e-08 *** income.log e-10 *** type *** type ** income.log:type ** income.log:type * EXERCISE: Sketch the fitted model and label with values of fitted coefficients Residual standard error: on 92 degrees of freedom Multiple R-Squared: , Adjusted R-squared: F-statistic: on 5 and 92 DF, p-value: < 2.2e-16
9 ANOVA output in R 9 anova( fit.int.s ) Analysis of Variance Table Df Sum Sq Mean Sq F value Pr(F) income.log < 2.2e-16 *** type < 2.2e-16 *** income.log:type *** Residuals anova( fit.add, fit.int.s ) Analysis of Variance Table EXERCISE: Fill in models Model 1: prestige ~ income.log + type Model 2: prestige ~ income.log * type Res.Df RSS Df Sum of Sq F Pr(F) EXERCISE: Fill in models *** EXERCISE: What happens if we use fit.add.s instead? anova( fit.int.s) Analysis of Variance Table Df Sum Sq Mean Sq F value Pr(F) income.log < 2.2e-16 *** type < 2.2e-16 *** income.log:type *** Residuals
10 ANOVA output in R 10 Anova Type II SS Anova( fit.int.s, type = "II") EXERCISE: Fill in models? Anova Table (Type II tests) Do you get same output as before? Sum Sq Df F value Pr(F) income.log e-11 *** type < 2.2e-16 *** income.log:type *** Residuals Anova Type III SS Anova( fit.int.s, type = "III") EXERCISE: Draw a sketch showing Anova Table (Type III tests) what is tested in this table. Sum Sq Df F value Pr(F) (Intercept) e-08 *** income.log e-10 *** type *** income.log:type *** Residuals Note: Type III SS using sum to 0 are quite popular in Psychology and Sociology. They are readily produced in SAS output. Two cautions: 1. Is the hypothesis meaningful in the presence of interaction? i.e. is the hypothesis that the average slope, averaging over levels of the factor, equals 0 of real interest? 2. If the number of observations in each category of the factor is quite variable, then the Type III SS will have low power in comparison with a hypothesis that uses a weighted average of slopes -- as the Type II SS.
11 ANOVA output in R 11 ## ## 3 way interaction ## summary( lm( prestige ~ income.log * type * women )) 3 way interaction EXERCISE: Draw 2 sketches, one for Call: women = 0 and another for lm(formula = prestige ~ income.log * type * women) women = 100. Label sketches appropriately Residuals: Min 1Q Median 3Q Max Coefficients: Estimate Std. Error t value Pr( t ) (Intercept) e-09 *** income.log e-10 *** typeprof e-05 *** typewc women income.log:typeprof e-05 *** income.log:typewc income.log:women typeprof:women typewc:women income.log:typeprof:women income.log:typewc:women Residual standard error: on 86 degrees of freedom Multiple R-Squared: , Adjusted R-squared: F-statistic: on 11 and 86 DF, p-value: < 2.2e-16 fit3.s <- lm( prestige ~ income.log * type * women, contrasts = list( type = contr.sum)) Sum to 0 coding
12 ANOVA output in R 12 ANOVA tables with treatment coding and with sum to 0 coding anova(fit3) Analysis of Variance Table Df Sum Sq Mean Sq F value Pr(F) income.log < 2.2e-16 *** type < 2.2e-16 *** women * income.log:type e-05 *** income.log:women type:women ** income.log:type:women Residuals Anova(fit3, type = "II") Anova Table (Type II tests) Sum Sq Df F value Pr(F) income.log e-11 *** type e-12 *** women *** income.log:type e-06 *** income.log:women type:women ** income.log:type:women Residuals Anova(fit3, type = "III") Anova Table (Type III tests) Sum Sq Df F value Pr(F) (Intercept) e-09 *** income.log e-10 *** type e-05 *** women
13 ANOVA output in R 13 income.log:type e-05 *** income.log:women type:women income.log:type:women Residuals anova(fit3.s) Analysis of Variance Table Df Sum Sq Mean Sq F value Pr(F) income.log < 2.2e-16 *** type < 2.2e-16 *** women * income.log:type e-05 *** income.log:women type:women ** income.log:type:women Residuals Anova(fit3.s, type = "II") Anova Table (Type II tests) Sum Sq Df F value Pr(F) income.log e-11 *** type e-12 *** women *** income.log:type e-06 *** income.log:women type:women ** income.log:type:women Residuals
14 ANOVA output in R 14 Anova(fit3.s, type = "III") Anova Table (Type III tests) Sum Sq Df F value Pr(F) (Intercept) e-07 *** income.log e-08 *** type e-05 *** women income.log:type e-05 *** income.log:women type:women income.log:type:women Residuals ## Refit without 3-way interaction: fit3.2 <- lm( prestige ~ (income.log + type + women)^2 ) summary(fit3.2) Call: lm(formula = prestige ~ (income.log + type + women)^2) Residuals: Min 1Q Median 3Q Max Coefficients: Estimate Std. Error t value Pr( t ) (Intercept) e-10 *** income.log e-11 *** typeprof e-06 *** typewc women income.log:typeprof e-06 *** income.log:typewc income.log:women
15 ANOVA output in R 15 typeprof:women * typewc:women Residual standard error: on 88 degrees of freedom Multiple R-Squared: , Adjusted R-squared: F-statistic: on 9 and 88 DF, p-value: < 2.2e-16 Anova( fit3.2, type = "II") Why would you prefer Type II here Anova Table (Type II tests) Sum Sq Df F value Pr(F) income.log e-12 *** type e-12 *** women *** income.log:type e-06 *** income.log:women type:women ** Residuals Drop income.log:women interaction fit3.22 <- lm( prestige ~ (income.log + women) * type ) summary(fit3.22) Call: lm(formula = prestige ~ (income.log + women) * type) Residuals: Min 1Q Median 3Q Max Coefficients: Estimate Std. Error t value Pr( t ) (Intercept) e-11 ***
16 ANOVA output in R 16 income.log e-12 *** women ** typeprof e-06 *** typewc income.log:typeprof e-06 *** income.log:typewc women:typeprof * women:typewc Residual standard error: on 89 degrees of freedom Multiple R-Squared: , Adjusted R-squared: F-statistic: 71.2 on 8 and 89 DF, p-value: < 2.2e-16 Anova( fit3.22, type = "II") Anova Table (Type II tests) Sum Sq Df F value Pr(F) income.log e-12 *** women *** type e-15 *** income.log:type e-07 *** women:type ** Residuals
17 ANOVA output in R 17 Graphical presentation of models with higher-order interactions summary(income.log) Min. 1st Qu. Median Mean 3rd Qu. Max pred <- expand.grid( income.log = seq(9,10.5,.1), + type = levels(type), + women = c(0,100)) pred$prestige <- predict( fit3.22, newdata = pred) library(lattice) td(new=t) xyplot( prestige ~ income.log women, pred, groups = type, type = 'l', + auto.key = T)
18 ANOVA output in R 18 bc prof wc women women 100 prestige income.log
19 ANOVA output in R 19 td( col = c('red','blue','black'), lwd = 1.5) xyplot( prestige ~ income.log + factor(paste("women =",women,"%")), pred, + groups = type, type = 'l', + auto.key = list(columns= 3, lines=t, points = F))
20 ANOVA output in R 20 bc prof wc women = 0 % women = 100 % prestige income.log
21 ANOVA output in R 21 xyplot( prestige ~ income.log type, + pred, + groups = factor(paste("women =",women,"%")), type = 'l', + auto.key = list(columns= 2, lines=t, points = F))
22 ANOVA output in R 22 women = 0 % women = 100 % wc prestige 150 bc prof \ income.log
STAT 5302 Applied Regression Analysis. Hawkins
Homework 3 sample solution 1. MinnLand data STAT 5302 Applied Regression Analysis. Hawkins newdata
More informationwine 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 informationThe 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 informationTo: 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 informationComparing 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 informationFrom 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 informationAJAE 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 informationGail 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 informationMissing 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 informationWine-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 informationINSTITUTE 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 informationMultiple 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 informationBusiness 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 informationActivity 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 informationFebruary 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 informationIT 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 informationCurtis 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 informationSummary 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 informationAppendix 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 informationFlexible 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 informationCredit 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 informationWine 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 informationOnline 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 informationMissing 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 informationDecision 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 informationThe 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 informationBORDEAUX 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 informationNovember K. J. Martijn Cremers Lubomir P. Litov Simone M. Sepe
ONLINE APPENDIX TABLES OF STAGGERED BOARDS AND LONG-TERM FIRM VALUE, REVISITED November 016 K. J. Martijn Cremers Lubomir P. Litov Simone M. Sepe The paper itself is available at https://papers.ssrn.com/sol3/papers.cfm?abstract-id=364165.
More informationStatistics 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 informationOnline Appendix to The Effect of Liquidity on Governance
Online Appendix to The Effect of Liquidity on Governance Table OA1: Conditional correlations of liquidity for the subsample of firms targeted by hedge funds This table reports Pearson and Spearman correlations
More informationMissing 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> library(sem) > cor.mat<-read.moments(names=c("ten1", "ten2", "ten3", "wor1", "wor2", + "wor3", "irthk1", "irthk2", "irthk3", "body1", "body2",
> library(sem) > cor.mat
More informationGasoline 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 informationRelationships 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 informationFinal 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 informationDietary Diversity in Urban and Rural China: An Endogenous Variety Approach
Dietary Diversity in Urban and Rural China: An Endogenous Variety Approach Jing Liu September 6, 2011 Road Map What is endogenous variety? Why is it? A structural framework illustrating this idea An application
More informationThe 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 informationProblem 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 informationClimate 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 informationPredicting 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 informationHandling 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 informationNotes 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 informationNapa 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 informationThe 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 informationPanel 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 informationMini Project 3: Fermentation, Due Monday, October 29. For this Mini Project, please make sure you hand in the following, and only the following:
Mini Project 3: Fermentation, Due Monday, October 29 For this Mini Project, please make sure you hand in the following, and only the following: A cover page, as described under the Homework Assignment
More informationMethod for the imputation of the earnings variable in the Belgian LFS
Method for the imputation of the earnings variable in the Belgian LFS Workshop on LFS methodology, Madrid 2012, May 10-11 Astrid Depickere, Anja Termote, Pieter Vermeulen Outline 1. Introduction 2. Imputation
More informationAWRI Refrigeration Demand Calculator
AWRI Refrigeration Demand Calculator Resources and expertise are readily available to wine producers to manage efficient refrigeration supply and plant capacity. However, efficient management of winery
More informationFair Trade and Free Entry: Can a Disequilibrium Market Serve as a Development Tool? Online Appendix September 2014
Fair Trade and Free Entry: Can a Disequilibrium Market Serve as a Development Tool? 1. Data Construction Online Appendix September 2014 The data consist of the Association s records on all coffee acquisitions
More informationA Comparison of Approximate Bayesian Bootstrap and Weighted Sequential Hot Deck for Multiple Imputation
A Comparison of Approximate Bayesian Bootstrap and Weighted Sequential Hot Deck for Multiple Imputation Darryl V. Creel RTI International 1 RTI International is a trade name of Research Triangle Institute.
More informationOnline 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 informationRelation 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 informationThis 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 informationFlexible 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 informationA 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 informationAppendix 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 informationThe premium for organic wines
Enometrics XV Collioure May 29-31, 2008 Estimating a hedonic price equation from the producer side Points of interest: - assessing whether there is a premium for organic wines, and which one - estimating
More informationNorthern 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 informationModeling 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 informationLollapalooza 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 information1. Identify environmental conditions (temperature) and nutritional factors (i.e. sugar and fat) that encourages the growth of bacteria.
Food Explorations Lab: Magnificent Microbes STUDENT LAB INVESTIGATIONS Name: Lab Overview In this investigation, you will use bacterial fermentation to produce yogurt. Fat content, sugar content (lactose),
More informationValuation 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 informationZeitschrift 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 informationReturn to wine: A comparison of the hedonic, repeat sales, and hybrid approaches
Return to wine: A comparison of the hedonic, repeat sales, and hybrid approaches James J. Fogarty a* and Callum Jones b a School of Agricultural and Resource Economics, The University of Western Australia,
More informationWhich 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 informationPoisson GLM, Cox PH, & degrees of freedom
Poisson GLM, Cox PH, & degrees of freedom Michael C. Donohue Alzheimer s Therapeutic Research Institute Keck School of Medicine University of Southern California December 13, 2017 1 Introduction We discuss
More informationAnalysis 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 informationBiologist 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 informationPreferred citation style
Preferred citation style Axhausen, K.W. (2016) How many cars are too many? A second attempt, distinguished transport lecture at the University of Hong Kong, Hong Kong, October 2016.. How many cars are
More informationCointegration 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 information1. Identify environmental conditions (temperature) and nutritional factors (i.e. sugar and fat) that encourages the growth of bacteria.
Food Explorations Lab II: Magnificent Microbes STUDENT LAB INVESTIGATIONS Name: Lab Overview In this investigation, you will use bacterial fermentation to produce yogurt. Fat content, sugar content (lactose),
More informationImputation 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 informationMissing Data: Part 2 Implementing Multiple Imputation in STATA and SPSS. Carol B. Thompson Johns Hopkins Biostatistics Center SON Brown Bag 4/24/13
Missing Data: Part 2 Implementing Multiple Imputation in STATA and SPSS Carol B. Thompson Johns Hopkins Biostatistics Center SON Brown Bag 4/24/13 Overview Reminder Steps in Multiple Imputation Implementation
More informationAppendix A. Table A1: Marginal effects and elasticities on the export probability
Appendix A Table A1: Marginal effects and elasticities on the export probability Variable PROP [1] PROP [2] PROP [3] PROP [4] Export Probability 0.207 0.148 0.206 0.141 Marg. Eff. Elasticity Marg. Eff.
More informationOnline 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 informationUpdate to A Comprehensive Look at the Empirical Performance of Equity Premium Prediction
Update to A Comprehensive Look at the Empirical Performance of Equity Premium Prediction Amit Goyal UNIL Ivo Welch UCLA September 17, 2014 Abstract This file contains updates, one correction, and links
More informationTHE 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 informationRheological and physicochemical studies on emulsions formulated with chitosan previously dispersed in aqueous solutions of lactic acid
SUPPLEMENTARY MATERIAL (SM) FOR Rheological and physicochemical studies on emulsions formulated with chitosan previously dispersed in aqueous solutions of lactic acid Lucas de Souza Soares a, Janaína Teles
More informationThe age of reproduction The effect of university tuition fees on enrolment in Quebec and Ontario,
The age of reproduction The effect of university tuition fees on enrolment in Quebec and Ontario, 1946 2011 Benoît Laplante, Centre UCS de l INRS Pierre Doray, CIRST-UQAM Nicolas Bastien, CIRST-UQAM Research
More informationManaging many models. February Hadley Chief Scientist, RStudio
Managing many models February 2016 Hadley Wickham @hadleywickham Chief Scientist, RStudio There are 7 key components of data science Import Visualise Communicate Tidy Transform Model Automate Understand
More informationPEEL 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 informationPitfalls for the Construction of a Welfare Indicator: An Experimental Analysis of the Better Life Index
Clemens Hetschko, Louisa von Reumont & Ronnie Schöb Pitfalls for the Construction of a Welfare Indicator: An Experimental Analysis of the Better Life Index University Alliance of Sustainability Spring
More information> 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 informationLong term impacts of facilitating temporary contracts: A comparative analysis of Italy and Spain using birth cohorts
Long term impacts of facilitating temporary contracts: A comparative analysis of Italy and Spain using birth cohorts Miguel Á. Malo (University of Salamanca, Spain) Dario Sciulli (University of Chietti
More informationProduction scale-up and adjustments in the specifications of a new type of biscuit
Production scale-up and adjustments in the specifications of a new type of biscuit Inês Daniela Graça ABSTRACT New products development is a constant in the food industry. The biscuit S is a new kind of
More informationNot 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 informationInvestment 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 informationTable 1: Number of patients by ICU hospital level and geographical locality.
Web-based supporting materials for Evaluating the performance of Australian and New Zealand intensive care units in 2009 and 2010, by J. Kasza, J. L. Moran and P. J. Solomon Table 1: Number of patients
More informationAssignment # 1: Answer key
INTERNATIONAL TRADE Autumn 2004 NAME: Section: Assignment # 1: Answer key Problem 1 The following Ricardian table shows the number of days of labor input needed to make one unit of output of each of the
More informationtutorial_archetypes_prototypes_siqd_ensembles.r michael Sat Oct 29 21:38:
Electronic Supplementary Material (ESI) for Nanoscale. This journal is The Royal Society of Chemistry 2016 tutorial_archetypes_prototypes_siqd_ensembles.r michael Sat Oct 29 21:38:31 2016 #R script to
More informationSTACKING 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 informationWhat makes a good muffin? Ivan Ivanov. CS229 Final Project
What makes a good muffin? Ivan Ivanov CS229 Final Project Introduction Today most cooking projects start off by consulting the Internet for recipes. A quick search for chocolate chip muffins returns a
More informationOptimization of Saccharomyces cerevisiae immobilization in bacterial cellulose by adsorption- incubation method
(009) Optimization of Saccharomyces cerevisiae immobilization in bacterial cellulose by adsorption- incubation method Nguyen, D. N., Ton, N. M. N. and * Le, V. V. M. Department of Food Technology, Ho Chi
More informationBags 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 informationEestimated coefficient. t-value
Table 1: Estimated wage curves for men, 1983 2009 Dependent variable: log (real wage rate) Dependent variable: log real wage rate Men 1983-2009 Men, 1983-2009 Rendom-effect Fixed-effect z-vae t-vae Men
More informationEffect of Inocucor on strawberry plants growth and production
Effect of Inocucor on strawberry plants growth and production Final report For Inocucor Technologies Inc. 20 Grove, Knowlton, Quebec, J0E 1V0 Jae Min Park, Dr. Soledad Saldías, Kristen Delaney and Dr.
More informationEvaluation of Alternative Imputation Methods for 2017 Economic Census Products 1 Jeremy Knutson and Jared Martin
Evaluation of Alternative Imputation Methods for 2017 Economic Census Products 1 Jeremy Knutson and Jared Martin Abstract In preparation for the 2017 change to the North American Product Classification
More informationStatewide Monthly Participation for Asian Non-Hispanic by Cultural Identities, Months of Prenatal Participation & Breastfeeding Amount
MN WIC PROGRAM Statewide Monthly Participation for Asian Non-Hispanic by Cultural Identities, Months of Prenatal Participation & Breastfeeding Amount COUNTS/MONTHLY PARTICIPATION (9.15.16) Report Overview
More informationIntroduction to Management Science Midterm Exam October 29, 2002
Answer 25 of the following 30 questions. Introduction to Management Science 61.252 Midterm Exam October 29, 2002 Graphical Solutions of Linear Programming Models 1. Which of the following is not a necessary
More information