diff --git a/.~lock.lab5_expression3.docx# b/.~lock.lab5_expression3.docx#
index 9c2cc5e..1269447 100644
--- a/.~lock.lab5_expression3.docx#
+++ b/.~lock.lab5_expression3.docx#
@@ -1 +1 @@
-,noah,NovaArchSys,13.10.2022 15:00,file:///home/noah/.config/libreoffice/4;
\ No newline at end of file
+,noah,NovaArchSys,13.10.2022 16:10,file:///home/noah/.config/libreoffice/4;
\ No newline at end of file
diff --git a/Schrick-Noah_CS-6643_Lab-5.R b/Schrick-Noah_CS-6643_Lab-5.R
index 35eb43f..48e295c 100644
--- a/Schrick-Noah_CS-6643_Lab-5.R
+++ b/Schrick-Noah_CS-6643_Lab-5.R
@@ -48,4 +48,24 @@ dim(GxS.covfilter)
 pheno.factor <- as.factor(colnames(GxS.covfilter))
 pheno.factor
 str(pheno.factor)
-levels(pheno.factor)
\ No newline at end of file
+levels(pheno.factor)
+
+
+#### Part A: Logistic Regression
+# make sure HC is the reference level
+pheno.factor.relevel <- relevel(pheno.factor,"HC")  
+levels(pheno.factor.relevel) 
+# also rename levels "0"/"1" from 1/2
+levels(pheno.factor.relevel)[levels(pheno.factor.relevel)=="MDD"] <- 1
+levels(pheno.factor.relevel)[levels(pheno.factor.relevel)=="HC"] <- 0
+
+# Fit logistic model of first gene to phenotype data
+gene.row <- 2
+gene.name <-rownames(GxS.covfilter)[gene.row]
+gene.expr <- GxS.covfilter[gene.row,] 
+gene.fit <- glm(pheno.factor.relevel~gene.expr, family=binomial)  
+summary(gene.fit)
+coeff.mat <- coef(summary(gene.fit))
+b0 <- coeff.mat[1,1]
+b1 <- coeff.mat[2,1]
+b1.pval <- coeff.mat[2,4]
\ No newline at end of file
diff --git a/lab5_expression3.docx b/lab5_expression3.docx
index a901f24..74dc797 100644
Binary files a/lab5_expression3.docx and b/lab5_expression3.docx differ