www = "http://www.biostatisticaumg.it/dataset/tossicologia.csv"
tossicologia = read.csv(www, header = TRUE)
# tossicologia = read.csv(file.choose(), header = TRUE)
attach(tossicologia)
head(tossicologia)
str(tossicologia)

library(xtable)
xtable(head(tossicologia))

boxplot(eta, horizontal = TRUE, col = "orange", main = "Age")

boxplot(eta ~ genere, horizontal = TRUE, col = c("magenta", "cyan"), varwidth = TRUE, main = "Age vs. Gender")

xtable(table(genere))

xtable(table(causa))

par(mfrow = c(1,2))
plot(genere  ~ tossicologico, col = c("magenta", "cyan"), main = "Gender vs. Toxicology")
boxplot(eta ~ tossicologico, horizontal = TRUE, col = c("green", "red"), varwidth = TRUE, main = "Age vs. Toxicology")

xtable(table(genere,  tossicologico))

summary(eta)

hist(alcol[alcolemia == "positivo"], col = "red")

modello = glm(alcolemia ~ genere + eta + causa + tossicologico, family = binomial)
summary(modello)
xtable(summary(modello))

modello = glm(alcolemia ~ genere + eta + tossicologico,
    family = binomial)
summary(modello)
xtable(summary(modello))

modello = glm(alcolemia ~ genere + tossicologico,
    family = binomial)
summary(modello)
xtable(summary(modello))

modello = glm(alcolemia ~ genere, family = binomial)
summary(modello)
xtable(summary(modello))

finale = glm(alcolemia ~ genere, family = quasibinomial)
summary(finale)
xtable(summary(finale))