options ls=75 ps=120 nodate;

data one;
    do fe=2 to 3;
        do j=1 to 3;
        if j=1 then dose=0.3;
        else if j=2 then dose=1.2;
        else dose=10.2;
        input y @;
        ly=log(y);
        y2=y*y;
        sy=sqrt(y);
        trt=10*fe+dose;
        output;
        end;
    end;
cards; 
            0.71  2.20  2.25  2.20  4.04  2.71
            1.66  2.93  3.93  2.69  4.16  5.43
            2.01  3.08  5.08  3.54  4.42  6.38
            2.16  3.49  5.82  3.75  4.93  6.38
            2.42  4.11  5.84  3.83  5.49  8.32
            2.42  4.95  6.89  4.08  5.77  9.04 
            2.56  5.16  8.50  4.27  5.86  9.56
            2.60  5.54  8.56  4.53  6.28  10.01
            3.31  5.68  9.44  5.32  6.97  10.08
            3.64  6.25  10.52 6.18  7.02  10.62
            3.74  7.25  13.46 6.22  7.78  13.80
            3.74  7.90  13.57 6.33  9.23  15.99
            4.39  8.85  14.76 6.97  9.34  17.90
            4.50  11.96 16.41 6.97  9.91  18.25
            5.07  15.54 16.96 7.52  13.46 19.32
            5.26  15.89 17.56 8.36  18.4  19.87
            8.15  18.3  22.82 11.65 23.89 21.60
            8.24  18.59 29.13 12.45 26.39 22.25
;
run;

proc print;
   var fe dose y;
run;