# clinical unpaired stats: DONE
# chi sq or similar test for binary outcome vars : DONE
# regression analysis
# redo obesity analysis without asthmatics: although you can account for this in reg
but just do it for completion!: DONE


=================================================================================
death ~ obesity + il12p70_npa1

Call:
glm(formula = p_form, family = binomial, data = my_reg_data_med)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-1.1117  -0.5249  -0.4427  -0.4302   2.2027  

Coefficients:
             Estimate Std. Error z value Pr(>|z|)    
(Intercept)  -2.36991    0.41816  -5.668 1.45e-08 ***
obesity       0.18899    0.63841   0.296   0.7672    
il12p70_npa1  0.12389    0.04843   2.558   0.0105 *  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 89.632  on 104  degrees of freedom
Residual deviance: 82.487  on 102  degrees of freedom
  (28 observations deleted due to missingness)
AIC: 88.487

Number of Fisher Scoring iterations: 4


##############################
=================================================================================
death ~ il12p70_npa1

Call:
glm(formula = p_form, family = binomial, data = my_reg_data_med)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-1.0765  -0.4935  -0.4376  -0.4376   2.1880  

Coefficients:
             Estimate Std. Error z value Pr(>|z|)    
(Intercept)  -2.33522    0.39856  -5.859 4.65e-09 ***
il12p70_npa1  0.12699    0.04741   2.679  0.00739 ** 
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 89.632  on 104  degrees of freedom
Residual deviance: 82.574  on 103  degrees of freedom
  (28 observations deleted due to missingness)
AIC: 86.574

Number of Fisher Scoring iterations: 4

Waiting for profiling to be done...
                     OR      2.5 %   97.5 %
(Intercept)  0.09678937 0.04094537 0.198725
il12p70_npa1 1.13540367 1.03471817 1.249931
[1] 105

=================================================================================
death ~ il1_npa1

Call:
glm(formula = p_form, family = binomial, data = my_reg_data_med)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-1.3883  -0.4893  -0.4793  -0.4792   2.1073  

Coefficients:
              Estimate Std. Error z value Pr(>|z|)    
(Intercept) -2.1068569  0.3347016  -6.295 3.08e-10 ***
il1_npa1     0.0010360  0.0003762   2.754  0.00588 ** 
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 89.632  on 104  degrees of freedom
Residual deviance: 82.215  on 103  degrees of freedom
  (28 observations deleted due to missingness)
AIC: 86.215

Number of Fisher Scoring iterations: 4

Waiting for profiling to be done...
                   OR      2.5 %    97.5 %
(Intercept) 0.1216196 0.05949241 0.2236722
il1_npa1    1.0010366 1.00030121 1.0018148
[1] 105
=================================================================================