diff --git a/finalproject b/finalproject index 3678e1f0d79d414eb912f4ba1b02237376ad1582..bb74efd964810a97c97f4ed9456f96cf2e1df3bc 100644 --- a/finalproject +++ b/finalproject @@ -199,10 +199,10 @@ fiveDrinks_mod <- glm(fiveDrinks_var ~ pride + truth + responsibility + friends summary(fiveDrinks_mod) fiveDrinks_mod_1 <- glm(fiveDrinks_var ~ Wpdrink + friends, data = data, family = "binomial") -summary(fiveDrinks_mod_1) +#summary(fiveDrinks_mod_1) fiveDrinks_mod_2 <- glm(fiveDrinks_var ~ friends, data = data, family = "binomial") -summary(fiveDrinks_mod_2) +#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 @@ -211,13 +211,19 @@ 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) +#summary(fiveDrinks_mod_3) fiveDrinks_mod_4 <- glm(fiveDrinks_var ~ Wpdrink, data = data, family = "binomial") -summary(fiveDrinks_mod_4) +#summary(fiveDrinks_mod_4) anova(fiveDrinks_mod_3, fiveDrinks_mod_4, test="LRT") +#drop in deviance test compares residual deivances from two models +# source: https://bookdown.org/roback/bookdown-BeyondMLR/ch-poissonreg.html#cs-philippines +# source: https://bookdown.org/roback/bookdown-BeyondMLR/ch-logreg.html +anova(fiveDrinks_mod_3, fiveDrinks_mod, test = "Chisq") + + ```