#remove outliers

cooksd <- cooks.distance(glm(consensus ~ .,
                             family = "binomial",
                             data = green))

plot(cooksd,
     pch="*",
     cex=2,
     main="Influential Obs by Cooks distance")
abline(h = 4*mean(cooksd, na.rm=T), col="red")

outliers <- rownames(green[cooksd > 4*mean(cooksd, na.rm=T), ])
print(outliers)