Hi all,
I have an RCBD experiment with a two-way factorial treatment design, which are fully randomized within each block. I also don't have subsampling.
My question is: In the random statement of the proc Glimmix, should I include only block, or block block*factor(a) and block*factor(b)?
Thank you in advance for the support,
Gustavo
I would start with the more specific random statement, and reduce as needed. In other words:
RANDOM intercept A B/subject=block;
would be the first I would try. If there is sufficient data and variability to estimate all of these components, then you are done. If however, you get the message that the G matrix isn't positive definite, or there are zero values for the block*A or block*B variance components, you can reduce the statement and remove the zero-valued variance component terms.
Steve Denham
I would start with the more specific random statement, and reduce as needed. In other words:
RANDOM intercept A B/subject=block;
would be the first I would try. If there is sufficient data and variability to estimate all of these components, then you are done. If however, you get the message that the G matrix isn't positive definite, or there are zero values for the block*A or block*B variance components, you can reduce the statement and remove the zero-valued variance component terms.
Steve Denham
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!
Learn the difference between classical and Bayesian statistical approaches and see a few PROC examples to perform Bayesian analysis in this video.
Find more tutorials on the SAS Users YouTube channel.