Logi Reg in-hosp mortality vs gender in POST- cabg MI patients using "where" option with RACE + income for non-s temi <50 T he LOGISTIC Procedure M od

Size: px
Start display at page:

Download "Logi Reg in-hosp mortality vs gender in POST- cabg MI patients using "where" option with RACE + income for non-s temi <50 T he LOGISTIC Procedure M od"

Transcription

1 Logi Reg in-hosp mortality vs gender in POST- cabg MI patients using "where" option with RACE + income for non-s temi <50 T he LOGISTIC Procedure M odel Information D ata Set L IBRARY.NISMICATHCABG4 R esponse Variable D IED D ied during hospitalization N umber of Response L evels 2 W eight Variable D ISCWT Weight to discharges in A HA universe M odel b inary logit Optimization T echnique F isher's scoring P robability O rdered V alue N umber of Observations Read 641 N umber of Observations Used 631 S um of Weights Read Sum of W eights Used R esponse Profile D IED T otal F requency T otal W eight 1 1 : Died in hospital : Did not die in hospital modeled is DIED='1: Died in hospital'. N ote: 10 observations were deleted due to missing values for the response or explanatory v ariables. Level Information V alue D esign Variables F EMALE 0 : Male 0 1 : Female 1 R ENLFAIL 0 0 1

2 Level Information V alue D esign Variables 1 1 D M D MCX H TN_C A IDS A LCOHOL A NEMDEF A RTH r ace Z IPINC_QRTL 1 : First quartile : Second quartile 0 3 : Third quartile 0 4 : Fourth quartile 0 H OSP_LOCATION 0 : Rural 1 1 : Urban 0 H _CONTRL 0 2 2

3 Level Information V alue D esign Variables 3 H OSP_TEACH 0 : Nonteaching 1 1 : Teaching 0 B LDLOSS C HF C HRNLUNG C OAG D EPRESS D RUG H YPOTHY L IVER L YMPH L YTES M ETS 0 N EURO O BESE 3

4 Level Information V alue D esign Variables P ARA P ERIVASC P SYCH P ULMCIRC T UMOR 0 U LCER 0 V ALVE 0 W GHTLOSS c ararrhythmia I teration M odel Convergence Status limit reached without convergence. W arning: C onvergence was not attained in 25 iterations. You may want to increase the maximum number of iterations (MAXITER= option) or change the convergence criteria (ABSFCONV=, FCONV=, GCONV=, XCONV= options) in the MODEL s tatement. W arning: The LOGISTIC procedure continues in spite of the above warning. Results s hown are based on the last maximum likelihood iteration. Validity of the model f it is questionable. M odel Fit Statistics C riterion I ntercept O nly I ntercept a nd C ovariates 4

5 M odel Fit Statistics C riterion I ntercept O nly I ntercept a nd C ovariates A IC S C Log L Testing Global Null Hypothe sis: BETA= 0 T est Chi-S quare D F L ikelihood Ratio S core W ald T ype E ffect D F 3 Analysis of Effects W ald Chi-S quare F EMALE D M D MCX H TN_C A IDS A LCOHOL A NEMDEF A RTH r ace Z IPINC_QRTL H OSP_LOCATION H _CONTRL H OSP_TEACH T OTAL_DISC B LDLOSS

6 T ype E ffect D F 3 Analysis of Effects W ald Chi-S quare C HF C HRNLUNG C OAG D EPRESS D RUG H YPOTHY L IVER L YMPH L YTES M ETS 0.. N EURO O BESE P ARA P ERIVASC P SYCH P ULMCIRC R ENLFAIL T UMOR 0.. U LCER 0.. V ALVE 0.. W GHTLOSS c ararrhythmia A nalysis P arameter D F of Maximum Likelihood Estimates S tandard E rror W ald Chi- S quare I ntercept

7 A nalysis P arameter D F of Maximum Likelihood Estimates S tandard E rror W ald Chi- S quare F EMALE 1 : Female D M D MCX H TN_C A IDS A LCOHOL A NEMDEF A RTH r ace r ace r ace r ace Z IPINC_QRTL 2: Second Z IPINC_QRTL 3 : Third Z IPINC_QRTL 4: Fourth H OSP_LOCATION 0 : Rural H _CONTRL H _CONTRL HO SP_TEACH 0: N onteaching T OTAL_DISC B LDLOSS C HF C HRNLUNG C OAG

8 A nalysis P arameter D F of Maximum Likelihood Estimates S tandard E rror W ald Chi- S quare D EPRESS D RUG H YPOTHY L IVER L YMPH L YTES N EURO O BESE P ARA P ERIVASC P SYCH P ULMCIRC RENLF AIL W GHTLOSS c ararrhythmia O dds Ratio Estimates E ffect Point 95% Wald Confidence L imits F EMALE 1: Female vs 0: Male D M 0 vs 1 D MCX 0 vs 1 H TN_C 0 vs A IDS 0 vs A LCOHOL 0 vs 1 A NEMDEF 0 vs 1 A RTH 0 vs 1 8

