# ==================================== # Exemple d'estimation du modèle logit # ==================================== # emploi de la procédure "glm" (pour Generalized Linear Model) # help(glm) si on veut lire l'aide de la commande "glm" # lecture de données d'abonnement abonnement = read.table("abon.dat", header=TRUE) attach(abonnement) abonnement # pour voir ou fix(abonnement) # pour éditer (refermer la fenêtre d'édition pour poursuivre) # variables # - abo (0/1) note l'abonnement à une certaine revue # - age # - sexe (0/1) # on estime le modèle logit expliquant l'abonnement par l'âge # et le sexe, variables "quantitatives" mlogit = glm(abo~age+sexe, family=binomial(link="logit")) summary(mlogit)