Homework 10 - Solution
|
|
- Hollie Hawkins
- 5 years ago
- Views:
Transcription
1 STAT Spring 2011 Homework 10 - Solution Olga Vitek Each part of the problems 5 points 1. Faraway Ch. 4 problem 1 (page 93) : The dataset parstum contains cross-classified data on marijuana usage by collge students as it relates to the alcohol and drug usage of the parents. Analyze the data as if both factors were nominal. Redo the analysis treating both factors as ordinal. Contrast the results. Answer : (a) As nominal : Since the calculated test statistic is with degree of freedom=4, there is strong evidence that the fitted model is insufficient with p-value= glm(formula = count ~ parent + student, family = poisson, data = parstum) Deviance Residuals: Estimate Std. Error z value Pr(> z ) (Intercept) < 2e-16 *** parentneither < 2e-16 *** parentone e-14 *** studentoccasional e-10 *** studentregular e-10 *** --- Signif. codes: 0 *** ** 0.01 * (Dispersion parameter for poisson family taken to be 1) Null deviance: on 8 degrees of freedom Residual deviance: on 4 degrees of freedom AIC: (b) As ordinal : Parent( Neither=0, One=1, Both=2), Student(Never=0, Occasional=1, regular=2) Since the calculated test statistic is with degree of freedom=3, there is not enough evidence that the fitted model is not sufficient to describe the data with p-value = glm(formula = count ~ parent + student + I(oparent * ostudent), family = poisson, data = parstum) Deviance Residuals: Estimate Std. Error z value Pr(> z ) (Intercept) e-08 *** parentneither e-16 *** 1
2 parentone e-15 *** studentoccasional e-14 *** studentregular e-13 *** I(oparent * ostudent) e-06 *** --- Signif. codes: 0 *** ** 0.01 * (Dispersion parameter for poisson family taken to be 1) Null deviance: on 8 degrees of freedom Residual deviance: on 3 degrees of freedom AIC: Number of Fisher Scoring iterations: 4 (c) Contrast the results A comparison of the two fitted models indicate that the model with ordinal factors is significantly better than the one with the nominal factors with p-value for the interaction is e-06 R code Analysis of Deviance Table Model 1: count ~ parent + student Model 2: count ~ parent + student + I(oparent * ostudent) Resid. Df Resid. Dev Df Deviance ## nominal nominal<-glm(count~parent+student, family=poisson, parstum) summary(nominal) pchisq(nominal$deviance, nominal$df.residual, lower=false) ## ordinal parstum$oparent<-unclass(parstum$parent) parstum$ostudent<-unclass(parstum$student) parstum$oparent[parstum$parent == "Neither"]<-0 parstum$oparent[parstum$parent == "One"]<-1 parstum$oparent[parstum$parent == "Both"]<-2 parstum ordinal<-glm(count~parent+student+i(oparent*ostudent), family=poisson, parstum) summary(ordinal2) pchisq(ordinal$deviance, ordinal$df.residual, lower=false) ### contrast anova(nominal, ordinal) pchisq(20.809, 1, lower=false) 2
3 2. [Methods qualifying exam, August 2010: use paper and pencil.] You conduct a survey of students in the university according to their book reading level and notes reading level, in order to study similarities and differences of these two types of study. You acquire the following data. Book Reading Total Notes Reading Very Good Good Bad Very Bad Very Good Good Bad Vary Bad Total (a) It is possible to consider up to 5 loglinear models, which express various levels of dependence between these two studies: the independence model, the linear-by-linear model, the row-effect model, the column-effect model, and the saturated model. Display all these models in the mathematical notation, define model terms and the assumptions, and provide the associated residual degrees of freedom. Let Y ij denote the count in row i and column j. We assume that Y ij P oisson(λ ij ), or, conditional on the total, Y ij Multinomial(π ij ). Let α i and β j denote row and column effects respectively. Let x i, y j denote the scores of rows and columns of the independent variables, defined on a continuous scale. The loglinear models are Models Formula Residual Df Independence log(λ ij ) = µ + α i + β j 9 Linear-by-linear log(λ ij ) = µ + α i + β j + ω(x i y j ) 8 Row-effect log(λ ij ) = µ + α i + β j + ω i y j 6 Column-effect log(λ ij ) = µ + α i + β j + ω j x i 6 Saturated log(λ ij ) = µ + α i + β j + (αβ) ij 0 subject to constraints α i = i=1 β j = j=1 ω i = i=1 ω j = j=1 (αβ) ij = i=1 (αβ) ij = 0. j=1 (b) Consider two extreme cases of the models above: the independence and the saturated models. Test both models for goodness of fit using deviance test, at the confidence level of 95%. State the null and the alternative hypotheses, the test statistic, and the conclusion. For the independence model, we test H 0 : The independence model is appropriate against H a : The saturated model is appropriate. The deviance statistic is G 2 = 2 i=1 j=1 n ij log n ij ˆn ij, where n ij is the observed count and ˆn ij is the predicted count. Based on the independence model, 4 i=1 ˆn ij = n n ij j=1 n ij 4 4 i=1 j=1 n ij 3
4 and we have the following predicted table Book Reading Notes Reading Very Good Good Bad Very Bad Very Good Good Bad Vary Bad Therefore, G 2 = > χ 2 9(0.95) = Therefore we reject H 0. For the saturated model, the same model is specified by both H 0 and H a. Therefore G 2 = 0, and we cannot conduct the test. (c) The residual deviance of the linear-by-linear model is Test whether the linear-by-linear model is more appropriate for the data than the independence model. Use the confidence level of 95%. The deviance of the linear-by-linear term is G 2 = = > χ 2 1(0.95) = Therefore, the linear-by-linear model is more appropriate for the data. 3. [Adapted from methods qualifying exam, January 2008: use paper and pencil.] A sample of 1083 voters were surveyed for their political ideology and political party affiliation in a presidential primary in Wisconsin. See the table below. Political Ideology Party Affiliation Liberal Moderate Conservative Total Democrat Independent Republican Total (a) The following commands were executed in R to check the independence between the political ideology and party affiliation. > mod1 <- glm(y ~ paff + polid, family = poisson, data = pidata) > summary(mod1) Estimate Std. Error z value P(> z ) (Intercept) < 2e-16 *** paffind * paffrep e-13 *** polidmod e-09 *** polidcon *** Null deviance: on 8 degrees of freedom Residual deviance: on 4 degrees of freedom AIC:
5 Please write down the fitted model and report a test of the independence hypothesis (include the hypotheses). Answer : Let Y ij denote the count in the ith row and j column of the table, with i = 1, 2, 3 = {Dep, Ind, Rep}, j = 1, 2, 3 = {Liberal, Moderate, Conservative}. We assume that Y ij P oisson(λ ij ), or equivalently Y j i Multinomial(π 1 i, π 2 i, π 3 i ) The fitted model assumes the independence of rows and columns, i.e. log(λ ij ) = µ + α i + β j, or equvalently Y j i Multinomial(π 1, π 2, π 3 ) log ˆλ ij = I Ind I Rep I moderate I conservative Hypothesis : H 0 : The independence model log(λ ij ) = µ+α i +β j, or equivalently Y j i Multinomial(π 1, π 2, π 3 ) H a : The saturated model log(λ ij ) = µ+α i +β j +(αβ) ij, or equivalently Y j i Multinomial(π 1 i, π 2 i, π 3 i ) Test statistic = Residual deviance = G 2 = χ 2 4 critical value = χ 2 4,0.95 = p-value = 6.126e 22 0 Reject H 0. The independence assumption is rejected. (b) As both political ideology and party affiliation may be considered ordinal, the following models are explored. > ct1 <- xtabs(y~polid+paff,pidata); > pid1<-data.frame(dem=ct1[,1],ind=ct1[,2],rep=ct1[,3],polid=c(1,2,3)) > mod2 <- multinom(cbind(dem,ind,rep)~polid,data=pid1); > summary(mod2) (Intercept) polid ind rep > ct2 <- xtabs(y~paff+polid,pidata); > pid2<-data.frame(lib=ct2[,1],mod=ct2[,2],con=ct2[,3],paff=c(1,2,3)) > mod3 <- multinom(cbind(lib,mod,con)~paff,data=pid2) ; > summary(mod3) (Intercept) paff mod con State the models that are fitted by mod2 and mod3. State the equivalent log-linear models. 5
6 i. mod2 : Political ideology is assigned a continuous score v j This is the row-effect model Y i j Multinomial(π 1 j, π 2 j, π 3 j ), log π i j π 1 j = α i + β i v j where π 1 j = P (Democrat v j ), π 2 j = P (Independent v j ), and π 3 j = P (Republican v j ). The fitted multinominal model is log π 2 j π 1 j = v j log π 3 j π 1 j = v j The equivalent log-linear model is log(λ ij ) = µ + α i + β j + γ i v j, where γ i : separate parameter of v j for each row α i : deviation from democrat party affiliation β j : deviation from liberal political ideology α 1 = β 1 = γ 3 = 0 ii. mod3 : Party affiliation is assigned a continuous score u i This is the column-effect model Y j i Multinomial(π 1 i, π 2 i, π 3 i ), log π j i π 1 i = α j + β j v i where π 1 i = P (Liberal u i ), π 2 i = P (Moderate u i ), and π 3 i = P (Conservative u i ). The fitted multinominal model is log π 2 i π 1 i = u i log π 3 i π 1 i = u i The equivalent log-linear model is log(λ ij ) = µ + α i + β j + γ j u i, where γ j : separate parameter of u i for each column α i : deviation from democrat party affiliation β j : deviation from liberal political ideology α 1 = β 1 = γ 3 = 0 (c) What other model(s) could be potentially appropriate for these data if the models above do not fit well? linear-by-linear model : Saturated model : log E(Y ij ) = µ + α i + β j + γ(u i v j ) log E(Y ij ) = µ + α i + β j + (αβ) ij 6
7 4. [Methods qualifying exam, August 2010: use paper and pencil.] As part of housing market research, you collected data on number or apartments sold for various price ranges (low, medium, high), in different locations (district A, B and C), building types (high-rise, low-rise, townhouse), and for three buyer incomes (low, medium, high). You d like to model the probability of selling an apartment in a specific price range, as function of the remaining characteristics. You consider two models: (1) price range is independent of all other variables, and (2) log of price range is determined by an additive contribution of location, building type, and buyer income. (a) You only have access to software which fits Poisson regression. Write models (1) and (2) that you d fit in terms of Poisson distribution. Clearly define the notation and state the assumptions. The surrogate linear model 1: Y ijkl ind P oisson(λ ijkl ), where log(λ ijkl ) = µ + loc i + build j + buy k + (loc build) ij + (loc buy) ik + (build buy) jk + (loc build buy) ijk + price l The surrogate linear model 2: Y ijkl ind P oisson(λ ijkl ), where log(λ ijkl ) = µ + loc i + build j + buy k + (loc build) ij + (loc buy) ik + (build buy) jk + (loc build buy) ijk + price l + (loc price) il + (build price) jl + (buy price) kl (b) What test statistic would you use to compare the two models, and what is its reference distribution of the test statistic under H 0 (including degrees of freedom)? We use the G 2 statistic which compares the deviances of the two models. Under H 0 it follows χ 2 distribution with (3 1)(3 1) + (3 1)(3 1) + (3 1)(3 1) = 12 degrees of freedom. 5. Download the dataset minn38 available in the package MASS, and check R documentation for further details. The goal of this problem is to use surrogate log-linear models to determine how phs is affected by the other three factors. (a) Analysis Use variable selection to find an appropriate surrogate log-linear model for the dataset. You may use step() to find a candidate model, then check locally the nearby models by drop1() and add1(), which would have larger AIC but may be worth attention due to other considerations, such as interpretation, log likelihood ratio tests, etc. Report the basic facts such as the minimum and maximum models you are considering, the meaning of the final model(s), etc., but be selective with the computer outputs you include. Answer : 7
8 i. initial model (minimum model) : no effect of covariate on phs. Surrograte log linear models was used for analysis. Therefore we model the response as a Poisson random variable, conditional on the total number of observations for each covariate pattern hs*sex*fol. This is done by including in the model a full set of linear and interaction terms involving hs*sex*fol. We start from the minimal model, which expresses a constant probability of phs status across all covariate patterns. In other words, probability of the phs is independent of other covariate patterns. > fit<-glm(f~hs*fol*sex + phs, family=poisson, data=minn38) ii. additive effect of covariates on phs. We consider the dependence of phs status on one covariate at a time, by adding interaction terms covariate:phs > fit.linear<-glm(f~hs*fol*sex + phs + phs:(hs+fol+sex), family=poisson, data=minn38) > add1(fit.surrogate, ~. + phs:(hs+sex+fol), test="chisq") Model: f ~ hs * sex * fol + phs Df Deviance AIC LRT Pr(Chi) <none> hs:phs < 2.2e-16 *** sex:phs < 2.2e-16 *** fol:phs < 2.2e-16 *** --- Signif. codes: 0 *** ** 0.01 * All the interaction terms are significantly different from 0, so we add all the three interaction terms into the model. iii. select variables using step > fullformula<-f~hs*fol*sex + phs + phs:(hs+fol+sex)^3 > stepresult<-step(fit, fullformula) > stepresult Call: glm(formula = f ~ hs + fol + sex + phs + hs:fol + hs:sex + fol:sex + fol:phs + hs:phs + sex:phs + hs:fol:sex + fol:sex:phs + hs:fol:phs, family = poisson, data = minn38) iv. non-additive effect of covariates on phs. We further attempt to express the dependence of the probability of phs status on covariates, while accounting for possible interactions between the covariates. > fit2<-update(fit1,.~.+ phs:(hs+fol+sex)^2) Compare four models : > anova(fit, fit.linear, stepresult, fit2) Analysis of Deviance Table Model 1: f ~ hs * fol * sex + phs Model 2: f ~ hs * fol * sex + phs + phs:(hs + fol + sex) Model 3: f ~ hs + fol + sex + phs + hs:fol + hs:sex + fol:sex + fol:phs + hs:phs + sex:phs + hs:fol:sex + fol:sex:phs + hs:fol:phs Model 4: f ~ hs + fol + sex + phs + hs:fol + hs:sex + fol:sex + hs:phs + fol:phs + sex:phs + hs:fol:sex + hs:fol:phs + hs:sex:phs + 8
9 fol:sex:phs Resid. Df Resid. Dev Df Deviance AIC Model 2, 3, 4 have quite similar AIC. Therefore, Since we also perform a Likelihood Ratio test comparing the additive model to the model with interaction and which interactions are good for fitting the model. The test rejects the hypothesis of no interaction between model 2 and model 3. Hence, finally we retains model 3 as final model. (b) Presentation Present your final fit in the form of estimated cell probabilities. Running the final model above to predict the probabilty of observations in each cell. R Code > mnames <- lapply(minn38[,-5],levels) > p <- predict(interactions, expand.grid(mnames), type="response") > p <- matrix( p, ncol=3, byrow=t, dimnames=list(null, mnames[[1]]) ) > pr <- p/drop(p%*%rep(1,3)) > cbind(expand.grid(mnames), prob=round(pr,2)) hs phs fol sex prob.l prob.m prob.u 1 L C F1 F M C F1 F U C F1 F L E F1 F M E F1 F U E F1 F L N F1 F M N F1 F U N F1 F L O F1 F M O F1 F U O F1 F L C F2 F M C F2 F U C F2 F L E F2 F M E F2 F U E F2 F L N F2 F M N F2 F M N F7 M U N F7 M L O F7 M M O F7 M U O F7 M (c) Discussion Briefly discuss in words your findings. The model expresses the probability of phs status as a function of the three covariates. The model contains the additive terms for all three covariates. hs:fol and sex:fol have non-additive effects on the probability of phs. In mathematical terms, the model can be written as where phs Multinomial(π 1, π 2, π 3, π 4 ) log(π j ) = β 0j + β 1j hs + β 2j fol + β 3j sex + β 4j hs : fol + β 5j sex : fol 9
Loglinear models. STAT 526 Professor Olga Vitek
Loglinear models STAT 526 Professor Olga Vitek April 19, 2011 8 Can Use Poisson Likelihood To Model Both Poisson and Multinomial Counts 8-1 Recall: Poisson Distribution Probability distribution: Y - number
More informationLog-linear Models for Contingency Tables
Log-linear Models for Contingency Tables Statistics 149 Spring 2006 Copyright 2006 by Mark E. Irwin Log-linear Models for Two-way Contingency Tables Example: Business Administration Majors and Gender A
More informationSample solutions. Stat 8051 Homework 8
Sample solutions Stat 8051 Homework 8 Problem 1: Faraway Exercise 3.1 A plot of the time series reveals kind of a fluctuating pattern: Trying to fit poisson regression models yields a quadratic model if
More informationSTA 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) T In 2 2 tables, statistical independence is equivalent to a population
More informationSTAT 510 Final Exam Spring 2015
STAT 510 Final Exam Spring 2015 Instructions: The is a closed-notes, closed-book exam No calculator or electronic device of any kind may be used Use nothing but a pen or pencil Please write your name and
More informationNATIONAL UNIVERSITY OF SINGAPORE EXAMINATION. ST3241 Categorical Data Analysis. (Semester II: ) April/May, 2011 Time Allowed : 2 Hours
NATIONAL UNIVERSITY OF SINGAPORE EXAMINATION Categorical Data Analysis (Semester II: 2010 2011) April/May, 2011 Time Allowed : 2 Hours Matriculation No: Seat No: Grade Table Question 1 2 3 4 5 6 Full marks
More informationST3241 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 informationPoisson Regression. The Training Data
The Training Data Poisson Regression Office workers at a large insurance company are randomly assigned to one of 3 computer use training programmes, and their number of calls to IT support during the following
More informationTwo Hours. Mathematical formula books and statistical tables are to be provided THE UNIVERSITY OF MANCHESTER. 26 May :00 16:00
Two Hours MATH38052 Mathematical formula books and statistical tables are to be provided THE UNIVERSITY OF MANCHESTER GENERALISED LINEAR MODELS 26 May 2016 14:00 16:00 Answer ALL TWO questions in Section
More informationNATIONAL UNIVERSITY OF SINGAPORE EXAMINATION (SOLUTIONS) ST3241 Categorical Data Analysis. (Semester II: )
NATIONAL UNIVERSITY OF SINGAPORE EXAMINATION (SOLUTIONS) Categorical Data Analysis (Semester II: 2010 2011) April/May, 2011 Time Allowed : 2 Hours Matriculation No: Seat No: Grade Table Question 1 2 3
More informationSTAT 526 Spring Midterm 1. Wednesday February 2, 2011
STAT 526 Spring 2011 Midterm 1 Wednesday February 2, 2011 Time: 2 hours Name (please print): Show all your work and calculations. Partial credit will be given for work that is partially correct. Points
More informationLogistic Regressions. Stat 430
Logistic Regressions Stat 430 Final Project Final Project is, again, team based You will decide on a project - only constraint is: you are supposed to use techniques for a solution that are related to
More information9 Generalized Linear Models
9 Generalized Linear Models The Generalized Linear Model (GLM) is a model which has been built to include a wide range of different models you already know, e.g. ANOVA and multiple linear regression models
More informationSTA 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 informationCh 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
More informationHomework 3 - Solution
STAT 526 - Spring 2011 Homework 3 - Solution Olga Vitek Each part of the problems 5 points 1. KNNL 25.17 (Note: you can choose either the restricted or the unrestricted version of the model. Please state
More informationGeneralized 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 informationPoisson Regression. Gelman & Hill Chapter 6. February 6, 2017
Poisson Regression Gelman & Hill Chapter 6 February 6, 2017 Military Coups Background: Sub-Sahara Africa has experienced a high proportion of regime changes due to military takeover of governments for
More informationUNIVERSITY OF TORONTO Faculty of Arts and Science
UNIVERSITY OF TORONTO Faculty of Arts and Science December 2013 Final Examination STA442H1F/2101HF Methods of Applied Statistics Jerry Brunner Duration - 3 hours Aids: Calculator Model(s): Any calculator
More informationExercise 5.4 Solution
Exercise 5.4 Solution Niels Richard Hansen University of Copenhagen May 7, 2010 1 5.4(a) > leukemia
More informationSTA 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 informationReview 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
More informationIntroduction to Statistical Data Analysis Lecture 7: The Chi-Square Distribution
Introduction to Statistical Data Analysis Lecture 7: The Chi-Square Distribution James V. Lambers Department of Mathematics The University of Southern Mississippi James V. Lambers Statistical Data Analysis
More informationPsych 230. Psychological Measurement and Statistics
Psych 230 Psychological Measurement and Statistics Pedro Wolf December 9, 2009 This Time. Non-Parametric statistics Chi-Square test One-way Two-way Statistical Testing 1. Decide which test to use 2. State
More informationStat 8053, Fall 2013: Multinomial Logistic Models
Stat 8053, Fall 2013: Multinomial Logistic Models Here is the example on page 269 of Agresti on food preference of alligators: s is size class, g is sex of the alligator, l is name of the lake, and f is
More informationLecture 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 informationModelling 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 informationMSH3 Generalized linear model
Contents MSH3 Generalized linear model 7 Log-Linear Model 231 7.1 Equivalence between GOF measures........... 231 7.2 Sampling distribution................... 234 7.3 Interpreting Log-Linear models..............
More informationPoisson Regression. James H. Steiger. Department of Psychology and Human Development Vanderbilt University
Poisson Regression James H. Steiger Department of Psychology and Human Development Vanderbilt University James H. Steiger (Vanderbilt University) Poisson Regression 1 / 49 Poisson Regression 1 Introduction
More informationVarious Issues in Fitting Contingency Tables
Various Issues in Fitting Contingency Tables Statistics 149 Spring 2006 Copyright 2006 by Mark E. Irwin Complete Tables with Zero Entries In contingency tables, it is possible to have zero entries in a
More informationR Hints for Chapter 10
R Hints for Chapter 10 The multiple logistic regression model assumes that the success probability p for a binomial random variable depends on independent variables or design variables x 1, x 2,, x k.
More informationExam Applied Statistical Regression. Good Luck!
Dr. M. Dettling Summer 2011 Exam Applied Statistical Regression Approved: Tables: Note: Any written material, calculator (without communication facility). Attached. All tests have to be done at the 5%-level.
More informationSTAT 526 Advanced Statistical Methodology
STAT 526 Advanced Statistical Methodology Fall 2017 Lecture Note 7 Contingency Table 0-0 Outline Introduction to Contingency Tables Testing Independence in Two-Way Contingency Tables Modeling Ordinal Associations
More informationSCHOOL OF MATHEMATICS AND STATISTICS. Linear and Generalised Linear Models
SCHOOL OF MATHEMATICS AND STATISTICS Linear and Generalised Linear Models Autumn Semester 2017 18 2 hours Attempt all the questions. The allocation of marks is shown in brackets. RESTRICTED OPEN BOOK EXAMINATION
More informationSTAC51: Categorical data Analysis
STAC51: Categorical data Analysis Mahinda Samarakoon April 6, 2016 Mahinda Samarakoon STAC51: Categorical data Analysis 1 / 25 Table of contents 1 Building and applying logistic regression models (Chap
More informationA Handbook of Statistical Analyses Using R. Brian S. Everitt and Torsten Hothorn
A Handbook of Statistical Analyses Using R Brian S. Everitt and Torsten Hothorn CHAPTER 6 Logistic Regression and Generalised Linear Models: Blood Screening, Women s Role in Society, and Colonic Polyps
More informationGeneralised linear models. Response variable can take a number of different formats
Generalised linear models Response variable can take a number of different formats Structure Limitations of linear models and GLM theory GLM for count data GLM for presence \ absence data GLM for proportion
More informationHomework 5 - Solution
STAT 526 - Spring 2011 Homework 5 - Solution Olga Vitek Each part of the problems 5 points 1. Agresti 10.1 (a) and (b). Let Patient Die Suicide Yes No sum Yes 1097 90 1187 No 203 435 638 sum 1300 525 1825
More informationR Output for Linear Models using functions lm(), gls() & glm()
LM 04 lm(), gls() &glm() 1 R Output for Linear Models using functions lm(), gls() & glm() Different kinds of output related to linear models can be obtained in R using function lm() {stats} in the base
More informationModeling Overdispersion
James H. Steiger Department of Psychology and Human Development Vanderbilt University Regression Modeling, 2009 1 Introduction 2 Introduction In this lecture we discuss the problem of overdispersion in
More informationMultiple Linear Regression
Multiple Linear Regression ST 430/514 Recall: a regression model describes how a dependent variable (or response) Y is affected, on average, by one or more independent variables (or factors, or covariates).
More informationDuration of Unemployment - Analysis of Deviance Table for Nested Models
Duration of Unemployment - Analysis of Deviance Table for Nested Models February 8, 2012 The data unemployment is included as a contingency table. The response is the duration of unemployment, gender and
More informationClassification. Chapter Introduction. 6.2 The Bayes classifier
Chapter 6 Classification 6.1 Introduction Often encountered in applications is the situation where the response variable Y takes values in a finite set of labels. For example, the response Y could encode
More informationUNIVERSITY 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 informationSTAT 526 Spring Final Exam. Thursday May 5, 2011
STAT 526 Spring 2011 Final Exam Thursday May 5, 2011 Time: 2 hours Name (please print): Show all your work and calculations. Partial credit will be given for work that is partially correct. Points will
More informationHomework 1 Solutions
36-720 Homework 1 Solutions Problem 3.4 (a) X 2 79.43 and G 2 90.33. We should compare each to a χ 2 distribution with (2 1)(3 1) 2 degrees of freedom. For each, the p-value is so small that S-plus reports
More informationLogistic Regression 21/05
Logistic Regression 21/05 Recall that we are trying to solve a classification problem in which features x i can be continuous or discrete (coded as 0/1) and the response y is discrete (0/1). Logistic regression
More information3. (a) (8 points) There is more than one way to correctly express the null hypothesis in matrix form. One way to state the null hypothesis is
Stat 501 Solutions and Comments on Exam 1 Spring 005-4 0-4 1. (a) (5 points) Y ~ N, -1-4 34 (b) (5 points) X (X,X ) = (5,8) ~ N ( 11.5, 0.9375 ) 3 1 (c) (10 points, for each part) (i), (ii), and (v) are
More informationGeneralized linear models for binary data. A better graphical exploratory data analysis. The simple linear logistic regression model
Stat 3302 (Spring 2017) Peter F. Craigmile Simple linear logistic regression (part 1) [Dobson and Barnett, 2008, Sections 7.1 7.3] Generalized linear models for binary data Beetles dose-response example
More informationBooklet 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 NBA attendance data........................ 2 2 Regression model for NBA attendances...............
More information8 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 informationSolution to Tutorial 7
1. (a) We first fit the independence model ST3241 Categorical Data Analysis I Semester II, 2012-2013 Solution to Tutorial 7 log µ ij = λ + λ X i + λ Y j, i = 1, 2, j = 1, 2. The parameter estimates are
More informationStatistical Methods III Statistics 212. Problem Set 2 - Answer Key
Statistical Methods III Statistics 212 Problem Set 2 - Answer Key 1. (Analysis to be turned in and discussed on Tuesday, April 24th) The data for this problem are taken from long-term followup of 1423
More informationLinear Regression Models P8111
Linear Regression Models P8111 Lecture 25 Jeff Goldsmith April 26, 2016 1 of 37 Today s Lecture Logistic regression / GLMs Model framework Interpretation Estimation 2 of 37 Linear regression Course started
More informationClinical Trials. Olli Saarela. September 18, Dalla Lana School of Public Health University of Toronto.
Introduction to Dalla Lana School of Public Health University of Toronto olli.saarela@utoronto.ca September 18, 2014 38-1 : a review 38-2 Evidence Ideal: to advance the knowledge-base of clinical medicine,
More informationLogistic Regression - problem 6.14
Logistic Regression - problem 6.14 Let x 1, x 2,, x m be given values of an input variable x and let Y 1,, Y m be independent binomial random variables whose distributions depend on the corresponding values
More informationFinal Exam. Name: Solution:
Final Exam. Name: Instructions. Answer all questions on the exam. Open books, open notes, but no electronic devices. The first 13 problems are worth 5 points each. The rest are worth 1 point each. HW1.
More informationSTAT 525 Fall Final exam. Tuesday December 14, 2010
STAT 525 Fall 2010 Final exam Tuesday December 14, 2010 Time: 2 hours Name (please print): Show all your work and calculations. Partial credit will be given for work that is partially correct. Points will
More informationStatistics 191 Introduction to Regression Analysis and Applied Statistics Practice Exam
Statistics 191 Introduction to Regression Analysis and Applied Statistics Practice Exam Prof. J. Taylor You may use your 4 single-sided pages of notes This exam is 14 pages long. There are 4 questions,
More informationssh 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 informationReview: what is a linear model. Y = β 0 + β 1 X 1 + β 2 X 2 + A model of the following form:
Outline for today What is a generalized linear model Linear predictors and link functions Example: fit a constant (the proportion) Analysis of deviance table Example: fit dose-response data using logistic
More information7/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
More information22s:152 Applied Linear Regression. RECALL: The Poisson distribution. Let Y be distributed as a Poisson random variable with the single parameter λ.
22s:152 Applied Linear Regression Chapter 15 Section 2: Poisson Regression RECALL: The Poisson distribution Let Y be distributed as a Poisson random variable with the single parameter λ. P (Y = y) = e
More informationSTA 450/4000 S: January
STA 450/4000 S: January 6 005 Notes Friday tutorial on R programming reminder office hours on - F; -4 R The book Modern Applied Statistics with S by Venables and Ripley is very useful. Make sure you have
More informationNon-Gaussian Response Variables
Non-Gaussian Response Variables What is the Generalized Model Doing? The fixed effects are like the factors in a traditional analysis of variance or linear model The random effects are different A generalized
More informationInteractions in Logistic Regression
Interactions in Logistic Regression > # UCBAdmissions is a 3-D table: Gender by Dept by Admit > # Same data in another format: > # One col for Yes counts, another for No counts. > Berkeley = read.table("http://www.utstat.toronto.edu/~brunner/312f12/
More informationPAPER 206 APPLIED STATISTICS
MATHEMATICAL TRIPOS Part III Thursday, 1 June, 2017 9:00 am to 12:00 pm PAPER 206 APPLIED STATISTICS Attempt no more than FOUR questions. There are SIX questions in total. The questions carry equal weight.
More informationChecking the Poisson assumption in the Poisson generalized linear model
Checking the Poisson assumption in the Poisson generalized linear model The Poisson regression model is a generalized linear model (glm) satisfying the following assumptions: The responses y i are independent
More informationIntroduction to General and Generalized Linear Models
Introduction to General and Generalized Linear Models Generalized Linear Models - part III Henrik Madsen Poul Thyregod Informatics and Mathematical Modelling Technical University of Denmark DK-2800 Kgs.
More informationRegression 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
More informationGeneralized Linear Models in R
Generalized Linear Models in R NO ORDER Kenneth K. Lopiano, Garvesh Raskutti, Dan Yang last modified 28 4 2013 1 Outline 1. Background and preliminaries 2. Data manipulation and exercises 3. Data structures
More informationMultinomial 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 informationChi-Square. Heibatollah Baghi, and Mastee Badii
1 Chi-Square Heibatollah Baghi, and Mastee Badii Different Scales, Different Measures of Association Scale of Both Variables Nominal Scale Measures of Association Pearson Chi-Square: χ 2 Ordinal Scale
More informationStat 579: Generalized Linear Models and Extensions
Stat 579: Generalized Linear Models and Extensions Yan Lu Jan, 2018, week 3 1 / 67 Hypothesis tests Likelihood ratio tests Wald tests Score tests 2 / 67 Generalized Likelihood ratio tests Let Y = (Y 1,
More informationMATH 644: Regression Analysis Methods
MATH 644: Regression Analysis Methods FINAL EXAM Fall, 2012 INSTRUCTIONS TO STUDENTS: 1. This test contains SIX questions. It comprises ELEVEN printed pages. 2. Answer ALL questions for a total of 100
More informationLecture 28 Chi-Square Analysis
Lecture 28 STAT 225 Introduction to Probability Models April 23, 2014 Whitney Huang Purdue University 28.1 χ 2 test for For a given contingency table, we want to test if two have a relationship or not
More informationMSH3 Generalized linear model Ch. 6 Count data models
Contents MSH3 Generalized linear model Ch. 6 Count data models 6 Count data model 208 6.1 Introduction: The Children Ever Born Data....... 208 6.2 The Poisson Distribution................. 210 6.3 Log-Linear
More informationCategorical Variables and Contingency Tables: Description and Inference
Categorical Variables and Contingency Tables: Description and Inference STAT 526 Professor Olga Vitek March 3, 2011 Reading: Agresti Ch. 1, 2 and 3 Faraway Ch. 4 3 Univariate Binomial and Multinomial Measurements
More informationMultiple comparisons - subsequent inferences for two-way ANOVA
1 Multiple comparisons - subsequent inferences for two-way ANOVA the kinds of inferences to be made after the F tests of a two-way ANOVA depend on the results if none of the F tests lead to rejection of
More informationSections 3.4, 3.5. Timothy Hanson. Department of Statistics, University of South Carolina. Stat 770: Categorical Data Analysis
Sections 3.4, 3.5 Timothy Hanson Department of Statistics, University of South Carolina Stat 770: Categorical Data Analysis 1 / 22 3.4 I J tables with ordinal outcomes Tests that take advantage of ordinal
More informationExperimental Design and Statistical Methods. Workshop LOGISTIC REGRESSION. Jesús Piedrafita Arilla.
Experimental Design and Statistical Methods Workshop LOGISTIC REGRESSION Jesús Piedrafita Arilla jesus.piedrafita@uab.cat Departament de Ciència Animal i dels Aliments Items Logistic regression model Logit
More informationMatched Pair Data. Stat 557 Heike Hofmann
Matched Pair Data Stat 557 Heike Hofmann Outline Marginal Homogeneity - review Binary Response with covariates Ordinal response Symmetric Models Subject-specific vs Marginal Model conditional logistic
More informationProportional Odds Logistic Regression. stat 557 Heike Hofmann
Proportional Odds Logistic Regression stat 557 Heike Hofmann Outline Proportional Odds Logistic Regression Model Definition Properties Latent Variables Intro to Loglinear Models Ordinal Response Y is categorical
More informationQuantitative Analysis and Empirical Methods
Hypothesis testing Sciences Po, Paris, CEE / LIEPP Introduction Hypotheses Procedure of hypothesis testing Two-tailed and one-tailed tests Statistical tests with categorical variables A hypothesis A testable
More informationReaction Days
Stat April 03 Week Fitting Individual Trajectories # Straight-line, constant rate of change fit > sdat = subset(sleepstudy, Subject == "37") > sdat Reaction Days Subject > lm.sdat = lm(reaction ~ Days)
More informationUsing R in 200D Luke Sonnet
Using R in 200D Luke Sonnet Contents Working with data frames 1 Working with variables........................................... 1 Analyzing data............................................... 3 Random
More informationSTATS216v Introduction to Statistical Learning Stanford University, Summer Midterm Exam (Solutions) Duration: 1 hours
Instructions: STATS216v Introduction to Statistical Learning Stanford University, Summer 2017 Remember the university honor code. Midterm Exam (Solutions) Duration: 1 hours Write your name and SUNet ID
More informationSTATISTICS 174: APPLIED STATISTICS FINAL EXAM DECEMBER 10, 2002
Time allowed: 3 HOURS. STATISTICS 174: APPLIED STATISTICS FINAL EXAM DECEMBER 10, 2002 This is an open book exam: all course notes and the text are allowed, and you are expected to use your own calculator.
More informationYou may use your calculator and a single page of notes.
LAST NAME (Please Print): KEY FIRST NAME (Please Print): HONOR PLEDGE (Please Sign): Statistics 111 Midterm 4 This is a closed book exam. You may use your calculator and a single page of notes. The room
More informationExercise 7.4 [16 points]
STATISTICS 226, Winter 1997, Homework 5 1 Exercise 7.4 [16 points] a. [3 points] (A: Age, G: Gestation, I: Infant Survival, S: Smoking.) Model G 2 d.f. (AGIS).008 0 0 (AGI, AIS, AGS, GIS).367 1 (AG, AI,
More informationInference for Regression
Inference for Regression Section 9.4 Cathy Poliak, Ph.D. cathy@math.uh.edu Office in Fleming 11c Department of Mathematics University of Houston Lecture 13b - 3339 Cathy Poliak, Ph.D. cathy@math.uh.edu
More informationUnit 9: Inferences for Proportions and Count Data
Unit 9: Inferences for Proportions and Count Data Statistics 571: Statistical Methods Ramón V. León 12/15/2008 Unit 9 - Stat 571 - Ramón V. León 1 Large Sample Confidence Interval for Proportion ( pˆ p)
More informationLecture 41 Sections Mon, Apr 7, 2008
Lecture 41 Sections 14.1-14.3 Hampden-Sydney College Mon, Apr 7, 2008 Outline 1 2 3 4 5 one-proportion test that we just studied allows us to test a hypothesis concerning one proportion, or two categories,
More informationPAPER 218 STATISTICAL LEARNING IN PRACTICE
MATHEMATICAL TRIPOS Part III Thursday, 7 June, 2018 9:00 am to 12:00 pm PAPER 218 STATISTICAL LEARNING IN PRACTICE Attempt no more than FOUR questions. There are SIX questions in total. The questions carry
More informationStat 135 Fall 2013 FINAL EXAM December 18, 2013
Stat 135 Fall 2013 FINAL EXAM December 18, 2013 Name: Person on right SID: Person on left There will be one, double sided, handwritten, 8.5in x 11in page of notes allowed during the exam. The exam is closed
More information12 Modelling Binomial Response Data
c 2005, Anthony C. Brooms Statistical Modelling and Data Analysis 12 Modelling Binomial Response Data 12.1 Examples of Binary Response Data Binary response data arise when an observation on an individual
More information1. Hypothesis testing through analysis of deviance. 3. Model & variable selection - stepwise aproaches
Sta 216, Lecture 4 Last Time: Logistic regression example, existence/uniqueness of MLEs Today s Class: 1. Hypothesis testing through analysis of deviance 2. Standard errors & confidence intervals 3. Model
More informationˆπ(x) = exp(ˆα + ˆβ T x) 1 + exp(ˆα + ˆβ T.
Exam 3 Review Suppose that X i = x =(x 1,, x k ) T is observed and that Y i X i = x i independent Binomial(n i,π(x i )) for i =1,, N where ˆπ(x) = exp(ˆα + ˆβ T x) 1 + exp(ˆα + ˆβ T x) This is called the
More informationSection 9.4. Notation. Requirements. Definition. Inferences About Two Means (Matched Pairs) Examples
Objective Section 9.4 Inferences About Two Means (Matched Pairs) Compare of two matched-paired means using two samples from each population. Hypothesis Tests and Confidence Intervals of two dependent means
More informationA 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