ST706: Spring One-Way Random effects example.
|
|
- Adrian Conley
- 5 years ago
- Views:
Transcription
1 ST706: Spring One-Way Random effects example. DATA IS FROM CONTROL ARM OF SKIN CANCER STUDY AT DARTMOUTH. DC0 TO DC5 ARE MEASURES OF DIETARY INTAKE OVER 6 YEARS. For illustration here we treat the six years as if they were replicates. Later we will allow year effects. The data has all of the reps on one line so we first create a file with one line for each rep (for each person). This analysis is for the 154 people all with six years of data. n = 6, m = 154. Since the data is balanced some exact confidence intervals are available. But, the analysis here assumes normality and constant within subject variance, which is clearly suspect here. See figure at end. Below is the analysis of variance table with MSA = and MSE = σ 2 = and σ a 2 = (MSA MSE)/6 = A 95% confidence interval for σ 2 based on the chi-square with 770 degrees of freedom is ( , ) An approximate 95% confidence interval for σa 2 is ( , ). This is based on Satterthwaite s method which uses here a chi-square with degrees of freedom. A Wald confidence interval for σa 2 is (2.6234, ) Estimate of population mean µ is Ȳ.. = with an estimated standard error of (MSA/nm) 1/2 = and a 95% confidence interval of ( , ). options linesize=90; /* DATA IS FROM CONTROL ARM OF SKIN CANCER STUDY AT DARTMOUTH. DC0 TO DC5 ARE MEASURES OF DIETARY INTAKE OVER 6 YEARS. HERE WE TREAT THE 6 YEARS AS REPLICATES FOR ILLUSTRATION. */ data a; infile f:\s506\dietrep.dat ; input id age sex smkestat quetelet labcar0 labcar1 labcar2 labcar3 labcar4 labcar5 dc0 dc1 dc2 dc3 dc4 dc5; nmiss = nmiss(of dc0-dc5); diet=dc0; rep=0;output; diet=dc1;rep=1;output; diet=dc2;rep=2;output; diet=dc3;rep=3;output; diet=dc4;rep=4;output; diet=dc5;rep=5;output; /* THIS CREATES A DATA SET WITH JUST THE ONES HAVING ALL 6 YEARS OF DATA*/ data balanced;set a; if nmiss=0; proc means data=balanced; var diet; The MEANS Procedure 1
2 Analysis Variable : diet N Mean Std Dev Minimum Maximum * title original data unbalanced ; * proc glm data=a; /* UNCOMMENT TWO LINES ABOVE AND COMMENT OUT THE NEXT TWO TO RUN WITH FULL UNBALANCED DATA */ title original data. Balanced ; proc glm data=balanced; class id; model diet=id; Sum of Source DF Squares Mean Square F Value Pr > F Model <.0001 Error Corrected Total id <.0001 USING VARCOMP proc varcomp; class id; model diet=id; proc varcomp method=reml; class id; model diet=id; Variance Components Estimation Procedure MIVQUE(0) Estimates Variance Component diet Var(id) Var(Error) Variance Component Estimate Var(id) Var(Error) Asymptotic Covariance Matrix of Estimates Var(id) Var(Error) Var(id) Var(Error) If you run varcomp with method=type1, it will display the expected mean squares for this example. Type 1 Analysis of Variance Sum of Source DF Squares Mean Square Expected Mean Square id Var(Error) + 6 Var(id) Error Var(Error) Corrected Total
3 /* USING PROC MIXED */ proc mixed asycov covtest cl nobound; /* default method in proc mixed is reml*/ class id; model diet=/ ddfm=satterth solution cl; random id/cl ; The Mixed Procedure Data Set WORK.BALANCED Dependent Variable diet Covariance Structure Variance Components Estimation Method REML Residual Variance Method Profile Fixed Effects SE Method Model-Based Degrees of Freedom Method Satterthwaite Covariance Parameter Estimates Standard Z Cov Parm Estimate Error Value Pr Z Alpha Lower Upper id < Residual < Asymptotic Covariance Matrix of Estimates Row Cov Parm CovP1 CovP2 1 id Residual Solution for Fixed Effects Standard Effect Estimate Error DF t Value Pr > t Alpha Lower Upper Intercept < Solution for Random Effects Std Err Effect id Estimate Pred DF t Value Pr > t Alpha Lower Upper id id id Additional comments about proc mixed. The random statement indicates there is a random effect associated with id. Note that the model is diet = since there is just an overall mean (i.e, intercept) in the fixed effects part of the model. In the model statement the ddfm = satterth is referring to the use of a Satterthwaite method associated with getting a degrees of freedom to go into the confidence interval for the fixed effect µ. Here it gives 3
4 the exact interval based on 153 degrees of freedom. You get the same answer if you use no ddfm= option here. With only a random statement (and no use of a repeated statement which we will discuss later), the program automatically incorporates the ǫ and the variance for this is referred to as residual variability. The results for σ 2 are the ones associated with the label Residual. Those for σa 2 are associated with the label id (since id is the grouping variable). The confidence interval given for σ 2 is based on the chi-square with m (n 1) = = 770 degrees of freedom. The proc mixed above had a nobound in the proc mixed statement. This allows for a negative estimate of a variance component. This option doesn t come into play here with respect to point estimation since the estimate of σ 2 a is positive. However, since nobound is used, the confidence interval for σ2 a is based on a Wald interval. With nobound, all intervals for all variance components associated with the random statement are computed using the Wald method. To get the Satterthwaite based confidence interval for σ 2 a, which we computed directly above, eliminate the nobound option as below. The interval associated with id is now the Satterthwaite interval for σ 2 a. proc mixed asycov covtest cl; class id; model diet=/ solution cl; random id; Standard Z Cov Parm Estimate Error Value Pr Z Alpha Lower Upper id < Residual < Doing some direct computing to double check what is getting done in varcomp and mixed. proc iml; n=6; m=154; dfe=r*(n-1); alpha=.05; MSA = ; MSE = ; sigma2hat = mse; sigmaa2hat = (msa-mse)/n; cvalue1 = cinv(1-(alpha/2),dfe); cvalue2 = cinv(alpha/2,dfe); lowsigma2 = (m*(n-1)*mse)/cvalue1; upsigma2= (m*(n-1)*mse)/cvalue2; print Estimation of sigma2 sigma2hat dfe lowsigma2 upsigma2; /* Computing the chi-square interval for sigma2mu using satterthwaite approximate degrees of freedom */ den = (msa**2/(m-1)) + (mse**2/(m*(n-1))); dfmu = (n*sigmaa2hat)**2/den; cvalmu1 = cinv(1-(alpha/2),dfmu); cvalmu2 = cinv(alpha/2,dfmu); 4
5 lowsigmaa2 = (dfmu*sigmaa2hat)/cvalmu1; upsigmaa2 = (dfmu*sigmaa2hat)/cvalmu2; print Estimation of sigmaa2 sigmaa2hat dfmu lowsigmaa2 upsigmaa2; /* Estimation of mu */ ybar = ; seybar = sqrt(msa/(m*n)); tval = tinv(1 - (alpha/2),r-1); lowermu = ybar -tval*seybar; uppermu = ybar + tval*seybar; print Estimate of population mean ybar seybar lowermu uppermu; quit; SIGMA2HAT DFE LOWSIGMA2 UPSIGMA2 Estimation of sigma SIGMAA2HAT DFMU LOWSIGMAA2 UPSIGMAA2 Estimation of sigmaa YBAR SEYBAR LOWERMU UPPERMU Estimate of population mean Figure 1: Measure of diet betacarotene intake. 5
6 Figure 2: Plot of id st.deviation versus id mean. 6 replicates 6
Repeated Measures Part 4: Blood Flow data
Repeated Measures Part 4: Blood Flow data /* bloodflow.sas */ options linesize=79 pagesize=100 noovp formdlim='_'; title 'Two within-subjecs factors: Blood flow data (NWK p. 1181)'; proc format; value
More informationSTAT 5200 Handout #25. R-Square & Design Matrix in Mixed Models
STAT 5200 Handout #25 R-Square & Design Matrix in Mixed Models I. R-Square in Mixed Models (with Example from Handout #20): For mixed models, the concept of R 2 is a little complicated (and neither PROC
More informationStat 500 lab notes c Philip M. Dixon, Week 10: Autocorrelated errors
Week 10: Autocorrelated errors This week, I have done one possible analysis and provided lots of output for you to consider. Case study: predicting body fat Body fat is an important health measure, but
More informationModule 3: SAS. 3.1 Initial explorative analysis 02429/MIXED LINEAR MODELS PREPARED BY THE STATISTICS GROUPS AT IMM, DTU AND KU-LIFE
St@tmaster 02429/MIXED LINEAR MODELS PREPARED BY THE STATISTICS GROUPS AT IMM, DTU AND KU-LIFE Module 3: SAS 3.1 Initial explorative analysis....................... 1 3.1.1 SAS JMP............................
More informationGeneralized Least Squares (GLS) and Estimated Generalized Least Squares (EGLS)
Generalized Least Squares (GLS) and Estimated Generalized Least Squares (EGLS) Linear Model in matrix notation for the population Y = Xβ + Var ( ) = In GLS, the error covariance matrix is known In EGLS
More informationAn introduction to SPSS
An introduction to SPSS To open the SPSS software using U of Iowa Virtual Desktop... Go to https://virtualdesktop.uiowa.edu and choose SPSS 24. Contents NOTE: Save data files in a drive that is accessible
More informationFactorial ANOVA with SAS
Factorial ANOVA with SAS /* potato305.sas */ options linesize=79 noovp formdlim='_' ; title 'Rotten potatoes'; title2 ''; proc format; value tfmt 1 = 'Cool' 2 = 'Warm'; data spud; infile 'potato2.data'
More informationTHIS IS NOT REPRESNTATIVE OF CURRENT CLASS MATERIAL. STOR 455 Midterm 1 September 28, 2010
THIS IS NOT REPRESNTATIVE OF CURRENT CLASS MATERIAL STOR 455 Midterm September 8, INSTRUCTIONS: BOTH THE EXAM AND THE BUBBLE SHEET WILL BE COLLECTED. YOU MUST PRINT YOUR NAME AND SIGN THE HONOR PLEDGE
More informationCentering and Interactions: The Training Data
Centering and Interactions: The Training Data A random sample of 150 technical support workers were first given a test of their technical skill and knowledge, and then randomly assigned to one of three
More informationby John Barnard and Walter T. Federer Cornell University
Genstat and SAS Programs for Recovering Interblock, Interrow-column, and Intergradient Information by John Barnard and Walter T. Federer Cornell University BU-1383-M 1 anuary 1997 0. Abstrast Genstat and
More informationConditional and Unconditional Regression with No Measurement Error
Conditional and with No Measurement Error /* reg2ways.sas */ %include 'readsenic.sas'; title2 ''; proc reg; title3 'Conditional Regression'; model infrisk = stay census; proc calis cov; /* Analyze the
More informationPLS205 Lab 1 January 9, Laboratory Topics 1 & 2
PLS205 Lab 1 January 9, 2014 Laboratory Topics 1 & 2 Welcome, introduction, logistics, and organizational matters Introduction to SAS Writing and running programs saving results checking for errors Different
More informationFactorial ANOVA. Skipping... Page 1 of 18
Factorial ANOVA The potato data: Batches of potatoes randomly assigned to to be stored at either cool or warm temperature, infected with one of three bacterial types. Then wait a set period. The dependent
More information1 The SAS System 23:01 Friday, November 9, 2012
2101f12HW9chickwts.log Saved: Wednesday, November 14, 2012 6:50:49 PM Page 1 of 3 1 The SAS System 23:01 Friday, November 9, 2012 NOTE: Copyright (c) 2002-2010 by SAS Institute Inc., Cary, NC, USA. NOTE:
More informationLab #9: ANOVA and TUKEY tests
Lab #9: ANOVA and TUKEY tests Objectives: 1. Column manipulation in SAS 2. Analysis of variance 3. Tukey test 4. Least Significant Difference test 5. Analysis of variance with PROC GLM 6. Levene test for
More informationLaboratory Topics 1 & 2
PLS205 Lab 1 January 12, 2012 Laboratory Topics 1 & 2 Welcome, introduction, logistics, and organizational matters Introduction to SAS Writing and running programs; saving results; checking for errors
More informationSAS data statements and data: /*Factor A: angle Factor B: geometry Factor C: speed*/
STAT:5201 Applied Statistic II (Factorial with 3 factors as 2 3 design) Three-way ANOVA (Factorial with three factors) with replication Factor A: angle (low=0/high=1) Factor B: geometry (shape A=0/shape
More informationWithin-Cases: Multivariate approach part one
Within-Cases: Multivariate approach part one /* sleep2.sas */ options linesize=79 noovp formdlim=' '; title "Student's Sleep data: Matched t-tests with proc reg"; data bedtime; infile 'studentsleep.data'
More informationLinear Mixed Models: Methodology and Algorithms
Linear Mixed Models: Methodology and Algorithms David M. Allen University of Kentucky April 10, 2017 18SAS Ddfm Options You have been encouraged to use the Satterthwaite option for degrees of freedom.
More informationIntroductory Guide to SAS:
Introductory Guide to SAS: For UVM Statistics Students By Richard Single Contents 1 Introduction and Preliminaries 2 2 Reading in Data: The DATA Step 2 2.1 The DATA Statement............................................
More informationConfirmatory Factor Analysis on the Twin Data: Try One
Confirmatory Factor Analysis on the Twin Data: Try One /************ twinfac2.sas ********************/ TITLE2 'Confirmatory Factor Analysis'; %include 'twinread.sas'; proc calis corr; /* Analyze the correlation
More informationSAS/STAT 14.2 User s Guide. The PLM Procedure
SAS/STAT 14.2 User s Guide The PLM Procedure This document is an individual chapter from SAS/STAT 14.2 User s Guide. The correct bibliographic citation for this manual is as follows: SAS Institute Inc.
More informationST512. Fall Quarter, Exam 1. Directions: Answer questions as directed. Please show work. For true/false questions, circle either true or false.
ST512 Fall Quarter, 2005 Exam 1 Name: Directions: Answer questions as directed. Please show work. For true/false questions, circle either true or false. 1. (42 points) A random sample of n = 30 NBA basketball
More informationOutline. Topic 16 - Other Remedies. Ridge Regression. Ridge Regression. Ridge Regression. Robust Regression. Regression Trees. Piecewise Linear Model
Topic 16 - Other Remedies Ridge Regression Robust Regression Regression Trees Outline - Fall 2013 Piecewise Linear Model Bootstrapping Topic 16 2 Ridge Regression Modification of least squares that addresses
More informationStat 302 Statistical Software and Its Applications SAS: Data I/O
Stat 302 Statistical Software and Its Applications SAS: Data I/O Yen-Chi Chen Department of Statistics, University of Washington Autumn 2016 1 / 33 Getting Data Files Get the following data sets from the
More informationWeek 4: Simple Linear Regression III
Week 4: Simple Linear Regression III Marcelo Coca Perraillon University of Colorado Anschutz Medical Campus Health Services Research Methods I HSMP 7607 2017 c 2017 PERRAILLON ARR 1 Outline Goodness of
More informationMixed Effects Models. Biljana Jonoska Stojkova Applied Statistics and Data Science Group (ASDa) Department of Statistics, UBC.
Mixed Effects Models Biljana Jonoska Stojkova Applied Statistics and Data Science Group (ASDa) Department of Statistics, UBC March 6, 2018 Resources for statistical assistance Department of Statistics
More informationSet up of the data is similar to the Randomized Block Design situation. A. Chang 1. 1) Setting up the data sheet
Repeated Measure Analysis (Univariate Mixed Effect Model Approach) (Treatment as the Fixed Effect and the Subject as the Random Effect) (This univariate approach can be used for randomized block design
More informationR-Square Coeff Var Root MSE y Mean
STAT:50 Applied Statistics II Exam - Practice 00 possible points. Consider a -factor study where each of the factors has 3 levels. The factors are Diet (,,3) and Drug (A,B,C) and there are n = 3 observations
More informationSAS Programs and Output for Alternative Experimentals Designs. Section 13-8 in Howell (2010)
SAS Programs and Output for Alternative Experimentals Designs. Section 13-8 in Howell (2010) In Statistical Methods for Psychology (7 th ed.) I discuss various alternative experimental designs involving
More informationThe SAS Denomininator Degrees of Freedom Option
The SAS Denomininator Degrees of Freedom Option David Allen University of Kentucky April 12, 2016 1 SAS Degrees of Freedom Options This discussion pertains to SAS proc mixed or proc glimmix when used for
More informationContrasts and Multiple Comparisons
Contrasts and Multiple Comparisons /* onewaymath.sas */ title2 'Oneway with contrasts and multiple comparisons (Exclude Other/DK)'; %include 'readmath.sas'; if ethnic ne 6; /* Otherwise, throw the case
More informationEXST3201 Mousefeed01 Page 1
EXST3201 Mousefeed01 Page 1 3 /* 4 Examine differences among the following 6 treatments 5 N/N85 fed normally before weaning and 85 kcal/wk after 6 N/R40 fed normally before weaning and 40 kcal/wk after
More informationSTA431 Handout 9 Double Measurement Regression on the BMI Data
STA431 Handout 9 Double Measurement Regression on the BMI Data /********************** bmi5.sas **************************/ options linesize=79 pagesize = 500 noovp formdlim='-'; title 'BMI and Health:
More informationIntermediate SAS: Statistics
Intermediate SAS: Statistics OIT TSS 293-4444 oithelp@mail.wvu.edu oit.wvu.edu/training/classmat/sas/ Table of Contents Procedures... 2 Two-sample t-test:... 2 Paired differences t-test:... 2 Chi Square
More informationStat 302 Statistical Software and Its Applications SAS: Data I/O & Descriptive Statistics
Stat 302 Statistical Software and Its Applications SAS: Data I/O & Descriptive Statistics Fritz Scholz Department of Statistics, University of Washington Winter Quarter 2015 February 19, 2015 2 Getting
More informationBiostat Methods STAT 5820/6910 Handout #9 Meta-Analysis Examples
Biostat Methods STAT 5820/6910 Handout #9 Meta-Analysis Examples Example 1 A RCT was conducted to consider whether steroid therapy for expectant mothers affects death rate of premature [less than 37 weeks]
More informationPoisson Regression and Model Checking
Poisson Regression and Model Checking Readings GH Chapter 6-8 September 27, 2017 HIV & Risk Behaviour Study The variables couples and women_alone code the intervention: control - no counselling (both 0)
More informationESTIMATING DENSITY DEPENDENCE, PROCESS NOISE, AND OBSERVATION ERROR
ESTIMATING DENSITY DEPENDENCE, PROCESS NOISE, AND OBSERVATION ERROR Coinvestigators: José Ponciano, University of Idaho Subhash Lele, University of Alberta Mark Taper, Montana State University David Staples,
More informationIntroduction to Mixed Models: Multivariate Regression
Introduction to Mixed Models: Multivariate Regression EPSY 905: Multivariate Analysis Spring 2016 Lecture #9 March 30, 2016 EPSY 905: Multivariate Regression via Path Analysis Today s Lecture Multivariate
More informationPage 1 of 8. Language Development Study
Page 1 of 8 Language Development Study /* cread3.sas Read Castilla's language development data */ options nodate linesize=79 noovp formdlim=' '; title "Castilla's Language development Study"; proc format;
More informationWeek 6, Week 7 and Week 8 Analyses of Variance
Week 6, Week 7 and Week 8 Analyses of Variance Robyn Crook - 2008 In the next few weeks we will look at analyses of variance. This is an information-heavy handout so take your time reading it, and don
More informationThe Kenton Study. (Applied Linear Statistical Models, 5th ed., pp , Kutner et al., 2005) Page 1 of 5
The Kenton Study The Kenton Food Company wished to test four different package designs for a new breakfast cereal. Twenty stores, with approximately equal sales volumes, were selected as the experimental
More informationIntroduction to SAS proc calis
Introduction to SAS proc calis /* path1.sas */ %include 'SenicRead.sas'; title2 ''; /************************************************************************ * * * Cases are hospitals * * * * stay Average
More informationBivariate (Simple) Regression Analysis
Revised July 2018 Bivariate (Simple) Regression Analysis This set of notes shows how to use Stata to estimate a simple (two-variable) regression equation. It assumes that you have set Stata up on your
More informationSTAT 5200 Handout #21. Split Split Plot Example (Ch. 16)
STAT 5200 Handout #21 Split Split Plot Example (Ch. 16) Build off the Split-Plot Design: (see Handouts #20 and #20a) We can add a third treatment factor and split the split-plot units into split-split-plot
More informationCSSS 510: Lab 2. Introduction to Maximum Likelihood Estimation
CSSS 510: Lab 2 Introduction to Maximum Likelihood Estimation 2018-10-12 0. Agenda 1. Housekeeping: simcf, tile 2. Questions about Homework 1 or lecture 3. Simulating heteroskedastic normal data 4. Fitting
More informationT-test og variansanalyse i SAS. T-test og variansanalyse i SAS p.1/18
T-test og variansanalyse i SAS T-test og variansanalyse i SAS p.1/18 T-test og variansanalyse i SAS T-test (Etstik, tostik, parrede observationer) Variansanalyse SAS-procedurer: PROC TTEST PROC GLM T-test
More informationXO = X +(JQ -Y)/ b I. y= y +~(X -X)+E. 174 Statistics, Data Analysis and Modeling. Proceedings of MWSUG '95
DOING INVERSE REGRESSION ESTIMATION VIA SAS/IML E. C. Drago and T. G. Filloon Abstract Inverse regression (also known as linear calibration) is used to determine when/where a linear regression line (±
More informationSTAT:5201 Applied Statistic II
STAT:5201 Applied Statistic II Two-Factor Experiment (one fixed blocking factor, one fixed factor of interest) Randomized complete block design (RCBD) Primary Factor: Day length (short or long) Blocking
More informationPredictive Checking. Readings GH Chapter 6-8. February 8, 2017
Predictive Checking Readings GH Chapter 6-8 February 8, 2017 Model Choice and Model Checking 2 Questions: 1. Is my Model good enough? (no alternative models in mind) 2. Which Model is best? (comparison
More informationSAS/STAT 13.1 User s Guide. The SURVEYFREQ Procedure
SAS/STAT 13.1 User s Guide The SURVEYFREQ Procedure This document is an individual chapter from SAS/STAT 13.1 User s Guide. The correct bibliographic citation for the complete manual is as follows: SAS
More informationCH5: CORR & SIMPLE LINEAR REFRESSION =======================================
STAT 430 SAS Examples SAS5 ===================== ssh xyz@glue.umd.edu, tap sas913 (old sas82), sas https://www.statlab.umd.edu/sasdoc/sashtml/onldoc.htm CH5: CORR & SIMPLE LINEAR REFRESSION =======================================
More informationOne Factor Experiments
One Factor Experiments 20-1 Overview Computation of Effects Estimating Experimental Errors Allocation of Variation ANOVA Table and F-Test Visual Diagnostic Tests Confidence Intervals For Effects Unequal
More informationRegression Analysis and Linear Regression Models
Regression Analysis and Linear Regression Models University of Trento - FBK 2 March, 2015 (UNITN-FBK) Regression Analysis and Linear Regression Models 2 March, 2015 1 / 33 Relationship between numerical
More informationBaruch College STA Senem Acet Coskun
Baruch College STA 9750 BOOK BUY A Predictive Mode Senem Acet Coskun Table of Contents Summary 3 Why this topic? 4 Data Sources 6 Variable Definitions 7 Descriptive Statistics 8 Univariate Analysis 9 Two-Sample
More informationEcon 3790: Business and Economics Statistics. Instructor: Yogesh Uppal
Econ 3790: Business and Economics Statistics Instructor: Yogesh Uppal Email: yuppal@ysu.edu Chapter 8: Interval Estimation Population Mean: Known Population Mean: Unknown Margin of Error and the Interval
More informationStat 5100 Handout #14.a SAS: Logistic Regression
Stat 5100 Handout #14.a SAS: Logistic Regression Example: (Text Table 14.3) Individuals were randomly sampled within two sectors of a city, and checked for presence of disease (here, spread by mosquitoes).
More informationNCSS Statistical Software. Design Generator
Chapter 268 Introduction This program generates factorial, repeated measures, and split-plots designs with up to ten factors. The design is placed in the current database. Crossed Factors Two factors are
More informationWeek 5: Multiple Linear Regression II
Week 5: Multiple Linear Regression II Marcelo Coca Perraillon University of Colorado Anschutz Medical Campus Health Services Research Methods I HSMP 7607 2017 c 2017 PERRAILLON ARR 1 Outline Adjusted R
More informationPackage simr. April 30, 2018
Type Package Package simr April 30, 2018 Title Power Analysis for Generalised Linear Mixed Models by Simulation Calculate power for generalised linear mixed models, using simulation. Designed to work with
More informationWeek 4: Simple Linear Regression II
Week 4: Simple Linear Regression II Marcelo Coca Perraillon University of Colorado Anschutz Medical Campus Health Services Research Methods I HSMP 7607 2017 c 2017 PERRAILLON ARR 1 Outline Algebraic properties
More informationIntroduction to SAS proc calis
Introduction to SAS proc calis /* path1.sas */ %include 'SenicRead.sas'; title2 'Path Analysis Example for 3 Observed Variables'; /************************************************************************
More informationSAS/STAT 14.3 User s Guide The SURVEYFREQ Procedure
SAS/STAT 14.3 User s Guide The SURVEYFREQ Procedure This document is an individual chapter from SAS/STAT 14.3 User s Guide. The correct bibliographic citation for this manual is as follows: SAS Institute
More information2.830J / 6.780J / ESD.63J Control of Manufacturing Processes (SMA 6303) Spring 2008
MIT OpenCourseWare http://ocw.mit.edu.83j / 6.78J / ESD.63J Control of Manufacturing Processes (SMA 633) Spring 8 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationCSC 328/428 Summer Session I 2002 Data Analysis for the Experimenter FINAL EXAM
options pagesize=53 linesize=76 pageno=1 nodate; proc format; value $stcktyp "1"="Growth" "2"="Combined" "3"="Income"; data invstmnt; input stcktyp $ perform; label stkctyp="type of Stock" perform="overall
More informationIQR = number. summary: largest. = 2. Upper half: Q3 =
Step by step box plot Height in centimeters of players on the 003 Women s Worldd Cup soccer team. 157 1611 163 163 164 165 165 165 168 168 168 170 170 170 171 173 173 175 180 180 Determine the 5 number
More informationChapter 15 Mixed Models. Chapter Table of Contents. Introduction Split Plot Experiment Clustered Data References...
Chapter 15 Mixed Models Chapter Table of Contents Introduction...309 Split Plot Experiment...311 Clustered Data...320 References...326 308 Chapter 15. Mixed Models Chapter 15 Mixed Models Introduction
More informationRegression on the trees data with R
> trees Girth Height Volume 1 8.3 70 10.3 2 8.6 65 10.3 3 8.8 63 10.2 4 10.5 72 16.4 5 10.7 81 18.8 6 10.8 83 19.7 7 11.0 66 15.6 8 11.0 75 18.2 9 11.1 80 22.6 10 11.2 75 19.9 11 11.3 79 24.2 12 11.4 76
More informationBeta-Regression with SPSS Michael Smithson School of Psychology, The Australian National University
9/1/2005 Beta-Regression with SPSS 1 Beta-Regression with SPSS Michael Smithson School of Psychology, The Australian National University (email: Michael.Smithson@anu.edu.au) SPSS Nonlinear Regression syntax
More informationSTATS PAD USER MANUAL
STATS PAD USER MANUAL For Version 2.0 Manual Version 2.0 1 Table of Contents Basic Navigation! 3 Settings! 7 Entering Data! 7 Sharing Data! 8 Managing Files! 10 Running Tests! 11 Interpreting Output! 11
More information5 Bootstrapping. 4.7 Extensions. 4.8 References. 5.1 Bootstrapping for random samples (the i.i.d. case) ST697F: Topics in Regression.
ST697F: Topics in Regression. Spring 2007 c 21 4.7 Extensions The approach above can be readily extended to the case where we are interested in inverse prediction or regulation about one predictor given
More informationSolution to Bonus Questions
Solution to Bonus Questions Q2: (a) The histogram of 1000 sample means and sample variances are plotted below. Both histogram are symmetrically centered around the true lambda value 20. But the sample
More informationy ij= µ + τ i + b j + ε ij, i=1,2,...,t, j=1,2,...,b where b j ~ iid N(0, σ b 2 ) and ε ij ~ iid N(0, σ ε 2 ).
Mixed Models Analysis Using JMP Software 4. George A. Milliken, Kansas State University, Manhattan, KS Annie L. Dudley and John P. Sall, SAS Institute Inc. Cary, NC ABSTRACT The analysis of most designed
More informationChapter 10: Extensions to the GLM
Chapter 10: Extensions to the GLM 10.1 Implement a GAM for the Swedish mortality data, for males, using smooth functions for age and year. Age and year are standardized as described in Section 4.11, for
More informationOrganizing data in R. Fitting Mixed-Effects Models Using the lme4 Package in R. R packages. Accessing documentation. The Dyestuff data set
Fitting Mixed-Effects Models Using the lme4 Package in R Deepayan Sarkar Fred Hutchinson Cancer Research Center 18 September 2008 Organizing data in R Standard rectangular data sets (columns are variables,
More informationLaboratory for Two-Way ANOVA: Interactions
Laboratory for Two-Way ANOVA: Interactions For the last lab, we focused on the basics of the Two-Way ANOVA. That is, you learned how to compute a Brown-Forsythe analysis for a Two-Way ANOVA, as well as
More informationPaper SDA-11. Logistic regression will be used for estimation of net error for the 2010 Census as outlined in Griffin (2005).
Paper SDA-11 Developing a Model for Person Estimation in Puerto Rico for the 2010 Census Coverage Measurement Program Colt S. Viehdorfer, U.S. Census Bureau, Washington, DC This report is released to inform
More informationCompare Linear Regression Lines for the HP-67
Compare Linear Regression Lines for the HP-67 by Namir Shammas This article presents an HP-67 program that calculates the linear regression statistics for two data sets and then compares their slopes and
More informationSTAT 5200 Handout #24: Power Calculation in Mixed Models
STAT 5200 Handout #24: Power Calculation in Mixed Models Statistical power is the probability of finding an effect (i.e., calling a model term significant), given that the effect is real. ( Effect here
More informationWeek 10: Heteroskedasticity II
Week 10: Heteroskedasticity II Marcelo Coca Perraillon University of Colorado Anschutz Medical Campus Health Services Research Methods I HSMP 7607 2017 c 2017 PERRAILLON ARR 1 Outline Dealing with heteroskedasticy
More informationSoci Statistics for Sociologists
University of North Carolina Chapel Hill Soci708-001 Statistics for Sociologists Fall 2009 Professor François Nielsen Stata Commands for Module 7 Inference for Distributions For further information on
More informationQuality Control-2 The ASTA team
Quality Control-2 The ASTA team Contents 0.1 OC curves............................................... 1 0.2 OC curves - example......................................... 2 0.3 CUSUM chart.............................................
More informationSTAT 503 Fall Introduction to SAS
Getting Started Introduction to SAS 1) Download all of the files, sas programs (.sas) and data files (.dat) into one of your directories. I would suggest using your H: drive if you are using a computer
More informationExample Using Missing Data 1
Ronald H. Heck and Lynn N. Tabata 1 Example Using Missing Data 1 Creating the Missing Data Variable (Miss) Here is a data set (achieve subset MANOVAmiss.sav) with the actual missing data on the outcomes.
More informationStata versions 12 & 13 Week 4 Practice Problems
Stata versions 12 & 13 Week 4 Practice Problems SOLUTIONS 1 Practice Screen Capture a Create a word document Name it using the convention lastname_lab1docx (eg bigelow_lab1docx) b Using your browser, go
More informationStandard Errors in OLS Luke Sonnet
Standard Errors in OLS Luke Sonnet Contents Variance-Covariance of ˆβ 1 Standard Estimation (Spherical Errors) 2 Robust Estimation (Heteroskedasticity Constistent Errors) 4 Cluster Robust Estimation 7
More informationSAS/STAT 13.1 User s Guide. The NESTED Procedure
SAS/STAT 13.1 User s Guide The NESTED Procedure This document is an individual chapter from SAS/STAT 13.1 User s Guide. The correct bibliographic citation for the complete manual is as follows: SAS Institute
More informationSimulating Multivariate Normal Data
Simulating Multivariate Normal Data You have a population correlation matrix and wish to simulate a set of data randomly sampled from a population with that structure. I shall present here code and examples
More informationMultivariate Normal Random Numbers
Multivariate Normal Random Numbers Revised: 10/11/2017 Summary... 1 Data Input... 3 Analysis Options... 4 Analysis Summary... 5 Matrix Plot... 6 Save Results... 8 Calculations... 9 Summary This procedure
More informationFor Additional Information...
For Additional Information... The materials in this handbook were developed by Master Black Belts at General Electric Medical Systems to assist Black Belts and Green Belts in completing Minitab Analyses.
More informationE-Campus Inferential Statistics - Part 2
E-Campus Inferential Statistics - Part 2 Group Members: James Jones Question 4-Isthere a significant difference in the mean prices of the stores? New Textbook Prices New Price Descriptives 95% Confidence
More informationNCSS Statistical Software
Chapter 327 Geometric Regression Introduction Geometric regression is a special case of negative binomial regression in which the dispersion parameter is set to one. It is similar to regular multiple regression
More informationThe SURVEYREG Procedure
SAS/STAT 9.2 User s Guide The SURVEYREG Procedure (Book Excerpt) SAS Documentation This document is an individual chapter from SAS/STAT 9.2 User s Guide. The correct bibliographic citation for the complete
More informationIntroduction to Hierarchical Linear Model. Hsueh-Sheng Wu CFDR Workshop Series January 30, 2017
Introduction to Hierarchical Linear Model Hsueh-Sheng Wu CFDR Workshop Series January 30, 2017 1 Outline What is Hierarchical Linear Model? Why do nested data create analytic problems? Graphic presentation
More informationPSY 9556B (Feb 5) Latent Growth Modeling
PSY 9556B (Feb 5) Latent Growth Modeling Fixed and random word confusion Simplest LGM knowing how to calculate dfs How many time points needed? Power, sample size Nonlinear growth quadratic Nonlinear growth
More informationUnit 8 SUPPLEMENT Normal, T, Chi Square, F, and Sums of Normals
BIOSTATS 540 Fall 017 8. SUPPLEMENT Normal, T, Chi Square, F and Sums of Normals Page 1 of Unit 8 SUPPLEMENT Normal, T, Chi Square, F, and Sums of Normals Topic 1. Normal Distribution.. a. Definition..
More informationStat 5303 (Oehlert): Unbalanced Factorial Examples 1
Stat 5303 (Oehlert): Unbalanced Factorial Examples 1 > section
More informationThe NESTED Procedure (Chapter)
SAS/STAT 9.3 User s Guide The NESTED Procedure (Chapter) SAS Documentation This document is an individual chapter from SAS/STAT 9.3 User s Guide. The correct bibliographic citation for the complete manual
More informationBIOMETRICS INFORMATION
BIOMETRICS INFORMATION (You re 95% likely to need this information) PAMPHLET NO. # 57 DATE: September 5, 1997 SUBJECT: Interpreting Main Effects when a Two-way Interaction is Present Interpreting the analysis
More information