Action Principale : bartProbit
| 1 | DATA casuser.simulated_data(caslib='casuser'); |
| 2 | DO i = 1 to 1000; |
| 3 | array x[10]; |
| 4 | DO j = 1 to 10; |
| 5 | x[j] = rannor(12345); |
| 6 | END; |
| 7 | mu = 10 * sin(3.14 * x[1] * x[2]) + 20 * (x[3] - 0.5)**2 + 10 * x[4] + 5 * x[5]; |
| 8 | pr = 1 / (1 + exp(-mu)); |
| 9 | y = rand('BERNOULLI', pr); |
| 10 | OUTPUT; |
| 11 | END; |
| 12 | RUN; |