Action Principale : countregFitModel
| 1 | DATA mycas.d_counts; |
| 2 | call streaminit(12345); |
| 3 | DO i = 1 to 1000; |
| 4 | Sex = ifn(rand('UNIFORM') > 0.5, 'F', 'M'); |
| 5 | Age = 20 + floor(rand('UNIFORM') * 40); |
| 6 | x1 = rand('UNIFORM'); |
| 7 | x2 = rand('NORMAL'); |
| 8 | lambda = exp(-1 + 0.5 * (Sex='M') + 0.02 * Age + 0.1 * x1 + 0.3 * x2); |
| 9 | NumInjur = rand('POISSON', lambda); |
| 10 | OUTPUT; |
| 11 | END; |
| 12 | RUN; |