Latin Square Design. Design of Experiments - Montgomery Section 4-2
|
|
- Barnaby Wheeler
- 6 years ago
- Views:
Transcription
1 Latin Square Design Design of Experiments - Montgomery Section 4-2
2 Latin Square Design Can be used when goal is to block on two nuisance factors Constructed so blocking factors orthogonal to treatment Must have same the number of blocks and treatments y ijk = µ + α i + τ j + β k + ǫ ijk α i - ith Block 1 effect (row) τ j - jth treatment effect β k - kth Block 2 effect (column) ǫ ijk N(0, σ 2 ) i = 1, 2,..., p j = 1, 2,..., p k = 1, 2,..., p Completely additive model (no interaction) Latin Square 2
3 Latin Square Design Design commonly represented in p p grid There are now two randomization restrictions One trt per row (row = blocking factor 1) One trt per column (column = blocking factor 2) Can randomly shuffle rows, columns, and treatments of standard square to get other variations of layout Standard square has treatment levels written alphabetically in the first row and first column. Other cells are filled using the same alphabetic order Latin Square 3
4 Layout Examples C B A B A C A C B D B C A B A C C D A B A C B B A D C C B A A C B D Latin Square 4
5 Proc PLAN title Latin Square Design ; proc plan seed=12; factors rows=4 ordered cols=4 ordered / NOPRINT; treatments tmts=4 cyclic; output out=g rows cvals=( Day 1 Day 2 Day 3 Day 4 ) random cols cvals=( Lab 1 Lab 2 Lab 3 Lab 4 ) random tmts nvals=( ) random; proc tabulate; class rows cols; var tmts; table rows, cols*(tmts*f=6.) / rts=8; run; Latin Square 5
6 Proc PLAN Output cols Lab 1 Lab 2 Lab 3 Lab 4 tmts tmts tmts tmts Sum Sum Sum Sum rows Day Day Day Day Latin Square 6
7 Partitioning the SS Rewrite observation as: y ijk = y... + (y i.. y... ) + (y.j. y... ) + (y..k y... ) + (y ijk y i.. y.j. y..k + 2y.. ) = ˆµ + ˆα i + ˆτ j + ˆβk + ˆǫ ijk p Partition SS T into (y i.. y... ) 2 + p (y.j. y... ) 2 + p (y..k y... ) 2 + ˆǫ 2 ijk SS Row + SS Treatment + SS Col + SS E Under H 0, all SS/σ 2 independent chi-squared RVs Usual F-test analysis Caution testing column and row effects Latin Square 7
8 Analysis of Variance Table Source of Sum of Degrees of Mean F 0 Variation Squares Freedom Square Rows SS Row p 1 MS Row Treatment SS Treatment p 1 MS Treatment F 0 Column SS Column p 1 MS Column Error SS E (p 2)(p 1) MS E Total SS T p 2 1 If F 0 > F α,p 1,(p 2)(p 1) then reject H 0 Latin Square 8
9 Comparisons of Treatments Multiple Comparisons/Contrasts Similar procedures as before with CRD Replace n with p in all standard error formulas Degrees of freedom error is (p 2)(p 1) y.j. = µ + ᾱ. + τ j + β. + ǫ.j. Var(y.j. ) = σ 2 /p y.j. y.j. = τ j τ j + ǫ.j. ǫ.j. Var(y.j. y.j.) = σ 2 /p + σ 2 /p Latin Square 9
10 Missing Values When missing (and not due to treatment) Design unbalanced Orthogonality lost Order of fit important Procedures 1 Regression approach Use Type III sum of squares 2 Estimate missing value Choose value to minimize SS E Take derivative and set equal to zero y ijk = p(y i.. +y.j. +y..k ) 2y... (p 2)(p 1) Latin Square 10
11 Using SAS Consider experiment to investigate the effect of 4 diets on milk production. There are 4 cows. Each lactation period the cows receive a different diet. Assume there is a washout period so previous diet does not affect future results. Will therefore block on lactation period and cow. data new; input cow period trt cards; ; proc glm; class cow trt period; model resp=trt period cow; means trt/ lines tukey; lsmeans trt / adjust=tukey lines; output out=new1 r=res p=pred; symbol1 v=circle; proc gplot; plot res*pred; plot res*cow; plot res*trt; plot res*period; proc univariate noprint; histogram res / normal (L=1 mu=0 sigma=est) kernel (L=2); run; Latin Square 11
12 ANOVA Table Dependent Variable: resp Sum of Source DF Squares Mean Square F Value Pr > F Model Error Corrected Total Source DF Type I SS Mean Square F Value Pr > F trt period <.0001 cow Source DF Type III SS Mean Square F Value Pr > F trt period <.0001 cow Latin Square 12
13 Multiple Comparisons Tukey s Studentized Range (HSD) Test for resp Alpha 0.05 Error Degrees of Freedom 6 Error Mean Square Critical Value of Studentized Range Minimum Significant Difference Mean N trt A A B B Latin Square 13
14 Latin Square 14
15 Latin Square 15
16 Latin Square 16
17 Using Proc Mixed Often at least one block is considered random Example: 4 cows randomly chosen from large herd Want the inference to extend to the herd Sometimes time is one of the blocking factors Example: Each cow measured over 4 lactation periods This is an example of a crossover design Could observations closer together be more correlated? Proc Mixed can incorporate both concepts into model Latin Square 17
18 Random Effects Similar approach and results as with the RCBD Must be careful if using Proc GLM Standard error for a mean: Proc GLM incorrect Standard error for a contrast: Proc GLM correct GLM very limited in terms of handling repeated measures Latin Square 18
19 Comparison of Approaches proc glm; class cow trt period; model resp=cow trt period; random cow; lsmeans trt / stderr tdiff lines; proc mixed; class cow trt period; model resp=trt period / ddfm=kr; ***recommended df adjustment option random cow; lsmeans trt/ diff; run; ****F test for treatment is the same for both models ***** Latin Square 19
20 GLM Selected Output Standard LSMEAN trt resp LSMEAN Error Pr > t Number < < < < resp LSMEAN LSMEAN trt Number A A A B B B Latin Square 20
21 MIXED Selected Output Least Squares Means Effect trt Estimate Std Error DF t Value Pr > t trt <.0001 trt <.0001 trt <.0001 trt <.0001 Differences of Least Squares Means Effect trt _trt Estimate Std Error DF t Value Pr > t trt trt trt trt trt trt Latin Square 21
22 Correlated Errors Random effects induce correlations among observations Errors from same EU may also be correlated Errors closer in time may be more similar Can incorporate various correlation structures Given p obs per EU, consider a p p covariance matrix Main diagonal contains the variances Off-diagonal elements represent covariances Latin Square 22
23 Covariance Matrices Uncorrelated errors (for p = 4 obs/cow) σ σ σ 2 0 = σ σ st order autoregressive (for p = 4 obs/cow) σ 2 1 ρ ρ 2 ρ 3 ρ 1 ρ ρ 2 ρ 2 ρ 1 ρ ρ 3 ρ 2 ρ 1 Latin Square 23
24 Covariance Matrices Compound symmetry (for p = 4 obs/cow) σ 2 + σ 1 σ 1 σ 1 σ 1 σ 1 σ 2 + σ 1 σ 1 σ 1 σ 1 σ 1 σ 2 + σ 1 σ 1 σ 1 σ 1 σ 1 σ 2 + σ 1 Toeplitz (for p = 4 obs/cow) σ 2 σ 1 σ 2 σ 3 σ 1 σ 2 σ 1 σ 2 σ 2 σ 1 σ 2 σ 1 σ 3 σ 2 σ 1 σ 2 Latin Square 24
25 Covariance Matrices Banded Toeplitz d = 2(for p = 4 obs/cow) σ 2 σ σ 1 σ 2 σ σ 1 σ 2 σ σ 1 σ 2 Unstructured (for p = 4 obs/cow) σ1 2 σ 12 σ 13 σ 14 σ 12 σ2 2 σ 23 σ 24 σ 13 σ 23 σ3 2 σ 34 σ 14 σ 24 σ 34 σ4 2 Latin Square 25
26 Mixed Model Can be expressed Y = Xβ + Zδ + e X and Z describe the treatment and design structure β is a vector of the fixed effect parameters δ is a vector of the random effect parameters Typically assume δ and e are multivariate Normal (MVN) Mean 0 Covariance matrices G and R uncorrelated Latin Square 26
27 Conditional/Marginal Models Conditional model of Y δ is MVN E(Y δ) = Xβ + Zδ Cov(Y δ) = R Marginal model of Y is MVN E(Y ) = Xβ Cov(Y ) = ZGZ + R Can specify either model in MIXED. Random statement needed for conditional model. When dists Normal, both models give identical estimates, unless covariances negative Latin Square 27
28 SAS Code /* Fit a simple correlation structure for residuals */ /* Assuming cow is fixed - matches glm output on Slide 20 */ proc mixed covtest cl; class cow trt period; model resp=trt period cow / ddfm=kr outp=diag; repeated period / subject=cow type=simple r=1; lsmeans trt / diff; run; /* Assuming cow is random - matches mixed output on Slide 21 */ proc mixed covtest cl; class cow trt period; model resp=trt period / ddfm=kr outp=diag; random intercept / subject=cow v=1; repeated period / subject=cow type=simple r=1; lsmeans trt / diff; run; Latin Square 28
29 Simple - Fixed effects Estimated R Matrix for cow 1 Row Col1 Col2 Col3 Col Covariance Parameter Estimates Cov Standard Z Parm Subject Estimate Error Value Pr > Z Alpha Lower Upper period cow Fit Statistics -2 Res Log Likelihood 26.9 AIC (Smaller is Better) 28.9 BIC (Smaller is Better) 28.3 Type 3 Tests of Fixed Effects Num Den Effect DF DF F Value Pr > F trt Effect trt Estimate Std Error DF t Value Pr > t trt <.0001 Effect trt _trt Estimate Std Error DF t Value Pr > t trt Latin Square 29
30 Simple - Random effects Estimated V Matrix for cow 1 Row Col1 Col2 Col3 Col Covariance Parameter Estimates Standard Z Cov Parm Subject Estimate Error Value Pr > Z Alpha Lower Upper Intercept cow period cow Fit Statistics -2 Res Log Likelihood 41.3 AIC (Smaller is Better) 45.3 BIC (Smaller is Better) 44.1 Type 3 Tests of Fixed Effects Num Den Effect DF DF F Value Pr > F trt Effect trt Estimate Std Error DF t Value Pr > t trt <.0001 Effect trt _trt Estimate Std Error DF t Value Pr > t trt Latin Square 30
31 SAS Code /* Fit Compound Symmetry - matches mixed output on Slide 21 */ proc mixed covtest cl; class cow trt period; model resp=trt period / ddfm=kr outp=diag; repeated period / subject=cow type=cs r=1; lsmeans trt / diff; run; /* Fit an Unstructured correlation structure */ proc mixed covtest cl; class cow trt period; model resp=trt period / ddfm=kr outp=diag; repeated period / type=un subject=cow; lsmeans trt / diff; run; Latin Square 31
32 Compound Symmetry Estimated R Matrix for cow 1 Row Col1 Col2 Col3 Col Covariance Parameter Estimates Standard Z Cov Parm Subject Estimate Error Value Pr > Z Alpha Lower Upper CS cow Residual Fit Statistics -2 Res Log Likelihood 41.3 AIC (Smaller is Better) 45.3 BIC (Smaller is Better) 44.1 Type 3 Tests of Fixed Effects Num Den Effect DF DF F Value Pr > F trt Effect trt Estimate Std Error DF t Value Pr > t trt <.0001 Effect trt _trt Estimate Std Error DF t Value Pr > t trt Latin Square 32
33 Unstructured WARNING: Stopped because of too many likelihood evaluations. Covariance Parameter Values At Last Iteration Cov Parm Subject Estimate UN(1,1) cow UN(2,1) cow UN(2,2) cow UN(3,1) cow UN(3,2) cow UN(3,3) cow UN(4,1) cow UN(4,2) cow UN(4,3) cow UN(4,4) cow **** SOME COVARIANCE STRUCTURES ARE TOO COMPLEX FOR THE DATA **** Latin Square 33
34 SAS Code /* Fit an AR(1) correlation structure for residuals*/ proc mixed covtest cl; class cow trt period; model resp=trt period / ddfm=kr outp=diag; repeated period / type=ar(1) subject=cow r=1; lsmeans trt / diff; run; /* Fit an AR(1) correlation structure for residuals*/ /* Add cow random effect */ proc mixed covtest cl; class cow trt period; model resp=trt period / ddfm=kr outp=diag; random intercept / subject=cow v=1; repeated period / type=ar(1) subject=cow r=1; lsmeans trt / diff; run; Latin Square 34
35 AR(1) Estimated R Matrix for cow 1 Row Col1 Col2 Col3 Col Covariance Parameter Estimates Standard Z Cov Parm Subject Estimate Error Value Pr Z Alpha Lower Upper AR(1) cow < Residual Fit Statistics -2 Res Log Likelihood 42.0 AIC (Smaller is Better) 46.0 BIC (Smaller is Better) 44.8 Type 3 Tests of Fixed Effects Num Den Effect DF DF F Value Pr > F trt Effect trt Estimate Std Error DF t Value Pr > t trt <.0001 Effect trt _trt Estimate Std Error DF t Value Pr > t trt trt Latin Square 35
36 AR(1) with random cow effect Estimated V Matrix for cow 1 Row Col1 Col2 Col3 Col Covariance Parameter Estimates Standard Z Cov Parm Subject Estimate Error Value Pr Z Alpha Lower Upper Intercept cow AR(1) cow Residual Fit Statistics -2 Res Log Likelihood 41.2 AIC (Smaller is Better) 47.2 BIC (Smaller is Better) 45.3 Type 3 Tests of Fixed Effects Num Den Effect DF DF F Value Pr > F trt Effect trt Estimate Std Error DF t Value Pr > t trt <.0001 Effect trt _trt Estimate Std Error DF t Value Pr > t trt trt Latin Square 36
37 SAS Code /* Fit an TOEP(2) correlation structure for residuals*/ proc mixed covtest cl; class cow trt period; model resp=trt period / ddfm=kr outp=diag; repeated period / type=toep(2) subject=cow r=1; lsmeans trt / diff; run; /* Fit an Heterogeneous CS correlation structure for residuals*/ proc mixed covtest cl; class cow trt period; model resp=trt period / ddfm=kr outp=diag; repeated period / type=csh subject=cow r=1; lsmeans trt / diff; run; Latin Square 37
38 Banded TOEP(2) Estimated R Matrix for cow 1 Row Col1 Col2 Col3 Col Covariance Parameter Estimates Standard Z Cov Parm Subject Estimate Error Value Pr Z Alpha Lower Upper TOEP(2) cow Residual Fit Statistics -2 Res Log Likelihood 45.0 AIC (Smaller is Better) 49.0 BIC (Smaller is Better) 47.8 Type 3 Tests of Fixed Effects Num Den Effect DF DF F Value Pr > F trt Effect trt Estimate Std Error DF t Value Pr > t trt <.0001 Effect trt _trt Estimate Std Error DF t Value Pr > t trt trt Latin Square 38
39 CSH Estimated R Matrix for cow 1 Row Col1 Col2 Col3 Col Covariance Parameter Estimates Cov Standard Z Parm Subject Estimate Error Value Pr Z Alpha Lower Upper Var(1) cow Var(2) cow Var(3) cow Var(4) cow CSH cow < Fit Statistics -2 Res Log Likelihood 34.8 AIC (Smaller is Better) 44.8 BIC (Smaller is Better) 41.7 Type 3 Tests of Fixed Effects Num Den Effect DF DF F Value Pr > F trt Effect trt Estimate Std Error DF t Value Pr > t trt <.0001 Effect trt _trt Estimate Std Error DF t Value Pr > t trt trt Latin Square 39
40 Summary When comparing models with a model selection criterion, REML fitting can only be used when fixed effects structure does not change. If different models with different fixed effects were compared, would have to use ML estimation, even though variance/covariance parameters are biased This is because REML uses a transformed set of data that takes into account loss of DF due to estimating fixed effects From the results here, heterogenous compound symmetry is the best choice provided residuals look ok Latin Square 40
Replicated Latin Square and Crossover Designs
Replicated Latin Square and Crossover Designs Replicated Latin Square Latin Square Design small df E, low power If 3 treatments 2 df error If 4 treatments 6 df error Can use replication to increase df
More informationSubject-specific observed profiles of change from baseline vs week trt=10000u
Mean of age 1 The MEANS Procedure Analysis Variable : age N Mean Std Dev Minimum Maximum ƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒƒ 109 55.5321101 12.1255537 26.0000000 83.0000000
More informationStatistical Consulting Topics. RCBD with a covariate
Statistical Consulting Topics RCBD with a covariate Goal: to determine the optimal level of feed additive to maximize the average daily gain of steers. VARIABLES Y = Average Daily Gain of steers for 160
More informationGLM Example: One-Way Analysis of Covariance
Understanding Design and Analysis of Research Experiments An animal scientist is interested in determining the effects of four different feed plans on hogs. Twenty four hogs of a breed were chosen and
More informationRANDOMIZED COMPLETE BLOCK DESIGN (RCBD) Probably the most used and useful of the experimental designs.
Description of the Design RANDOMIZED COMPLETE BLOCK DESIGN (RCBD) Probably the most used and useful of the experimental designs. Takes advantage of grouping similar experimental units into blocks or replicates.
More informationMixed Models Lecture Notes By Dr. Hanford page 151 More Statistics& SAS Tutorial at Type 3 Tests of Fixed Effects
Assessing fixed effects Mixed Models Lecture Notes By Dr. Hanford page 151 In our example so far, we have been concentrating on determining the covariance pattern. Now we ll look at the treatment effects
More informationPROC GLM AND PROC MIXED CODES FOR TREND ANALYSES FOR ROW-COLUMN DESIGNED EXPERIMENTS
PROC GLM AND PROC MIXED CODES FOR TREND ANALYSES FOR ROW-COLUMN DESIGNED EXPERIMENTS BU-1491-M June,2000 Walter T. Federer Dept. of Biometrics Cornell University Ithaca, NY 14853 wtfl@cornell.edu and Russell
More informationModelling Intervention Effects in Clustered Randomized Pretest/Posttest Studies. Ed Stanek
Modelling Intervention Effects in Clustered Randomized Pretest/Posttest Studies Introduction Ed Stanek We consider a study design similar to the design for the Well Women Project, and discuss analyses
More informationProblem Points Score USE YOUR TIME WISELY USE CLOSEST DF AVAILABLE IN TABLE SHOW YOUR WORK TO RECEIVE PARTIAL CREDIT
Stat 514 EXAM I Stat 514 Name (6 pts) Problem Points Score 1 32 2 30 3 32 USE YOUR TIME WISELY USE CLOSEST DF AVAILABLE IN TABLE SHOW YOUR WORK TO RECEIVE PARTIAL CREDIT WRITE LEGIBLY. ANYTHING UNREADABLE
More informationMixed Effects Models Yan Wang, Bristol-Myers Squibb, Wallingford, CT
PharmaSUG 2016 - Paper PO06 Mixed Effects Models Yan Wang, Bristol-Myers Squibb, Wallingford, CT ABSTRACT The MIXED procedure has been commonly used at the Bristol-Myers Squibb Company for quality of life
More informationRCBD with Sampling Pooling Experimental and Sampling Error
RCBD with Sampling Pooling Experimental and Sampling Error As we had with the CRD with sampling, we will have a source of variation for sampling error. Calculation of the Experimental Error df is done
More informationMixed models in R using the lme4 package Part 2: Longitudinal data, modeling interactions
Mixed models in R using the lme4 package Part 2: Longitudinal data, modeling interactions Douglas Bates 2011-03-16 Contents 1 sleepstudy 1 2 Random slopes 3 3 Conditional means 6 4 Conclusions 9 5 Other
More information1'-tq/? BU-- _-M August 2000 Technical Report Series of the Department of Biometrics, Cornell University, Ithaca, New York 14853
SAS/GLM AND SAS/MIXED FOR TREND ANALYSES US:ING FOURIER AND POLYNOMIAL REGRESSION FOR CENTERED AND NON-CENTERED VARIATES BY Walter T. Federer, Murari Singh, and Russell D. Wolfinger ABSTRACT Spatial variation
More informationLinear mixed models and when implied assumptions not appropriate
Mixed Models Lecture Notes By Dr. Hanford page 94 Generalized Linear Mixed Models (GLMM) GLMMs are based on GLM, extended to include random effects, random coefficients and covariance patterns. GLMMs are
More informationBlock Block Block
Advanced Biostatistics Quiz 3 Name March 16, 2005 9 or 10 Total Points Directions: Thoroughly, clearly and neatly answer the following two problems in the space given, showing all relevant calculations.
More informationExercises. ASReml Tutorial: B4 Bivariate Analysis p. 55
Exercises Coopworth data set - see Reference manual Five traits with varying amounts of data. No depth of pedigree (dams not linked to sires) Do univariate analyses Do bivariate analyses. Use COOP data
More informationI. Model. Q29a. I love the options at my fingertips today, watching videos on my phone, texting, and streaming films. Main Effect X1: Gender
1 Hopewell, Sonoyta & Walker, Krista COM 631/731 Multivariate Statistical Methods Dr. Kim Neuendorf Film & TV National Survey dataset (2014) by Jeffres & Neuendorf MANOVA Class Presentation I. Model INDEPENDENT
More informationModel II ANOVA: Variance Components
Model II ANOVA: Variance Components Model II MS A = s 2 + ns 2 A MS A MS W = ns 2 A (MS A MS W )/n = ns 2 A /n = s2 A Usually Expressed: s 2 A /(s2 A + s2 W ) x 100 Assumptions of ANOVA Random Sampling
More informationMANOVA/MANCOVA Paul and Kaila
I. Model MANOVA/MANCOVA Paul and Kaila From the Music and Film Experiment (Neuendorf et al.) Covariates (ONLY IN MANCOVA) X1 Music Condition Y1 E20 Contempt Y2 E21 Anticipation X2 Instrument Interaction
More informationSupplementary Figures Supplementary Figure 1 Comparison of among-replicate variance in invasion dynamics
1 Supplementary Figures Supplementary Figure 1 Comparison of among-replicate variance in invasion dynamics Scaled posterior probability densities for among-replicate variances in invasion speed (nine replicates
More informationChapter 27. Inferences for Regression. Remembering Regression. An Example: Body Fat and Waist Size. Remembering Regression (cont.)
Chapter 27 Inferences for Regression Copyright 2007 Pearson Education, Inc. Publishing as Pearson Addison-Wesley Slide 27-1 Copyright 2007 Pearson Education, Inc. Publishing as Pearson Addison-Wesley An
More informationTWO-FACTOR ANOVA Kim Neuendorf 4/9/18 COM 631/731 I. MODEL
1 TWO-FACTOR ANOVA Kim Neuendorf 4/9/18 COM 631/731 I. MODEL Using the Humor and Public Opinion Data, a two-factor ANOVA was run, using the full factorial model: MAIN EFFECT: Political Philosophy (3 groups)
More informationDo delay tactics affect silking date and yield of maize inbreds? Stephen Zimmerman Creative Component November 2015
Do delay tactics affect silking date and yield of maize inbreds? Stephen Zimmerman Creative Component November 2015 Overview Acknowledgements My Background Introduction Materials and Methods Results and
More informationResampling Statistics. Conventional Statistics. Resampling Statistics
Resampling Statistics Introduction to Resampling Probability Modeling Resample add-in Bootstrapping values, vectors, matrices R boot package Conclusions Conventional Statistics Assumptions of conventional
More informationMore About Regression
Regression Line for the Sample Chapter 14 More About Regression is spoken as y-hat, and it is also referred to either as predicted y or estimated y. b 0 is the intercept of the straight line. The intercept
More informationPaired plot designs experience and recommendations for in field product evaluation at Syngenta
Paired plot designs experience and recommendations for in field product evaluation at Syngenta 1. What are paired plot designs? 2. Analysis and reporting of paired plot designs 3. Case study 1 : analysis
More informationMANOVA COM 631/731 Spring 2017 M. DANIELS. From Jeffres & Neuendorf (2015) Film and TV Usage National Survey
1 MANOVA COM 631/731 Spring 2017 M. DANIELS I. MODEL From Jeffres & Neuendorf (2015) Film and TV Usage National Survey INDEPENDENT VARIABLES DEPENDENT VARIABLES X1: GENDER Q23a. I often watch a favorite
More informationMaster's thesis FACULTY OF SCIENCES Master of Statistics
2013 2014 FACULTY OF SCIENCES Master of Statistics Master's thesis Power linear
More informationK ABC Mplus CFA Model. Syntax file (kabc-mplus.inp) Data file (kabc-mplus.dat)
K ABC Mplus CFA Model Syntax file (kabc-mplus.inp) title: principles and practice of sem (4th ed.), rex kline two-factor model of the kabc-i, figure 9.7, table 13.1 data: file is "kabc-mplus.dat"; type
More informationDV: Liking Cartoon Comedy
1 Stepwise Multiple Regression Model Rikki Price Com 631/731 March 24, 2016 I. MODEL Block 1 Block 2 DV: Liking Cartoon Comedy 2 Block Stepwise Block 1 = Demographics: Item: Age (G2) Item: Political Philosophy
More informationBootstrap Methods in Regression Questions Have you had a chance to try any of this? Any of the review questions?
ICPSR Blalock Lectures, 2003 Bootstrap Resampling Robert Stine Lecture 3 Bootstrap Methods in Regression Questions Have you had a chance to try any of this? Any of the review questions? Getting class notes
More informationMoving on from MSTAT. March The University of Reading Statistical Services Centre Biometrics Advisory and Support Service to DFID
Moving on from MSTAT March 2000 The University of Reading Statistical Services Centre Biometrics Advisory and Support Service to DFID Contents 1. Introduction 3 2. Moving from MSTAT to Genstat 4 2.1 Analysis
More informationMultiple-point simulation of multiple categories Part 1. Testing against multiple truncation of a Gaussian field
Multiple-point simulation of multiple categories Part 1. Testing against multiple truncation of a Gaussian field Tuanfeng Zhang November, 2001 Abstract Multiple-point simulation of multiple categories
More informationAlgebra I Module 2 Lessons 1 19
Eureka Math 2015 2016 Algebra I Module 2 Lessons 1 19 Eureka Math, Published by the non-profit Great Minds. Copyright 2015 Great Minds. No part of this work may be reproduced, distributed, modified, sold,
More informationsubplots (30-m by 33-m) without space between potential subplots. Depending on the size of the
REM-S-13-00090 Online Supplemental Information Pyke et al. Appendix I Subplot Selection within Arid SageSTEP whole plots Each of the four whole plots (fuel reduction treatments) was gridded into potential
More informationSECTION I. THE MODEL. Discriminant Analysis Presentation~ REVISION Marcy Saxton and Jenn Stoneking DF1 DF2 DF3
Discriminant Analysis Presentation~ REVISION Marcy Saxton and Jenn Stoneking COM 631/731--Multivariate Statistical Methods Instructor: Prof. Kim Neuendorf (k.neuendorf@csuohio.edu) Cleveland State University,
More informationCOMP Test on Psychology 320 Check on Mastery of Prerequisites
COMP Test on Psychology 320 Check on Mastery of Prerequisites This test is designed to provide you and your instructor with information on your mastery of the basic content of Psychology 320. The results
More informationFor these items, -1=opposed to my values, 0= neutral and 7=of supreme importance.
1 Factor Analysis Jeff Spicer F1 F2 F3 F4 F9 F12 F17 F23 F24 F25 F26 F27 F29 F30 F35 F37 F42 F50 Factor 1 Factor 2 Factor 3 Factor 4 For these items, -1=opposed to my values, 0= neutral and 7=of supreme
More informationTI-Inspire manual 1. Real old version. This version works well but is not as convenient entering letter
TI-Inspire manual 1 Newest version Older version Real old version This version works well but is not as convenient entering letter Instructions TI-Inspire manual 1 General Introduction Ti-Inspire for statistics
More informationSociology 7704: Regression Models for Categorical Data Instructor: Natasha Sarkisian
OLS Regression Assumptions Sociology 7704: Regression Models for Categorical Data Instructor: Natasha Sarkisian A1. All independent variables are quantitative or dichotomous, and the dependent variable
More informationNormalization Methods for Two-Color Microarray Data
Normalization Methods for Two-Color Microarray Data 1/13/2009 Copyright 2009 Dan Nettleton What is Normalization? Normalization describes the process of removing (or minimizing) non-biological variation
More informationReproducibility Assessment of Independent Component Analysis of Expression Ratios from DNA microarrays.
Reproducibility Assessment of Independent Component Analysis of Expression Ratios from DNA microarrays. David Philip Kreil David J. C. MacKay Technical Report Revision 1., compiled 16th October 22 Department
More informationLecture 5: Clustering and Segmentation Part 1
Lecture 5: Clustering and Segmentation Part 1 Professor Fei Fei Li Stanford Vision Lab 1 What we will learn today Segmentation and grouping Gestalt principles Segmentation as clustering K means Feature
More informationDiscriminant Analysis. DFs
Discriminant Analysis Chichang Xiong Kelly Kinahan COM 631 March 27, 2013 I. Model Using the Humor and Public Opinion Data Set (Neuendorf & Skalski, 2010) IVs: C44 reverse coded C17 C22 C23 C27 reverse
More informationCS229 Project Report Polyphonic Piano Transcription
CS229 Project Report Polyphonic Piano Transcription Mohammad Sadegh Ebrahimi Stanford University Jean-Baptiste Boin Stanford University sadegh@stanford.edu jbboin@stanford.edu 1. Introduction In this project
More informationASReml Tutorial: C1 Variance structures p. 1. ASReml tutorial. C1 Variance structures. Arthur Gilmour
ASReml tutorial C1 Variance structures Arthur Gilmour ASReml Tutorial: C1 Variance structures p. 1 ASReml tutorial C1 Variance structures Arthur Gilmour ASReml Tutorial: C1 Variance structures p. 2 Overview
More informationHIGH-DIMENSIONAL CHANGEPOINT DETECTION
HIGH-DIMENSIONAL CHANGEPOINT DETECTION VIA SPARSE PROJECTION 3 6 8 11 14 16 19 22 26 28 31 33 35 39 43 47 48 52 53 56 60 63 67 71 73 77 80 83 86 88 91 93 96 98 101 105 109 113 114 118 120 121 125 126 129
More informationunbiased , is zero. Yï) + iab Fuller and Burmeister [4] suggested the estimator: N =Na +Nb + Nab Na +NB =Nb +NA.
RELTIVE EFFICIENCY OF SOME TWO -FRME ESTIMTORS H. Huang, Minnesota State Department of Education 1. Introduction In sample surveys, a complete frame is often unavailable or too expensive to construct.
More informationK3. Why did the certain ethnic mother put her baby in a crib with 20-foot high legs? So she could hear it if it fell out of bed.
Factor Analysis 1 COM 531, Spring 2008 K. Neuendorf MODEL: From Group Humor Data Set-- Responses to jokes: K1 K2 F1. F2. F3. F4. F5 K29 F6 K30 K31 For all items K1-K31, 0=not funny at all, 10=extremely
More informationK3. Why did the certain ethnic mother put her baby in a crib with 20-foot high legs? So she could hear it if it fell out of bed.
Factor Analysis 1 COM 531, Spring 2009 K. Neuendorf MODEL: From Group Humor Data Set-- Responses to jokes: K1 K2 F1. F2. F3. F4. F5 K29 F6 K30 K31 For all items K1-K31, 0=not funny at all, 10=extremely
More informationK-Pop Idol Industry Minhyung Lee
K-Pop Idol Industry 20100663 Minhyung Lee 1. K-Pop Idol History 2. Idol Industry Factor 3. Regression Analysis 4. Result & Interpretation K-Pop Idol History (1990s) Turning point of Korean Music history
More informationSensors, Measurement systems Signal processing and Inverse problems Exercises
Files: http://djafari.free.fr/cours/master_mne/cours/cours_mne_2014_01.pdf A. Mohammad-Djafari, Sensors, Measurement systems, Signal processing and Inverse problems, Master MNE 2014, 1/17. Sensors, Measurement
More informationDetecting Medicaid Data Anomalies Using Data Mining Techniques Shenjun Zhu, Qiling Shi, Aran Canes, AdvanceMed Corporation, Nashville, TN
Paper SDA-04 Detecting Medicaid Data Anomalies Using Data Mining Techniques Shenjun Zhu, Qiling Shi, Aran Canes, AdvanceMed Corporation, Nashville, TN ABSTRACT The purpose of this study is to use statistical
More informationa user's guide to Probit Or LOgit analysis
United States Department of Agriculture Forest Service Pacific Southwest Forest and Range Experiment Station General Technical Report PSW-38 a user's guide to Probit Or LOgit analysis Jacqueline L. Robertson
More informationComparison of Mixed-Effects Model, Pattern-Mixture Model, and Selection Model in Estimating Treatment Effect Using PRO Data in Clinical Trials
Comparison of Mixed-Effects Model, Pattern-Mixture Model, and Selection Model in Estimating Treatment Effect Using PRO Data in Clinical Trials Xiaolei Zhou, 1,2 Jianmin Wang, 1 Jessica Zhang, 1 Hongtu
More informationUPDATED STANDARDIZED CATCH RATES OF BLUEFIN TUNA (THUNNUS THYNNUS) FROM THE TRAP FISHERY IN TUNISIA
SCRS/2004/083 Col. Vol. Sci. Pap. ICCAT, 58(2): 596-605 (2005) UPDATED STANDARDIZED CATCH RATES OF BLUEFIN TUNA (THUNNUS THYNNUS) FROM THE TRAP FISHERY IN TUNISIA A. Hattour 1, J.M. de la Serna 2 and J.M
More informationBlueline, Linefree, Accuracy Ratio, & Moving Absolute Mean Ratio Charts
INTRODUCTION This instruction manual describes for users of the Excel Standard Celeration Template(s) the features of each page or worksheet in the template, allowing the user to set up and generate charts
More information1. Model. Discriminant Analysis COM 631. Spring Devin Kelly. Dataset: Film and TV Usage National Survey 2015 (Jeffres & Neuendorf) Q23a. Q23b.
1 Discriminant Analysis COM 631 Spring 2016 Devin Kelly 1. Model Dataset: Film and TV Usage National Survey 2015 (Jeffres & Neuendorf) Q23a. Q23b. Q23c. DF1 Q23d. Q23e. Q23f. Q23g. Q23h. DF2 DF3 CultClass
More informationBER margin of COM 3dB
BER margin of COM 3dB Yasuo Hidaka Fujitsu Laboratories of America, Inc. September 9, 2015 IEEE P802.3by 25 Gb/s Ethernet Task Force Abstract I was curious how much actual margin we have with COM 3dB So,
More informationAUTOREGRESSIVE MFCC MODELS FOR GENRE CLASSIFICATION IMPROVED BY HARMONIC-PERCUSSION SEPARATION
AUTOREGRESSIVE MFCC MODELS FOR GENRE CLASSIFICATION IMPROVED BY HARMONIC-PERCUSSION SEPARATION Halfdan Rump, Shigeki Miyabe, Emiru Tsunoo, Nobukata Ono, Shigeki Sagama The University of Tokyo, Graduate
More informationThe Time Series Forecasting System Charles Hallahan, Economic Research Service/USDA, Washington, DC
INTRODUCTION The Time Series Forecasting System Charles Hallahan, Economic Research Service/USDA, Washington, DC The Time Series Forecasting System (TSFS) is a component of SAS/ETS that provides a menu-based
More informationHIGH-DIMENSIONAL CHANGEPOINT ESTIMATION
HIGH-DIMENSIONAL CHANGEPOINT ESTIMATION VIA SPARSE PROJECTION 3 6 8 11 14 16 19 22 26 28 31 33 35 39 43 47 48 52 53 56 6 63 67 71 73 77 8 83 86 88 91 93 96 98 11 15 19 113 114 118 12 121 125 126 129 133
More informationRestoration of Hyperspectral Push-Broom Scanner Data
Restoration of Hyperspectral Push-Broom Scanner Data Rasmus Larsen, Allan Aasbjerg Nielsen & Knut Conradsen Department of Mathematical Modelling, Technical University of Denmark ABSTRACT: Several effects
More informationMixed Linear Models. Case studies on speech rate modulations in spontaneous speech. LSA Summer Institute 2009, UC Berkeley
Mixed Linear Models Case studies on speech rate modulations in spontaneous speech LSA Summer Institute 2009, UC Berkeley Florian Jaeger University of Rochester Managing speech rate How do speakers determine
More informationSTAT 503 Case Study: Supervised classification of music clips
STAT 503 Case Study: Supervised classification of music clips 1 Data Description This data was collected by Dr Cook from her own CDs. Using a Mac she read the track into the music editing software Amadeus
More informationMATH 214 (NOTES) Math 214 Al Nosedal. Department of Mathematics Indiana University of Pennsylvania. MATH 214 (NOTES) p. 1/3
MATH 214 (NOTES) Math 214 Al Nosedal Department of Mathematics Indiana University of Pennsylvania MATH 214 (NOTES) p. 1/3 CHAPTER 1 DATA AND STATISTICS MATH 214 (NOTES) p. 2/3 Definitions. Statistics is
More informationAnalysis of Packet Loss for Compressed Video: Does Burst-Length Matter?
Analysis of Packet Loss for Compressed Video: Does Burst-Length Matter? Yi J. Liang 1, John G. Apostolopoulos, Bernd Girod 1 Mobile and Media Systems Laboratory HP Laboratories Palo Alto HPL-22-331 November
More informationAP Statistics Sampling. Sampling Exercise (adapted from a document from the NCSSM Leadership Institute, July 2000).
AP Statistics Sampling Name Sampling Exercise (adapted from a document from the NCSSM Leadership Institute, July 2000). Problem: A farmer has just cleared a field for corn that can be divided into 100
More informationBox-Jenkins Methodology: Linear Time Series Analysis Using R
Box-Jenkins Methodology: Linear Time Series Analysis Using R Melody Ghahramani Mathematics & Statistics January 29, 2014 Melody Ghahramani (U of Winnipeg) R Seminar Series January 29, 2014 1 / 67 Outline
More informationHigh Speed Optical Networking: Task 3 FEC Coding, Channel Models, and Evaluations
1 Sponsored High Speed Optical Networking: Task 3 FEC Coding, Channel Models, and Evaluations Joel M. Morris, PhD Communications and Signal Processing Laboratory (CSPL) UMBC/CSEE Department 1000 Hilltop
More informationSTAT 113: Statistics and Society Ellen Gundlach, Purdue University. (Chapters refer to Moore and Notz, Statistics: Concepts and Controversies, 8e)
STAT 113: Statistics and Society Ellen Gundlach, Purdue University (Chapters refer to Moore and Notz, Statistics: Concepts and Controversies, 8e) Learning Objectives for Exam 1: Unit 1, Part 1: Population
More informationSociology 704: Topics in Multivariate Statistics Instructor: Natasha Sarkisian
Sociology 704: Topics in Multivariate Statistics Instructor: Natasha Sarkisian OLS Regression in Stata To run an OLS regression:. reg agekdbrn educ born sex mapres80 Source SS df MS Number of obs = 1091
More information8 Nonparametric test. Question 1: Are (expected) value of x and y the same?
Econometrics A: Tokyo International University 2017 autumn Satoshi OHIRA 26 8 Nonparametric test Question 1: Are (expected) value of x and y the same? One of the simplest way to answer the question is
More informationChapter 5. Describing Distributions Numerically. Finding the Center: The Median. Spread: Home on the Range. Finding the Center: The Median (cont.
Chapter 5 Describing Distributions Numerically Copyright 2007 Pearson Education, Inc. Publishing as Pearson Addison-Wesley Copyright 2007 Pearson Education, Inc. Publishing as Pearson Addison-Wesley Slide
More informationMATH& 146 Lesson 11. Section 1.6 Categorical Data
MATH& 146 Lesson 11 Section 1.6 Categorical Data 1 Frequency The first step to organizing categorical data is to count the number of data values there are in each category of interest. We can organize
More informationA STATISTICAL VIEW ON THE EXPRESSIVE TIMING OF PIANO ROLLED CHORDS
A STATISTICAL VIEW ON THE EXPRESSIVE TIMING OF PIANO ROLLED CHORDS Mutian Fu 1 Guangyu Xia 2 Roger Dannenberg 2 Larry Wasserman 2 1 School of Music, Carnegie Mellon University, USA 2 School of Computer
More informationRepeated measures ANOVA
Repeated measures ANOVA Pronoun interpretation in direct and indirect speech 07-05-2013 1 Franziska Köder Seminar in Methodology and Statistics, May 23, 2013 24-10-2012 2 Overview 1. Experimental design
More informationSupplemental Material: Color Compatibility From Large Datasets
Supplemental Material: Color Compatibility From Large Datasets Peter O Donovan, Aseem Agarwala, and Aaron Hertzmann Project URL: www.dgp.toronto.edu/ donovan/color/ 1 Unmixing color preferences In the
More informationMultiple-Window Spectrogram of Peaks due to Transients in the Electroencephalogram
284 IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING, VOL. 48, NO. 3, MARCH 2001 Multiple-Window Spectrogram of Peaks due to Transients in the Electroencephalogram Maria Hansson*, Member, IEEE, and Magnus Lindgren
More informationDetecting Musical Key with Supervised Learning
Detecting Musical Key with Supervised Learning Robert Mahieu Department of Electrical Engineering Stanford University rmahieu@stanford.edu Abstract This paper proposes and tests performance of two different
More information%CHCKFRQS A Macro Application for Generating Frequencies for QC and Simple Reports
Paper PO10 %CHCKFRQS A Macro Application for Generating Frequencies for QC and Simple Reports John Iwaniszek, MSc, Stat-Tech Services, LLC, Chapel Hill, NC Natalie Walker, Stat-Tech Services, LLC, Chapel
More informationOrthogonal rotation in PCAMIX
Orthogonal rotation in PCAMIX Marie Chavent 1,2, Vanessa Kuentz 3 and Jérôme Saracco 2,4 1 Université de Bordeaux, IMB, CNRS, UMR 5251, France 2 INRIA Bordeaux Sud-Ouest, CQFD team, France 3 CEMAGREF,
More informationECONOMICS 351* -- INTRODUCTORY ECONOMETRICS. Queen's University Department of Economics. ECONOMICS 351* -- Winter Term 2005 INTRODUCTORY ECONOMETRICS
Queen's University Department of Economics ECONOMICS 351* -- Winter Term 2005 INTRODUCTORY ECONOMETRICS Winter Term 2005 Instructor: Web Site: Mike Abbott Office: Room A521 Mackintosh-Corry Hall or Room
More informationRelease Year Prediction for Songs
Release Year Prediction for Songs [CSE 258 Assignment 2] Ruyu Tan University of California San Diego PID: A53099216 rut003@ucsd.edu Jiaying Liu University of California San Diego PID: A53107720 jil672@ucsd.edu
More informationOptimum Frame Synchronization for Preamble-less Packet Transmission of Turbo Codes
! Optimum Frame Synchronization for Preamble-less Packet Transmission of Turbo Codes Jian Sun and Matthew C. Valenti Wireless Communications Research Laboratory Lane Dept. of Comp. Sci. & Elect. Eng. West
More informationJournal of Statistical Software
JSS Journal of Statistical Software MMMMMM YYYY, Volume VV, Issue II. doi: 10.18637/jss.v000.i00 The canonical approach to subset selection problem is to choose k out of p predictors for each k {0, 1,
More informationWEB APPENDIX. Managing Innovation Sequences Over Iterated Offerings: Developing and Testing a Relative Innovation, Comfort, and Stimulation
WEB APPENDIX Managing Innovation Sequences Over Iterated Offerings: Developing and Testing a Relative Innovation, Comfort, and Stimulation Framework of Consumer Responses Timothy B. Heath Subimal Chatterjee
More informationPredictability of Music Descriptor Time Series and its Application to Cover Song Detection
Predictability of Music Descriptor Time Series and its Application to Cover Song Detection Joan Serrà, Holger Kantz, Xavier Serra and Ralph G. Andrzejak Abstract Intuitively, music has both predictable
More informationAP Statistics Sec 5.1: An Exercise in Sampling: The Corn Field
AP Statistics Sec.: An Exercise in Sampling: The Corn Field Name: A farmer has planted a new field for corn. It is a rectangular plot of land with a river that runs along the right side of the field. The
More informationModeling television viewership
Modeling television viewership The Nielsen ratings are the best known measures of viewership of television shows. These ratings form the basis for the setting of advertising rates, and are thus crucial
More informationLOCOCODE versus PCA and ICA. Jurgen Schmidhuber. IDSIA, Corso Elvezia 36. CH-6900-Lugano, Switzerland. Abstract
LOCOCODE versus PCA and ICA Sepp Hochreiter Technische Universitat Munchen 80290 Munchen, Germany Jurgen Schmidhuber IDSIA, Corso Elvezia 36 CH-6900-Lugano, Switzerland Abstract We compare the performance
More informationTutorial 0: Uncertainty in Power and Sample Size Estimation. Acknowledgements:
Tutorial 0: Uncertainty in Power and Sample Size Estimation Anna E. Barón, Keith E. Muller, Sarah M. Kreidler, and Deborah H. Glueck Acknowledgements: The project was supported in large part by the National
More informationDesign and Analysis of New Methods on Passive Image Forensics. Advisor: Fernando Pérez-González. Signal Theory and Communications Department
Design and Analysis of New Methods on Passive Image Forensics Advisor: Fernando Pérez-González GPSC Signal Processing and Communications Group Vigo. November 8, 3. Why do we need Image Forensics? Because...
More informationHidden Markov Model based dance recognition
Hidden Markov Model based dance recognition Dragutin Hrenek, Nenad Mikša, Robert Perica, Pavle Prentašić and Boris Trubić University of Zagreb, Faculty of Electrical Engineering and Computing Unska 3,
More informationOpen Access Determinants and the Effect on Article Performance
International Journal of Business and Economics Research 2017; 6(6): 145-152 http://www.sciencepublishinggroup.com/j/ijber doi: 10.11648/j.ijber.20170606.11 ISSN: 2328-7543 (Print); ISSN: 2328-756X (Online)
More informationStudy of the Clocking Effect in the TRT Alignment
Draft version 2. Study of the Clocking Effect in the TRT Alignment John Alison Aart Heijboer Joel Heinrich Joe Kroll University of Pennsylvania Andrea Bocci Duke University April 15, 28 Abstract This paper
More informationSelling the Premium in the Freemium: Impact of Product Line Extensions
Selling the Premium in the Freemium: Impact of Product Line Extensions Xian Gu 1 P. K. Kannan Liye Ma August 2017 1 Xian Gu is Doctoral Candidate in Marketing, P. K. Kannan is Dean s Chair in Marketing
More informationarxiv: v1 [cs.dl] 9 May 2017
Understanding the Impact of Early Citers on Long-Term Scientific Impact Mayank Singh Dept. of Computer Science and Engg. IIT Kharagpur, India mayank.singh@cse.iitkgp.ernet.in Ajay Jaiswal Dept. of Computer
More informationPlacement Rent Exponent Calculation Methods, Temporal Behaviour, and FPGA Architecture Evaluation. Joachim Pistorius and Mike Hutton
Placement Rent Exponent Calculation Methods, Temporal Behaviour, and FPGA Architecture Evaluation Joachim Pistorius and Mike Hutton Some Questions How best to calculate placement Rent? Are there biases
More informationThe Great Beauty: Public Subsidies in the Italian Movie Industry
The Great Beauty: Public Subsidies in the Italian Movie Industry G. Meloni, D. Paolini,M.Pulina April 20, 2015 Abstract The aim of this paper to examine the impact of public subsidies on the Italian movie
More information