9 O dds Ratio Estimates E ffect Point 95% Wald Confidence L imits r ace1 2 vs 1 r ace1 3 vs 1 r ace1 4 vs 1 r ace1 5 vs 1 ZIPINC_QRTL 2: Second quartile vs 1: First Z IPINC_QRTL 3: Third quartile vs 1: First ZIPINC_QRTL 4: Fourth quartile vs 1: First H OSP_LOCATION 0: Rural vs 1: Urban H _CONTRL 2 vs 1 H _CONTRL 3 vs 1 < H OSP_TEACH 0: Nonteaching vs 1: Teaching T OTAL_DISC B LDLOSS 0 vs 1 C HF 0 vs 1 C HRNLUNG 0 vs 1 C OAG 0 vs 1 D EPRESS 0 vs 1 D RUG 0 vs H YPOTHY 0 vs 1 L IVER 0 vs 1 L YMPH 0 vs 1 L YTES 0 vs 1 N EURO 0 vs 1 O BESE 0 vs 1 P ARA 0 vs 1 9

10 O dds Ratio Estimates E ffect Point 95% Wald Confidence L imits P ERIVASC 0 vs 1 P SYCH 0 vs 1 P ULMCIRC 0 vs 1 R ENLFAIL 1 vs 0 WGHT LOSS 0 vs 1 c ararrhythmia 0 vs 1 A ssociation of Predicted Probabilities and O bserved Responses P ercent Concordant S omers' D P ercent Discordant 0. 0 G amma P ercent Tied 0. 0 Tau-a P airs 4984 c

11 Logi Reg in-hosp mortality vs gender in POST- cabg MI patients using "where" option with RACE + income for non-s temi <75 T he LOGISTIC Procedure M odel Information D ata Set L IBRARY.NISMICATHCABG4 R esponse Variable D IED Died durin g hospitalization Number of Response L evels 2 W eight Variable D ISCWT Weight to discharges in A HA universe M odel b inary logit Optimization T echnique F isher's scoring P robability O rdered V alue N umber of Observations Read 4626 N umber of Observations Used 4529 Sum of Weig hts Read S um of Weights Used R esponse Profile D IED T otal F requency T otal W eight 1 1 : Died in hospital : Did not die in hospital modeled is DIED='1: Died in hospital'. N ote: 97 observations were deleted due to missing values for the response or explanatory v ariables. Level Information V alue D esign Variables F EMALE 0 : Male 0 11

12 Level Information V alue D esign Variables 1 : Female 1 R ENLFAIL D M DMC X H TN_C A IDS A LCOHOL A NEMDEF A RTH r ace ZIPINC_QR TL 1 : First quartile : Second quartile 0 3 : Third quartile 0 4 : Fourth quartile 0 H OSP_LOCATION 0 : Rural 1 1 : Urban 0 12

13 Level Information V alue D esign Variables H _CONTRL H OSP_TEACH 0 : Nonteaching 1 1: Teaching 0 B LDLOSS C HF C HRNLUNG C OAG D EPRESS D RUG H YPOTHY L IVER L YMPH L YTES M ETS N EURO 13

14 Level Information V alue D esign Variables O BESE P ARA P ERIVASC P SYCH PULMCIR C T UMOR U LCER V ALVE W GHTLOSS c ararrhythmia M odel Convergence Status Convergence criterion (GCONV=1E-8) satisfi ed. M odel Fit Statistics C riterion I ntercept O nly I ntercept a nd C ovariates 14

15 M odel Fit Statistics C riterion I ntercept O nly I ntercept a nd C ovariates A IC S C Log L T esting Global Null Hypothesis: BETA=0 T est Chi-S quare D F L ikelihood Ratio < S core W ald. 42. W arning: T he information matrix is singular and thus the convergence is questionable. T ype E ffect D F 3 Analysis of Effects W ald Chi-S quare F EMALE D M D MCX H TN_C A IDS A LCOHOL A NEMDEF A RTH r ace Z IPINC_QRTL H OSP_LOCATION H _CONTRL H OSP_TEACH

16 T ype E ffect D F 3 Analysis of Effects W ald Chi-S quare T OTAL_DISC B LDLOSS C HF C HRNLUNG C OAG D EPRESS D RUG H YPOTHY L IVER L YMPH LY TES M ETS N EURO O BESE P ARA P ERIVASC P SYCH P ULMCIRC R ENLFAIL T UMOR U LCER 0.. VAL VE W GHTLOSS c ararrhythmia A nalysis of Maximum Likelihood Estimates 16

17 P arameter D F S tandard E rror W ald Chi- S quare I ntercept F EMALE 1 : Female D M D MCX H TN_C A IDS A LCOHOL A NEMDEF A RTH r ace r ace r ace r ace Z IPINC_QRTL 2: Second Z IPINC_QRTL 3: Third Z IPINC_QRTL 4: Fourth H OSP_LOCATION 0 : Rural H _CONTRL H _CONTRL H OSP_TEACH 0: N onteaching T OTAL_DISC E E B LDLOSS C HF C HRNLUNG C OAG

