Régression logistique [RESOLU]

Bonjour,

j'essai de taper sous matlab l'algorithme de la régression logistique (maximisation de vraisemblance + Newton-Raphson).

En fait mon objectif est dans un premier temps l'écrire puis de comparer mes résultats à ceux qu'on a via SAS.

Mais voilà en faisant l'analogie à la proc logistic justement, je me suis rendu compte que je voyais pas l'influence des références des modalités des variables, enfin je veux dire, où dans l'algorithme on a une modification des calculs en fonction d'elles.

Si quelqu'un pouvait m'expliquer, ou si quelqu'un connait un site ou je peux télécharger l'algo sous matlab pour pouvoir l'étudier.

Merci d'avance.

en fait c'est dans l'établissement de la matrice design que l'influence de la modalité de référence va se retrouver, parce que selon la modalité de référence choisie ta matrice sera différente :

Code:

# 1 exemple :
            y x
1  -1.3447118 1
2  -0.0648961 1
3  -1.6578271 1
4  -0.3853041 2
5  -0.4927110 2
6  -0.2557316 2
7  -0.9612893 3
8  -0.8368667 3
9  1.3924498 3
10 -0.8621528 1


Code:

# si 1 est la modalité de référence de x : 
(Intercept) x2 x3
1            1  0  0
2            1  0  0
3            1  0  0
4            1  1  0
5            1  1  0
6            1  1  0
7            1  0  1
8            1  0  1
9            1  0  1
10          1  0  0


Code:

# si 2 est la modalité de référence de x
  (Intercept) x1 x3
1            1  1  0
2            1  1  0
3            1  1  0
4            1  0  0
5            1  0  0
6            1  0  0
7            1  0  1
8            1  0  1
9            1  0  1
10          1  1  0

etc...