#################################################
#                                               #
# exemple de régression logistique "ratée"      #
#                                               #
#         (gH)  gilles.hunault@univ-angers.fr   #
#                                               #
#################################################



# 1. lecture des données pour lesquels on connait le groupe
# ---------------------------------------------------------

pg <- read.table("rlf.dar",head=TRUE,row.names=1)

print(pg)


# 2. tracé des données
# ---------------------

png(file="reglogi1.png",width=800,height=600)
plot(pg$TAILLE,pg$GROUPE,col=1+pg$GROUPE,pch=16,xlim=c(130,200))
dev.off()


# 3. Essai d'utilisation de glm
# -----------------------------

glm(pg$GROUPE~pg$TAILLE,family=binomial)