Proc Mixed for Random Blocks

suebayran — Tue, 03 May 2016 06:16:57 GMT

Hello,

I have a dataset with a continuous outcome variable and 3 treatments (1 of which is placebo). I have about 60 observations in my dataset; 60 subjects each of which got 1 treatment and was measured once - the dataset is balanced but measurements were taken in blocks of days but only 2 measurements were taken each day (incomplete blocks?). I need to include the block in my model as a random effect. I am trying to get the effect of the drug and im running a proc mixed with a random block effect, but i am not getting an anova table and am getting an error that the "Estimated G matrix is not positive definite." Can anyone tell me what I am doing wrong? Thank you.

Re: Proc Mixed for Random Blocks

SteveDenham — Thu, 05 May 2016 17:54:49 GMT

Please share your PROC MIXED code and enough of your data to enable some trouble shooting. I suspect that you may have misspecified the random effect, as the message 'Estimated G matrix is not positive definite' means that at least one variance parameter is being set to zero. There are a couple of ways around this, but it depends on design, data structure and code.

Steve Denham

topic Re: Proc Mixed for Random Blocks in SAS Procedures

Proc Mixed for Random Blocks

Re: Proc Mixed for Random Blocks