18 A nalysis P arameter D F of Maximum Likelihood Estimates S tandard E rror W ald Chi- S quare D EPRESS D RUG H YPOTHY L IVER L YMPH L YTES M ETS N EURO O BESE P ARA P ERIVASC P SYCH P ULMCIRC R ENLFAIL T UMOR U LCER V ALVE W GHTLOSS c ararrhythmia O dds Ratio Estimates E ffect Point 95% Wald Confidence L imits F EMALE 1: Female vs 0: Male D M 0 vs D MCX 0 vs H TN_C 0 vs

19 Z IPINC_QRTL 2: O dds Ratio Estimates E ffect Point 95% Wald Confidence L imits A IDS 0 vs 1 < A LCOHOL 0 vs A NEMDEF 0 vs A RTH 0 vs r ace1 2 vs r ace1 3 vs r ace1 4 vs r ace1 5 vs Second quartile vs 1: First ZIPINC_QRTL 3: Third quartile vs 1: First ZIPINC_QRTL 4: Fourth quartile vs 1: First H OSP_LOCATION 0: Rural vs 1: Urban H_CONT RL 2 vs H _CONTRL 3 vs H OSP_TEACH 0: Nonteaching vs 1: Teaching T OTAL_DISC B LDLOSS 0 vs C HF 0 vs C HRNLUNG 0 vs COA G 0 vs D EPRESS 0 vs D RUG 0 vs H YPOTHY 0 vs L IVER 0 vs L YMPH 0 vs 1 19

20 O dds Ratio Estimates E ffect Point 95% Wald Confidence L imits L YTES 0 vs M ETS 0 vs N EURO 0 vs O BESE 0 vs P ARA 0 vs P ERIVASC 0 vs P SYCH 0 vs P ULMCIRC 0 vs 1 R ENLFAIL 1 vs TUMOR 0 vs V ALVE 0 vs W GHTLOSS 0 vs c ararrhythmia 0 vs A ssociation of Predicted Probabilities and O bserved Responses P ercent Concordant S omers' D P ercent Discordant G amma P ercent Tied 2. 2 Tau-a P airs c

Logi Reg interaction in- h osp mortality vs gender in overall patients using "where" option with RACE T he GLIMMIX Procedure Model Informati on D ata

Logi Reg interaction in- h osp mortality vs gender in overall patients using where option with RACE T he GLIMMIX Procedure Model Informati on D ata p roc g limmix d ata = library.nismicathcabg4 ; */ridging=absolute not needed*/ ; c lass female dm dmcx htn_c aids alcohol ANEMDEF arth race1 ZIPINC_QRTL hosp_location h_contrl h osp_teach b ldloss chf

More information

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, 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

More information

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 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

More information

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.

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

More information

Section 9c. Propensity scores. Controlling for bias & confounding in observational studies

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

More information

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 ; 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

More information

Count data page 1. Count data. 1. Estimating, testing proportions

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

More information

