The MIANALYZE Procedure

Example 77.8 Reading Mixed Model Results with Classification Covariates

(View the complete code for this example.)

This example creates data sets that contains parameter estimates and corresponding covariance matrices with classification variables computed by a mixed regression model analysis for a set of imputed data sets. These estimates are then combined to generate valid statistical inferences about the model parameters.

The following statements use PROC MIXED to generate the parameter estimates and covariance matrix for each imputed data set:

ods select none;
proc mixed data=outfish2;
   class Species;
   model Length= Species Width/ solution;
   by _Imputation_;
   ods output SolutionF=mxparms;
run;
ods select all;

Because of the ODS SELECT statements, no output is displayed. The following statements display (in Output 77.8.1) the output mixed model coefficients from PROC MIXED for the first two imputed data sets:

proc print data=mxparms (obs=8);
   var _Imputation_ Effect Species Estimate StdErr;
   title 'MIXED Model Coefficients (First Two Imputations)';
run;

Output 77.8.1: PROC MIXED Model Coefficients

MIXED Model Coefficients (First Two Imputations)

Obs_Imputation_EffectSpeciesEstimateStdErr
11Intercept 4.51060.8244
21SpeciesParkki1.57740.7020
31SpeciesPerch0.
41Width 5.25850.1599
52Intercept 4.52500.8771
62SpeciesParkki1.48850.7693
72SpeciesPerch0.
82Width 5.23890.1701


The following statements use the MIANALYZE procedure with an input PARMS= data set:

proc mianalyze parms(classvar=full)=mxparms;
   class Species;
   modeleffects Intercept Species Width;
run;

The "Variance Information" table in Output 77.8.2 displays the between-imputation, within-imputation, and total variances for combining complete-data inferences.

Output 77.8.2: Variance Information

The MIANALYZE Procedure

Variance Information (25 Imputations)
ParameterSpeciesVarianceDFRelative
Increase
in Variance
Fraction
Missing
Information
Relative
Efficiency
BetweenWithinTotal
Intercept 0.0656650.6686670.7369592794.90.1021310.0933160.996281
SpeciesParkki0.0772910.5256400.6060231364.10.1529240.1339090.994672
SpeciesPerch0......
Width 0.0022760.0258760.02824334160.0914880.0843560.996637


The "Parameter Estimates" table in Output 77.8.3 displays the combined parameter estimates with associated standard errors.

Output 77.8.3: Parameter Estimates

Parameter Estimates (25 Imputations)
ParameterSpeciesEstimateStd Error95% Confidence LimitsDFMinimumMaximumTheta0t for H0:
Parameter=Theta0
Pr > |t|
Intercept 4.5192590.8584632.835976.2025452794.94.0136245.16246305.26<.0001
SpeciesParkki1.2779020.778475-0.249242.8050391364.10.6478981.71055901.640.1009
SpeciesPerch0....000..
Width 5.2842850.1680564.954785.61378634165.2065165.369231031.44<.0001