I have run this code in Proc Mixed for repeated measurements on subjects=Flask:
proc mixed data=FRF method=reml cl ic covtest;
class CCI4 CHCI3 Time Flask ;
model Leakage=CCI4 Time CCI4*Time ;
random CHCI3 CCI4*CHCI3 CHCI3*Time CCI4*CHCI3*Time /s;
repeated / subject=Flask(CCI4*CHCI3) type=arh(1) r;
lsmeans CCI4 / pdiff cl adjust=tukey;
run;
And part of the output shows the G Matrix below. How do we interpret the ARH(1) subject effect Flask(CCI4*CHCHI3) and its Estimate 0.9730?
What you are asking for might not be applicable for a complex covariance model like yours.
Jill
0.9730 in the Covariance Parameter Estimates table is not an element in the G matrix. It is a parameter in the R matrix. It is the correlation in the residuals between two adjacent observations. See the documentation below for more information on the ARH(1) structure. The estimated value for rho is 0.9730.
Also, your G matrix is not positive definite. You might want to take out some of the random effects in your RANDOM statement.
Hope this helps,
Jill
@jiltao Thank you for the material. You are very correct. Does the non-positive definite matrix impact negatively on the covariance parameters estimated by the model, or it is just a warning that something in not right? I relised that the more I include some random parameters the more i get problems with convergence as well.
G matrix not positive definite might not be devastating, but it is always a good idea to take care of this so you have a clean convergence. If you take out the random effects that have 0 estimated variance, this message is likely going away. The estimates for other effects might remain the same.
Hi @jiltao - removing the random effects when there is a non-positive G matrix seems to me to be a problem for randomized complete block designs, where the block estimate is zero. Removal shifts the denominator degrees of freedom from containment to residual, so the analysis no longer reflects the blocking in the design. Is the answer then to use a ddf= option to set the denominator degrees of freedom to what would be expected for the design?
SteveDenham
Actually, as far as I know, there is no consensus on what to do when the G matrix is not positive definite. While some advice removing the "offending" random effects, others argue the modeling of the random effects is dictated by the design (like a RCB design), and therefore should not change the model effect specifications. If there is only one random effect with 0 estimated variance, one approach is to re-specify that random effect with the statement -- repeated / subject=<random_effect> type=cs;
If the default ddf value becomes inappropriate when a random effect is removed, then you might use ddf= option to specify your own denominator degrees of freedom value, or use other estimation method such as ddfm=kr if that helps.
Jill
@jiltao @SteveDenham Thank you for the insights, I am following. However, based on the posted G and R matrices, I have further questions:
1. How do I calculate the percentage contribution of each estimated variance component?
2. How do we get the overall error variance when heterogeneous error variance for each subject is displayed?
What you are asking for might not be applicable for a complex covariance model like yours.
Jill
Are you ready for the spotlight? We're accepting content ideas for SAS Innovate 2025 to be held May 6-9 in Orlando, FL. The call is open until September 25. Read more here about why you should contribute and what is in it for you!
ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.
Find more tutorials on the SAS Users YouTube channel.