DATA bustilfr;
INPUT fors mutil util bog tilf mtilf ialt;
Utilfr = mutil + util ;
andel = Utilfr/ialt ;
logit = log(andel/(1-andel));
vegt = ialt*andel*(1-andel);
CARDS;
0  234 559 1157 5826 2553 10329
2  41  100  145  602  237  1125
5  42   76   89  254   72   533
7  35   48   39   95   27   244
;
RUN;
QUIT;