COMPLEMENTARY LOG-LOG MODEL

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) =Φ ( α

More information

CHAPTER 1: BINARY LOGIT MODEL

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

More information

ssh tap sas913, sas https://www.statlab.umd.edu/sasdoc/sashtml/onldoc.htm

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.

More information

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 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

More information

Simple logistic regression

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

More information

Class Notes: Week 8. Probit versus Logit Link Functions and Count Data

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

More information

Investigating Models with Two or Three Categories

Investigating Models with Two or Three Categories Ronald H. Heck and Lynn N. Tabata 1 Investigating Models with Two or Three Categories For the past few weeks we have been working with discriminant analysis. Let s now see what the same sort of model might

More information

Section IX. Introduction to Logistic Regression for binary outcomes. Poisson regression

Section IX. Introduction to Logistic Regression for binary outcomes. Poisson regression Section IX Introduction to Logistic Regression for binary outcomes Poisson regression 0 Sec 9 - Logistic regression In linear regression, we studied models where Y is a continuous variable. What about

More information

In Class Review Exercises Vartanian: SW 540

In Class Review Exercises Vartanian: SW 540 In Class Review Exercises Vartanian: SW 540 1. Given the following output from an OLS model looking at income, what is the slope and intercept for those who are black and those who are not black? b SE

More information

STA6938-Logistic Regression Model

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

More information

SAS Analysis Examples Replication C11

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\"

More information

Models for Binary Outcomes

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.

More information

BIOS 625 Fall 2015 Homework Set 3 Solutions

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

More information

Department of Mathematics The University of Toledo. Master of Science Degree Comprehensive Examination Applied Statistics.

Department of Mathematics The University of Toledo. Master of Science Degree Comprehensive Examination Applied Statistics. Department of Mathematics The University of Toledo Master of Science Degree Comprehensive Examination Applied Statistics April 8, 205 nstructions Do all problems. Show all of your computations. Prove all

More information

ADVANCED STATISTICAL ANALYSIS OF EPIDEMIOLOGICAL STUDIES. Cox s regression analysis Time dependent explanatory variables

ADVANCED STATISTICAL ANALYSIS OF EPIDEMIOLOGICAL STUDIES. Cox s regression analysis Time dependent explanatory variables ADVANCED STATISTICAL ANALYSIS OF EPIDEMIOLOGICAL STUDIES Cox s regression analysis Time dependent explanatory variables Henrik Ravn Bandim Health Project, Statens Serum Institut 4 November 2011 1 / 53

More information

Lecture 7 Time-dependent Covariates in Cox Regression

Lecture 7 Time-dependent Covariates in Cox Regression Lecture 7 Time-dependent Covariates in Cox Regression So far, we ve been considering the following Cox PH model: λ(t Z) = λ 0 (t) exp(β Z) = λ 0 (t) exp( β j Z j ) where β j is the parameter for the the

More information

Generalized logit models for nominal multinomial responses. Local odds ratios

Generalized logit models for nominal multinomial responses. Local odds ratios Generalized logit models for nominal multinomial responses Categorical Data Analysis, Summer 2015 1/17 Local odds ratios Y 1 2 3 4 1 π 11 π 12 π 13 π 14 π 1+ X 2 π 21 π 22 π 23 π 24 π 2+ 3 π 31 π 32 π

More information

PROC LOGISTIC: Traps for the unwary Peter L. Flom, Independent statistical consultant, New York, NY

PROC LOGISTIC: Traps for the unwary Peter L. Flom, Independent statistical consultant, New York, NY Paper SD174 PROC LOGISTIC: Traps for the unwary Peter L. Flom, Independent statistical consultant, New York, NY ABSTRACT Keywords: Logistic. INTRODUCTION This paper covers some gotchas in SAS R PROC LOGISTIC.

More information

Binary Logistic Regression

Binary Logistic Regression The coefficients of the multiple regression model are estimated using sample data with k independent variables Estimated (or predicted) value of Y Estimated intercept Estimated slope coefficients Ŷ = b

More information

Introduction to mtm: An R Package for Marginalized Transition Models

Introduction to mtm: An R Package for Marginalized Transition Models Introduction to mtm: An R Package for Marginalized Transition Models Bryan A. Comstock and Patrick J. Heagerty Department of Biostatistics University of Washington 1 Introduction Marginalized transition

More information

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: 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.

More information

STAT 7030: Categorical Data Analysis

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

More information

4/6/16. Non-parametric Test. Overview. Stephen Opiyo. Distinguish Parametric and Nonparametric Test Procedures

4/6/16. Non-parametric Test. Overview. Stephen Opiyo. Distinguish Parametric and Nonparametric Test Procedures Non-parametric Test Stephen Opiyo Overview Distinguish Parametric and Nonparametric Test Procedures Explain commonly used Nonparametric Test Procedures Perform Hypothesis Tests Using Nonparametric Procedures

More information

ECLT 5810 Linear Regression and Logistic Regression for Classification. Prof. Wai Lam

ECLT 5810 Linear Regression and Logistic Regression for Classification. Prof. Wai Lam ECLT 5810 Linear Regression and Logistic Regression for Classification Prof. Wai Lam Linear Regression Models Least Squares Input vectors is an attribute / feature / predictor (independent variable) The

More information

Logistic Regression Analysis

Logistic Regression Analysis Logistic Regression Analysis Predicting whether an event will or will not occur, as well as identifying the variables useful in making the prediction, is important in most academic disciplines as well

More information

Lecture 8: Summary Measures

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:

More information

Ron Heck, Fall Week 3: Notes Building a Two-Level Model

Ron Heck, Fall Week 3: Notes Building a Two-Level Model Ron Heck, Fall 2011 1 EDEP 768E: Seminar on Multilevel Modeling rev. 9/6/2011@11:27pm Week 3: Notes Building a Two-Level Model We will build a model to explain student math achievement using student-level

More information

ECLT 5810 Linear Regression and Logistic Regression for Classification. Prof. Wai Lam

ECLT 5810 Linear Regression and Logistic Regression for Classification. Prof. Wai Lam ECLT 5810 Linear Regression and Logistic Regression for Classification Prof. Wai Lam Linear Regression Models Least Squares Input vectors is an attribute / feature / predictor (independent variable) The

More information

(Where does Ch. 7 on comparing 2 means or 2 proportions fit into this?)

(Where does Ch. 7 on comparing 2 means or 2 proportions fit into this?) 12. Comparing Groups: Analysis of Variance (ANOVA) Methods Response y Explanatory x var s Method Categorical Categorical Contingency tables (Ch. 8) (chi-squared, etc.) Quantitative Quantitative Regression

More information

ESP 178 Applied Research Methods. 2/23: Quantitative Analysis

ESP 178 Applied Research Methods. 2/23: Quantitative Analysis ESP 178 Applied Research Methods 2/23: Quantitative Analysis Data Preparation Data coding create codebook that defines each variable, its response scale, how it was coded Data entry for mail surveys and

More information

Linear Regression With Special Variables

Linear Regression With Special Variables Linear Regression With Special Variables Junhui Qian December 21, 2014 Outline Standardized Scores Quadratic Terms Interaction Terms Binary Explanatory Variables Binary Choice Models Standardized Scores:

More information

Introduction to lnmle: An R Package for Marginally Specified Logistic-Normal Models for Longitudinal Binary Data

Introduction to lnmle: An R Package for Marginally Specified Logistic-Normal Models for Longitudinal Binary Data Introduction to lnmle: An R Package for Marginally Specified Logistic-Normal Models for Longitudinal Binary Data Bryan A. Comstock and Patrick J. Heagerty Department of Biostatistics University of Washington

More information

Advanced Quantitative Data Analysis

Advanced Quantitative Data Analysis Chapter 24 Advanced Quantitative Data Analysis Daniel Muijs Doing Regression Analysis in SPSS When we want to do regression analysis in SPSS, we have to go through the following steps: 1 As usual, we choose

More information

Binary Dependent Variables

Binary Dependent Variables Binary Dependent Variables In some cases the outcome of interest rather than one of the right hand side variables - is discrete rather than continuous Binary Dependent Variables In some cases the outcome

More information

Exercise Sheet 6: Solutions

Exercise Sheet 6: Solutions Exercise Sheet 6: Solutions R.G. Pierse 1. (a) Regression yields: Dependent Variable: LC Date: 10/29/02 Time: 18:37 Sample(adjusted): 1950 1985 Included observations: 36 after adjusting endpoints C 0.244716

More information

Nonparametric statistic methods. Waraphon Phimpraphai DVM, PhD Department of Veterinary Public Health

Nonparametric statistic methods. Waraphon Phimpraphai DVM, PhD Department of Veterinary Public Health Nonparametric statistic methods Waraphon Phimpraphai DVM, PhD Department of Veterinary Public Health Measurement What are the 4 levels of measurement discussed? 1. Nominal or Classificatory Scale Gender,

More information

Adjusting Analyses of Survey Results using a Predicted Probability of Response Rob Gately, INNOVUS: Epidemiology, Waltham, MA

Adjusting Analyses of Survey Results using a Predicted Probability of Response Rob Gately, INNOVUS: Epidemiology, Waltham, MA PharmaSUG2011 Paper HS07 Adjusting Analyses of Survey Results using a Predicted Probability of Response Rob Gately, INNOVUS: Epidemiology, Waltham, MA ABSTRACT When conducting surveys, an issue that is

More information

Two sample Hypothesis tests in R.

Two sample Hypothesis tests in R. Example. (Dependent samples) Two sample Hypothesis tests in R. A Calculus professor gives their students a 10 question algebra pretest on the first day of class, and a similar test towards the end of the

More information

Models for Count and Binary Data. Poisson and Logistic GWR Models. 24/07/2008 GWR Workshop 1

Models for Count and Binary Data. Poisson and Logistic GWR Models. 24/07/2008 GWR Workshop 1 Models for Count and Binary Data Poisson and Logistic GWR Models 24/07/2008 GWR Workshop 1 Outline I: Modelling counts Poisson regression II: Modelling binary events Logistic Regression III: Poisson Regression

More information

2.1 Linear regression with matrices

2.1 Linear regression with matrices 21 Linear regression with matrices The values of the independent variables are united into the matrix X (design matrix), the values of the outcome and the coefficient are represented by the vectors Y and

More information

Multinomial Logistic Regression Models

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

More information

Cohen s s Kappa and Log-linear Models

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

More information

LOGISTIC REGRESSION. Lalmohan Bhar Indian Agricultural Statistics Research Institute, New Delhi

LOGISTIC REGRESSION. Lalmohan Bhar Indian Agricultural Statistics Research Institute, New Delhi LOGISTIC REGRESSION Lalmohan Bhar Indian Agricultural Statistics Research Institute, New Delhi- lmbhar@gmail.com. Introduction Regression analysis is a method for investigating functional relationships

More information

STAT 501 EXAM I NAME Spring 1999

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

More information

Ron Heck, Fall Week 8: Introducing Generalized Linear Models: Logistic Regression 1 (Replaces prior revision dated October 20, 2011)

Ron Heck, Fall Week 8: Introducing Generalized Linear Models: Logistic Regression 1 (Replaces prior revision dated October 20, 2011) Ron Heck, Fall 2011 1 EDEP 768E: Seminar in Multilevel Modeling rev. January 3, 2012 (see footnote) Week 8: Introducing Generalized Linear Models: Logistic Regression 1 (Replaces prior revision dated October

More information

Basic Medical Statistics Course

Basic Medical Statistics Course Basic Medical Statistics Course S7 Logistic Regression November 2015 Wilma Heemsbergen w.heemsbergen@nki.nl Logistic Regression The concept of a relationship between the distribution of a dependent variable

More information

Q30b Moyale Observed counts. The FREQ Procedure. Table 1 of type by response. Controlling for site=moyale. Improved (1+2) Same (3) Group only

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

More information

Booklet of Code and Output for STAD29/STA 1007 Midterm Exam

Booklet of Code and Output for STAD29/STA 1007 Midterm Exam Booklet of Code and Output for STAD29/STA 1007 Midterm Exam List of Figures in this document by page: List of Figures 1 Packages................................ 2 2 Hospital infection risk data (some).................

More information

Linear Regression. y» F; Ey = + x Vary = ¾ 2. ) y = + x + u. Eu = 0 Varu = ¾ 2 Exu = 0:

Linear Regression. y» F; Ey = + x Vary = ¾ 2. ) y = + x + u. Eu = 0 Varu = ¾ 2 Exu = 0: Linear Regression 1 Single Explanatory Variable Assume (y is not necessarily normal) where Examples: y» F; Ey = + x Vary = ¾ 2 ) y = + x + u Eu = 0 Varu = ¾ 2 Exu = 0: 1. School performance as a function

