BE test in S-plus

["He Sun 301-827-5919 FAX t-" <SUNH@cder.fda.gov>]

    Dear S-Plus users:

    A practical SAS code for conducting bioequivalence analysis (BE) on
    maximum concentration (cmax) for a two-way crossover trial using SAS
    PROC GLM could be:

    Proc glm data=data1;
    class subj seq per trt;
    model cmax=seq subj(seq) per trt/ss1 ss2 ss3 ss4;
    test h=seq e=subj(seq);
    estimate 'trt 1 vs trt 2' trt 1 -1;
    lsmeans trt/s;

    Could anyone on the list please advise me how to perform the same
    analysis using S-plus? Or please provide me the proper reading resource
    to study.  When I use <-glm(cmax~seq+subj%in%seq+per+trt) function, I
    cannot get the same results as using SAS.

    Thanks a million for your help and Happy Holidays!

    Henry Sun
    Pharmacometrics
    FDA
    e-mail  SUNH@cder.fda.gov