The MIXED Procedure

requests that a t-type confidence interval be constructed for each of the random-effect estimates with confidence level 1 – number. The value of number must be between 0 and 1; the default is 0.05.

requests that t-type confidence limits be constructed for each of the random-effect estimates. The confidence level is 0.95 by default; this can be changed with the ALPHA= option.

requests that the estimated matrix be displayed. PROC MIXED displays blanks for values that are 0. If you specify the SUBJECT= option, then the block of the matrix corresponding to the first subject is displayed. The ODS name of the table is G.

displays the lower-triangular Cholesky root of the estimated matrix according to the rules listed under the G option. The ODS name of the table is CholG.

displays the inverse Cholesky root of the estimated matrix according to the rules listed under the G option. The ODS name of the table is InvCholG.

displays the correlation matrix corresponding to the estimated matrix according to the rules listed under the G option. The ODS name of the table is GCorr.

requests that the matrix be read in from a SAS data set. This matrix is assumed to be known; therefore, only -side parameters from effects in the REPEATED statement are included in the Newton-Raphson iterations. If no REPEATED statement is specified, then only a residual variance is estimated.

The information in the GDATA= data set can appear in one of two ways. The first is a sparse representation for which you include Row, Col, and Value variables to indicate the row, column, and value of , respectively. All unspecified locations are assumed to be 0. The second representation is for dense matrices. In it you include Row and Col1–Coln variables to indicate, respectively, the row and columns of , which is a symmetric matrix of order n. For both representations, you must specify effects in the RANDOM statement that generate a matrix that contains n columns. (See Example 84.4.)

If you have more than one RANDOM statement, only one GDATA= option is required in any one of them, and the data set you specify must contain the entire matrix defined by all of the RANDOM statements.

If the GDATA= data set contains variance ratios instead of the variances themselves, then use the RATIOS option.

Known parameters of can also be input by using the PARMS statement with the HOLD= option.

displays the inverse of the estimated matrix according to the rules listed under the G option. The ODS name of the table is InvG.

defines an effect specifying heterogeneity in the covariance structure of . All observations having the same level of the group effect have the same covariance parameters. Each new level of the group effect produces a new set of covariance parameters with the same structure as the original group. You should exercise caution in defining the group effect, because strange covariance patterns can result from its misuse. Also, the group effect can greatly increase the number of estimated covariance parameters, which can adversely affect the optimization process.

Continuous variables are permitted as arguments to the GROUP= option. PROC MIXED does not sort by the values of the continuous variable; rather, it considers the data to be from a new subject or group whenever the value of the continuous variable changes from the previous observation. Using a continuous variable decreases execution time for models with a large number of subjects or groups and also prevents the production of a large "Class Level Information" table.

reads the coefficient matrices associated with the TYPE=LIN(number) option. The data set must contain the variables Parm, Row, Col1–Coln or Parm, Row, Col, Value. The Parm variable denotes which of the number coefficient matrices is currently being constructed, and the Row, Col1–Coln, or Row, Col, Value variables specify the matrix values, as they do with the GDATA= option. Unspecified values of these matrices are set equal to 0.

eliminates the columns in corresponding to missing levels of random effects involving CLASS variables. By default, these columns are included in .

indicates that ratios with the residual variance are specified in the GDATA= data set instead of the covariance parameters themselves. The default GDATA= data set contains the individual covariance parameters.

requests that the solution for the random-effects parameters be produced. Using notation from the section Mixed Models Theory, these estimates are the empirical best linear unbiased predictors (EBLUPs) . They can be useful for comparing the random effects from different experimental units and can also be treated as residuals in performing diagnostics for your mixed model.

The numbers displayed in the SE Pred column of the "Solution for Random Effects" table are not the standard errors of the displayed in the Estimate column; rather, they are the standard errors of predictions , where is the ith EBLUP and is the ith random-effect parameter.

identifies the subjects in your mixed model. Complete independence is assumed across subjects; thus, for the RANDOM statement, the SUBJECT= option produces a block-diagonal structure in with identical blocks. The matrix is modified to accommodate this block diagonality. In fact, specifying a subject effect is equivalent to nesting all other effects in the RANDOM statement within the subject effect.

Continuous variables are permitted as arguments to the SUBJECT= option. PROC MIXED does not sort by the values of the continuous variable; rather, it considers the data to be from a new subject or group whenever the value of the continuous variable changes from the previous observation. Using a continuous variable decreases execution time for models with a large number of subjects or groups and also prevents the production of a large "Class Level Information" table.

When you specify the SUBJECT= option and a classification random effect, computations are usually much quicker if the levels of the random effect are duplicated within each level of the SUBJECT= effect.

specifies the covariance structure of . Valid values for covariance-structure and their descriptions are listed in Table 17 and Table 18. Although a variety of structures are available, most applications call for either TYPE=VC or TYPE=UN. The TYPE=VC (variance components) option is the default structure, and it models a different variance component for each random effect.

The TYPE=UN (unstructured) option is useful for correlated random coefficient models. For example, the following statement specifies a random intercept-slope model that has different variances for the intercept and slope and a covariance between them:

random intercept age / type=un subject=person;

You can also use TYPE=FA0(2) here to request a estimate that is constrained to be nonnegative definite.

If you are constructing your own columns of with continuous variables, you can use the TYPE=TOEP(1) structure to group them together to have a common variance component. If you want to have different covariance structures in different parts of , you must use multiple RANDOM statements with different TYPE= options.

requests that blocks of the estimated matrix be displayed. The first block determined by the SUBJECT= effect is the default displayed block. PROC MIXED displays entries that are 0 as blanks in the table.

You can optionally use the value-list specification, which indicates the subjects for which blocks of are to be displayed. For example, the following statement displays block matrices for the first, third, and seventh persons:

random int time / type=un subject=person v=1,3,7;

The ODS table name is V.

displays the Cholesky root of the blocks of the estimated matrix. The value-list specification is the same as in the V option. The ODS table name is CholV.

displays the inverse of the Cholesky root of the blocks of the estimated matrix. The value-list specification is the same as in the V option. The ODS table name is InvCholV.

displays the correlation matrix corresponding to the blocks of the estimated matrix. The value-list specification is the same as in the V option. The ODS table name is VCorr.

displays the inverse of the blocks of the estimated matrix. The value-list specification is the same as in the V option. The ODS table name is InvV.