More information

Module 4 Practice problem and Homework answers

Module 4 Practice problem and Homework answers Module 4 Practice problem and Homework answers Practice Problems, pages 2 3 OBSERVED cat no cat total over 3 7 10 under 9 3 12 total 12 10 22 The observed table simply summarizes the data from page 1.

More information

Lecture 14: Introduction to Poisson Regression

Lecture 14: Introduction to Poisson Regression Lecture 14: Introduction to Poisson Regression Ani Manichaikul amanicha@jhsph.edu 8 May 2007 1 / 52 Overview Modelling counts Contingency tables Poisson regression models 2 / 52 Modelling counts I Why

More information

Modelling counts. Lecture 14: Introduction to Poisson Regression. Overview

Modelling counts. Lecture 14: Introduction to Poisson Regression. Overview Modelling counts I Lecture 14: Introduction to Poisson Regression Ani Manichaikul amanicha@jhsph.edu Why count data? Number of traffic accidents per day Mortality counts in a given neighborhood, per week

More information

Answers to Problem Set #4

Answers to Problem Set #4 Answers to Problem Set #4 Problems. Suppose that, from a sample of 63 observations, the least squares estimates and the corresponding estimated variance covariance matrix are given by: bβ bβ 2 bβ 3 = 2

More information

STA6938-Logistic Regression Model

