Tell me more ×
Cross Validated is a question and answer site for statisticians, data analysts, data miners and data visualization experts. It's 100% free, no registration required.

PROC MIXED in SAS sometimes estimates a non-positive definite matrix of the random effects. However PROC MIXED nonetheless provides the estimates and the standard errors of the fixed effects in such a case. Can we consider that these estimates and standard errors are correct ?

share|improve this question

1 Answer

up vote 2 down vote accepted

You can make valid inferences on the fixed effects as long as the matrix $\mathbf{V}$ is positive definite, where $\mathbf{V}$ is the variance/covariance matrix of the implied marginal model.

Remember to specify the option nobound in the proc mixed statement.

share|improve this answer
1  
Thank you, I will check the $V$ matrix. Do you have a reference for your claim ? And why would I specify the nobound option ? – Stéphane Laurent Jun 21 '12 at 13:03
1  
Yes, this is my reference: Linear Mixed Models - A Practical Guide Using Statistical Software (B.West, K.B.Welch and A.T.Galecki) Chapman & Hall/CRC. If the matrix is not positive definite because, for example, one variance is estimated to be negative and you don't specify the nobound option, that variance will be forced to be equal to 0. – andrea Jun 21 '12 at 13:52
Ok, thanks. But in my example there is no zero estimated variance. – Stéphane Laurent Jun 21 '12 at 14:09
andrea, do I have to use my hands to derive the $V$ matrix or is it possible to get it from SAS ? – Stéphane Laurent Jun 21 '12 at 14:46
1  
Take a look at the option V of the random statement. – andrea Jun 21 '12 at 14:51
show 5 more comments

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.