diff --git a/logistic_regression.R b/logistic_regression.R
index 932935b..e0d7a17 100644
--- a/logistic_regression.R
+++ b/logistic_regression.R
@@ -261,23 +261,25 @@ for(i in my_reg_params_meds1) {
   cat("=================================================================================\n")
 }
 
-
 #-----------------------------
-# outcome: death 
+# outcome: t1_resp_recoded ~ obesity + asthma
 # data: fp adults
 #-----------------------------
-my_reg_params_meds1v2 = my_reg_params_meds[!my_reg_params_meds%in%c("mosaic"
+my_reg_params_meds2v2 = my_reg_params_meds[!my_reg_params_meds%in%c("mosaic"
                                                                     , "flustat"
                                                                     , "onset_2_initial"
                                                                     , "onsfindeath"
                                                                     , "ia_exac_copd"
-                                                                    , "death")]
+                                                                    , "death"
+                                                                    , "t1_resp_recoded"
+                                                                    , "T1_resp_score"
+                                                                    , "asthma")]
 
 
 
-for(i in my_reg_params_meds1v2) { 
+for(i in my_reg_params_meds2v2) { 
   #print (i)
-  p_form = as.formula(paste("death ~ ", i ,sep = ""))
+  p_form = as.formula(paste("t1_resp_recoded ~ obesity + asthma + ", i ,sep = ""))
   print(p_form)
   model_reg = glm(p_form , family = binomial, data = my_reg_data_med)
   print(summary(model_reg))