STA6938-Logistic Regression Model Dr. Ying Zhang STA6938-Logistic Regression Model Topic 6-Logistic Regression for Case-Control Studies Outlines: 1. Biomedical Designs 2. Logistic Regression Models for Case-Control Studies 3. Logistic

More information

Relative-risk regression and model diagnostics. 16 November, 2015

Relative-risk regression and model diagnostics. 16 November, 2015 Relative-risk regression and model diagnostics 16 November, 2015 Relative risk regression More general multiplicative intensity model: Intensity for individual i at time t is i(t) =Y i (t)r(x i, ; t) 0

More information

Logistic Regression. Interpretation of linear regression. Other types of outcomes. 0-1 response variable: Wound infection. Usual linear regression

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

More information

Regression so far... Lecture 21 - Logistic Regression. Odds. Recap of what you should know how to do... At this point we have covered: Sta102 / BME102

Regression so far... Lecture 21 - Logistic Regression. Odds. Recap of what you should know how to do... At this point we have covered: Sta102 / BME102 Background Regression so far... Lecture 21 - Sta102 / BME102 Colin Rundel November 18, 2014 At this point we have covered: Simple linear regression Relationship between numerical response and a numerical

More information

Sociology 362 Data Exercise 6 Logistic Regression 2

Sociology 362 Data Exercise 6 Logistic Regression 2 Sociology 362 Data Exercise 6 Logistic Regression 2 The questions below refer to the data and output beginning on the next page. Although the raw data are given there, you do not have to do any Stata runs

More information

Biostat Methods STAT 5820/6910 Handout #5a: Misc. Issues in Logistic Regression

Biostat Methods STAT 5820/6910 Handout #5a: Misc. Issues in Logistic Regression Biostat Methods STAT 5820/6910 Handout #5a: Misc. Issues in Logistic Regression Recall general χ 2 test setu: Y 0 1 Trt 0 a b Trt 1 c d I. Basic logistic regression Previously (Handout 4a): χ 2 test of

More information

BIOSTATS Intermediate Biostatistics Spring 2017 Exam 2 (Units 3, 4 & 5) Practice Problems SOLUTIONS

BIOSTATS Intermediate Biostatistics Spring 2017 Exam 2 (Units 3, 4 & 5) Practice Problems SOLUTIONS BIOSTATS 640 - Intermediate Biostatistics Spring 2017 Exam 2 (Units 3, 4 & 5) Practice Problems SOLUTIONS Practice Question 1 Both the Binomial and Poisson distributions have been used to model the quantal

More information

