################################################################################

PROC MIXED data=df0;
CLASSES subject sequence period formulation;
MODEL  var = sequence period formulation/ DDFM=SATTERTH LComponents E solution covb covbi;
CONTRAST 'CONTR' formulation    1 -1;
RANDOM  formulation/TYPE=CSH SUB=subject G V;
REPEATED/GRP=formulation SUB=subject R;
ESTIMATE 'T vs. R' formulation 1 -1/CL ALPHA=0.1 E;
RUN;

REML: 10.06523862

################################################################################

PROC MIXED data=df0;
CLASSES subject sequence period formulation;
MODEL  var = sequence period formulation/ DDFM=SATTERTH LComponents E solution covb covbi;
CONTRAST 'CONTR' formulation    1 -1;
RANDOM  formulation/TYPE=VC SUB=subject G V;
REPEATED/GRP=formulation SUB=subject R;
ESTIMATE 'T vs. R' formulation 1 -1/CL ALPHA=0.1 E;
RUN;

REML: 16.06148160

################################################################################

PROC MIXED data=df0;
CLASSES subject sequence period formulation;
MODEL  var = sequence period formulation/ DDFM=SATTERTH LComponents E solution covb covbi;
CONTRAST 'CONTR' formulation    1 -1;
RANDOM  subject/TYPE=VC G V;
ESTIMATE 'T vs. R' formulation 1 -1/CL ALPHA=0.1 E;
RUN;

REML: 10.86212458

################################################################################

PROC MIXED data=df0;
CLASSES subject sequence period formulation;
MODEL  var = sequence period formulation/ DDFM=SATTERTH LComponents E solution covb covbi;
CONTRAST 'CONTR' formulation    1 -1;
RANDOM  period/TYPE=VC G V;
RANDOM  formulation/TYPE=VC G V;
ESTIMATE 'T vs. R' formulation 1 -1/CL ALPHA=0.1 E;
RUN;

REML: 25.12948063

################################################################################

PROC MIXED data=df0;
CLASSES subject sequence period formulation;
MODEL  var = sequence period formulation/ DDFM=SATTERTH LComponents E solution covb covbi;
CONTRAST 'CONTR' formulation    1 -1;
RANDOM  formulation/TYPE=UN(1) SUB=subject G V;
ESTIMATE 'T vs. R' formulation 1 -1/CL ALPHA=0.1 E;
RUN;

REML: 16.24111264

################################################################################
