## Documentation Example 9 for PROC MIXED

```/****************************************************************/
/*          S A S   S A M P L E   L I B R A R Y                 */
/*                                                              */
/*    NAME: MIXEX9                                              */
/*   TITLE: Documentation Example 9 for PROC MIXED              */
/*          Examining Individual Test Components                */
/* PRODUCT: STAT                                                */
/*  SYSTEM: ALL                                                 */
/*    KEYS: Mixed Linear Models                                 */
/*   PROCS: MIXED                                               */
/*    DATA:                                                     */
/*                                                              */
/* SUPPORT: Tianlin Wang                                        */
/*     REF:                                                     */
/*    MISC:                                                     */
/*                                                              */
/****************************************************************/

data sp;
input Block A B Y @@;
datalines;
1 1 1  56  1 1 2  41
1 2 1  50  1 2 2  36
1 3 1  39  1 3 2  35
2 1 1  30  2 1 2  25
2 2 1  36  2 2 2  28
2 3 1  33  2 3 2  30
3 1 1  32  3 1 2  24
3 2 1  31  3 2 2  27
3 3 1  15  3 3 2  19
4 1 1  30  4 1 2  25
4 2 1  35  4 2 2  30
4 3 1  17  4 3 2  18
;
proc mixed data=sp;
class a b block;
model y = a b a*b / LComponents e3;
random block a*block;
run;
proc mixed data=sp;
class a b block ;
model y = a b a*b;
random block a*block;
estimate 'a    1' a 1 0 -1;
estimate 'a    2' a 0 1 -1;
estimate 'b    1' b   1 -1;
estimate 'a*b  1' a*b 1 -1 0  0 -1 1;
estimate 'a*b  2' a*b 0  0 1 -1 -1 1;
ods select Estimates;
run;
data polynomial;
do x=1 to 20; input y@@; output; end;
datalines;
1.092   1.758   1.997   3.154   3.880
3.810   4.921   4.573   6.029   6.032
6.291   7.151   7.154   6.469   7.137
6.374   5.860   4.866   4.155   2.711
;
proc mixed data=polynomial;
model y = x x*x x*x*x / s lcomponents htype=1,3;
run;

```