NR | CHOICE | TIME | COST | INCOME
1  | 1      | 30   | 4    | 0
2  | 0      | 20   | 6    | 24000
3  | 1      | 23   | 3    | 0
4  | 0      | 34   | 7    | 19500
...| ...    | ...  | ...  | ...

probit <- glm (CHOICE ~ TIME + COST + INCOME + SEX,
family = binomial(link = "probit"), data = dataset)