ESTIMATE PROP. IMPAIRED PRE- AND POST-INTERVENTION FOR THIN LIQUID SWALLOW TASKS. The SURVEYFREQ Procedure
|
|
|
- Scott Adams
- 5 years ago
- Views:
Transcription
1 ESTIMATE PROP. IMPAIRED PRE- AND POST-INTERVENTION FOR THIN LIQUID SWALLOW TASKS 18:58 Sunday, July 26, The SURVEYFREQ Procedure Data Summary Number of Clusters 30 Number of Observations 360 time_cat pas_ind Frequency Percent pre not impaired Std Err of Percent Table of time_cat by pas_ind 95% Limits for Percent Row Percent Std Err of Row Percent 95% Limits for Row Percent impaired Total post not impaired impaired Total Total not impaired impaired Total Frequency Missing = 5 Rao-Scott Chi-Square Test Pearson Chi-Square Design Correction Rao-Scott Chi-Square DF 1 Pr > ChiSq <.0001 F Value Num DF 1 Den DF 29 Pr > F <.0001 Sample Size = 355
2 ESTIMATE PROP. IMPAIRED PRE- AND POST-INTERVENTION AT EACH THIN LIQUID VOLUME 18:58 Sunday, July 26, The SURVEYFREQ Procedure Data Summary Number of Clusters 30 Number of Observations 360 pas_ind volume Frequency Percent Std Err of Percent Table of pas_ind by volume Controlling for time_cat=pre 95% Limits for Percent Row Percent Std Err of Row Percent 95% Limits for Row Percent not impaired 5ml ml cup sip Total impaired 5ml ml cup sip Total Total 5ml ml cup sip Total Rao-Scott Chi-Square Test Pearson Chi-Square Design Correction Rao-Scott Chi-Square DF 2 Pr > ChiSq <.0001 F Value Num DF 2 Den DF 58 Pr > F <.0001 Sample Size = 355
3 ESTIMATE PROP. IMPAIRED PRE- AND POST-INTERVENTION AT EACH THIN LIQUID VOLUME 18:58 Sunday, July 26, pas_ind volume Frequency Percent The SURVEYFREQ Procedure Std Err of Percent Table of pas_ind by volume Controlling for time_cat=post 95% Limits for Percent Row Percent Std Err of Row Percent 95% Limits for Row Percent not impaired 5ml ml cup sip Total impaired 5ml ml cup sip Total Total 5ml ml cup sip Total Rao-Scott Chi-Square Test Pearson Chi-Square Design Correction Rao-Scott Chi-Square DF 2 Pr > ChiSq <.0001 F Value Num DF 2 Den DF 58 Pr > F <.0001 Sample Size = 355
4 GEE ANALYSIS FOR THIN LIQUID TASKS 18:58 Sunday, July 26, Model Information Data Set WORK.THINLIQUID Distribution Binomial Link Function Logit Dependent Variable pas_ind Penetration-aspiration impairment indicator Number of Observations Read 360 Number of Observations Used 355 Number of Events 216 Number of Trials 355 Missing Values 5 Class Class Level Information Value time_cat post 1 pre 0 Design Variables volume 15ml 1 0 5ml 0 0 cup sip 0 1 inst MUSC 0 VA 1 Ordered Value Response Profile pas_ind Total Frequency 1 impaired not impaired 139 PROC GENMOD is modeling the probability that pas_ind='impaired'. One way to change this to model the probability that pas_ind='not impaired' is to specify the DESCENDING option in the PROC statement. Parameter Information Parameter Effect time_cat volume inst Prm1 Intercept Prm2 time_cat post Prm3 volume 15ml Prm4 volume cup sip Prm5 inst VA Algorithm converged.
5 GEE ANALYSIS FOR THIN LIQUID TASKS 18:58 Sunday, July 26, d Covariance Matrix Prm1 Prm2 Prm3 Prm4 Prm5 Prm Prm Prm Prm Prm GEE Model Information Correlation Structure Exchangeable Subject Effect subid (30 levels) Number of Clusters 30 Clusters With Missing Values 5 Correlation Matrix Dimension 12 Maximum Cluster Size 12 Minimum Cluster Size 11 Covariance Matrix (Model-Based) Prm1 Prm2 Prm3 Prm4 Prm5 Prm Prm Prm Prm Prm Algorithm converged. Working Correlation Matrix Col1 Col2 Col3 Col4 Col5 Col6 Col7 Col8 Col9 Col10 Col11 Col12 Row Row Row Row Row Row Row Row Row Row
6 GEE ANALYSIS FOR THIN LIQUID TASKS 18:58 Sunday, July 26, Working Correlation Matrix Col1 Col2 Col3 Col4 Col5 Col6 Col7 Col8 Col9 Col10 Col11 Col12 Row Row Exchangeable Working Correlation Correlation GEE Fit Criteria QIC QICu Parameter Analysis Of GEE Parameter s Empirical Standard Error s Standard Error 95% Limits Z Pr > Z Intercept time_cat post <.0001 volume 15ml <.0001 volume cup sip <.0001 inst VA Score Statistics For Type 3 GEE Analysis Source DF Chi-Square Pr > ChiSq time_cat <.0001 volume <.0001 inst Label Limits Contrast Results Standard Error Alpha Limits Chi-Square Pr > ChiSq treatment effect <.0001 Exp(treatment effect) ml vs 5ml <.0001 Exp(15ml vs 5ml) cupsip vs 5ml <.0001 Exp(cupsip vs 5ml)
7 GEE ANALYSIS FOR THIN LIQUID TASKS 18:58 Sunday, July 26, Label Limits Contrast Results Standard Error Alpha Limits Chi-Square Pr > ChiSq cupsip vs 15ml Exp(cupsip vs 15ml)
8 GEE ANALYSIS FOR FOR ALL SWALLOW TASKS 18:58 Sunday, July 26, Model Information Data Set WORK.SWALLOW Distribution Binomial Link Function Logit Dependent Variable pas_ind Penetration-aspiration impairment indicator Number of Observations Read 1080 Number of Observations Used 1064 Number of Events 618 Number of Trials 1064 Missing Values 16 Class Class Level Information Value time_cat post 1 pre 0 Design Variables viscosity honey 0 0 nectar 1 0 thin liquid 0 1 volume 15ml 1 0 5ml 0 0 cup sip 0 1 inst MUSC 0 VA 1 Ordered Value Response Profile pas_ind Total Frequency 1 impaired not impaired 446 PROC GENMOD is modeling the probability that pas_ind='impaired'. One way to change this to model the probability that pas_ind='not impaired' is to specify the DESCENDING option in the PROC statement. Parameter Information Parameter Effect time_cat viscosity volume inst Prm1 Intercept Prm2 time_cat post Prm3 viscosity nectar Prm4 viscosity thin liquid
9 GEE ANALYSIS FOR FOR ALL SWALLOW TASKS 18:58 Sunday, July 26, Parameter Information Parameter Effect time_cat viscosity volume inst Prm5 volume 15ml Prm6 volume cup sip Prm7 viscosity*volume nectar 15ml Prm8 viscosity*volume nectar cup sip Prm9 viscosity*volume thin liquid 15ml Prm10 viscosity*volume thin liquid cup sip Prm11 inst VA Algorithm converged. d Covariance Matrix Prm1 Prm2 Prm3 Prm4 Prm5 Prm6 Prm7 Prm8 Prm Prm Prm Prm Prm Prm Prm Prm Prm Prm Prm d Covariance Matrix Prm9 Prm10 Prm11 Prm Prm Prm Prm Prm Prm Prm Prm Prm E-6 Prm E-6 Prm E E
10 GEE ANALYSIS FOR FOR ALL SWALLOW TASKS 18:58 Sunday, July 26, GEE Model Information Correlation Structure Exchangeable Subject Effect subid (30 levels) Number of Clusters 30 Clusters With Missing Values 8 Correlation Matrix Dimension 36 Maximum Cluster Size 36 Minimum Cluster Size 29 Covariance Matrix (Model-Based) Prm1 Prm2 Prm3 Prm4 Prm5 Prm6 Prm7 Prm Prm Prm Prm Prm Prm Prm Prm Prm Prm Prm Covariance Matrix (Model-Based) Prm8 Prm9 Prm10 Prm11 Prm Prm Prm Prm Prm Prm Prm Prm Prm Prm Prm Algorithm converged.
11 GEE ANALYSIS FOR FOR ALL SWALLOW TASKS 18:58 Sunday, July 26, Exchangeable Working Correlation Correlation GEE Fit Criteria QIC QICu Parameter Analysis Of GEE Parameter s Empirical Standard Error s Standard Error 95% Limits Z Pr > Z Intercept time_cat post <.0001 viscosity nectar viscosity thin liquid volume 15ml volume cup sip viscosity*volume nectar 15ml viscosity*volume nectar cup sip viscosity*volume thin liquid 15ml <.0001 viscosity*volume thin liquid cup sip <.0001 inst VA Score Statistics For Joint Tests For GEE Source DF Chi-Square Pr > ChiSq time_cat <.0001 viscosity volume viscosity*volume inst Note: Under full-rank parameterizations, Type 3 effect tests are replaced by joint tests. The joint test for an effect is a test that all the parameters associated with that effect are zero. Such joint tests might not be equivalent to Type 3 effect tests under GLM parameterization.
12 GEE ANALYSIS FOR FOR ALL SWALLOW TASKS 18:58 Sunday, July 26, Label Limits Contrast Results Standard Error Alpha Limits Chi-Square Pr > ChiSq treatment effect <.0001 Exp(treatment effect) nectar: 15ml vs 5ml Exp(nectar: 15ml vs 5ml) honey: cupsip vs 15ml Exp(honey: cupsip vs 15ml)
13 18:58 Sunday, July 26, GEE ANALYSIS FOR FOR THIN LIQUID TASKS CLUSTERED ORDINAL RESPONSE PROPORTIONAL ODDS MODEL Model Information Data Set WORK.THINLIQUID Distribution Multinomial Link Function Cumulative Logit Dependent Variable pas_class Number of Observations Read 360 Number of Observations Used 355 Missing Values 5 Class Class Level Information Value time_cat post 1 pre 0 Design Variables volume 15ml 1 0 5ml 0 0 cup sip 0 1 inst MUSC 0 VA 1 Ordered Value Response Profile pas_class Total Frequency PROC GENMOD is modeling the probabilities of levels of pas_class having LOWER Ordered Values in the response profile table. One way to change this to model the probabilities of HIGHER Ordered Values is to specify the DESCENDING option in the PROC statement. Parameter Information Parameter Effect time_cat volume inst Prm1 time_cat post Prm2 volume 15ml Prm3 volume cup sip Prm4 inst VA Algorithm converged.
14 18:58 Sunday, July 26, GEE ANALYSIS FOR FOR THIN LIQUID TASKS CLUSTERED ORDINAL RESPONSE PROPORTIONAL ODDS MODEL d Covariance Matrix Intercept1 Intercept2 Prm1 Prm2 Prm3 Prm4 Intercept Intercept Prm Prm Prm Prm GEE Model Information Correlation Structure Independent Subject Effect subid (30 levels) Number of Clusters 30 Clusters With Missing Values 5 Correlation Matrix Dimension 12 Maximum Cluster Size 12 Minimum Cluster Size 11 Algorithm converged. GEE Fit Criteria QIC QICu Parameter Analysis Of GEE Parameter s Empirical Standard Error s Standard Error 95% Limits Z Pr > Z Intercept Intercept <.0001 time_cat post volume 15ml <.0001 volume cup sip <.0001 inst VA
15 18:58 Sunday, July 26, GEE ANALYSIS FOR FOR THIN LIQUID TASKS CLUSTERED ORDINAL RESPONSE PROPORTIONAL ODDS MODEL Score Statistics For Type 3 GEE Analysis Source DF Chi-Square Pr > ChiSq time_cat volume <.0001 inst Label Limits Contrast Results Standard Error Alpha Limits Chi-Square Pr > ChiSq treatment effect Exp(treatment effect) ml vs 5ml <.0001 Exp(15ml vs 5ml) cupsip vs 5ml <.0001 Exp(cupsip vs 5ml) cupsip vs 15ml Exp(cupsip vs 15ml)
Q30b Moyale Observed counts. The FREQ Procedure. Table 1 of type by response. Controlling for site=moyale. Improved (1+2) Same (3) Group only
Moyale Observed counts 12:28 Thursday, December 01, 2011 1 The FREQ Procedure Table 1 of by Controlling for site=moyale Row Pct Improved (1+2) Same () Worsened (4+5) Group only 16 51.61 1.2 14 45.16 1
You can specify the response in the form of a single variable or in the form of a ratio of two variables denoted events/trials.
The GENMOD Procedure MODEL Statement MODEL response = < effects > < /options > ; MODEL events/trials = < effects > < /options > ; You can specify the response in the form of a single variable or in the
Contrasting Marginal and Mixed Effects Models Recall: two approaches to handling dependence in Generalized Linear Models:
Contrasting Marginal and Mixed Effects Models Recall: two approaches to handling dependence in Generalized Linear Models: Marginal models: based on the consequences of dependence on estimating model parameters.
Overdispersion Workshop in generalized linear models Uppsala, June 11-12, Outline. Overdispersion
Biostokastikum Overdispersion is not uncommon in practice. In fact, some would maintain that overdispersion is the norm in practice and nominal dispersion the exception McCullagh and Nelder (1989) Overdispersion
SAS Analysis Examples Replication C8. * SAS Analysis Examples Replication for ASDA 2nd Edition * Berglund April 2017 * Chapter 8 ;
SAS Analysis Examples Replication C8 * SAS Analysis Examples Replication for ASDA 2nd Edition * Berglund April 2017 * Chapter 8 ; libname ncsr "P:\ASDA 2\Data sets\ncsr\" ; data c8_ncsr ; set ncsr.ncsr_sub_13nov2015
UNIVERSITY OF TORONTO. Faculty of Arts and Science APRIL 2010 EXAMINATIONS STA 303 H1S / STA 1002 HS. Duration - 3 hours. Aids Allowed: Calculator
UNIVERSITY OF TORONTO Faculty of Arts and Science APRIL 2010 EXAMINATIONS STA 303 H1S / STA 1002 HS Duration - 3 hours Aids Allowed: Calculator LAST NAME: FIRST NAME: STUDENT NUMBER: There are 27 pages
The GENMOD Procedure. Overview. Getting Started. Syntax. Details. Examples. References. SAS/STAT User's Guide. Book Contents Previous Next
Book Contents Previous Next SAS/STAT User's Guide Overview Getting Started Syntax Details Examples References Book Contents Previous Next Top http://v8doc.sas.com/sashtml/stat/chap29/index.htm29/10/2004
Poisson Data. Handout #4
Poisson Data The other response variable of interest records the number of blue spots observed after incubation. This type of data, i.e. count data, is often skewed showing numerous small values with occasional
Models for Binary Outcomes
Models for Binary Outcomes Introduction The simple or binary response (for example, success or failure) analysis models the relationship between a binary response variable and one or more explanatory variables.
COMPLEMENTARY LOG-LOG MODEL
COMPLEMENTARY LOG-LOG MODEL Under the assumption of binary response, there are two alternatives to logit model: probit model and complementary-log-log model. They all follow the same form π ( x) =Φ ( α
Model Based Statistics in Biology. Part V. The Generalized Linear Model. Chapter 18.1 Logistic Regression (Dose - Response)
Model Based Statistics in Biology. Part V. The Generalized Linear Model. Logistic Regression ( - Response) ReCap. Part I (Chapters 1,2,3,4), Part II (Ch 5, 6, 7) ReCap Part III (Ch 9, 10, 11), Part IV
CHL 5225 H Crossover Trials. CHL 5225 H Crossover Trials
CHL 55 H Crossover Trials The Two-sequence, Two-Treatment, Two-period Crossover Trial Definition A trial in which patients are randomly allocated to one of two sequences of treatments (either 1 then, or
Logistic Regression. Interpretation of linear regression. Other types of outcomes. 0-1 response variable: Wound infection. Usual linear regression
Logistic Regression Usual linear regression (repetition) y i = b 0 + b 1 x 1i + b 2 x 2i + e i, e i N(0,σ 2 ) or: y i N(b 0 + b 1 x 1i + b 2 x 2i,σ 2 ) Example (DGA, p. 336): E(PEmax) = 47.355 + 1.024
Subject-specific observed profiles of log(fev1) vs age First 50 subjects in Six Cities Study
Subject-specific observed profiles of log(fev1) vs age First 50 subjects in Six Cities Study 1.4 0.0-6 7 8 9 10 11 12 13 14 15 16 17 18 19 age Model 1: A simple broken stick model with knot at 14 fit with
Longitudinal Modeling with Logistic Regression
Newsom 1 Longitudinal Modeling with Logistic Regression Longitudinal designs involve repeated measurements of the same individuals over time There are two general classes of analyses that correspond to
Simple logistic regression
Simple logistic regression Biometry 755 Spring 2009 Simple logistic regression p. 1/47 Model assumptions 1. The observed data are independent realizations of a binary response variable Y that follows a
ST3241 Categorical Data Analysis I Multicategory Logit Models. Logit Models For Nominal Responses
ST3241 Categorical Data Analysis I Multicategory Logit Models Logit Models For Nominal Responses 1 Models For Nominal Responses Y is nominal with J categories. Let {π 1,, π J } denote the response probabilities
Appendix: Computer Programs for Logistic Regression
Appendix: Computer Programs for Logistic Regression In this appendix, we provide examples of computer programs to carry out unconditional logistic regression, conditional logistic regression, polytomous
Review. Timothy Hanson. Department of Statistics, University of South Carolina. Stat 770: Categorical Data Analysis
Review Timothy Hanson Department of Statistics, University of South Carolina Stat 770: Categorical Data Analysis 1 / 22 Chapter 1: background Nominal, ordinal, interval data. Distributions: Poisson, binomial,
Example 7b: Generalized Models for Ordinal Longitudinal Data using SAS GLIMMIX, STATA MEOLOGIT, and MPLUS (last proportional odds model only)
CLDP945 Example 7b page 1 Example 7b: Generalized Models for Ordinal Longitudinal Data using SAS GLIMMIX, STATA MEOLOGIT, and MPLUS (last proportional odds model only) This example comes from real data
Outline. Linear OLS Models vs: Linear Marginal Models Linear Conditional Models. Random Intercepts Random Intercepts & Slopes
Lecture 2.1 Basic Linear LDA 1 Outline Linear OLS Models vs: Linear Marginal Models Linear Conditional Models Random Intercepts Random Intercepts & Slopes Cond l & Marginal Connections Empirical Bayes
Testing Independence
Testing Independence Dipankar Bandyopadhyay Department of Biostatistics, Virginia Commonwealth University BIOS 625: Categorical Data & GLM 1/50 Testing Independence Previously, we looked at RR = OR = 1
Count data page 1. Count data. 1. Estimating, testing proportions
Count data page 1 Count data 1. Estimating, testing proportions 100 seeds, 45 germinate. We estimate probability p that a plant will germinate to be 0.45 for this population. Is a 50% germination rate
Stat 642, Lecture notes for 04/12/05 96
Stat 642, Lecture notes for 04/12/05 96 Hosmer-Lemeshow Statistic The Hosmer-Lemeshow Statistic is another measure of lack of fit. Hosmer and Lemeshow recommend partitioning the observations into 10 equal
STAT 5200 Handout #23. Repeated Measures Example (Ch. 16)
Motivating Example: Glucose STAT 500 Handout #3 Repeated Measures Example (Ch. 16) An experiment is conducted to evaluate the effects of three diets on the serum glucose levels of human subjects. Twelve
Using PROC GENMOD to Analyse Ratio to Placebo in Change of Dactylitis. Irmgard Hollweck / Meike Best 13.OCT.2013
Using PROC GENMOD to Analyse Ratio to Placebo in Change of Dactylitis Irmgard Hollweck / Meike Best 13.OCT.2013 Agenda 2 Introduction to Dactylitis Background Definitions: Trial Definitions:Terms Statistics:
Multinomial Logistic Regression Models
Stat 544, Lecture 19 1 Multinomial Logistic Regression Models Polytomous responses. Logistic regression can be extended to handle responses that are polytomous, i.e. taking r>2 categories. (Note: The word
The GENMOD Procedure (Book Excerpt)
SAS/STAT 9.22 User s Guide The GENMOD Procedure (Book Excerpt) SAS Documentation This document is an individual chapter from SAS/STAT 9.22 User s Guide. The correct bibliographic citation for the complete
Using PROC GENMOD to Analyse Ratio to Placebo in Change of Dactylitis
Paper SP03 Using PROC GENMOD to Analyse Ratio to Placebo in Change of Dactylitis Irmgard Hollweck, UCB Biosciences GmbH, Monheim, Germany Meike Best, UCB Biosciences GmbH, Monheim, Germany ABSTRACT A common
STA 303 H1S / 1002 HS Winter 2011 Test March 7, ab 1cde 2abcde 2fghij 3
STA 303 H1S / 1002 HS Winter 2011 Test March 7, 2011 LAST NAME: FIRST NAME: STUDENT NUMBER: ENROLLED IN: (circle one) STA 303 STA 1002 INSTRUCTIONS: Time: 90 minutes Aids allowed: calculator. Some formulae
Repeated Measures Part 2: Cartoon data
Repeated Measures Part 2: Cartoon data /*********************** cartoonglm.sas ******************/ options linesize=79 noovp formdlim='_'; title 'Cartoon Data: STA442/1008 F 2005'; proc format; /* value
Inference for Binomial Parameters
Inference for Binomial Parameters Dipankar Bandyopadhyay, Ph.D. Department of Biostatistics, Virginia Commonwealth University D. Bandyopadhyay (VCU) BIOS 625: Categorical Data & GLM 1 / 58 Inference for
Models for binary data
Faculty of Health Sciences Models for binary data Analysis of repeated measurements 2015 Julie Lyng Forman & Lene Theil Skovgaard Department of Biostatistics, University of Copenhagen 1 / 63 Program for
3 Joint Distributions 71
2.2.3 The Normal Distribution 54 2.2.4 The Beta Density 58 2.3 Functions of a Random Variable 58 2.4 Concluding Remarks 64 2.5 Problems 64 3 Joint Distributions 71 3.1 Introduction 71 3.2 Discrete Random
SAS/STAT 14.2 User s Guide. The GENMOD Procedure
SAS/STAT 14.2 User s Guide The GENMOD 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.
Analysis of Longitudinal Data: Comparison Between PROC GLM and PROC MIXED. Maribeth Johnson Medical College of Georgia Augusta, GA
Analysis of Longitudinal Data: Comparison Between PROC GLM and PROC MIXED Maribeth Johnson Medical College of Georgia Augusta, GA Overview Introduction to longitudinal data Describe the data for examples
Chapter 5: Logistic Regression-I
: Logistic Regression-I Dipankar Bandyopadhyay Department of Biostatistics, Virginia Commonwealth University BIOS 625: Categorical Data & GLM [Acknowledgements to Tim Hanson and Haitao Chu] D. Bandyopadhyay
IP WEIGHTING AND MARGINAL STRUCTURAL MODELS (CHAPTER 12) BIOS IPW and MSM
IP WEIGHTING AND MARGINAL STRUCTURAL MODELS (CHAPTER 12) BIOS 776 1 12 IPW and MSM IP weighting and marginal structural models ( 12) Outline 12.1 The causal question 12.2 Estimating IP weights via modeling
EPSY 905: Fundamentals of Multivariate Modeling Online Lecture #7
Introduction to Generalized Univariate Models: Models for Binary Outcomes EPSY 905: Fundamentals of Multivariate Modeling Online Lecture #7 EPSY 905: Intro to Generalized In This Lecture A short review
ST3241 Categorical Data Analysis I Two-way Contingency Tables. 2 2 Tables, Relative Risks and Odds Ratios
ST3241 Categorical Data Analysis I Two-way Contingency Tables 2 2 Tables, Relative Risks and Odds Ratios 1 What Is A Contingency Table (p.16) Suppose X and Y are two categorical variables X has I categories
Model Estimation Example
Ronald H. Heck 1 EDEP 606: Multivariate Methods (S2013) April 7, 2013 Model Estimation Example As we have moved through the course this semester, we have encountered the concept of model estimation. Discussions
BIOS 625 Fall 2015 Homework Set 3 Solutions
BIOS 65 Fall 015 Homework Set 3 Solutions 1. Agresti.0 Table.1 is from an early study on the death penalty in Florida. Analyze these data and show that Simpson s Paradox occurs. Death Penalty Victim's
Generalized Models: Part 1
Generalized Models: Part 1 Topics: Introduction to generalized models Introduction to maximum likelihood estimation Models for binary outcomes Models for proportion outcomes Models for categorical outcomes
STAT 7030: Categorical Data Analysis
STAT 7030: Categorical Data Analysis 5. Logistic Regression Peng Zeng Department of Mathematics and Statistics Auburn University Fall 2012 Peng Zeng (Auburn University) STAT 7030 Lecture Notes Fall 2012
Sections 4.1, 4.2, 4.3
Sections 4.1, 4.2, 4.3 Timothy Hanson Department of Statistics, University of South Carolina Stat 770: Categorical Data Analysis 1/ 32 Chapter 4: Introduction to Generalized Linear Models Generalized linear
Class Notes: Week 8. Probit versus Logit Link Functions and Count Data
Ronald Heck Class Notes: Week 8 1 Class Notes: Week 8 Probit versus Logit Link Functions and Count Data This week we ll take up a couple of issues. The first is working with a probit link function. While
UNIVERSITY OF MASSACHUSETTS Department of Mathematics and Statistics Applied Statistics Friday, January 15, 2016
UNIVERSITY OF MASSACHUSETTS Department of Mathematics and Statistics Applied Statistics Friday, January 15, 2016 Work all problems. 60 points are needed to pass at the Masters Level and 75 to pass at the
Correlated data. Non-normal outcomes. Reminder on binary data. Non-normal data. Faculty of Health Sciences. Non-normal outcomes
Faculty of Health Sciences Non-normal outcomes Correlated data Non-normal outcomes Lene Theil Skovgaard December 5, 2014 Generalized linear models Generalized linear mixed models Population average models
Section 9c. Propensity scores. Controlling for bias & confounding in observational studies
Section 9c Propensity scores Controlling for bias & confounding in observational studies 1 Logistic regression and propensity scores Consider comparing an outcome in two treatment groups: A vs B. In a
Cohen s s Kappa and Log-linear Models
Cohen s s Kappa and Log-linear Models HRP 261 03/03/03 10-11 11 am 1. Cohen s Kappa Actual agreement = sum of the proportions found on the diagonals. π ii Cohen: Compare the actual agreement with the chance
Lecture 3.1 Basic Logistic LDA
y Lecture.1 Basic Logistic LDA 0.2.4.6.8 1 Outline Quick Refresher on Ordinary Logistic Regression and Stata Women s employment example Cross-Over Trial LDA Example -100-50 0 50 100 -- Longitudinal Data
Introduction to Generalized Models
Introduction to Generalized Models Today s topics: The big picture of generalized models Review of maximum likelihood estimation Models for binary outcomes Models for proportion outcomes Models for categorical
Generalized linear models
Generalized linear models Douglas Bates November 01, 2010 Contents 1 Definition 1 2 Links 2 3 Estimating parameters 5 4 Example 6 5 Model building 8 6 Conclusions 8 7 Summary 9 1 Generalized Linear Models
Generalized Linear Model under the Extended Negative Multinomial Model and Cancer Incidence
Generalized Linear Model under the Extended Negative Multinomial Model and Cancer Incidence Sunil Kumar Dhar Center for Applied Mathematics and Statistics, Department of Mathematical Sciences, New Jersey
ANALYSING BINARY DATA IN A REPEATED MEASUREMENTS SETTING USING SAS
Libraries 1997-9th Annual Conference Proceedings ANALYSING BINARY DATA IN A REPEATED MEASUREMENTS SETTING USING SAS Eleanor F. Allan Follow this and additional works at: http://newprairiepress.org/agstatconference
T-test: means of Spock's judge versus all other judges 1 12:10 Wednesday, January 5, judge1 N Mean Std Dev Std Err Minimum Maximum
T-test: means of Spock's judge versus all other judges 1 The TTEST Procedure Variable: pcwomen judge1 N Mean Std Dev Std Err Minimum Maximum OTHER 37 29.4919 7.4308 1.2216 16.5000 48.9000 SPOCKS 9 14.6222
STA6938-Logistic Regression Model
Dr. Ying Zhang STA6938-Logistic Regression Model Topic 2-Multiple Logistic Regression Model Outlines:. Model Fitting 2. Statistical Inference for Multiple Logistic Regression Model 3. Interpretation of
SAS Analysis Examples Replication C11
SAS Analysis Examples Replication C11 * SAS Analysis Examples Replication for ASDA 2nd Edition * Berglund April 2017 * Chapter 11 ; libname d "P:\ASDA 2\Data sets\hrs 2012\HRS 2006_2012 Longitudinal File\"
wa.lst wa.lst Printed by Alison Gibbs Oct 14, 12 12:53 Page 2/14 ALL DATA 2 12:53 Sunday, October 14, 2012 The FREQ Procedure
Oct 14, 12 12:53 Page 1/14 ALL DATA 1 The MEANS Procedure Variable N Mean Std Dev Minimum subjage 60 59.6833333 6.2504011 45.0000000 aftermean 60 51.7800000 7.1625864 37.2000000 beforemean 60 57.9566667
dm'log;clear;output;clear'; options ps=512 ls=99 nocenter nodate nonumber nolabel FORMCHAR=" = -/\<>*"; ODS LISTING;
dm'log;clear;output;clear'; options ps=512 ls=99 nocenter nodate nonumber nolabel FORMCHAR=" ---- + ---+= -/\*"; ODS LISTING; *** Table 23.2 ********************************************; *** Moore, David
MIXED MODELS FOR REPEATED (LONGITUDINAL) DATA PART 2 DAVID C. HOWELL 4/1/2010
MIXED MODELS FOR REPEATED (LONGITUDINAL) DATA PART 2 DAVID C. HOWELL 4/1/2010 Part 1 of this document can be found at http://www.uvm.edu/~dhowell/methods/supplements/mixed Models for Repeated Measures1.pdf
Covariance Structure Approach to Within-Cases
Covariance Structure Approach to Within-Cases Remember how the data file grapefruit1.data looks: Store sales1 sales2 sales3 1 62.1 61.3 60.8 2 58.2 57.9 55.1 3 51.6 49.2 46.2 4 53.7 51.5 48.3 5 61.4 58.7
EDF 7405 Advanced Quantitative Methods in Educational Research. Data are available on IQ of the child and seven potential predictors.
EDF 7405 Advanced Quantitative Methods in Educational Research Data are available on IQ of the child and seven potential predictors. Four are medical variables available at the birth of the child: Birthweight
Topic 20: Single Factor Analysis of Variance
Topic 20: Single Factor Analysis of Variance Outline Single factor Analysis of Variance One set of treatments Cell means model Factor effects model Link to linear regression using indicator explanatory
Ch 6: Multicategory Logit Models
293 Ch 6: Multicategory Logit Models Y has J categories, J>2. Extensions of logistic regression for nominal and ordinal Y assume a multinomial distribution for Y. In R, we will fit these models using the
Semiparametric Generalized Linear Models
Semiparametric Generalized Linear Models North American Stata Users Group Meeting Chicago, Illinois Paul Rathouz Department of Health Studies University of Chicago prathouz@uchicago.edu Liping Gao MS Student
Changes Report 2: Examples from the Australian Longitudinal Study on Women s Health for Analysing Longitudinal Data
ChangesReport: ExamplesfromtheAustralianLongitudinal StudyonWomen shealthforanalysing LongitudinalData June005 AustralianLongitudinalStudyonWomen shealth ReporttotheDepartmentofHealthandAgeing ThisreportisbasedonthecollectiveworkoftheStatisticsGroupoftheAustralianLongitudinal
7/28/15. Review Homework. Overview. Lecture 6: Logistic Regression Analysis
Lecture 6: Logistic Regression Analysis Christopher S. Hollenbeak, PhD Jane R. Schubart, PhD The Outcomes Research Toolbox Review Homework 2 Overview Logistic regression model conceptually Logistic regression
Three Factor Completely Randomized Design with One Continuous Factor: Using SPSS GLM UNIVARIATE R. C. Gardner Department of Psychology
Data_Analysis.calm Three Factor Completely Randomized Design with One Continuous Factor: Using SPSS GLM UNIVARIATE R. C. Gardner Department of Psychology This article considers a three factor completely
Introduction to SAS proc mixed
Faculty of Health Sciences Introduction to SAS proc mixed Analysis of repeated measurements, 2017 Julie Forman Department of Biostatistics, University of Copenhagen 2 / 28 Preparing data for analysis The
R code and output of examples in text. Contents. De Jong and Heller GLMs for Insurance Data R code and output. 1 Poisson regression 2
R code and output of examples in text Contents 1 Poisson regression 2 2 Negative binomial regression 5 3 Quasi likelihood regression 6 4 Logistic regression 6 5 Ordinal regression 10 6 Nominal regression
Lecture 8: Summary Measures
Lecture 8: Summary Measures Dipankar Bandyopadhyay, Ph.D. BMTRY 711: Analysis of Categorical Data Spring 2011 Division of Biostatistics and Epidemiology Medical University of South Carolina Lecture 8:
MULTIVARIATE HOMEWORK #5
MULTIVARIATE HOMEWORK #5 Fisher s dataset on differentiating species of Iris based on measurements on four morphological characters (i.e. sepal length, sepal width, petal length, and petal width) was subjected
unadjusted model for baseline cholesterol 22:31 Monday, April 19,
unadjusted model for baseline cholesterol 22:31 Monday, April 19, 2004 1 Class Level Information Class Levels Values TRETGRP 3 3 4 5 SEX 2 0 1 Number of observations 916 unadjusted model for baseline cholesterol
Introduction to SAS proc mixed
Faculty of Health Sciences Introduction to SAS proc mixed Analysis of repeated measurements, 2017 Julie Forman Department of Biostatistics, University of Copenhagen Outline Data in wide and long format
4.1 Computing section Example: Bivariate measurements on plants Post hoc analysis... 7
Master of Applied Statistics ST116: Chemometrics and Multivariate Statistical data Analysis Per Bruun Brockhoff Module 4: Computing 4.1 Computing section.................................. 1 4.1.1 Example:
Model and Working Correlation Structure Selection in GEE Analyses of Longitudinal Data
The 3rd Australian and New Zealand Stata Users Group Meeting, Sydney, 5 November 2009 1 Model and Working Correlation Structure Selection in GEE Analyses of Longitudinal Data Dr Jisheng Cui Public Health
More Accurately Analyze Complex Relationships
SPSS Advanced Statistics 17.0 Specifications More Accurately Analyze Complex Relationships Make your analysis more accurate and reach more dependable conclusions with statistics designed to fit the inherent
Regression Models for Risks(Proportions) and Rates. Proportions. E.g. [Changes in] Sex Ratio: Canadian Births in last 60 years
![Regression Models for Risks(Proportions) and Rates. Proportions. E.g. [Changes in] Sex Ratio: Canadian Births in last 60 years](/thumbs/83/87443434.jpg "Regression Models for Risks(Proportions) and Rates. Proportions. E.g. [Changes in] Sex Ratio: Canadian Births in last 60 years")
Regression Models for Risks(Proportions) and Rates Proportions E.g. [Changes in] Sex Ratio: Canadian Births in last 60 years Parameter of Interest: (male)... just above 0.5 or (male)/ (female)... "ODDS"...
Where have all the puffins gone 1. An Analysis of Atlantic Puffin (Fratercula arctica) Banding Returns in Newfoundland
Where have all the puffins gone 1 An Analysis of Atlantic Puffin (Fratercula arctica) Banding Returns in Newfoundland Dave Fifield 1 With apologies to Pete Seeger. Introduction Seabirds are the most visible
,..., θ(2),..., θ(n)
Likelihoods for Multivariate Binary Data Log-Linear Model We have 2 n 1 distinct probabilities, but we wish to consider formulations that allow more parsimonious descriptions as a function of covariates.
ssh tap sas913, sas https://www.statlab.umd.edu/sasdoc/sashtml/onldoc.htm
Kedem, STAT 430 SAS Examples: Logistic Regression ==================================== ssh abc@glue.umd.edu, tap sas913, sas https://www.statlab.umd.edu/sasdoc/sashtml/onldoc.htm a. Logistic regression.
Answer to exercise: Blood pressure lowering drugs
Answer to exercise: Blood pressure lowering drugs The data set bloodpressure.txt contains data from a cross-over trial, involving three different formulations of a drug for lowering of blood pressure:
ANALYSES OF NCGS DATA FOR ALCOHOL STATUS CATEGORIES 1 22:46 Sunday, March 2, 2003
ANALYSES OF NCGS DATA FOR ALCOHOL STATUS CATEGORIES 1 22:46 Sunday, March 2, 2003 The MEANS Procedure DRINKING STATUS=1 Analysis Variable : TRIGL N Mean Std Dev Minimum Maximum 164 151.6219512 95.3801744
Review of Multinomial Distribution If n trials are performed: in each trial there are J > 2 possible outcomes (categories) Multicategory Logit Models
Chapter 6 Multicategory Logit Models Response Y has J > 2 categories. Extensions of logistic regression for nominal and ordinal Y assume a multinomial distribution for Y. 6.1 Logit Models for Nominal Responses
General Linear Model (Chapter 4)
General Linear Model (Chapter 4) Outcome variable is considered continuous Simple linear regression Scatterplots OLS is BLUE under basic assumptions MSE estimates residual variance testing regression coefficients
2. We care about proportion for categorical variable, but average for numerical one.
Probit Model 1. We apply Probit model to Bank data. The dependent variable is deny, a dummy variable equaling one if a mortgage application is denied, and equaling zero if accepted. The key regressor is
A course in statistical modelling. session 06b: Modelling count data
A Course in Statistical Modelling University of Glasgow 29 and 30 January, 2015 session 06b: Modelling count data Graeme Hutcheson 1 Luiz Moutinho 2 1 Manchester Institute of Education Manchester university
The SAS System 18:28 Saturday, March 10, Plot of Canonical Variables Identified by Cluster
The SAS System 18:28 Saturday, March 10, 2018 1 The FASTCLUS Procedure Replace=FULL Radius=0 Maxclusters=2 Maxiter=10 Converge=0.02 Initial Seeds Cluster SepalLength SepalWidth PetalLength PetalWidth 1
Chapter 14 Logistic and Poisson Regressions
STAT 525 SPRING 2018 Chapter 14 Logistic and Poisson Regressions Professor Min Zhang Logistic Regression Background In many situations, the response variable has only two possible outcomes Disease (Y =
Some comments on Partitioning
Some comments on Partitioning Dipankar Bandyopadhyay Department of Biostatistics, Virginia Commonwealth University BIOS 625: Categorical Data & GLM 1/30 Partitioning Chi-Squares We have developed tests
BMI 541/699 Lecture 22
BMI 541/699 Lecture 22 Where we are: 1. Introduction and Experimental Design 2. Exploratory Data Analysis 3. Probability 4. T-based methods for continous variables 5. Power and sample size for t-based
CHAPTER 1: BINARY LOGIT MODEL
CHAPTER 1: BINARY LOGIT MODEL Prof. Alan Wan 1 / 44 Table of contents 1. Introduction 1.1 Dichotomous dependent variables 1.2 Problems with OLS 3.3.1 SAS codes and basic outputs 3.3.2 Wald test for individual
ZERO INFLATED POISSON REGRESSION
STAT 6500 ZERO INFLATED POISSON REGRESSION FINAL PROJECT DEC 6 th, 2013 SUN JEON DEPARTMENT OF SOCIOLOGY UTAH STATE UNIVERSITY POISSON REGRESSION REVIEW INTRODUCING - ZERO-INFLATED POISSON REGRESSION SAS
SAS Software to Fit the Generalized Linear Model
SAS Software to Fit the Generalized Linear Model Gordon Johnston, SAS Institute Inc., Cary, NC Abstract In recent years, the class of generalized linear models has gained popularity as a statistical modeling
Chapter 1. Modeling Basics
Chapter 1. Modeling Basics What is a model? Model equation and probability distribution Types of model effects Writing models in matrix form Summary 1 What is a statistical model? A model is a mathematical
Regression models. Generalized linear models in R. Normal regression models are not always appropriate. Generalized linear models. Examples.
Regression models Generalized linear models in R Dr Peter K Dunn http://www.usq.edu.au Department of Mathematics and Computing University of Southern Queensland ASC, July 00 The usual linear regression
Descriptive Statistics
*following creates z scores for the ydacl statedp traitdp and rads vars. *specifically adding the /SAVE subcommand to descriptives will create z. *scores for whatever variables are in the command. DESCRIPTIVES
STAT 705: Analysis of Contingency Tables
STAT 705: Analysis of Contingency Tables Timothy Hanson Department of Statistics, University of South Carolina Stat 705: Analysis of Contingency Tables 1 / 45 Outline of Part I: models and parameters Basic
STAT 501 EXAM I NAME Spring 1999
STAT 501 EXAM I NAME Spring 1999 Instructions: You may use only your calculator and the attached tables and formula sheet. You can detach the tables and formula sheet from the rest of this exam. Show your