Logit regression, binaire

Bonjour à tous,

Je suis débutant sur Matlab et je souhaite modéliser la probabilité d'un évènement en utilisant le logit model. J'ai pour variable dépendante:
y=[0 0 0 0 1 1 1 1 0 0 0 0 0 0 0 0 0 0 1]

et pour quatre variables explicatives x1 x2 x3 et x4, toutes de même dimension.

Voici mes questions:
1) Dois je premièrement lancer une régression ayant pour variable dépendante Y (étant la série de donnée à partir de laquelle est dérivé y et qui n'est donc pas sous une forme binaire) et pour variable explicative x1 x2 x3 et x4 pour tester leur significance. Ou dois directement régresser ces variables explicatives contre y (forme binaire)?

Si la seconde option est la bonne, quelle fonction dois je utiliser avec matlab et avec quels paramètres?

2) Je ne suis pas certain des étapes à suivre pour obtenir la probabilité de mon évènement:

Si l'on assume que x1 et x2 sont "significant", les lignes qui suivent sont elles correctes?

X=[x1 x2]
[b,dev,stats] = glmfit(X,y,'binomial','link','logit')
yfit = glmval(b, X,'logit');

3) Pour finir, si les étapes précédentes sont bonnes, parmi les outputs, lequel me donne la probabilité de mon évènement?

Je m'excuse pour le possible manque de clarté des termes employés.

Merci d'avance,

Victor