Update finalproject

finalproject

@@ -194,13 +194,29 @@ round(cor(data[, -c(1, 2, 3, 4, 6)], 3))
 # making models
 fiveDrinks_threshold <- mean(data$X30.5drinks)
 data$fiveDrinks_var <- ifelse(data$X30.5drinks >= fiveDrinks_threshold, 1, 0)
-#fiveDrinks_mod <- glm(fiveDrinks_var ~ pride + truth + responsibility + friends + fix.problems + decision + excite + hard.work + safe + best.school + talk.adult + grades + Wpdrink + N.safe, data = data, family = "binomial")
 
-fiveDrinks_mod <- glm(fiveDrinks_var ~ Wpdrink, data = data, family = "binomial")
+fiveDrinks_mod <- glm(fiveDrinks_var ~ pride + truth + responsibility + friends + fix.problems + decision + excite + hard.work + safe + best.school + talk.adult + grades + Wpdrink + N.safe, data = data, family = "binomial")
 summary(fiveDrinks_mod)
 
-fiveDrinks_mod <- glm(fiveDrinks_var ~ friends, data = data, family = "binomial")
-summary(fiveDrinks_mod)
+fiveDrinks_mod_1 <- glm(fiveDrinks_var ~ Wpdrink + friends, data = data, family = "binomial")
+summary(fiveDrinks_mod_1)
+
+fiveDrinks_mod_2 <- glm(fiveDrinks_var ~ friends, data = data, family = "binomial")
+summary(fiveDrinks_mod_2)
+
+#likelihood ratio test to test whether the observed difference in model fits is statistically significant
+# source: https://www.listendata.com/2016/07/insignificant-levels-of-categorical-variable.html
+
+anova(fiveDrinks_mod_1, fiveDrinks_mod_2, test="LRT")
+
+
+fiveDrinks_mod_3 <- glm(fiveDrinks_var ~ Wpdrink + friends, data = data, family = "binomial")
+summary(fiveDrinks_mod_3)
+
+fiveDrinks_mod_4 <- glm(fiveDrinks_var ~ Wpdrink, data = data, family = "binomial")
+summary(fiveDrinks_mod_4)
+
+anova(fiveDrinks_mod_3, fiveDrinks_mod_4, test="LRT")
 
 ```