turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

Find a Community

- Home
- /
- SAS Programming
- /
- SAS Procedures
- /
- Proc Glimmix random statement

Topic Options

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

04-29-2016 01:57 PM

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

Accepted Solutions

Solution

05-05-2016
06:16 PM

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

Posted in reply to Gubardella

05-02-2016 09:16 AM

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

All Replies

Solution

05-05-2016
06:16 PM

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

Posted in reply to Gubardella

05-02-2016 09:16 AM

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

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

Posted in reply to SteveDenham

05-02-2016 09:50 PM

Thank you for the help, Steve.