GROUPED DATA E.G. FOR SAMPLE OF RAW DATA (E.G. 4, 12, 7, 5, MEAN G x / n STANDARD DEVIATION MEDIAN AND QUARTILES STANDARD DEVIATION

GROUPED DATA E.G. FOR SAMPLE OF RAW DATA (E.G. 4, 12, 7, 5, MEAN G x / n STANDARD DEVIATION MEDIAN AND QUARTILES STANDARD DEVIATION FOR SAMPLE OF RAW DATA (E.G. 4, 1, 7, 5, 11, 6, 9, 7, 11, 5, 4, 7) BE ABLE TO COMPUTE MEAN G / STANDARD DEVIATION MEDIAN AND QUARTILES Σ ( Σ) / 1 GROUPED DATA E.G. AGE FREQ. 0-9 53 10-19 4...... 80-89

More information

36-309/749 Experimental Design for Behavioral and Social Sciences. Dec 1, 2015 Lecture 11: Mixed Models (HLMs)

36-309/749 Experimental Design for Behavioral and Social Sciences. Dec 1, 2015 Lecture 11: Mixed Models (HLMs) 36-309/749 Experimental Design for Behavioral and Social Sciences Dec 1, 2015 Lecture 11: Mixed Models (HLMs) Independent Errors Assumption An error is the deviation of an individual observed outcome (DV)

More information

Chapter 19: Logistic regression

Chapter 19: Logistic regression Chapter 19: Logistic regression Self-test answers SELF-TEST Rerun this analysis using a stepwise method (Forward: LR) entry method of analysis. The main analysis To open the main Logistic Regression dialog

More information

Sociology Exam 1 Answer Key Revised February 26, 2007

Sociology Exam 1 Answer Key Revised February 26, 2007 Sociology 63993 Exam 1 Answer Key Revised February 26, 2007 I. True-False. (20 points) Indicate whether the following statements are true or false. If false, briefly explain why. 1. An outlier on Y will

More information

MANOVA MANOVA,$/,,# ANOVA ##$%'*!# 1. $!;' *$,$!;' (''

MANOVA MANOVA,$/,,# ANOVA ##$%'*!# 1. $!;' *$,$!;' ('' 14 3! "#!$%# $# $&'('$)!! (Analysis of Variance : ANOVA) *& & "#!# +, ANOVA -& $ $ (+,$ ''$) *$#'$)!!#! (Multivariate Analysis of Variance : MANOVA).*& ANOVA *+,'$)$/*! $#/#-, $(,!0'%1)!', #($!#$ # *&,

More information

Didacticiel - Études de cas. In this tutorial, we show how to use TANAGRA ( and higher) for measuring the association between ordinal variables.

Didacticiel - Études de cas. In this tutorial, we show how to use TANAGRA ( and higher) for measuring the association between ordinal variables. Subject Association measures for ordinal variables. In this tutorial, we show how to use TANAGRA (1.4.19 and higher) for measuring the association between ordinal variables. All the measures that we present

More information

MMWS Software Program Manual

MMWS Software Program Manual MMWS Software Program Manual 1 Software Development The MMWS program is regularly updated. The latest beta version can be downloaded from http://hlmsoft.net/ghong/ MMWS Click here to get MMWS. For a sample

More information

Sections 4.1, 4.2, 4.3

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

More information

Biostat Methods STAT 5500/6500 Handout #12: Methods and Issues in (Binary Response) Logistic Regression

Biostat Methods STAT 5500/6500 Handout #12: Methods and Issues in (Binary Response) Logistic Regression Biostat Methods STAT 5500/6500 Handout #12: Methods and Issues in (Binary Resonse) Logistic Regression Recall general χ 2 test setu: Y 0 1 Trt 0 a b Trt 1 c d I. Basic logistic regression Previously (Handout

More information

Logistic Regression Analyses in the Water Level Study

Logistic Regression Analyses in the Water Level Study Logistic Regression Analyses in the Water Level Study A. Introduction. 166 students participated in the Water level Study. 70 passed and 96 failed to correctly draw the water level in the glass. There

More information

β j = coefficient of x j in the model; β = ( β1, β2,

β j = coefficient of x j in the model; β = ( β1, β2, Regression Modeling of Survival Time Data Why regression models? Groups similar except for the treatment under study use the nonparametric methods discussed earlier. Groups differ in variables (covariates)

More information

8 Nominal and Ordinal Logistic Regression

8 Nominal and Ordinal Logistic Regression 8 Nominal and Ordinal Logistic Regression 8.1 Introduction If the response variable is categorical, with more then two categories, then there are two options for generalized linear models. One relies on

More information

STA 4504/5503 Sample Exam 1 Spring 2011 Categorical Data Analysis. 1. Indicate whether each of the following is true (T) or false (F).

STA 4504/5503 Sample Exam 1 Spring 2011 Categorical Data Analysis. 1. Indicate whether each of the following is true (T) or false (F). STA 4504/5503 Sample Exam 1 Spring 2011 Categorical Data Analysis 1. Indicate whether each of the following is true (T) or false (F). (a) (b) (c) (d) (e) In 2 2 tables, statistical independence is equivalent

More information

Immigration attitudes (opposes immigration or supports it) it may seriously misestimate the magnitude of the effects of IVs

Immigration attitudes (opposes immigration or supports it) it may seriously misestimate the magnitude of the effects of IVs Logistic Regression, Part I: Problems with the Linear Probability Model (LPM) Richard Williams, University of Notre Dame, https://www3.nd.edu/~rwilliam/ Last revised February 22, 2015 This handout steals

More information

Introducing Generalized Linear Models: Logistic Regression

Introducing Generalized Linear Models: Logistic Regression Ron Heck, Summer 2012 Seminars 1 Multilevel Regression Models and Their Applications Seminar Introducing Generalized Linear Models: Logistic Regression The generalized linear model (GLM) represents and

More information

Jun Tu. Department of Geography and Anthropology Kennesaw State University

Jun Tu. Department of Geography and Anthropology Kennesaw State University Examining Spatially Varying Relationships between Preterm Births and Ambient Air Pollution in Georgia using Geographically Weighted Logistic Regression Jun Tu Department of Geography and Anthropology Kennesaw

More information

Claudia Spix. Medizinisches Institut fur Umwelthygiene an der Universitat Dusseldorf

Claudia Spix. Medizinisches Institut fur Umwelthygiene an der Universitat Dusseldorf AN ITERATIVE WEIGHTING PROCEDURE FOR REGRESSION OF AREA MORTALITY RATES ON EXPLANATORY VARIABLES AFTER AN IDEA BY S.J.POCOCK ET AL. 1981. PROGR"AMMED WITH STANDARD SAS PROCEDURES WITHOUT MACROS Problem

More information

Statistics 135 Fall 2008 Final Exam

Statistics 135 Fall 2008 Final Exam Name: SID: Statistics 135 Fall 2008 Final Exam Show your work. The number of points each question is worth is shown at the beginning of the question. There are 10 problems. 1. [2] The normal equations

More information

Logistic Regression. Some slides from Craig Burkett. STA303/STA1002: Methods of Data Analysis II, Summer 2016 Michael Guerzhoy

Logistic Regression. Some slides from Craig Burkett. STA303/STA1002: Methods of Data Analysis II, Summer 2016 Michael Guerzhoy Logistic Regression Some slides from Craig Burkett STA303/STA1002: Methods of Data Analysis II, Summer 2016 Michael Guerzhoy Titanic Survival Case Study The RMS Titanic A British passenger liner Collided

More information

Truck prices - linear model? Truck prices - log transform of the response variable. Interpreting models with log transformation

Truck prices - linear model? Truck prices - log transform of the response variable. Interpreting models with log transformation Background Regression so far... Lecture 23 - Sta 111 Colin Rundel June 17, 2014 At this point we have covered: Simple linear regression Relationship between numerical response and a numerical or categorical

More information

Study Guide #3: OneWay ANALYSIS OF VARIANCE (ANOVA)

Study Guide #3: OneWay ANALYSIS OF VARIANCE (ANOVA) Study Guide #3: OneWay ANALYSIS OF VARIANCE (ANOVA) About the ANOVA Test In educational research, we are most often involved finding out whether there are differences between groups. For example, is there

More information

Introduction to Regression Analysis. Dr. Devlina Chatterjee 11 th August, 2017

Introduction to Regression Analysis. Dr. Devlina Chatterjee 11 th August, 2017 Introduction to Regression Analysis Dr. Devlina Chatterjee 11 th August, 2017 What is regression analysis? Regression analysis is a statistical technique for studying linear relationships. One dependent

More information

Chapter 20: Logistic regression for binary response variables

Chapter 20: Logistic regression for binary response variables Chapter 20: Logistic regression for binary response variables In 1846, the Donner and Reed families left Illinois for California by covered wagon (87 people, 20 wagons). They attempted a new and untried

More information

emerge Network: CERC Survey Survey Sampling Data Preparation

emerge Network: CERC Survey Survey Sampling Data Preparation emerge Network: CERC Survey Survey Sampling Data Preparation Overview The entire patient population does not use inpatient and outpatient clinic services at the same rate, nor are racial and ethnic subpopulations

More information

More Accurately Analyze Complex Relationships

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

More information

Survival Regression Models

Survival Regression Models Survival Regression Models David M. Rocke May 18, 2017 David M. Rocke Survival Regression Models May 18, 2017 1 / 32 Background on the Proportional Hazards Model The exponential distribution has constant

More information

Addition and subtraction: element by element, and dimensions must match.

Addition and subtraction: element by element, and dimensions must match. Matrix Essentials review: ) Matrix: Rectangular array of numbers. ) ranspose: Rows become columns and vice-versa ) single row or column is called a row or column) Vector ) R ddition and subtraction: element

