Pregunta que me han hecho hoy: cómo hacer un test de bondad de ajuste con SAS. Aquí la respuesta que he dado:

data datos_aleatorios;
    do i = 1 to 200000;
        * GENERAMOS UNAS VARIABLES ALEATORIAS;
        variable_gamma = rangam(89, 450);
        variable_exponencial = ranexp(23) * 100 + 0.17045;
        output;
    end;
run;

* ods select ParameterEstimates GoodnessOfFit;
proc univariate data=datos_aleatorios;
   var variable_:;
   histogram / gamma;
run;

Mucho cuidado con estos test de hipótesis. Yo suelo conformarme con ver la tabla de cuantiles. Saludos.