Documentation Example 7 for PROC MIXED

/****************************************************************/
/*          S A S   S A M P L E   L I B R A R Y                 */
/*                                                              */
/*    NAME: MIXEX7                                              */
/*   TITLE: Documentation Example 7 for PROC MIXED              */
/*          Influence in Heterogeneous Variance Model           */
/* PRODUCT: STAT                                                */
/*  SYSTEM: ALL                                                 */
/*    KEYS: Mixed Linear Models ODS Graphics                    */
/*   PROCS: MIXED, PRINT                                        */
/*    DATA:                                                     */
/*                                                              */
/* SUPPORT: Tianlin Wang                                        */
/*     REF:                                                     */
/*    MISC: Influence diagnostics                               */
/*                                                              */
/****************************************************************/

*---Influence in Heterogeneous Variance Model (Experimental)---*
| A one-way classification model with heterogeneous            |
| variances is fit. Data from Snedecor and Cochran (1980),p216 |
*--------------------------------------------------------------*;

data absorb;
   input FatType Absorbed @@;
   datalines;
 1 164  1 172  1 168  1 177  1 156  1 195
 2 178  2 191  2 197  2 182  2 185  2 177
 3 175  3 193  3 178  3 171  3 163  3 176
 4 155  4 166  4 149  4 164  4 170  4 168
;

ods graphics on;

proc mixed data=absorb asycov;
   class FatType;
   model Absorbed = FatType / s
                    influence(iter=10 estimates);
   repeated / group=FatType;
   ods output Influence=inf;
run;

ods graphics off;

proc print data=inf label;
   var parm1-parm5 covp1-covp4;
run;

proc print data=inf label;
   var observed predicted residual pressres
       student Rstudent;
run;

proc print data=inf label;
   var leverage observed CookD DFFITS CovRatio RLD;
run;

proc print data=inf label;
   var iter CookDCP CovRatioCP;
run;