More information

Testing and Model Selection

Testing and Model Selection Testing and Model Selection This is another digression on general statistics: see PE App C.8.4. The EViews output for least squares, probit and logit includes some statistics relevant to testing hypotheses

More information

A Generalized Linear Model for Binomial Response Data. Copyright c 2017 Dan Nettleton (Iowa State University) Statistics / 46

A Generalized Linear Model for Binomial Response Data. Copyright c 2017 Dan Nettleton (Iowa State University) Statistics / 46 A Generalized Linear Model for Binomial Response Data Copyright c 2017 Dan Nettleton (Iowa State University) Statistics 510 1 / 46 Now suppose that instead of a Bernoulli response, we have a binomial response

More information

Applied Survival Analysis Lab 10: Analysis of multiple failures

Applied Survival Analysis Lab 10: Analysis of multiple failures Applied Survival Analysis Lab 10: Analysis of multiple failures We will analyze the bladder data set (Wei et al., 1989). A listing of the dataset is given below: list if id in 1/9 +---------------------------------------------------------+

More information

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 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

More information

Generalized linear models

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

More information

Description Syntax for predict Menu for predict Options for predict Remarks and examples Methods and formulas References Also see

Description Syntax for predict Menu for predict Options for predict Remarks and examples Methods and formulas References Also see Title stata.com logistic postestimation Postestimation tools for logistic Description Syntax for predict Menu for predict Options for predict Remarks and examples Methods and formulas References Also see

More information