diff --git a/logistic_regression.R b/logistic_regression.R
index d610cd4..e47289b 100644
--- a/logistic_regression.R
+++ b/logistic_regression.R
@@ -288,3 +288,17 @@ sink()
 
 cat("=================================================================================\n")
 
+model = glm(death ~ obesity, family = binomial, data = my_reg_data)
+summary(model); nobs(model)
+
+model2 = glm(t1_resp_recoded ~ obesity, family = binomial, data = my_reg_data)
+summary(model2); nobs(model2)
+
+model3 = glm(t1_resp_recoded ~ obesity + ia_cxr+ com_noasthma
+             , family = binomial, data = my_reg_data)
+summary(model3); nobs(model3)
+
+
+model4 = glm(death ~ obesity + ia_cxr + com_noasthma
+             , family = binomial, data = my_reg_data)
+summary(model4); nobs(model4)