Programming the statistical procedures from SAS

Need Help with Proc MIXED, nested random factor removal

Accepted Solution Solved
Reply
Occasional Contributor
Posts: 16
Accepted Solution

Need Help with Proc MIXED, nested random factor removal

[ Edited ]

Hello,
I run an analysis by using PROC MIXED as implemented in the EG 7.1 software. I have the code below:


PROC MIXED DATA = WORK.SORTTempTableSorted covtest
PLOTS(ONLY)=ALL
METHOD=REML;
CLASS IsolStr    IsolRepeatStr    IsolVialStr    InfStaStr;
MODEL TotFec=InfStaStr;
RANDOM   IsolStr     IsolRepeatStr(IsolStr)     IsolVialStr(IsolRepeatStr   IsolStr ) / TYPE=VC;
RUN; QUIT;

 

As can be seen in the code, I have a fixed factor, "InfStaStr," and several nested random factors. The random factors "IsolStr" and "IsolVialStr(IsolRepeatStr IsolStr)" had a varainace zero when the full model was run, but the nested term "IsolRepearStr(IsolStr)" did not a have a variance of zero (see attached file), and I get the "NOTE: Estimated G matrix is not positive definite."


I read in the SAS manual that in a case like this, random factors with variance zero can be romoved from the model. If I follow that advise, I would be left with a model with a single nested factor, "IsolRepearStr(IsolStr)," where the factor it is nested whitin, "IsolStr," is not present.

 

So, my question is: is it statistically sound to remove "IsolStr" in this case? Thank you for any help on this matter.
Igforek


Accepted Solutions
Solution
‎07-16-2016 02:13 AM
Occasional Contributor
Posts: 16

Re: Need Help with Proc MIXED, nested random factor removal

OK, it seems reasonable. I checked and there is not much change in the output when removing the random nested factors with value zero from the model. Thank you for your answer.

View solution in original post


All Replies
Valued Guide
Valued Guide
Posts: 684

Re: Need Help with Proc MIXED, nested random factor removal

For nested variance component (VC) effects, you can leave the 0-variance terms in the model. The results will be the same with or without them. Also, there is nothing wrong with leaving in a single random effects term as you described.

Solution
‎07-16-2016 02:13 AM
Occasional Contributor
Posts: 16

Re: Need Help with Proc MIXED, nested random factor removal

OK, it seems reasonable. I checked and there is not much change in the output when removing the random nested factors with value zero from the model. Thank you for your answer.

☑ This topic is solved.

Need further help from the community? Please ask a new question.

Discussion stats
  • 2 replies
  • 241 views
  • 0 likes
  • 2 in conversation