Bioequivalence and Bioavailability Forum

Hi all, especially Detlew

Do I need to do anything out of the ordinary in SAS to analyze a 3 treatment study dataset having one of the treatment arms removed?

I use PROC GLM to analyze a BE dataset from a 3-period, 3-treatment, 6 sequence crossover study. For some reason when I removed the data from one treatment arm (I just put "." for the data, leaving the "subject, period, treatment, sequence" info intact), SAS failed to generate results for T/R Ratio, 90% CI for the treatment that I kept.

SAS errror message (with respect to PROC GLM):
NOTE: TYPE I EMS not available without the E1 option.
WARNING: More coefficients than levels specified for effect TREATMENT. Some coefficients will be ignored.
NOTE: Test2 - Reference is not estimable.
NOTE: The data set WORK.L_S_MEANS has 5 observations and 9 variables.

WARNING: Output 'Estimates' was not created. Make sure that the output object name, label, or path is spelled correctly. Also, verify that the appropriate procedure options are used to produce the requested output object. For example, verify that the NOPRINT option is not used.

My SAS code:

PROC GLM data=&glmdata order=internal outstat=tanova;
  CLASS TREATMENT PERIOD SEQUENCE SUBJECT;
  MODEL &pk= SEQUENCE SUBJECT(SEQUENCE) PERIOD TREATMENT/SS1 SS3 CLPARM alpha=0.1;
  output out=pred p=yhat r=resid student=stresid;
  RANDOM SUBJECT(SEQUENCE) /TEST;
  LSMEANS TREATMENT PERIOD /stderr pdiff out=lsmeanl;
  estimat = ESTIMATE "Test2 - Reference" TREATMENT 0 1 -1;
run;

Thanks
John