Action Principale : frontierCost
| 1 | DATA mycas.cost_data; |
| 2 | call streaminit(12345); |
| 3 | DO i = 1 to 500; |
| 4 | log_q = rand('Normal', 5, 1); /* Log de la production */ |
| 5 | log_p1 = rand('Normal', 2, 0.5); /* Log prix intrant 1 */ |
| 6 | log_p2 = rand('Normal', 3, 0.5); /* Log prix intrant 2 */ |
| 7 | v = rand('Normal', 0, 0.1); /* Erreur aléatoire */ |
| 8 | u = rand('Exponential'); /* Inefficacité de coût (positive) */ |
| 9 | /* Modèle de coût: C = f(Q, P) + v + u */ |
| 10 | log_cost = 1 + 0.7*log_q + 0.4*log_p1 + 0.6*log_p2 + v + u; |
| 11 | OUTPUT; |
| 12 | END; |
| 13 | RUN; |