Consider a dataset consisting of 610 males involved in a study of coronary heart disease. The outcome variable is CHD status (1 = case, 0 = non case), the exposure variable of interest is CAT which is a dichotomous variable that indicates high (coded 1) or normal (coded 0) catecholamine level. The only other variables recorded in the data set are AGE (1 = age > 55, 0 = age ≤ 55) and ECG (1 = abnormal, 0 = normal). The dataset involving the above variables is given as follows: data data1; input cases total CAT AGE ECG; cards; 17 275 0 0 0 15 121 0 1 0 7 59 0 0 1 5 32 0 1 1 1 8 1 0 0 10 40 1 1 0 4 17 1 0 1 14 58 1 1 1 ; run; We are interested in the following logistic model: As given in the SAS code above, the model is a full model with all main effects and interactions (both two way and three way interactions). A main effect model can be obtained from SAS by this model statement – model cases/total = AGE CAT ECG / cl; Perform a hypothesis test to see if the interactions (including all of the two-way and the three-way Interactions) help with the model using a likelihood ratio test (LRT) to compare the full model and the main effect model using alpha of 0.05.
Consider a dataset consisting of 610 males involved in a study of coronary heart disease. The outcome variable is CHD status (1 = case, 0 = non case), the exposure variable of interest is CAT which is a dichotomous variable that indicates high (coded 1) or normal (coded 0) catecholamine level. The only other variables recorded in the data set are AGE (1 = age > 55, 0 = age ≤ 55) and ECG (1 = abnormal, 0 = normal). The dataset involving the above variables is given as follows:
data data1;
input cases total CAT AGE ECG;
cards;
17 275 0 0 0
15 121 0 1 0
7 59 0 0 1
5 32 0 1 1
1 8 1 0 0
10 40 1 1 0
4 17 1 0 1
14 58 1 1 1
;
run;
We are interested in the following logistic model:
As given in the SAS code above, the model is a full model with all main effects and interactions (both two way and three way interactions). A main effect model can be obtained from SAS by this model statement –
model cases/total = AGE CAT ECG / cl;
Perform a hypothesis test to see if the interactions (including all of the two-way and the three-way Interactions) help with the model using a likelihood ratio test (LRT) to compare the full model and the main effect model using alpha of 0.05.
![+
Criterion
AIC
SC
-2 Log L
Intercept Only
448.851
453.264
446.851
Criterion
AIC
448.851
SC
453.264
-2 Log L 446.851
Model Fit Statistics
Intercept Only
Log Likelihood
438.648
473.955
422.648
Output from SAS - reduced model (without interaction)
Intercept and Covariates
Model Fit Statistics
432.352
450.006
424.352
Log Likelihood
Full Log Likelihood
44.990
80.297
28.990
Intercept and Covariates
Full Log Likelihood
38.694
56.348
30.694
0](/v2/_next/image?url=https%3A%2F%2Fcontent.bartleby.com%2Fqna-images%2Fquestion%2F45415dcf-fee2-485b-baaf-08b1e98d4ec0%2F47a90aea-6d96-49ee-af5e-734fc48864bf%2Fo18wpyc_processed.png&w=3840&q=75)
![T₁ = P(CHD; = 1| AGE,, CAT, ECG;) logit() =B₁ + B₁AGE; + B₂CAT; + B₂ECG;](/v2/_next/image?url=https%3A%2F%2Fcontent.bartleby.com%2Fqna-images%2Fquestion%2F45415dcf-fee2-485b-baaf-08b1e98d4ec0%2F47a90aea-6d96-49ee-af5e-734fc48864bf%2Foacd2ja_processed.png&w=3840&q=75)
![](/static/compass_v2/shared-icons/check-mark.png)
Trending now
This is a popular solution!
Step by step
Solved in 5 steps with 22 images
![Blurred answer](/static/compass_v2/solution-images/blurred-answer.jpg)
![Glencoe Algebra 1, Student Edition, 9780079039897…](https://www.bartleby.com/isbn_cover_images/9780079039897/9780079039897_smallCoverImage.jpg)
![Calculus For The Life Sciences](https://www.bartleby.com/isbn_cover_images/9780321964038/9780321964038_smallCoverImage.gif)
![Glencoe Algebra 1, Student Edition, 9780079039897…](https://www.bartleby.com/isbn_cover_images/9780079039897/9780079039897_smallCoverImage.jpg)
![Calculus For The Life Sciences](https://www.bartleby.com/isbn_cover_images/9780321964038/9780321964038_smallCoverImage.gif)