# Essais_loglin.r

load("PoidsBebe.rda")
summary(PoidsBebe)

PB <- PoidsBebe

# Utilisation de loglm
library(MASS)

Ta <- xtabs(~PB$LOW+PB$HT+PB$RACE)
Ta

M1 <- loglm(~PB$LOW+PB$HT+PB$RACE,data=Ta)
M1
summary(M1)

M2 <- loglm(~(PB$LOW+PB$HT+PB$RACE)^2,data=Ta)
M2

M3 <- loglm(~(PB$LOW+PB$HT+PB$RACE)^2-PB$LOW:PB$HT,data=Ta)
M3


# Utilisation de glm
DF <- as.data.frame(Ta)
DF

M1b <- glm(Freq~PB.LOW+PB.HT+PB.RACE,data=DF,family=poisson)
M1b
summary(M1b)

M1b2 <- glm(Freq~relevel(PB.LOW,2)+PB.HT+PB.RACE,data=DF,family=poisson)
summary(M1b2)


M2b <- glm(Freq~(PB.LOW+PB.HT+PB.RACE)^2,data=DF,family=poisson)
M2b
summary(M2b)

M3b <- glm(Freq~(PB.LOW+PB.HT+PB.RACE)^2-PB.LOW:PB.HT,data=DF,family=poisson)
M3b
summary(M3b)

anova(M2b)
anova(M1,M3,M2)
anova(M1b,M3b,M2b)