Stat 431 Assignment 3 solutions
1. [ 12 marks]
(a) We start by fitting the interaction model
logit(πi) = β0 + β1xi1 + β2xi2 + β3xi1xi2
The Wald-test for H0 : β3 (coefficient of interaction term) gives a p−value << 0.05,
hence we can not drop the interaction term. Since we can not...
po
The Wald-test for H0 : β3 (coefficient of interaction term) gives a p − value << 0.05,
hence we can not drop the interaction term. Since we can not further simplify this
model, it is the best logit model we can fit. [2]
t
no
cells<-read.table("cells.txt", header=T)
cells$resp<-cbind(cells$Cell, 200-cells$Cell)
D
Call:
glm(formula = resp ~ tnf + ifn + tnf * ifn, family = binomial,
-
data = cells)
ht
Deviance Residuals:
Min 1Q Median 3Q Max
-5.2331 -3.4359 -0.7136 2.8511 4.8272
ig
yr
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -1.730e+00 7.055e-02 -24.524 < 2e-16 ***
op
(Dispersion parameter for binomial family taken to be 1)
th
Null deviance: 1309.96 on 15 degrees of freedom
Residual deviance: 192.36 on 12 degrees of freedom
ns
AIC: 281.84
ow
However, the deviance test indicate that this model is not adequate (compared to
the saturated model).
r
Interpretation [4]
to
• exp(βb0 ) = 0.1772844: odds of response when both TNF dose xi1 = 0 and IFN
c
dose xi2 = 0. This odds is < 1, which means the probability of cell differentiation
ru
is very low (compared to no cell differentiation).
• β2 is the change in the log odds of cell differentiation associated with 1 unit
st
of increase in IFN, when TNF dose xi1 = 0. As βb2 = 1.068, the log odds (or
In
the probability) of cell differentiation increase as IFN dose increases, when TNF
dose is zero.
• When the dose of TNF (xi1 ) increase by 1 unit, the log odds of response change
by β1 + β3 xi2 , which is a linear function of IFN dose (xi2 ).
– βb1 = 2.521 is the amount of increase in log odds due to 1 unit of increase in
TNF, when IFN dose xi2 = 0.
1
, – We see that βb3 = 3.839 (the slope of the change in log odds), it means
that the higher the dose of IFN, the bigger the increase in log odds of cell
differentiation due to the increase of TNF dose.
In summary, the significant positive effect of interaction term β3 implies
that, the higher the TNF dose, the stronger the impact of IFN dose on the
probability of cell differentiation, and vice versa.
!
st
• Similarly when the dose of IFN dose (xi2 ) increase by 1 unit, the log odds of
po
response change by β2 + β3 xi1 , which is a linear function of TNF dose (xi1 ).
– βb2 = 1.068 is the amount of increase in log odds due to 1 unit of increase in
t
IFN, when TNF dose xi1 = 0.
no
(b) Deviance residuals plots demonstrate the poor fit of the selected logistic model in
(a), which agrees with the Deviance test result. [2]
plot(cells$ifn, rd,ylim=c(-8,8), xlab="IFN dose",ylab="DEVIANCE RESIDUALS")
abline(h=-2); abline(h= 2)
plot(fv, rd,ylim=c(-8,8), xlab="FITTED VALUE",ylab="DEVIANCE RESIDUALS")
abline(h=-2); abline(h= 2)
r
to
(c) The best fitted cloglog model is the main effects model [2]
c
log(− log(1 − πi )) = β0 + β1 xi1 + β2 xi2
ru
However, changing to cloglog link seems do not improve the fit at all, and the residual
st
plots are equally bad as the logit model. [2]
In
2
The benefits of buying summaries with Stuvia:
Guaranteed quality through customer reviews
Stuvia customers have reviewed more than 700,000 summaries. This how you know that you are buying the best documents.
Quick and easy check-out
You can quickly pay through credit card or Stuvia-credit for the summaries. There is no membership needed.
Focus on what matters
Your fellow students write the study notes themselves, which is why the documents are always reliable and up-to-date. This ensures you quickly get to the core!
Frequently asked questions
What do I get when I buy this document?
You get a PDF, available immediately after your purchase. The purchased document is accessible anytime, anywhere and indefinitely through your profile.
Satisfaction guarantee: how does it work?
Our satisfaction guarantee ensures that you always find a study document that suits you well. You fill out a form, and our customer service team takes care of the rest.
Who am I buying these notes from?
Stuvia is a marketplace, so you are not buying this document from us, but from seller Abbyy01. Stuvia facilitates payment to the seller.
Will I be stuck with a subscription?
No, you only buy these notes for $8.49. You're not tied to anything after your purchase.
