data sim;
    do i = 1 to 100;
        x1 = rand("Normal");
        x2 = rand("Binomial", 0.5, 100);
        output;
    end;
run;

proc means data=sashelp.class;
    class sex;
    var height weight;
    output out = mean_by_sex;
run;